简单的PDE求解偏微分方程

简单的PDE求解偏微分方程

clear;
T=20;
L=10;
D=0.8;
w=pi/40;
d=（80/pi）^0.5;
M=200;
N=200;
t=linspace（0TM）;
x=linspace（0LN）;
t0=0.5+0.5*exp（-x/d）.*sin（-x/d）;
tin=0.5+0.5*sin（w*t）;
options=odeset;
sol=pdepe（0@pdefun@tic@tbcxtoptions）;
[xt] = meshgrid（xt）;
k=0.5+0.5*exp（-x/d）.*sin（w*t-x/d）;
u=sol（::1）; 
surf（xtu）
hold on;
mesh（xtk）
figure;
imagesc（sol‘）;
colorbar;
%imagesc（k‘）;
t=linspace（035）; 
sol=pdepe（0@pdefun@tic@tbcxtoptions）;     
u=sol（::1）;
M=length（t）;  xout=linspace（010200）; 
[uoutdudx]=pdeval（mxu（M:）xout）; 
plot（xoutuout）; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        579  2018-06-04 23:28  homework529\homework529.m

     文件       1271  2018-06-04 23:32  homework529\homework5292.m

     文件         65  2018-06-04 22:49  homework529\pdefun.m

     文件        105  2018-06-04 22:41  homework529\tbc.m

     文件         79  2018-06-04 22:28  homework529\tic.m

     目录          0  2018-06-04 23:33  homework529

----------- ---------  ---------- -----  ----

                 2099                    6