24小时热门版块排行榜    

查看: 2933  |  回复: 7

凌风7388

铁虫 (正式写手)

[交流] matlab偏微分方程的参数估计

求大神告知如何用matlab对偏微分方程进行参数估计?或者有什么好的例子或书籍推荐下也行。
我现在看的黄花江那本《实用化工计算机模拟》,关于偏微分参数估计这一块讲的不是特别详细,希望能有大神给指导下,谢谢!!
附件是要求未知的参数D,希望懂的大神可以帮忙编下代码。
关于x=0时的那个微分边界条件,如果采用MOL法,该如何如何处理?

matlab偏微分方程的参数估计
12.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
凌风7388: 金币+5 2015-09-10 13:53:19
Mol法就是把一个方向离散掉,比如时间t 方向离散掉,比如t 方向离散成100个点,这样原偏微分方程转换为一个含有100个方程的常微分方程组,然后用经典的龙格-库塔法求解该常微分方程组。

关于偏微分方程参数拟合,个人觉得很难实现,即便转化为常微分方程组,也是问题重重。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2015-09-10 09:39:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaochong8693

木虫之王 (文坛精英)

★ ★
凌风7388: 金币+2 2015-09-10 13:55:48
为楼主留下满满的、真诚的祝福
祝福楼主和各位虫友事事如意,梦想成真
祝福与祈祷中...
3楼2015-09-10 09:55:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

凌风7388

铁虫 (正式写手)

引用回帖:
2楼: Originally posted by 月只蓝 at 2015-09-10 09:39:56
Mol法就是把一个方向离散掉,比如时间t 方向离散掉,比如t 方向离散成100个点,这样原偏微分方程转换为一个含有100个方程的常微分方程组,然后用经典的龙格-库塔法求解该常微分方程组。

关于偏微分方程参数拟合, ...

嗯,我也是从书上看到的,先对其离散化转化为ODE再参数拟合,可是边界条件不知道如何处理了。书上也没例子,所以希望看看其他什么书上有讲的。
4楼2015-09-10 13:55:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wyj210

禁虫 (小有名气)

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
凌风7388: 金币+5 2015-09-11 12:06:40
本帖内容被屏蔽

5楼2015-09-11 00:03:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿Q~~

至尊木虫 (文坛精英)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
凌风7388: 金币+2 2015-09-11 12:06:59
路过的看了一下,帮顶,祝福你好运!~~~~~~~~~~~~~~~~~~
自强不息,厚德载物;独立精神,自由思想。
6楼2015-09-11 06:50:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

凌风7388

铁虫 (正式写手)

引用回帖:
5楼: Originally posted by wyj210 at 2015-09-11 00:03:15
看这本新的吧《化工数值计算与MATLAB》...

好的,我看看。谢谢了
7楼2015-09-11 12:06:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr__Right

专家顾问 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
终于有时间了,尝试回答下这个问题。

这个问题直观看跟MOL没有关系。
这个D如果是待定的常数,则偏微分方程式线性的,可以得到精确的解析解。
解析解可以通过分离变量得到。


有3+1=4个待定参数C1,C2,C3和D,初值和边值条件共3个
得到三个等式(约束条件)

一般来说是一个“欠定”问题,可能存在一个以上的解(非线性情况还可能无解)。

得到三个等式之后求带非线性方程组的解(或最小二乘解),随机优化方法即可。

所以,跟MOL没有什么关系。
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
8楼2016-02-11 20:40:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 凌风7388 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 9zw1gfyi7x 2026-07-31 10/500 2026-08-01 20:23 by Kddn1pB6pAzX
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 10/500 2026-08-01 20:22 by Kddn1pB6pAzX
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 20:10 by Kddn1pB6pAzX
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 jzg3xk1po8 2026-07-31 8/400 2026-08-01 20:03 by Kddn1pB6pAzX
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 jzg3xk1po8 2026-07-31 11/550 2026-08-01 20:03 by Kddn1pB6pAzX
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 jzg3xk1po8 2026-07-31 13/650 2026-08-01 20:02 by Kddn1pB6pAzX
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 9zw1gfyi7x 2026-08-01 3/150 2026-08-01 18:57 by Kddn1pB6pAzX
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 9/450 2026-08-01 18:43 by Kddn1pB6pAzX
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 9zw1gfyi7x 2026-08-01 7/350 2026-08-01 18:43 by Kddn1pB6pAzX
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 9zw1gfyi7x 2026-08-01 7/350 2026-08-01 18:42 by Kddn1pB6pAzX
[考博] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-08-01 8/400 2026-08-01 18:31 by Kddn1pB6pAzX
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 8/400 2026-08-01 18:22 by Kddn1pB6pAzX
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-07-31 10/500 2026-08-01 18:20 by Kddn1pB6pAzX
[博后之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:02 by Kddn1pB6pAzX
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jzg3xk1po8 2026-07-31 7/350 2026-08-01 17:42 by Kddn1pB6pAzX
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了? +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
信息提示
请填处理意见