偏微分方程组解法.docVIP

偏微分方程组解法 某厚度为10cm平壁原温度为20，现其两侧面分别维持在20和120，试求经过8秒后平壁内温度分布，并分析温度分布随时间的变化直至温度分布稳定为止。 式中为导温系数，；。 解： 模型转化为标准形式： 初始条件为： 边界条件为： ， 函数： 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; 命令： x=linspace(0,10,20)*1e-2; t=linspace(0,15,16); sol=pdepe(0,@pdefun,@icfun,@bcfun,x,t); mesh(x,t,sol(:,:,1)) %温度与时间和空间位置的关系图 %画1、2、4、6、8、15s时刻温度分布图 plot(x,sol(2,:,1)) 1s时刻，（因为本题sol第一行为0时刻） hold on plot(x,sol(3,:,1))