|
1) 设x为自变量,y为因变量,两者满足方程
9 b7 I! p$ ^: K3 T9 Y, L- E, y(x=0) = 1
: l, I1 |; D7 m/ D6 B用数值方法求在[0,10]区间,步长为1的x点所对应的y值,即x=0:101 { R$ |0 v+ D- L1 I
并用这些数据点生成插值多项式公式,求[0,10]区间任意一点的值。* O- e: z( d9 L% x* J: [
* P0 N7 E$ O% }2 R+ t7 l要求:
' A; H1 D" r4 R5 ]1 |编写常微分方程的四阶Runge-Kutta求解函数,和Matlab内建的ode45对比
2 ?( Z! U9 h3 V( c编写Lagrange插值函数,要求支持任意多的输入点
4 X8 K% d( C( K! R2 M- j- w. [( p - f1 ]( a0 J1 h6 j4 m# j
2) 计算 在区间[-5, 5]上的值,x的步长为1,对求出的数据点(x,y)用上述Lagrange函数生成插值多项式。在全区间上比较通过插值多项式和原函数计算得到的结果的差异,并设法改进 |
|