该用户从未签到
|
1) 设x为自变量,y为因变量,两者满足方程 _# N7 q! B) y
, y(x=0) = 12 t- b `8 [% k1 e
用数值方法求在[0,10]区间,步长为1的x点所对应的y值,即x=0:10
8 e* s% @! |. e( J, ~, v4 a6 t并用这些数据点生成插值多项式公式,求[0,10]区间任意一点的值。
$ @/ t! L+ ` @: R% p ' @. M! f5 F, @% t$ B
要求:
: F+ ]" Q, l* X1 a5 X编写常微分方程的四阶Runge-Kutta求解函数,和Matlab内建的ode45对比) y, u! p) c) \' m0 a! r1 ~
编写Lagrange插值函数,要求支持任意多的输入点) @5 L; q# o6 ^8 G0 N, t* W( B# {2 \# c: c
6 c4 G o% K- Z1 ?8 d
2) 计算 在区间[-5, 5]上的值,x的步长为1,对求出的数据点(x,y)用上述Lagrange函数生成插值多项式。在全区间上比较通过插值多项式和原函数计算得到的结果的差异,并设法改进 |
|