论文部分内容阅读
本文主要基于肖小涛博士开发的一个经过检验的适用于长时间尺度的导心轨道模拟程序,运用离散傅立叶变换的方法,得到了性能更优越的计算程序。
为了满足对称性以及求高阶导数的要求,本文在原有的程序的础上,对磁场等物理量的离散化信息进行了离散傅立叶变换,用三角多项式逼近了文中的所有二元函数,理论上可以求到任意阶导数。新的计算方法,和传统的三次样条插值和B样条插值进行对比,得到了完全相同的轨道模拟结果。并且可以清晰的看到各个物理量在各级次的上的变化规律。而且保证了角向上的对称性,还有滤波的效果,计算稳定性更好。
本文得出结论,用傅立叶变换来产生磁场等物理量的方法是完全可行的。同时,针对采样效应带来的拟合误差的来源进行分析,最后对快速傅立叶变换带来的一些问题进行了讨论,给出了FFT应用到托卡马克中物理量的计算中要注意的问题。