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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求帮忙运行一下下面的1stopt的程序，微分方程-参数估计
71/1返回列表
查看: 779  |  回复: 6
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

hua99dushu

铁虫 (小有名气)

[求助] 求帮忙运行一下下面的1stopt的程序，微分方程-参数估计已有1人参与

Parameters a1=0.001[0,1],a2=0.001[0,1],a3=0.001[0,1],b1=0.001[0,1],b2=0.001[0,1],b3=0.001[0,1],c1=0.001[0,1],c2=0.001[0,1],c3=0.001[0,1],M=[0,1],C=[0,1],E=[0,1],N=[0,1];
Sharedmodel;
Variable t, x,y,z;
Function x'=a1*x*(y/M-1)-a2*y+a3*z;
y'=-b1*x+b2*y*(1-y/C)+b3*z*(1-z/E);
z'=c1*x*(x/N-1)-c2*y-c3*z;
RowData：
x=[1.0000,1.098,1.226,1.426,1.709,2.102,2.442,2.906,3.438,3.796,4.532,5.189];
y=[1.0000,1.017,1.081,1.219,1.470,1.838,2.064,2.247,2.397,2.449,2.742,3.209];
z=[1.0000,1.106,1.240,1.455,1.754,2.174,2.542,3.042,3.622,4.028,4.843,5.741];
最好能给出误差图和相图呢。。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2014-03-22 11:36:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua99dushu

铁虫 (小有名气)
★
fegg7502: 金币+1, 3ks 2014-03-24 08:32:54
fegg7502: 鼓励交流 2014-03-24 08:33:13
a1=0.001[0,1]这种格式应该是表示初值是0.001，范围是0~1之间哟
我也是从教程里看到的，如果不对，请帮忙改正呢。
赞一下 回复此楼
2楼2014-03-22 11:38:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua99dushu

铁虫 (小有名气)
额，忘了t了
t=[0,1,2,3,4,5,6,7,8,9,10,11];
回复此楼
3楼2014-03-22 11:39:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua99dushu

铁虫 (小有名气)
★
fegg7502: 金币+1, 鼓励交流 2014-03-24 08:33:10
上面的错误挺多的，程序更正确定如下：
Parameters a1=0.001[0,1],a2=0.001[0,1],a3=0.001[0,1],b1=0.001[0,1],b2=0.001[0,1],b3=0.001[0,1],
c1=0.001[0,1],c2=0.001[0,1],c3=0.001[0,1],M=[0,1],C=[0,1],E=[0,1],N=[0,1];
Sharedmodel;
Variable t,x,y,z;
Function x'=a1*x*(y/M-1)-a2*y+a3*z;
y'=-b1*x+b2*y*(1-y/C)+b3*z*(1-z/E);
z'=c1*x*(x/N-1)-c2*y-c3*z;
RowData：
t=[0,1,2,3,4,5,6,7,8,9,10,11];
x=[1.0000,1.098,1.226,1.426,1.709,2.102,2.442,2.906,3.438,3.796,4.532,5.189];
y=[1.0000,1.017,1.081,1.219,1.470,1.838,2.064,2.247,2.397,2.449,2.742,3.209];
z=[1.0000,1.106,1.240,1.455,1.754,2.174,2.542,3.042,3.622,4.028,4.843,5.741];
希望能够得到估计的参数。。最好能有相关的图
赞一下 回复此楼
4楼2014-03-22 14:24:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua99dushu

铁虫 (小有名气)
请版主看此楼。。。。。快点帮忙运行一下吧。。
Parameters a1=0.001[0,1],a2=0.001[0,1],a3=0.001[0,1],b1=0.001[0,1],b2=0.001[0,1],b3=0.001[0,1],
c1=0.001[0,1],c2=0.001[0,1],c3=0.001[0,1],M=[0,5],C=[0,10],E=[0,10],N=[0,5];
Sharedmodel;
Variable t,x,y,z;
Function x'=a1*x*(y/M-1)-a2*y+a3*z;
y'=-b1*x+b2*y*(1-y/C)+b3*z*(1-z/E);
z'=c1*x*(x/N-1)-c2*y-c3*z;
RowData：
t=[0,1,2,3,4,5,6,7,8,9,10,11];
x=[1.0000,1.0777,1.1443,1.3497,1.5881,1.9814,2.1905,2.5799,2.9763,3.2445,3.5765,3.9315];
y=[1.0000,1.0597,1.0982,1.2035,1.3566,1.4610,1.4771,1.4862,1.6052,1.6094,1.8892,2.0334];
z=[1.0000,1.0920,1.2033,1.4031,1.6920,1.9382,2.2159,2.6186,2.9489,3.1536,3.5979,4.0233];
赞一下 回复此楼
5楼2014-03-23 20:18:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
fegg7502: 金币+2, 应助指数+1, 3ks 2014-03-24 08:33:25
只大概算了下，供参考：
均方差(RMSE):0.227821606236476
残差平方和(SSE):1.71278858084954
相关系数(R): 0.968347715467601
相关系数之平方(R^2): 0.937697298051321
确定系数(DC): 0.810472079088207
F统计(F-Statistic): -6.59531364009876

参数                  最佳估算
--------------------        -------------
a1        3.45198672383651E-11
a2        5.83733964128734E-8
a3        0.138472304076658
b1        1.44029340345734E-6
b2        0.0518147236689287
b3        0.0406931861755383
c1        3.70541984860404E-6
c2        1.37133289578244E-7
c3        8.59791920071737E-10
m        3.16205888404511
c        2.08529659127593
e        3.70210463181138
n        6.20898311572885E-5

====== 结果输出 ======

文件: 数据文件-1
No        t        目标 x        计算 x        目标 y        计算 y        目标 z        计算 z
1        1        1.0777        1.14300503145253        1.0597        1.05729662499709        1.092        1.06851434515071
2        2        1.1443        1.29685305954089        1.0982        1.11588828762999        1.2033        1.1572932896264
3        3        1.3497        1.46465905303486        1.2035        1.17578665477138        1.4031        1.2709880836245
4        4        1.5881        1.65027065718425        1.3566        1.23690130994621        1.692        1.41560520068449
5        5        1.9814        1.85849139764935        1.461        1.29894832240681        1.9382        1.59905260863551
6        6        2.1905        2.09539536339523        1.4771        1.36129294495917        2.2159        1.8319321926032
7        7        2.5799        2.36877609710511        1.4862        1.42267433005819        2.6186        2.12871429124705
8        8        2.9763        2.68879681981474        1.6052        1.48071370354723        2.9489        2.50951876009061
9        9        3.2445        3.06894995110712        1.6094        1.53101337771447        3.1536        3.00288297139041
10        10        3.5765        3.52750422345532        1.8892        1.56545537275185        3.5979        3.65018117220115
11        11        3.9315        4.08974278597234        2.0334        1.56886778417041        4.0233        4.51289544111373
赞一下 回复此楼
6楼2014-03-23 21:52:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua99dushu

铁虫 (小有名气)
引用回帖:
6楼: Originally posted by dingd at 2014-03-23 21:52:53
只大概算了下，供参考：
均方差(RMSE):0.227821606236476
残差平方和(SSE):1.71278858084954
相关系数(R): 0.968347715467601
相关系数之平方(R^2): 0.937697298051321
确定系数(DC): 0.810472079088207
F统计 ...

请问版主，能否把n调整为和之前三个参数一个数量级啊？
赞一下 回复此楼
7楼2014-03-24 20:49:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hua99dushu 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭