该用户从未签到
|
1) 设x为自变量,y为因变量,两者满足方程
( V# S. v5 @6 z+ Z1 {! n/ h, y(x=0) = 1/ i; v# `( p! z) {* b
用数值方法求在[0,10]区间,步长为1的x点所对应的y值,即x=0:10
3 X+ B& k* n( z; o: c) s并用这些数据点生成插值多项式公式,求[0,10]区间任意一点的值。
. U: ]5 e9 }" S" [8 Z5 q " h- O" X/ u) d5 {8 E, [2 {
要求:
; h3 _/ j$ Y4 f& q( l3 t( U5 g6 P2 r编写常微分方程的四阶Runge-Kutta求解函数,和Matlab内建的ode45对比6 ]; W- `) l% N
编写Lagrange插值函数,要求支持任意多的输入点( o* h: ~8 C. W% {
! u5 ]5 r" C8 b
2) 计算 在区间[-5, 5]上的值,x的步长为1,对求出的数据点(x,y)用上述Lagrange函数生成插值多项式。在全区间上比较通过插值多项式和原函数计算得到的结果的差异,并设法改进 |
|