该用户从未签到
|
1) 设x为自变量,y为因变量,两者满足方程. O$ T" x1 k8 b/ J3 [4 i Y, ^; [
, y(x=0) = 10 `/ s" I% B/ ]( z5 ^1 M2 C4 u
用数值方法求在[0,10]区间,步长为1的x点所对应的y值,即x=0:10
4 {% U- r9 U$ j$ n9 g2 X: ?并用这些数据点生成插值多项式公式,求[0,10]区间任意一点的值。) _6 E r' H8 p1 k
4 ~3 Z) Z. r; h2 `9 h' D
要求:/ z" c* w" g+ y1 r
编写常微分方程的四阶Runge-Kutta求解函数,和Matlab内建的ode45对比
# [- u b' H' P! K+ z编写Lagrange插值函数,要求支持任意多的输入点4 h+ Z& Q$ [/ F7 p% |
5 E, O& D+ Q0 h1 C! T
2) 计算 在区间[-5, 5]上的值,x的步长为1,对求出的数据点(x,y)用上述Lagrange函数生成插值多项式。在全区间上比较通过插值多项式和原函数计算得到的结果的差异,并设法改进 |
|