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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 【求助】用matlab最优化方法进行参数拟合
182/2返回列表上一页12
查看: 3953  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 数学EPI ，点击这里进行查看

xiaohsu2006

木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
7楼: Originally posted by dingd at 2011-12-05 10:57:31:
这种动力学微分方程拟合，1stOpt很简单好使啊！
[code]
Variable T, yA, yB, yC, yD, yE;
ODEFunction yA'=-(k1+k2+k3+k4)*yA;
            yB'=k1*yA-(k5+k6+k7)*yB;
            yC'=k2*yA+k5*yB-(k8+k9)* ...

t/min        y1        y2        y3        y4
2.5                3.8                2.5                0.5                0.5
5                5.8                3.3                1.0                0.6
7.5                5.8                2.7                1.2                0.7
10                5.4                2.3                1.2                0.8
15                5.7                2.4                1.2                1.3
30                8.5                2.6                2.3                2.4
模型方程
Y1’=x*k4*y5*y7+k6*y5+k10*y2*y7-0.1942*k10*y1*y4;
Y2’=x*k3*y5*y7+k5*y5-k10*y2*y7+0.1942*k10*y1*y4-k11*y2*y4+2.7624E-6*k11*y3*y7;
Y3’=k7*y5+k11*y2*y4-2.7624E-6*k11*y3*y7;
Y4’=(x-z+y/2)*k3*y5*y7+(2x-z+y/2)*k4*y5*y7+k8*y5+k10*y2*y7-0.1942*k10*y1*y4-3*k11*y2*y4+8.2872E-6*k11*y3*y7;
Y5’=k2*y7-k3*y5*y7-k4*y5*y7-k9*y5;
Y6’=-k2*y6;
Y7’=-(2x-z)*k4*y5*y7-k10*y2*y7+0.1942*k10*y1*y4+k11*y2*y4-2.7624E-6*k11*y3*y7;
其中，自变量为t，因变量为y1,y2,y3,y4,y5,y6,y7,参数为x,y,z,k2,k3,k4,k5,k6,k7,k8,k9,k10,k11
想回归出其中的模型参数，为什么我的在1stopt上运行不出结果啊，能帮试一下吗？非常感谢！下面是我的程序：
Parameter x,y,z,k2,k3,k4,k5,k6,k7,k8,k9,k10,k11;
SharedModel;
Variable t,y(1:7);
Function   Y1'=x*k4*y5*y7+k6*y5+k10*y2*y7-0.1942*k10*y1*y4;
Y2'=x*k3*y5*y7+k5*y5-k10*y2*y7+0.1942*k10*y1*y4-k11*y2*y4+2.7624E-6*k11*y3*y7;
Y3'=k7*y5+k11*y2*y4-2.7624E-6*k11*y3*y7;
Y4'=(x-z+y/2)*k3*y5*y7+(2x-z+y/2)*k4*y5*y7+k8*y5+k10*y2*y7-0.1942*k10*y1*y4-3*k11*y2*y4+8.2872E-6*k11*y3*y7;
Y5'=k2*y7-k3*y5*y7-k4*y5*y7-k9*y5;
Y6'=-k2*y6;
Y7'=-(2x-z)*k4*y5*y7-k10*y2*y7+0.1942*k10*y1*y4+k11*y2*y4-2.7624E-6*k11*y3*y7;
Data;
//t, y1, y2, y3, y4
2.5                3.8                2.5                0.5                0.5
5                5.8                3.3                1.0                0.6
7.5                5.8                2.7                1.2                0.7
10                5.4                2.3                1.2                0.8
15                5.7                2.4                1.2                1.3
30                8.5                2.6                2.3                2.4

[ Last edited by xiaohsu2006 on 2011-12-12 at 20:14 ]
赞一下(1人) 回复此楼
11楼2011-12-12 20:04:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
★
小木虫(金币+0.5):给个红包，谢谢回帖
下面是正确代码，待求参数数目比数据行数还多，解估计不是唯一的，你自己用1stOpt慢慢跑下吧：
CODE:
Parameter x,y,z,k2,k3,k4,k5,k6,k7,k8,k9,k10,k11;
Variable t,y(1:4);
ODEFunction   Y1'=x*k4*y5*y7+k6*y5+k10*y2*y7-0.1942*k10*y1*y4;
Y2'=x*k3*y5*y7+k5*y5-k10*y2*y7+0.1942*k10*y1*y4-k11*y2*y4+2.7624E-6*k11*y3*y7;
Y3'=k7*y5+k11*y2*y4-2.7624E-6*k11*y3*y7;
Y4'=(x-z+y/2)*k3*y5*y7+(2x-z+y/2)*k4*y5*y7+k8*y5+k10*y2*y7-0.1942*k10*y1*y4-3*k11*y2*y4+8.2872E-6*k11*y3*y7;
Y5'=k2*y7-k3*y5*y7-k4*y5*y7-k9*y5;
Y6'=-k2*y6;
Y7'=-(2x-z)*k4*y5*y7-k10*y2*y7+0.1942*k10*y1*y4+k11*y2*y4-2.7624E-6*k11*y3*y7;
Data;
//t, y1, y2, y3, y4
2.5                3.8                2.5                0.5                0.5
5                5.8                3.3                1.0                0.6
7.5                5.8                2.7                1.2                0.7
10                5.4                2.3                1.2                0.8
15                5.7                2.4                1.2                1.3
30                8.5                2.6                2.3                2.4

PJ版就别跑了，没这些功能。
赞一下(1人) 回复此楼
12楼2011-12-12 21:21:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaohsu2006

木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
12楼: Originally posted by dingd at 2011-12-12 21:21:52:
下面是正确代码，待求参数数目比数据行数还多，解估计不是唯一的，你自己用1stOpt慢慢跑下吧：
[code]
Parameter x,y,z,k2,k3,k4,k5,k6,k7,k8,k9,k10,k11;
Variable t,y(1:4);
ODEFunction   Y1'=x*k4*y5*y7+ ...

迭代数: 0
计算用时(时:分:秒:微秒): 21:14:09:533
均方差(RMSE): 1.23760521438246
残差平方和(SSE): 9.19
相关系数(R): 0
相关系数之平方(R^2): 0
决定系数(DC): -2.56893203883495
卡方系数(Chi-Square): 6.3
F统计(F-Statistic): 0

参数        最佳估算
----------        ---------
我的结果显示这样，R显示0，这是不没运行完，要那么久吗？太谢谢你了。
我用的破解版的，能算吗?

[ Last edited by xiaohsu2006 on 2011-12-12 at 21:49 ]
赞一下(1人) 回复此楼
13楼2011-12-12 21:47:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
★
小木虫(金币+0.5):给个红包，谢谢回帖
这是一组结果图：
赞一下(1人) 回复此楼
14楼2011-12-12 21:51:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaohsu2006

木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
内容已删除
回复此楼
15楼2011-12-12 22:03:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ffjtouchlife
送红花一朵
16楼2013-12-27 10:20:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

春春小样

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
2楼: Originally posted by leedobb at 2011-03-21 15:36:02
% 哦，你的题目是这个啊，
% 这只能依靠数值解，matlab肯定没有现成的程序吧。
%
% 你这样做吧。构造一个能量函数（其实就是最小二乘法啦）
%
% F =  (yA(t2)- YYA(t2,k1,k2,k3...,k9) ) ^ 2
% ..嗯，算 ...

为什么我复制后运行不出来啊，求大神帮忙
赞一下 回复此楼
17楼2014-07-01 20:51:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leedobb

金虫 (正式写手)
什么问题啊
回复此楼
有一天，我打了个瞌睡就到了这里，但我知道我掉入了时光的循环中，虽得以永生，但只有第一个循环有意义。
18楼2014-07-01 22:41:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wodaifei 的主题更新
182/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +3 素年祭语 2026-03-15 3/150 2026-03-16 09:24 by 哦哦123
[考研] 东南大学364求调剂 +4 JasonYuiui 2026-03-15 4/200 2026-03-16 08:36 by Linda Hu
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 304求调剂 +5 小熊joy 2026-03-14 5/250 2026-03-14 21:07 by peike
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 求调剂 +5 鹤遨予卿 2026-03-09 5/250 2026-03-14 02:44 by JourneyLucky
[考研] 求调剂，药学 +3 归零lbm 2026-03-09 5/250 2026-03-14 02:21 by JourneyLucky
[考研] 333求调剂 +3 球球古力 2026-03-09 3/150 2026-03-14 01:57 by JourneyLucky
[考研] 云南财经大学信息学院计算机学硕专硕学位点 +3 zjptai 2026-03-10 5/250 2026-03-14 01:23 by 飞行琦
[考研] 材料工程,326分，求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 0703，333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 【0856】化学工程（085602）313 分，本科学科评估A类院校化学工程与工艺，诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 0856化工原理 +6 z2839474511 2026-03-10 6/300 2026-03-13 10:41 by houyaoxu
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 一志愿江南大学085701环境工程专硕总分287求调剂 +5 18266118446 2026-03-09 5/250 2026-03-11 16:51 by 2020015
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
信息提示
关闭
请填处理意见
关闭