该用户从未签到
|
1) 设x为自变量,y为因变量,两者满足方程
/ Q$ K# ]7 b/ @( o, y(x=0) = 1
) ~& o) z( A+ W: w$ Z用数值方法求在[0,10]区间,步长为1的x点所对应的y值,即x=0:104 i M. n/ f' Z& y
并用这些数据点生成插值多项式公式,求[0,10]区间任意一点的值。
& J: P5 Y9 z4 V0 N8 V. S, W
. e' u- l4 H6 ]+ H要求:
' A1 v/ G; N$ u! f3 k* \编写常微分方程的四阶Runge-Kutta求解函数,和Matlab内建的ode45对比
( g% t" Q' B+ w' I' m编写Lagrange插值函数,要求支持任意多的输入点
- P) P, h7 f0 V
0 p: @% m7 r+ f( `; p) f2) 计算 在区间[-5, 5]上的值,x的步长为1,对求出的数据点(x,y)用上述Lagrange函数生成插值多项式。在全区间上比较通过插值多项式和原函数计算得到的结果的差异,并设法改进 |
|