该用户从未签到
|
1) 设x为自变量,y为因变量,两者满足方程1 @0 ^# H) c. |, L+ P" C
, y(x=0) = 1
4 Q2 r+ n9 g- R: N1 U4 K用数值方法求在[0,10]区间,步长为1的x点所对应的y值,即x=0:10! q& t1 m: X% h
并用这些数据点生成插值多项式公式,求[0,10]区间任意一点的值。
3 L! K4 C+ R. V3 p; A# X9 _0 j2 U - p' N B" Y W, w! p. {+ i
要求:
9 x, a1 Q( I, s m9 D# B. K; I5 l9 c; C编写常微分方程的四阶Runge-Kutta求解函数,和Matlab内建的ode45对比( s# ~7 M4 B1 W
编写Lagrange插值函数,要求支持任意多的输入点
' J: s6 T0 ^3 s! k9 s- W6 e% t 5 F/ x' b0 X2 n6 k
2) 计算 在区间[-5, 5]上的值,x的步长为1,对求出的数据点(x,y)用上述Lagrange函数生成插值多项式。在全区间上比较通过插值多项式和原函数计算得到的结果的差异,并设法改进 |
|