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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 【求助】用matlab最优化方法进行参数拟合
182/2返回列表上一页12
查看: 4409  |  回复: 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
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +3 启萌科技 2026-08-12 11/550 2026-08-13 22:57 by 启萌科技
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 28/1400 2026-08-13 17:51 by 医学老男孩
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +9 家与远方 2026-08-10 14/700 2026-08-13 09:41 by 臭臭不臭01
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 是这周出结果还是下周出结果？ +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
关闭
请填处理意见
关闭