微分方程数值解报告三.doc

微分方程数值解报告三 迎风格式和Lax-Friedrichs格式求解微分方程一 、引言双曲型方程差分格式的性质和定解问题解析解的性质之间有着密切的联系，由于其对初值的局部依赖关系和特征关系是其他两类方程所没有的，其初值函数的一些性质也会沿特征线传播，从而使解不具有光滑性质，在构造双曲型方程的差分格式时应充分考虑这些特性，下面江将用两种常见格式求解一个简单的偏微分方程。二、方法原理本题求解首先用到一阶迎风格式（Upwind Scheme）其中定义：另外，为了克服中心差分格式的不稳定性，还有Lax-Friedrichs格式：三、求解问题四、编程设计和方法分析考虑网格参数：接下来我们按照迎风格式和Lax-Friedrichs格式分别对方程进行离散处理迎风格式：Lax-Fredrichs格式为：迎风格式代码：%迎风格式clear allclose allclc%设置计算网格的参数L=1.0;%计算区域的长度T=0.5;%计算的时间长度dx=0.01;dt=0.005;x=-L:dx:L;t=0:dt:T;xn=size(x,2)-1;%网格段的个数tn=size(t,2)-1;%时间段的个数u=zeros(tn+1,xn+1);%设置初始条件for xi=1:xn+1if(x(xi)=0)u(1,xi)=1;%x0时,uo(x)=0endend%在计算开始前，将中间的网格点的值设为零for i=2:xnu(i,1)=1;%将中间的网格点的值设置为0endfor ti=1:tnfor j=2:xnu(ti+1,j)=dt*u(ti,j-1)/dx+(1.0-dt/dx)*u(ti,j);endend%绘制计算结果for i=2:30:92subplot(2,2,1+(i-2)/30);plot(x,u(i,:);str=num2str(t(i),T=%f);xlabel(str);str=num2str(i,u(%d,x);ylabel(str);endLax-Friedrichs格式代码2%Lax-Friedrichs格式clear allclose allclc%设置计算网格的参数L=1.0;%计算区域的长度T=0.8;%计算的时间长度dx=0.01;%网格尺寸dt=0.005;%时间步长miu=1;x=-L:dx:L;t=0:dt:T;xn=size(x,2)-1;%x的分段数tn=size(t,2)-1;%时间的分段数u=zeros(tn+1,xn+1);%u的上标为时间标号,下标为坐标标号%设置初始条件for xi=1:xn+1if(x(xi)=0)u(1,xi)=1;%x0时,uo(x)=0endendfor i=2:xnu(i,1)=1;endfor ti=1:tnfor j=2:xnu(ti+1,j)=0.5*(1-dt/dx)*u(ti,j+1)+0.5*(1+dt/dx)*u(ti,j-1);endend%绘制计算结果for i=2:30:92subplot(2,2,1+(i-2)/30);plot(x,u(i,:);str=nu