作者zupo (外派越南工作者)
标题Re: [讨论] 网格法求解问题
时间Tue May 31 22:00:42 2016
※ 引述《danny555361 ()》之铭言：
: 小弟的热传作业需要用网格法解每个node上的温度,我现在手上有28条方程式＆28个变数
: 方成组长的全部都想这样 △x △y 是已知 T(m,n)是我想要解的
: (△y/△x)*T(1,2)+(△x/△y)*T(2,1)-((△y/△x)+(△x/△y))*T(1,1)=0
: 我把每个网格的温度令成T(1,1) T(1,2).....T(6,5)
: 想说最後用solve来解28条方程式
: 最後想要解方程组的时候发现
: syms T(1,1) 无法编译
: 所以有想到办法是不是把T(m,n)都换成Tmn这样似乎就行了
: 可是这样我方程式全部都要重写
: 想问大家有没有其他好办法
: 先谢谢各位的建议
用差分法的话是可以解热传方程式
不过建议还是把原始算式放上来比较好帮忙
另外边界条件是? 没有这些也没有办法解
你这是 dT/dx + dT/dy = f(T,x,y) 的形式吧?
总之还是先把原始方程式放上来看看
--
数值方法 : 求根值 . ODE系统 . DDE系统. DAE系统. PDE系统
非线性联立系统. 化工程序控制系统 . 力学系统
RLC 系统 . MCK 系统 . Laplace 转换 系统
其他 : 工程数学 .微积分.化工程序.输送现象.化工程序线性控制
Simulator . MATLAB tool box . Simulink . FORTRAN
不会影像处理.....囧>
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 14.163.255.245
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/MATLAB/M.1464703247.A.3FC.html
1F：→ rex0707: 热传导是 dT/dt = k_x*(d^2T/dx^2) + k_y*(d^2T/d^2) 吧 06/01 00:50
2F：→ rex0707: 後面可能再多一项source 06/01 00:51
3F：→ zupo: 有一种是方向热传，解反矩阵即可求解 06/02 23:37
4F：→ rex0707: 恩...我只听过热传导、热对流、热辐射... 06/03 00:23
5F：→ rex0707: 解矩阵是隐式方法 不用解矩阵是显式方法 看原原PO用哪种 06/03 00:24