应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
51/1返回列表
查看: 10321  |  回复: 22
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 数学EPI ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

tmw1987

铁杆木虫 (正式写手)

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

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

» 猜你喜欢

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

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

查看全部散金贴
1楼 2010-11-30 06:44:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tmw1987

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

请大侠帮忙看看我的解题思路如何？
非常感谢！！！
赞一下(1人) 回复此楼
4楼2010-12-07 14:44:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答

saladin983

铁杆木虫 (正式写手)
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是什么，什么已知，什么未知？边值条件，初值条件？
赞一下(1人) 回复此楼
2楼2010-12-01 05:30:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saladin983

铁杆木虫 (正式写手)
tmw1987(金币+15):我把两个时间层组成矩阵求解，请见回帖，望指教！！！ 2010-12-07 08:33:48
风云箭(数学EPI+1):谢谢你详细的解答 2010-12-14 11:42:16
需要知道用于时间的离散格式，比如说，你用的显式格式，那么每一个时间点上的所有方程堆在一起构成一个线性方程组，规模不大的话直接用"\"运算求解，按时间逐层计算就可以了。不过你的方程写成这个格式，是不是意味着一个方程里同时有两个未知的时间层需要求解？这个似乎条件就不够用了。
赞一下(1人) 回复此楼
3楼2010-12-02 03:41:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saladin983

铁杆木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
cqsmath(金币+2):谢谢参与！ 2010-12-10 11:12:13
中心差分和欧拉用来离散基本上是最简洁明了的做法了，看起来没有什么问题。不过这样你一个方程里有三个时间层，在初始步算上已知的初值条件，还有两个未知的时间层。我对这样的问题不熟悉，但是感觉上还是少了个条件，或者子问题的规模要放大一倍。你可以仔细写清楚线性方程组，看看能不能归并处理掉一些。如果剩下的是可解的，自然就没问题了。当然，算得快慢依赖于算法，但是首先得把能不能解的问题弄清楚。
赞一下(3人) 回复此楼
5楼2010-12-08 03:56:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭