首页 > TAG信息列表 > 龙格
数值分析：研究高次插值的龙格现象
1 实验目的 研究高次插值的龙格现象。考虑函数在某区间范围内，构造拉格朗日插值多项式L(x)，分别画出不同n值下的拉格朗日插值函数。 2 实验内容 2.1 实验1.1 研究高次插值的龙格现象。考虑函数在[-1,1]上取n+1个等距节点。构造拉格朗日插值多项L(x)。 （1）分别画出n=2,4,6,8,10,14阶经典龙格库塔数值计算matlab程序
%用途：4阶经典龙格库塔格式解常微分方程组y’=f(x,y),y(x0)=y0 %格式：[x,y]=marunge4s(dyfun,xspan,y0,h) %dyfun为向量函数f(x,y),xspan为求解区间[x0,xn], %init为初值向量，N为步数，x返回节点，y返回数值解向量 function[x,y]=marunge4s(dyfun,xspan,init,h) % h=(xspan(1)-xspanMATLAB中ode23函数，龙格库塔函数
今天说一说MATLAB中ode23函数的原理，在网上看了好多，但是不知道是怎么计算的，就知道是那么用的，但是最后结果咋回事不知道，今天来讲一讲是怎么计算的。 首先来个程序： function f=eg6fun(t,y) f=-y^3-2; end 上面是我定义的一个函数，看着挺简单的哈！不多说了。 [t,y]=ode23(@eg龙格-库塔法求解连续系统输出
import control import numpy as np import matplotlib.pyplot as plt G1 = control.tf([1,1],[1,4]) G2 = control.tf([1,2],[1,1,1.25]) G = control.series(2.5*G1,G2) sys = control.tf2ss(G) h = 0.3 #步长 L = np.arange(0,15,h) # 仿真步数 u = np.ones(len(L)) # 输入