8.1 欧接法.ppt

步长与误差 离散误差 设y=y(t)是初值问题的唯一解, 是离散近似解集 →从tk到tk+1步的计算误差 注 : ,每步忽略 那么在末端进行M步之后,误差积累到 欧拉方法的精度 设y=y(t)是初值问题 的解,若 是欧拉方法计算的近似值序列,则 注: 区间末端的误差称为最终全局误差. 实现 通过计算 求[a,b]上的初值问题 的近似解. Matlab Code f=inline(‘x’); a=0;b=10; Ya=1;M=20; h=(b-a)/M; T=zeros(1,M+1); Y=zeros(1,M+1); T=a:h:b; Y(1)=ya; for k=1:M Y(k+1)=Y(k)+h*f(T(k),Y(k)); end E=[T’ Y’]; Plot(E) 欧拉公式及其改进 tn tn+1 欧拉公式 欧拉公式的改进办法 后退的欧拉公式 差商代替导数即得欧拉公式 若对于点tn+1列出方程, 若以向后差商 替代导数 代入离散化得 ---------后退欧拉公式 隐式与显式欧拉公式 注: 隐式前后均有yn+1,无法直接求解yn+1. 迭代法,逐步显示化. 欧拉公式得出yn+1的初值,隐式欧拉公式逐步迭代. 隐式收敛性 只hL<1,则迭代法收敛 误差分析 欧拉公式的误差为: 而后退欧拉公式的误差为: 如果我们将这两种方法进行算术平均,即可消除误差的主要部分 从而获得更高的精度 梯形公式 几何意义 tn tn+1 由于后退Euler公式是隐形公式,计算例1将很麻烦 事实上大多数情况下用后退Euler公式都较困难 就可得到新的Euler公式 此方法称为预测—校正系统 Matlab Code f=inline(‘x’); a=0;b=10; Ya=1;M=20; h=(b-a)/M; T=zeros(l,M+l); Y=zeros(l,M+l); T=a:h:b; Y(l)=ya; for j=l:M kl=feval(f,T(j),Y(j)); k2=feval(f,T(j+l),Y(j)+h*kl); Y(j+l)=Y(j)+(h/2)*(kl+k2); end H=[T'Y']; Plot(H) Largran Code X=[]; %X is a vector that contains a list of abscissas Y=[]; % Y is a vector that contains a list of ordinates w=length(X); n=w-1; L=zeros(w,w); %Form the Lagrange coefficient polynomials for k=1:n+1 V=1; for j=1:n+1 if k~=j V=conv(V,poly(X(j)))/(X(k)-X(j)); end end L(k,: )=V; end %Determine the coefficients of the Lagrange interpolating polynomial C=Y*L; * * 24 * 24 * 24 * 24 * 24 第一节 拉格朗日插值 一般插值问题 拉格朗日线性插值 拉格朗日插值公式 算法设计 一、微分方程导论 是一个微分方程，因为它包含“未知函数”y=y( t )的 导数dy/dt . 由于只有独立变量t出现在上式的右端，因此 的 不定积分就是方程的一个解。可由积分公式求解 y( t ) 其中C为积分常数。任意C，y(t)均满足上式，因此微 分方程的解为一曲线族. 微分方程的曲线族 积分方法可用于求解式中函数的显式公式，在这样的解中有一个自由度C. 通过改变C的值可以向上或向下“移动解曲线”，从而可以找到过任意需要的曲线。 然而科学与工程问题中极少表现为显式公式，通常只能考查一个变量的变化如何影响另一个变量，得到的就是包含未知函数的变化率及自变量和应变量的方程。即 dy/dt=f (x,y) 曲线簇 以冷却物体的温