sunminmin 发表于 2012-5-15 21:45

插值——微分方程求解




ptfile:///C:/DOCUME~1/ADMINI~1/LOCALS~1/Temp/ksohtml/wps_clip_image-3445.png
x0y0
300025167.79
518.940.4773443535.23
835.880.9546770124.16
1246.41.432104563.8
17031.9093142869.1
2114.82.3867177416.1
23902.864200503.4
2490.33.3414208917.8
2442.53.8187204907.7

由于膛内时期内弹道给出的平均压力曲线p-t为表格函数，因而需对file:///C:/DOCUME~1/ADMINI~1/LOCALS~1/Temp/ksohtml/wps_clip_image-18442.pngfile:///C:/DOCUME~1/ADMINI~1/LOCALS~1/Temp/ksohtml/wps_clip_image-2785.png这一离散函数进行插值，以得到对应任意时刻的炮膛合力，便于微分方程的求解。插值计算用拉格朗日三点插值公式进行一元三点不等距插值，计算公式为file:///C:/DOCUME~1/ADMINI~1/LOCALS~1/Temp/ksohtml/wps_clip_image-16392.png其中file:///C:/DOCUME~1/ADMINI~1/LOCALS~1/Temp/ksohtml/wps_clip_image-26432.pngfile:///C:/DOCUME~1/ADMINI~1/LOCALS~1/Temp/ksohtml/wps_clip_image-7812.png   j=2,3,……,n-2

sunminmin 发表于 2012-5-17 11:36

我写了个程序，但不正确。请各位指点下！

sl0412 发表于 2012-5-17 12:22

学习学习
页: [1]
查看完整版本: 插值——微分方程求解