Question

Matlab 这段代码怎么成功运行

我打到matlab里面出错，求大神帮忙搞出来，代码如下： 函数：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 onplot(x,sol(3,:,1))plot(x,sol(5,:,1))plot(x,sol(7,:,1))plot(x,sol(9,:,1))plot(x,sol(16,:,1))结果如下：（看图片）搞的出结果的大神有加分啊啊 然后告诉我怎么搞的
（前两张图是结果 后两张图是整洁版问题描述）





Answer 1

将下面的代码保存为m文件，在命令行输入：htt

function  htt

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;