该用户从未签到
|
1) 设x为自变量,y为因变量,两者满足方程8 p' G0 ?! A% C2 [3 l# ^! i% C2 g
, y(x=0) = 1
: u4 d9 E9 @6 u用数值方法求在[0,10]区间,步长为1的x点所对应的y值,即x=0:10
2 u9 U- y2 M. U/ O( F9 @: M并用这些数据点生成插值多项式公式,求[0,10]区间任意一点的值。. ]! h; |; }/ c
1 j) S5 E5 x) v8 _/ }
要求:0 A( T8 C0 S& Y2 c- C1 ] N, r
编写常微分方程的四阶Runge-Kutta求解函数,和Matlab内建的ode45对比
& ~9 U b$ W2 c ~3 O+ P% d编写Lagrange插值函数,要求支持任意多的输入点$ k" \; c7 k& M6 F
. q3 }( ^9 ?* b
2) 计算 在区间[-5, 5]上的值,x的步长为1,对求出的数据点(x,y)用上述Lagrange函数生成插值多项式。在全区间上比较通过插值多项式和原函数计算得到的结果的差异,并设法改进 |
|