|
1) 设x为自变量,y为因变量,两者满足方程
" m6 W& A) f2 f$ N7 @& L, y(x=0) = 1+ Y Y. Q9 F, t% N% i& `
用数值方法求在[0,10]区间,步长为1的x点所对应的y值,即x=0:10+ A4 x! E+ r. F! W
并用这些数据点生成插值多项式公式,求[0,10]区间任意一点的值。
! W" D# A4 q8 |7 S% N- K; ~2 c
F9 F+ M5 l% J要求:! v1 c4 W; x* `
编写常微分方程的四阶Runge-Kutta求解函数,和Matlab内建的ode45对比/ }- _; Q. z" m- s0 S3 _
编写Lagrange插值函数,要求支持任意多的输入点8 r; c J# ^7 `; P" U$ o! x
0 \2 c& J& s1 ?2 C: R* e) H5 c
2) 计算 在区间[-5, 5]上的值,x的步长为1,对求出的数据点(x,y)用上述Lagrange函数生成插值多项式。在全区间上比较通过插值多项式和原函数计算得到的结果的差异,并设法改进 |
|