matlab求一维热传导方程数值解代码,一维热传导方程数值解法及matlab实现.DOC

一维热传导方程的Matlab解法分离变量法和有限差分法2010/12/20

问题描述

实验原理

分离变量法实验原理

有限差分法

实验目的

利用分离变量法和有限差分法解热传导方程问题 利用matlab进行建模构建图形 研究不同的情况下采用何种方法 从更深层次上理解热量分布与时间、空间分布关系。

模拟与仿真作业

分离变量法(代码)：

x=0:0.1*pi:pi;

y=0:0.04:1;

[x,t]=meshgrid(x,y);

s=0;

m=length(j);%matlab可计算的最大数 相当于无穷

for i=1:m

s=s+(200*(1-(-1)^i))/(i*pi)*(sin(i*x).*exp(-i^2*t));

end;

surf(x,t,s);

xlabel('x'),ylabel('t'),zlabel('T');

title(' 分离变量法(无穷)');

axis([0 pi 0 1 0 100]);

所得到的三维热传导图形为：

有限差分法：

u=zeros(10,25); %t=1 x=pi 构造一个1025列的矩阵(初始化为0)用于存放时间t和变量x 横坐标为x 纵坐标为t

s=(1/25)/(pi/10)^2;

fprintf('稳定性系数S为:\n');

disp(s);

for i=2:9

u(i,1)=100;

end;