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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 【求助】用matlab最优化方法进行参数拟合
北京石油化工学院2026年研究生招生接收调剂公告
182/2返回列表上一页12
查看: 4007  |  回复: 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
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿郑州大学085600求调剂 +15 吃的不少 2026-04-05 18/900 2026-04-05 23:15 by 来看流星雨10
[考研] 一志愿河北工业大学材料工程，初试344求专硕调剂 +3 15933906766 2026-04-05 3/150 2026-04-05 22:17 by dongzh2009
[考研] 262求调剂 +7 天下第一文 2026-04-04 8/400 2026-04-05 21:31 by 激流勇渡
[考研] 348求调剂 +3 车厘子zzz 2026-04-05 3/150 2026-04-05 20:30 by 啵啵啵0119
[考研] 327求调剂 +4 拾光任染 2026-04-05 4/200 2026-04-05 20:16 by 南航~万老师
[考研] 计算机11408，286分求调剂 +7 木子念晞 2026-04-05 7/350 2026-04-05 19:02 by chy09050039
[考研] 304求调剂 +3 luoye0105 2026-04-05 3/150 2026-04-05 18:16 by 土木硕士招生
[考研] 一志愿 南京航空航天大学 ，080500材料科学与工程学硕 +10 @taotao 2026-03-30 10/500 2026-04-05 17:57 by jj987
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 346分的生物与医药08600求调剂 +5 常雨阳上岸 2026-04-05 6/300 2026-04-05 13:42 by imissbao
[考研] 358求调剂 +7 秋gk 2026-04-04 7/350 2026-04-05 13:29 by huangmoli
[考研] 367求调剂 +3 谢28 2026-03-30 3/150 2026-04-05 13:27 by huangmoli
[考研] 求调剂 +3 小沢 2026-04-03 3/150 2026-04-05 09:10 by sihailian3
[考研] 338求调剂 +4 zzz，，r 2026-04-03 4/200 2026-04-03 16:39 by lijunpoly
[考研] 316求调剂 +14 舟自梗 2026-04-01 18/900 2026-04-03 10:28 by linyelide
[考研] 279求调剂 +6 qazplm0852 2026-04-02 6/300 2026-04-03 10:03 by 蓝云思雨
[考研] 一志愿山东大学化学与化工学院材料与化工专硕，360分求调剂 +4 不愿透露姓名的� 2026-04-02 4/200 2026-04-03 09:29 by 遗忘消失的灆
[考研] 化学070300-总分378-求调剂 +5 挪椅子的泡泡糖 2026-04-02 5/250 2026-04-02 22:20 by ZXlzxl0425
[考研] 材料调剂 +12 一样YWY 2026-04-01 12/600 2026-04-02 09:15 by olim
[考研] 食品学硕362求调剂 +3 xuanxianxian 2026-04-01 3/150 2026-04-01 21:05 by 啊李999
信息提示
关闭
请填处理意见
关闭