有限差分法热传导matlab,热传导方程有限差分法的MATLAB实现

【实例简介】

热传导方程有限差分法的MATLAB实现

适用于解决热传导方面的偏微分问题

MATLAB

29

end

0.45

数值解

ezzeros(M-1, 1)

数值解

0.35

for ii=2. M

0.3

0.25

ez(ii-1)=abs(u(ii, 1)-u(i1, 2))/u(ii, 2)

0.2

end

0.15

E=max(ez)

0

25

fprintf

1):di

00

0.40.6

(u);

fprintf

图2数值解与准确解的比较

);

disp([num2str(E*100)%);

5

plot(x0,u(;1),r:,x0,u(:2).b-);

MATLAB

xlabel(x'), ylabel(u(x, t ),

MATLAB

title(

2005,19(4):76-80

0,1

0.001

),2007,20(3):12-14

n

90

3 StephenJ Chapman MATLAB [ MI

MATLAB

LI

,2006,20(4):

12-13

.5723%。

MATLAB

),2003,22(4):21-27

2008,24(2):3