偏微分实验报告1

青岛农业大学微分方程数值解法实验报告姓名：学院：专业：班级：学号：完成时间：2012年11月4日微分方程数值解法实验课：实验一实验题目利用五点差分格式求解椭圆型方程Dirichlet边值问题实验目的掌握五点差分格式的建立及其求法实验地点及时间机房2452012年10月23日实验内容：利用五点差分格式求解椭圆型方程Dirichlet边值问题实验习题用五点差分格式求解椭圆型方程边值问题222222(1)sin(),02,01,(0,)sin(),(2,)sin(),01,(,0)0,(,1)0,02.xuueyxyxyuyyuyeyyuxuxx该问题的精确解为(,)sin().xuxyey要求：1、写出数值求解格式；hjijijijijijijiyxfuhuhuhhuhuh),(),,(11)11(211-1,22,121,2221,1211,22),(),,()(22jiyxfuujiijyijx),(),,(jiyxujiij2、给出所用程序代码；functionvarargout=hu(varargin)a=0;b=2;c=0;d=1;h1=1/10;h2=1/20;f=inline('(pi^2-1)*exp(x)*sin(pi*y)','x','y');g1x=inline('0');g2x=inline('0');g1y=inline('sin(pi*y)');g2y=inline('exp(2)*sin(pi*y)');[X,Y,Z]=chfenmethed(f,g1x,g2x,g1y,g2y,a,b,c,d,h1,h2);function[X,Y,Z]=chfenmethed(f,g1x,g2x,g1y,g2y,a,b,c,d,h1,h2)N=10000;x=a:h1:by=c:h2:dm=length(x);n=length(y);ee=;[Y,X]=meshgrid(y,x);Z=zeros(m,n);U=zeros(m,n);fori=2:m-1U(i,1)=feval(g1x,x(i));U(i,n)=feval(g2x,x(i));endforj=1:nU(1,j)=feval(g1y,y(j));U(m,j)=feval(g2y,y(j));endfori=2:m-1forj=2:n-1qu=f0(x(i),y(j));Z(i,j)=abs(qu-U(i,j));endendfunctionz=f0(x,y)z=exp(x)*sin(pi*y)3、给出实验结果。x=0y=0z=………实验总结：通过本次实验，我学会了运用五点差分法解决椭圆型方程Dirichlet边值问题，并且学会了运用Matlab进行程序的编写并且调试。实验成绩评阅时间评阅教师微分方程数值解法实验课：实验二实验题目利用Crank-Nicolson差分格式计算抛物型方程定解问题实验目的掌握Crank-Nicolson差分格式的建立及其求法实验地点及时间机房1272012年10月31日实验内容：利用Crank-Nicolson差分格式计算抛物型方程定解问题实验习题：用Crank-Nicolson差分格式求解抛物型方程定解问题：2210,01,01,(,0),01,(0,),(1,),01.xttuuxttxuxexuteutet该问题的精确解为：(,).xtuxte要求：1、写出数值求解格式；10,11),,(2)1(22)1(2-2112121121112nkmitxfuhuruhuhuruhkikikikikikiki2、给出所用程序代码；functionvarargout=liu(varargin)T=1;a=1;h=1/20;k=1/10;f=inline('0','x','t');fx=inline('exp(x)');ft1=inline('exp(t)');ft2=inline('exp(1+t)');[X,T,Z]=chfenmethed(f,fx,ft1,ft2,a,T,h,k);function[X,T,Z]=chfenmethed(f,fx,ft1,ft2,a,T,h,k)x=0:h:at=0:k:Tm=length(x);n=length(t);r=a*k/h^2;[T,X]=meshgrid(t,x);Z=zeros(m,n);U=zeros(m,n);fori=1:mU(i,1)=feval(fx,x(i));endforj=2:nU(1,j)=feval(ft1,t(j));U(m,j)=feval(ft2,t(j));endA=-r/2*ones(1,m-2);B=(1+r)*ones(1,m-2);C=A;UU=zeros(1,m-2);f1=UU;forj=2:ntk1=*(t(j-1)+t(j));fori=2:m-1UU(i-1)=f0(x(i),t(j));f1(i-1)=r/2*U(i-1,j-1)+(1-r)*U(i,j-1)+r/2*U(i+1,j-1)+k*feval(f,x(i),tk1);endf1(1)=f1(1)+r/2*U(1,j);f1(m-2)=f1(m-2)+r/2*U(m,j);U(2:m-1,j)=zgf(A,B,C,f1);Z(2:m-1,j)=abs(U(2:m-1,j)-UU');endfunctionx=zgf(A,B,C,f)n=length(B);B1=zeros(1,n-1);Y=zeros(1,n);x1=zeros(1,n);B1(1)=C(1)/B(1);fori=2:n-1B1(i)=C(i)/(B(i)-A(i)*B1(i-1));endY(1)=f(1)/B(1);fori=2:nY(i)=(f(i)-A(i)*Y(i-1))/(B(i)-A(i)*B1(i-1));endx1(n)=Y(n);fori=n-1:-1:1x1(i)=Y(i)-B1(i)*x1(i+1);endx=x1;functionz=f0(x,t)z=exp(x+t)3、给出实验结果。xtz00实验总结：本次实验我学会了利用Crank-Nicolson差分格式计算抛物型方程定解问题，掌握了用Matlab为其编写程序。实验成绩评阅时间评阅教师微分方程数值解法实验课：实验三实验题目利用显式差分格式计算二阶双曲型方程定解问题实验目的掌握双曲型方程显式差分格式的建立及其求法实验地点及时间机房2452012年11月2日实验内容：利用显式差分格式计算二阶双曲型方程定解问题实验习题：用显式差分格式求解双曲型方程定解问题222210,01,01,(,0),(,0),01,(0,),(1,),01.xxttuuxttxuuxexextuteutet该问题的精确解为：(,).xtuxte要求：1．写出数值求解格式；.11,11),,()1(2211222211nkmitxfuuhuhuhukikikikikiki2．给出所用程序代码；functionvarargout=wang(varargin)a=1;T=1;b=;h=1/5;k=1/20;f=inline('0','x','t');fx1=inline('exp(x)');fx2=inline('exp(x)');ft1=inline('exp(t)');ft2=inline('exp(1+t)');[X,Y,Z]=chfenmethed(f,fx1,fx2,ft1,ft2,a,T,h,k);function[X,T,Z]=chfenmethed(f,fx1,fx2,ft1,ft2,a,T,h,k)x=0:h:at=0:k:Tm=length(x);n=length(t);s=a*k/h;[X,T]=meshgrid(x,t);Z=zeros(n,m);U=zeros(n,m);fori=2:m-1U(1,i)=feval(fx1,x(i));U(2,i)=U(1,i)+k*feval(fx2,x(i))+k^2/2*(a^2/h^2*...(feval(fx1,x(i+1))-2*feval(fx1,x(i))+feval(fx1,x(i-1))+feval(f,x(i),0)));endforj=1:nU(j,1)=feval(ft1,t(j));U(j,m)=feval(ft2,t(j));endUU=f0(X,T);fori=2:n-1forj=2:m-1U(i+1,j)=s^2*U(i,j-1)+2*(1-s^2)*U(i,j)+s^2*U(i,j+1)-U(i-1,j)+k^2*feval(f,x(j),t(i));Z(i+1,j)=abs(U(i+1,j)-f0(x(j),t(i+1)));endendfunctionz=f0(x,t)z=exp(x+t)3．给出实验结果。x=0t=Columns1through180Columns19through21z=…实验总结：本次实验我学会了利用显式差分格式计算二阶双曲型方程定解问题，掌握双曲型方程显式差分格式的建立及其求法，学会了其对应的Matlab语言的编写。实验成绩评阅时间评阅教师