matlab圆柱内导热分离变量法,一维热传导方程数值解法及matlab实现分离变量法和有限差分法.doc...
一维热传导方程的Matlab解法分离变量法和有限差分法
问题描述
实验原理
分离变量法实验原理
有限差分法
实验目的
利用分离变量法和有限差分法解热传导方程问题 利用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;
for j=1:25
u(1,j)=0;
u(10,j)=0;
end;
for j=1:24
for i=2:9
u(i,j+1)=s*u(i+1,j)+(1-2*s)*u(i,j)+s*u(i-1,j);
end
end
disp(u);
[x,t]=meshgrid(1:25,1:10);
surf(x,t,u);
xlabel('t'),ylabel('x'),zlabel('T');
title(' 有限差分法解');
所得到的热传导图形为:
(2)
i分离变量法(取前100项和)
x=0:0.1*pi:pi;
y=0:0.04:1;
[x,t]=meshgrid(x,y);
s=0;
for i=1:100
s=s+(200*(1-(-1)^i))/(i*pi)*(sin(i*x).*exp(-i^2*t));
end;
surf(x,t,u);
xlabel('x'),ylabel('t'),zlabel('T');
title(' 分离变量法');
axis([0 pi 0 1 0 100]);
所得到的热传导图形为:
Ii有限差分法
根据(1)我们有如下图
结论:
比较可得这两幅图基本相同,有限差分法和分离变量法对本题都适应
(3)
第一种情况(取无穷项):
在原来程序代码的基础上加上 disp(s(:,6)); 可得出第六列(即x=pi/2)处温度随时间的变化情况
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]);
disp(s(:,6));
我们得到如下一组数据:
当温度低于50度是 时间为 t=23.5*0.04=0.94
第二种情况(取前100项和)
在原来程序代码的基础上加上 disp(s(:,6)); 可得出第六列(即x=pi/2)处温度随时间的变化情况
x=0:0.1*pi:pi;
y=0:0.04:1;
[x,t]=meshgrid(x,y);
r=0.04/(0.1*pi)^2;
fprintf('稳定性系数S为:')
disp(r);
s=0;
for i=1:100
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]);
disp(s(:,6));
当温度低于50度是 时间为 t=23.5*0.04=0.94
第三种情况(有限差分法)
在原来程序代码的基础上加上 disp(u(5,:));可得出第五行(即x=pi/2)处温度随时间的变化情况
u=
matlab圆柱内导热分离变量法,一维热传导方程数值解法及matlab实现分离变量法和有限差分法.doc...相关推荐
- matlab圆柱内导热分离变量法,一维热传导方程数值解法及matlab实现分离变量法和有限差分法...
一维热传导方程数值解法及matlab实现分离变量法和有限差分法 一维热传导方程的Matlab解法分离变量法和有限差分法问题描述实验原理分离变量法实验原理有限差分法实验目的利用分离变量法和有限差分法解热 ...
- 用matlab求解线性代数方程组,线性代数方程组数值解法与MATLAB实现综述
线性代数方程组数值解法及MATLAB 实现综述 廖淑芳 20122090 数计学院 12计算机科学与技术1班(职教本科) 一.分析课题 随着科学技术的发展,提出了大量复杂的数值计算问题,在建立电子计算 ...
- matlab定积分上界求解,定积分问题的数值求解及Matlab实现.pdf
定积分问题的数值求解及Matlab实现 第28卷第5期 哈 尔滨 商 业 大 学 学报 (自然科学版) Vo1.28No.5 2012年 10月 JournalofHarbinUniversityof ...
- 玩转matlab之一维 gauss 数值积分公式及matlab源代码
目录 标准区间 一般区间 数值实验 实验一 实验二 总结 下节预告 matlab代码 在数值分析中,尤其是有限元刚度矩阵.质量矩阵等的计算中,必然要求如下定积分: \[ I=\int_a^b f(x) ...
- matlab一维矩形积分,玩转matlab之一维 gauss 数值积分公式及matlab源代码
在数值分析中,尤其是有限元刚度矩阵.质量矩阵等的计算中,必然要求如下定积分: $$ I=\int_a^b f(x)dx $$学好gauss积分也是学好有限元的重要基础,学过高等数学的都知道,手动积分能 ...
- matlab的积分公式,玩转matlab之一维 gauss 数值积分公式及matlab源代码
释放双眼,带上耳机,听听看~! 目录 在数值分析中,尤其是有限元刚度矩阵.质量矩阵等的计算中,必然要求如下定积分: \\[ I=\\int_a^b f(x)dx \\]学好gauss积分也是学好有限元 ...
- 玩转 matlab 之一维 gauss 数值积分公式及matlab源代码
文章目录 标准区间 一般区间 数值实验 实验一 实验二 总结 下节预告 matlab代码 在数值分析中,尤其是 有限元刚度矩阵.质量矩阵等的计算中,必然要求如下定积分: I=∫abf(x)dxI=\i ...
- MATLAB模拟导热过程,一维热传导MATLAB模拟.doc
PAGE 昆 明 学 院 2015 届毕业设计(论文) 设计(论文)题目 一维热传导问题的数值解法及其MATLAB模拟 子课题题目 无 姓 名 伍有超 学 号 201117030225 所 属 系 物 ...
- matlab圆柱饶流,有限元法解二维圆柱绕流问题.pdf
目 录 目 录 1 1. 问题描述 1 2. 相关的有限元理论基础 1 2.1 二次泛函极值原理和里兹解法 1 2.2 伽辽金加权余数法 2 3. 二维圆柱绕流的有限元解法 ...
最新文章
- c语言链表创建递归,递归创建二叉树c语言实现+详细解释
- Quartz表达式触发
- 学python需要记得的单词_学习Python必背的初级单词有哪些?
- Mysql 存储引擎中InnoDB与Myisam的主要区别
- 面向初学者的JSF 2.0教程
- Linux 字符设备驱动开发基础(二)—— 编写简单 PWM 设备驱动
- LeetCode62 Unique Paths
- make_heap,pop_heap,push_heap
- FMS直播,在PC和手机PAD同时接收直播(使用RTMP和HLS协议)
- 【LeetCode】求众数(四种方法)
- 有什么能测试安卓硬件的软件吗,手机硬件检测工具有哪些 总有一款适合你
- 防火墙开放21端口linux,linux防火墙开放80,3306,21,443端口
- 代码从svn到工作空间,Myeclipse中java项目转成Web项目
- 马云登陆雅虎首次发表演讲稿全文。
- 怎么设置微信公众号自动回复蓝色字体,点击蓝色字体自动发送出去
- 封装尺寸与封装的对应关系
- python求绝对值_python中取绝对值简单方法总结
- AE基础教程(15)——第15章 质量图标和效果开关
- 2019-2-14 if条件判断语句+车费计算器
- linux种子文件制作工具,Centos 制作BT种子并获取BT种子信息