该用户从未签到
|
1) 设x为自变量,y为因变量,两者满足方程
8 x! Y( E# V! i' M* `, y(x=0) = 1
* c. R' [! U0 i0 n1 p用数值方法求在[0,10]区间,步长为1的x点所对应的y值,即x=0:10) R6 Q: G7 G* {0 B8 |$ Y" c
并用这些数据点生成插值多项式公式,求[0,10]区间任意一点的值。+ s1 M% i, O8 V7 v+ M4 d2 z, a
8 e+ R; w. \% N1 F要求:
% `9 X) M; s% \3 X# B2 G( Y编写常微分方程的四阶Runge-Kutta求解函数,和Matlab内建的ode45对比; X( {, l2 S9 e$ R
编写Lagrange插值函数,要求支持任意多的输入点1 D; l y+ q. [( u
0 F3 s$ k. }! ~/ X& }2 i8 s; w
2) 计算 在区间[-5, 5]上的值,x的步长为1,对求出的数据点(x,y)用上述Lagrange函数生成插值多项式。在全区间上比较通过插值多项式和原函数计算得到的结果的差异,并设法改进 |
|