|
1) 设x为自变量,y为因变量,两者满足方程
/ R" \7 ]( _6 d# X, y(x=0) = 1+ ~; g+ |8 y; A) ]
用数值方法求在[0,10]区间,步长为1的x点所对应的y值,即x=0:10
( @4 P {) H; h5 b9 P并用这些数据点生成插值多项式公式,求[0,10]区间任意一点的值。" h6 L) q" m" ]0 B' @- L8 \0 I" P
& v! n' A4 K2 o" F# _$ t F9 F要求:
8 N- b2 v2 p. j( @) ]5 C编写常微分方程的四阶Runge-Kutta求解函数,和Matlab内建的ode45对比
) `# Y% D. l9 B; i编写Lagrange插值函数,要求支持任意多的输入点
/ @+ Z+ N7 K8 m ' I" _- A- Z6 p9 e
2) 计算 在区间[-5, 5]上的值,x的步长为1,对求出的数据点(x,y)用上述Lagrange函数生成插值多项式。在全区间上比较通过插值多项式和原函数计算得到的结果的差异,并设法改进 |
|