|
|

楼主 |
发表于 2010-3-30 20:25:25
|
显示全部楼层
计算差分格式系数的小程序(第3讲的课件)
一个Fortran小程序,用来求差分格式中的系数。
使用方法:
按照提示,输入基架点的个数m,基架点最左边点的位置k及求导数的阶数m_order。
例如想计算一阶导数差分格式的系数,则输入m_order=1. 想使用从j-3到 ,j+2 这6个点,则输入m=6, k=3
即:
f'(j)=[a(1)*f(j-k)+a2*f(j-k+1)+....+a(m)*f(j-k+m-1)]/Dx
F(j+1/2)=[b1*f(j-k+1)+b2*f(j-k+2)+.....+b(m-1)*f(j-k+m-1)]/Dx
程序运行后 输出a(i) 及 b(i)
程序还可以把小数化成分数。这样格式的系数更加好记。
[ 本帖最后由 lixl-imech 于 2010-3-30 12:29 编辑 ] |
|