24小时热门版块排行榜    

查看: 1419  |  回复: 8
【悬赏金币】回答本帖问题,作者东林小王霸将赠送您 20 个金币

东林小王霸

新虫 (小有名气)

[求助] 一个两参数的微分方程数值拟合 已有2人参与

我是个matlab新手,刚csdn上看到一些一阶线性微分方程的数据拟合,感觉和我这个问题相近,但又有点不一样。大佬能给个思路吗

一个两参数的微分方程数值拟合


发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

东林小王霸

新虫 (小有名气)

2楼2024-05-09 12:03:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

将第一式代入第二式 得到 微分方程组。将各个sigma代入方程,实际上就是3组微分方程组 对应三组数据。目标函数可设置为三组微分方程残差平方和

发自小木虫Android客户端

» 本帖已获得的红花(最新10朵)

数值计算
3楼2024-05-09 14:12:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

东林小王霸

新虫 (小有名气)

送红花一朵
引用回帖:
3楼: Originally posted by 独孤神宇 at 2024-05-09 14:12:49
将第一式代入第二式 得到 微分方程组。将各个sigma代入方程,实际上就是3组微分方程组 对应三组数据。目标函数可设置为三组微分方程残差平方和

非常感谢,用ode计算能到H随时间的关系,带入一式可以得到不同时刻应变率,这样怎么应变求出来

发自小木虫Android客户端
4楼2024-05-10 12:11:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
4楼: Originally posted by 东林小王霸 at 2024-05-10 12:11:28
非常感谢,用ode计算能到H随时间的关系,带入一式可以得到不同时刻应变率,这样怎么应变求出来
...

组成微分方程组求解, 两个因变量不就直接求解出来了么

发自小木虫Android客户端

» 本帖已获得的红花(最新10朵)

数值计算
5楼2024-05-10 12:35:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

东林小王霸

新虫 (小有名气)

送红花一朵
引用回帖:
5楼: Originally posted by 独孤神宇 at 2024-05-10 12:35:54
组成微分方程组求解, 两个因变量不就直接求解出来了么
...

我联立求解,解得的值不变啊
一个两参数的微分方程数值拟合-2


一个两参数的微分方程数值拟合-3


一个两参数的微分方程数值拟合-4



发自小木虫Android客户端
6楼2024-05-10 14:19:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

实在不行 对应变率的结果积分

发自小木虫Android客户端
数值计算
7楼2024-05-10 15:30:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

johnliu1983

至尊木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
你是不是ode求解的时候就写了一个方程?是可以一起把这几个变量同时求出来, 然后一起求方差。如果时间点不是很符合就用插值计算一下就行。你可以到help 里面看看ode45这里面举的例子就行。然后把这些变成一个目标函数,再用最小二乘或者求极值的函数都可以计算收敛的参数。
8楼2024-05-11 16:38:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

东林小王霸

新虫 (小有名气)

引用回帖:
8楼: Originally posted by johnliu1983 at 2024-05-11 16:38:16
你是不是ode求解的时候就写了一个方程?是可以一起把这几个变量同时求出来, 然后一起求方差。如果时间点不是很符合就用插值计算一下就行。你可以到help 里面看看ode45这里面举的例子就行。然后把这些变成一个目标函 ...

谢谢您,我求解函数的问题解决了。我想问一下,我用差值计算出时间点的值和实验数据的差方作为目标函数,这个目标函数似乎和输入参数没有明确的表达式,具体怎么去优化呢

发自小木虫Android客户端
9楼2024-05-13 13:43:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 东林小王霸 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-01 15:08 by BCdvBAEpt3zm
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 15:04 by BCdvBAEpt3zm
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 14:55 by BCdvBAEpt3zm
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 14:50 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 14:50 by BCdvBAEpt3zm
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 14:29 by BCdvBAEpt3zm
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 14:24 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-01 12:53 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 12:19 by BCdvBAEpt3zm
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 12:11 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 10:27 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 07:01 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见