sovamatlab測試代碼

Ⅰ Matlab 這段代碼怎麼成功運行

functionhtt

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))

holdon

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,dx)
c=1/2e-4;f=dx;s=0;

%icbun.m%偏微分方程初始條件(一維動態傳熱)
functionu0=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;