|
1) 设x为自变量,y为因变量,两者满足方程
; R# A) f2 t# E3 q4 S; O8 @$ f, y(x=0) = 1, @3 P/ Y5 S6 M0 w
用数值方法求在[0,10]区间,步长为1的x点所对应的y值,即x=0:10
/ q- O' B1 K" p% G4 S2 Q7 c并用这些数据点生成插值多项式公式,求[0,10]区间任意一点的值。; s1 i# P" d& U
7 u( K& j2 C: d/ H/ I. P( B$ Q8 e2 j
要求:+ \8 \/ N4 Z* w2 H* v8 d
编写常微分方程的四阶Runge-Kutta求解函数,和Matlab内建的ode45对比+ V. ?$ Q9 L' `: L
编写Lagrange插值函数,要求支持任意多的输入点
% @& I- S- y- y2 G* t# X 8 N& t7 d5 T- q
2) 计算 在区间[-5, 5]上的值,x的步长为1,对求出的数据点(x,y)用上述Lagrange函数生成插值多项式。在全区间上比较通过插值多项式和原函数计算得到的结果的差异,并设法改进 |
|