matlab求一维热传导方程数值解代码,一维热传导方程数值解法及matlab实现
【实例简介】
含matlab程序,个人感觉很有帮助,在研究传热学的可以下来看看
能呈守恒定律:
因为内部无热源,净流入的热量应该等于介质在此时
间内温度升高所需要的热量。
cdmdu=dQ=[q(x, t)g(x+ dx, t)]dt
g(x, t)dxdt
∵ comdt= cpdd∴ perdu=-q,dxdt
cpm=-9,即cPm2=-9(2)
q(x, t) q(x+dx, t)
xxIx
X
:-k
COL
由(1)、(2)得cp2=-q2=ku
,其中
cp
介质内存在热源时
如果在介质内有热量产生(例如,有化学反应发生,
或者通有电流,…),单位时间内单位体积介质产生
的热量为F(x,t)
因为热传导的Foue定律没有变化,所以仍然有
q=-k
对于能量守恒定律,有
du
(x.tddt+ Fdxdt即
qx +F
.u,=a'u+--F=a'u+f(x,t)
实验原理
分离变量法实验原理
有界长杆的热传导问题
、考察齐次热传导方程的混合问题(边界条
件都是第一类的情形)
u(00)
a(0,t)=0,a(l,t)=0
(x0)=g(x)
其中q(x)为给定的已知函数
下面用分离变量法(或称驻波法)来求解定解问
题(17)。
首先令
(x,t)=X(x)r()
将其代入方程
并分离变量得两个常微分方程
r(t)+a7()=0
X"(x)+AX(x)=0,
由边界条件(0,)=0.u(,)=0.可得
X(O)=0.X()=0
求边值问题
X"(x)+AX(x)=0.X(0)=X()=0
的非0解。
(1)当<0时,该问题没有非平凡解
(2)当=0时,该问题也没有非平凡解
(3)当入>0时,该问题有非平凡解。
此时
n=12…
x (x)=B, sin
(n=123…)
现在考虑T()+ia2()=0.
将特征值
=Z=
n、2
n=1.2.…)
代入上方程得
T()+()2T()=0
其通解为
n=C.
1.2.…
于是可得定解问题(17)中的一维热传导方程且
且满足齐次边界条件的具有变量分离形式的特解
n(x,1)=∑ae
sIn
(18)
其中an=BCn:是任意常数。
再利用初值条件(x,0)=o(x),可得
∑a,=∞x
an=[o∞(x)sin"t,
(19)
(18)
o(x)sin -dx.
(19)
(18)(19)合在一起就是所求定解问题(17)的特解。
1=a2na(00)
a(02t)=0.a(,)=0,
(x,0)=(x)
7)
若问题中的边界条件出现第二类或第三类齐次边界
条件,解法类似
有限差分法
、有限差分法的特点
有限差分方法(FD)是计算机数值模拟最早采用的方法,
至今仍被广泛运用。该方法将求解域划分为差分网格,
用有限个网格节点代替连续的求解域。有限差分法以
Taylor级数展开等方法,把控制方程中的导数用网格
节点上的函数值的差商代替进行离散,从而建立以网
格节点上的值为未知数的代数方程组。
该方法是一种直接将微分可题变为代数可题的近似数
值解法,数学概念直观,表达简单,是发展较早且比
较成熟的数值方法。
有限差分法的缺点是必雩进行整个区域的剖分,并且
要求网格比较规则,空间网格最好为直角网格。
、热传导方程(抛物方程)
1.热传导方程的介绍
0,t)=a(Z,)=0
x10)=f(x)
2.离散光u=u(O,大k)=0a3y=a(L)=0
=(ih20)=f(ih)=f
1)向前差分格式
L1-2n;+
+1
k
h2
计算
=sl1+(1-25)x1+sa1s
h2
这是一个星式格式(四点格式
F+1
i+1
t =f
可以证明:当0
式是稳定的。所以x的步长h和的步长k取法要恰
当
(2)向后差分格式
+1
2u:'+u
j+1
k
sul+(1+29)1-u1=a
J+1
1+1
实验目的
利用分离变量法和有限差分法解热传导方程问题利
用 matlab进行建模构建图形研究不同的情况下采用
何种方法从更深层次上理解热量分布与时间、空间分
布关系
【实例截图】
【核心代码】
matlab求一维热传导方程数值解代码,一维热传导方程数值解法及matlab实现相关推荐
- matlab圆柱内导热分离变量法,一维热传导方程数值解法及matlab实现分离变量法和有限差分法...
一维热传导方程数值解法及matlab实现分离变量法和有限差分法 一维热传导方程的Matlab解法分离变量法和有限差分法问题描述实验原理分离变量法实验原理有限差分法实验目的利用分离变量法和有限差分法解热 ...
- 用matlab求解线性代数方程组,线性代数方程组数值解法与MATLAB实现综述
线性代数方程组数值解法及MATLAB 实现综述 廖淑芳 20122090 数计学院 12计算机科学与技术1班(职教本科) 一.分析课题 随着科学技术的发展,提出了大量复杂的数值计算问题,在建立电子计算 ...
- Matlab求微分方程的数值解
注:首先计算微分方程的解析解,如果发现没有解析解,再用数值解 一.Matlab中求微分方程的数值解函数 [x,y]=solver('f',ts,x0,options) 1)x代表自变量 2)y代表函数 ...
- 一维谐振子定态 Schrödinger 方程的数值解法
本文链接:个人站 | 简书 | CSDN 版权声明:除特别声明外,本博客文章均采用 BY-NC-SA 许可协议.转载请注明出处. 前几天整理电脑的时候发现了本科上量子力学讨论班时做的一个 Slide, ...
- 用matlab求不动点迭代,科学网—数值分析--非线性方程组不动点迭代法matlab程序 - 殷春武的博文...
%%%程序编写者 西北工业大学自动化学院 Email: yincwxa2013@mail.nwpu.edu.cn %% All rights reserved clear clc x1=in ...
- matlab中牛顿下山法实例,非线性方程的数值解法牛顿下山法matlab
非线性方程的数值解法牛顿下山法matlab 1 非线性方程的数值解法 --计算物理实验作业九 陈万 物理学2013级 130******** ● 题目: 用下列方法求0133=--=x x f(x)在 ...
- matlab中牛顿下山法实例,非线性方程的数值解法牛顿下山法matlab.docx
非线性方程的数值解法牛顿下山法matlab.docx 1 非线性方程的数值解法 --计算物理实验作业九 陈万 物理学2013级 13020011006 题目: 用下列方法求 在 附近的根.根的准确 ...
- matlab圆柱内导热分离变量法,一维热传导方程数值解法及matlab实现分离变量法和有限差分法.doc...
一维热传导方程的Matlab解法分离变量法和有限差分法 问题描述 实验原理 分离变量法实验原理 有限差分法 实验目的 利用分离变量法和有限差分法解热传导方程问题 利用matlab进行建模构建图形 研究 ...
- matlab求系统根轨迹代码_广州数控GSK RH06六轴焊接机器人的轨迹规划
技术参数 运动范围 任务点位置 问题:如实现上图中的 运动轨迹(以基坐标系为参照),请给出相关的轨迹规划算法及其实现方式(结合运动学及动力学方程进行相关的轨迹点选取及工程实现),其中各个点的坐标为A= ...
- matlab求系统根轨迹代码_根轨迹法、PID参数整定和matlab指令计算
收获 (1)理解根轨迹的概念及其在控制系统设计中的作用: (2)手绘根轨迹草图,以及如何使用极端及绘制根轨迹: (3)熟悉在反馈控制系统中应用广泛的关键部件:PID控制器: (4)理解根轨迹在参数设计 ...
最新文章
- 你不知道的shell用法和脚本,简单强大
- 【PHP】php生成一个不重复的数字(订单号、会员号)
- (三)html5的结构
- 一些移动端开发的细节记录
- 数据结构题(莫队算法)
- AJAX Control Toolkit 控件学习(转自csdn山巅)
- Python笔记-类的单例模式及重载打印
- 去掉ecshop后台版权
- Hadoop 高可用集群(HA)
- javascript鼠标滚轮滚动实现模块的伸缩功能
- 3月11日Linux课程笔记
- php图片异步上传,怎么在php中使用JavaScript实现图片异步上传功能
- vscode 格式化不加分号_大前端时代你的 VSCode 插件
- 将多个Excel工作簿合并到一个Excel工作簿中
- linux中使用rpm命令详解
- 如何进行SPSS的重复测量方差分析
- 4.Nginx优化,谁用谁说好
- docker健康检查
- 怎么主动发起话题_聊天怎么开启话题?学会这四点,再也不担心尴尬开场
- “SqlSession[xxx] was not registered for synchronization because synchronization is not active”问题成功解决
热门文章
- python pytz下载_使用 pytz 处理 Python 中的时区问题
- java可达性_可达性分析详解
- 面向 5G 的新型多载波传输技术比较
- 2011年12月13日 timeout 与 refused windows clipbrd
- 无线网络信号测试软件WirelessMon 免费赠送
- 破解无线路由 backtrack BT5 Linux R3 安装 SpoonWep2.deb
- sketchup生成面域插件_什么插件这么神奇,SketchUp一秒搞定99%异形建模
- 必须了解的EXCEL常用函数都有哪些?【附动图详解】数据分析
- MATLAB实现香农编码
- 计算机组成原理——总线标准