《表1 Lorenz方程的与MSRM (h=0.01) [6]相比较的数值模拟结果Tab.1 Numerical solution for the Lorenz system compared wit
利用多级移位Chebyshev配置方法,在[0,10]求解Lorenz方程组(1),在Intel(R)Core(TM)CPU 3.40GHz,Matlab R 2014a环境下进行计算.对于初值为(1,5,10)时,取MSSM步长h=0.05,在每个子区间,取N=10,计算结果如表1所示.文献[6]中MSRM的步长h=0.01,在精度基本保持一致的前提下,MSSM的计算时间仅为MSRM的45.8%,3个分量的计算结果如图1所示.实际图像分别在xoy平面与xoz平面的投影见图2.从表1可以看出,在数值解的有效数字位数基本保持不变时,MSSM方法的计算时间明显优于MSRM方法,这是因为MSSM方法的计算过程只利用了右端函数的Chebyshev级数的系数导出待求函数的Chebyshev级数的系数,进而用Chebyshev级数近似待求函数,整个过程避免了方程组的求解.
图表编号 | XD0026870800 严禁用于非法目的 |
---|---|
绘制时间 | 2018.11.28 |
作者 | 杨录峰、马宁 |
绘制单位 | 北方民族大学数学与信息科学学院、吴忠市气象局 |
更多格式 | 高清、无水印(增值服务) |