小弟的热传作业需要用网格法解每个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这样似乎就行了
可是这样我方程式全部都要重写
想问大家有没有其他好办法
先谢谢各位的建议

写成矩阵的形式 x=A\b你是用有限差分吗? 怎么觉得离散式写的怪怪的

没错 是有限差分 可能是第一个点两边绝热 看起来怪怪感谢大大的建议想问有办法可以用方程式直接求解吗因为每个变量都可能会出现在好几条方程式中

我懂你的意思 正常的作法就是写成矩阵的形式再解线性系统

用矩阵的话我天生眼残 蛮可能漏看的

可以用Gauss-Seidel方法求解不过28条式子写成矩阵 直接法求解就很好了

ok 感谢 原来一般就是用换成矩阵来解直接法是x=A\b吗

我不知道在matlab里x=A\b是用什么方法解