[Asm] 纯文本查看 复制代码
function ht
m = 0;
x=linspace(0,10,20)*1e-2;
t=linspace(0,15,16);
sol = pdepe(m,@pdefun,@icbun,@bcfun,x,t);
figure
mesh(x,t,sol(:,:,1)) %温度与时间和空间位置的关系图
figure
%画1、2、4、6、8、15s时刻温度分布图
%1s时刻，（因为本题sol第一行为0时刻）
plot(x,sol(2,:,1))
hold on
plot(x,sol(3,:,1))
plot(x,sol(5,:,1))
plot(x,sol(7,:,1))
plot(x,sol(9,:,1))
plot(x,sol(16,:,1))
% pdefun.m% 方程
function [c,f,s]=pdefun(x,t,u,dudx)
c=1/2e-4;f=dudx;s=0;
% icbun.m% 初始条件
function u0=icbun(x)
u0=20;
% bcfun.m% 边界条件
function [pl,ql,pr,qr] = bcfun(xl,ul,xr,ur,t)
pl=ul-120;ql=0;pr=ur-20;qr=0;
[复制链接]
发表于 2021-5-16 08:37
