微分方程数值解迎风格式.docVIP

宋玉霞 座机电话号码16 统数 信科0802 实验四：迎风格式 一、实验题目： 其中a1 1,b1 2,。 二、实验要求：用迎风格式求解双曲型对流扩散方程 迎风格式为 若令则迎风格式可整理为：稳定条件： 三、程序运行结果： e 7.9402e-004 四、误差图如下图所示： 五、程序附录： 1. function z ft x %求下边界 z exp x/2 ; 2.function z fx1 t %求左边界 z exp -t ; 3.function z fx2 t %求右边界 z exp 1/2-t ; 4.function z f x,t %求右端函数 z -exp x/2-t ; 5 .function z fu x,t %求真解 z exp x/2-t ; 6. function [X,T,z] upwindL a1,b1,a,b,h,k %用迎风格式求解 x 0:h:a;t 0:k:b; [T,X] meshgrid t,x ; m length x ;n length t ; r1 a1*k/h;r2 b1*k/h^2; uu zeros m,n ;%储存数值解 z uu;%储存误差 for i 1:m%求下边界 uu i,1 ft x i ; end for j 2:n%求左右边界 uu 1,j fx1 t j ; uu m,j fx2 t j ; end for j 2:n%迎风格式求内点，从下往上 for i 2:m-1%从左往右 uu i,j 1-r1-2*r2 *uu i,j-1 + r1+r2 *uu i-1,j-1 +r2*uu i+1,j-1 +k*f x i ,t j-1 ;%求数值解 z i,j abs uu i,j -fu x i ,t j ;%求误差 end end 7. a1 1;b1 2;a 1;b 1;h 0.1;k 0.001;%主函数，用于输出 [X,T,z] upwindL a1,b1,a,b,h,k ; mesh T,X,z e max max z title 误差图 xlabel x轴 ylabel t轴 zlabel z轴