24小时热门版块排行榜    

查看: 1676  |  回复: 10

苜蓿情缘

新虫 (小有名气)

[求助] 哪位大神指导一下,如何用matlab求解偏微分方程组 已有2人参与

图片中t是时间变量,z是空间变量。初始条件为Rliq(0,z)=-0.01,Rice(0,z)=-0.01。phi=0.4,gamma=0.68,aeq=1.0031。我用了向量方法,就是把R=(Rliq,Rice)',然后变成了向量型的一阶一维双曲方程,最后向前差分离散化,然后处理。可结果出了问题,哪位大神可以知道一下,求助

哪位大神指导一下,如何用matlab求解偏微分方程组


发自小木虫Android客户端

[ Last edited by 苜蓿情缘 on 2016-6-8 at 00:26 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

★ ★
感谢参与,应助指数 +1
苜蓿情缘: 金币+2 2016-06-08 19:45:29
此方程组属线性常系数偏微分方程组,但现在仅有初始条件,尚缺边界条件。因此必须补上以后才能求解。并且如果边界条件不复杂的话甚至可用拉氏变换法求出解析解。
2楼2016-06-08 06:19:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woshiqiujun

新虫 (正式写手)

尝试先用拉普拉斯变换处理 看看边界怎么样 然后用逆变换导出边界条件 没有非线性项应该可以做出来 然后初边值问题写好算法最后代码实现 祝楼主万事如意 给个金币呗

发自小木虫Android客户端
9楼2016-06-08 16:31:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

苜蓿情缘

新虫 (小有名气)

3楼2016-06-08 08:10:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

苜蓿情缘

新虫 (小有名气)

这是我的做法,只利用了初始条件!以上偏微分方程组只有初始条件,也就是因为缺少边界条件,所以求起来很麻烦!
哪位大神指导一下,如何用matlab求解偏微分方程组-2


哪位大神指导一下,如何用matlab求解偏微分方程组-3



发自小木虫Android客户端
4楼2016-06-08 08:12:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

苜蓿情缘

新虫 (小有名气)

引用回帖:
2楼: Originally posted by peterflyer at 2016-06-08 06:19:01
此方程组属线性常系数偏微分方程组,但现在仅有初始条件,尚缺边界条件。因此必须补上以后才能求解。并且如果边界条件不复杂的话甚至可用拉氏变换法求出解析解。

这个方程组只有初始条件,所以解析解求不出,数值解也求起来比较麻烦!我已经把我的做法发到下面去了,请您给看一看,哪里有错误

发自小木虫Android客户端
5楼2016-06-08 08:14:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

苜蓿情缘

新虫 (小有名气)

这是我编写的matlab程序,请大神指导下
哪位大神指导一下,如何用matlab求解偏微分方程组-4



发自小木虫Android客户端
6楼2016-06-08 08:51:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woshiqiujun

新虫 (正式写手)

稍微大一点区域上用零边界条件吧

发自小木虫Android客户端
7楼2016-06-08 16:25:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woshiqiujun

新虫 (正式写手)

引用回帖:
7楼: Originally posted by woshiqiujun at 2016-06-08 16:25:35
稍微大一点区域上用零边界条件吧

不行 初值还不是不具有紧支集。。。

发自小木虫Android客户端
8楼2016-06-08 16:27:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

铸剑人

金虫 (正式写手)

【答案】应助回帖

★ ★ ★
苜蓿情缘: 金币+3 2016-06-11 15:06:10
使用matlab 偏微分方程求解工具箱先搞一搞,画个图看看。然后在研究基础代码。
绝尘。
10楼2016-06-11 11:42:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 苜蓿情缘 的主题更新
信息提示
请填处理意见