24小时热门版块排行榜    

查看: 9230  |  回复: 22
本帖产生 1 个 数学EPI ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

tmw1987

铁杆木虫 (正式写手)


[交流] 【求助】MATLAB 有限差分法(FDM)求解偏微分方程

大侠:
  
    我最近计算某模型扩散模型,离散化结果如下所示,请高手指点如何求解!!!
   
   
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫(金币+0.5):给个红包,谢谢回帖交流
步骤写得很清楚,赞一个。不过还是这个问题,系数矩阵[ML1, ML2]是L*2L的,这个线性方程组怎么求解?这个系数矩阵连方阵都不是,怎么能用求逆的运算呢?如果强行计算的话,得到的应该只是一个特解,没什么意义吧。

另外你Matlab求解的这个函数是内建的?不认识,能简单说说么?
7楼2010-12-10 05:30:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答
tmw1987(金币+5):n是时间,j是坐标,a1-a3,b1-b3,c1-c3是已知,求解T,初值条件为T(x,t)=0,边值条件T(0,t)=0,T(L,t)=q。请大侠指正!!! 2010-12-01 06:57:47
n是什么,j是什么,什么已知,什么未知?边值条件,初值条件?
2楼2010-12-01 05:30:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
tmw1987(金币+15):我把两个时间层组成矩阵求解,请见回帖,望指教!!! 2010-12-07 08:33:48
风云箭(数学EPI+1):谢谢你详细的解答 2010-12-14 11:42:16
需要知道用于时间的离散格式,比如说,你用的显式格式,那么每一个时间点上的所有方程堆在一起构成一个线性方程组,规模不大的话直接用"\"运算求解,按时间逐层计算就可以了。不过你的方程写成这个格式,是不是意味着一个方程里同时有两个未知的时间层需要求解?这个似乎条件就不够用了。
3楼2010-12-02 03:41:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tmw1987

铁杆木虫 (正式写手)


★ ★
cqsmath(金币+2):谢谢参与! 2010-12-10 11:11:55
引用回帖:
Originally posted by saladin983 at 2010-12-02 03:41:58:
需要知道用于时间的离散格式,比如说,你用的显式格式,那么每一个时间点上的所有方程堆在一起构成一个线性方程组,规模不大的话直接用"\"运算求解,按时间逐层计算就可以了。不过你的方程写成这个格式 ...

请大侠帮忙看看我的解题思路如何?
非常感谢!!!
4楼2010-12-07 14:44:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 3/150 2026-02-08 03:05 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 3/150 2026-02-08 02:52 by vs90ilomwc
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 4/200 2026-02-08 02:45 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 3/150 2026-02-08 02:32 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 5/250 2026-02-08 02:27 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 5/250 2026-02-08 02:25 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 5/250 2026-02-08 02:12 by vs90ilomwc
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 02:07 by vs90ilomwc
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 02:05 by vs90ilomwc
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 6/300 2026-02-08 01:46 by vs90ilomwc
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 01:26 by vs90ilomwc
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 01:12 by vs90ilomwc
[有机交流] 酰胺脱乙酰基 10+5 chibby 2026-02-03 12/600 2026-02-07 19:29 by 江东闲人
[基金申请] 同年申请2项不同项目,第1个项目里不写第2个项目的信息,可以吗 +4 hitsdu 2026-02-06 4/200 2026-02-07 13:07 by jurkat.1640
[基金申请] 有时候真觉得大城市人没有县城人甚至个体户幸福 +9 苏东坡二世 2026-02-04 10/500 2026-02-07 12:37 by 小毛球
[考博] 天津大学招2026.09的博士生,欢迎大家推荐交流(博导是本人) +4 a793625982 2026-02-05 5/250 2026-02-07 10:57 by a793625982
[公派出国] CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生|MSCA经费充足、排名优 +4 雨念 2026-02-01 6/300 2026-02-06 23:32 by MelissaPon
[基金申请] 面上项目申报 +3 Tide man 2026-02-01 3/150 2026-02-05 22:56 by god_tian
[硕博家园] 博士延得我,科研能力直往上蹿 +7 偏振片 2026-02-02 7/350 2026-02-04 17:36 by 陈氏帝国
[教师之家] 遇见不省心的家人很难过 +18 otani 2026-02-03 22/1100 2026-02-04 11:06 by tangmnt
信息提示
请填处理意见