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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求版主帮忙跑下面的1stopt的程序，参数估计。
61/1返回列表
查看: 561  |  回复: 6
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

hua99dushu

铁虫 (小有名气)

[求助] 求版主帮忙跑下面的1stopt的程序，参数估计。已有1人参与

[/code]Title "test 1";
Parameters a1[0,],a2[0,],a3[0,],b1[0,],b2[0,],b3[0,],c1[0,],c2[0,],c3[0,],M[1.5,],C[y,],E[z,],N[1.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：
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];

感谢小木虫论坛，感谢斑竹们的耐心帮忙。
回复此楼

» 猜你喜欢

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

hua99dushu

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2014-02-20 20:40:16
参考：

Parameters a1=,a2=,a3=,b1=,b2=,b3=,c1=,c2=,c3=,M=,C=,E=,N=;
Variable t,x,y,z;
ODEFunction  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;
Dat ...

有几个参数是这样设置的：
M[1.5,],C[1,],E[1,],N[1.5,]。
求帮着再跑一下啊。。
赞一下 回复此楼
3楼2014-02-21 14:54:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua99dushu

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by dingd at 2014-02-20 20:40:16
参考：

Parameters a1=,a2=,a3=,b1=,b2=,b3=,c1=,c2=,c3=,M=,C=,E=,N=;
Variable t,x,y,z;
ODEFunction  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;
Dat ...

之前程序里参数的设置我没有讲清楚。
是这样的：M[1.5,],C[1,],E[1,],N[1.5,].
帮忙再来一遍。。急着对比结果啊
赞一下 回复此楼
4楼2014-02-24 10:09:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua99dushu

铁虫 (小有名气)
引用回帖:
5楼: Originally posted by dingd at 2014-02-24 11:52:03
参考下：

均方差(RMSE):0.0743428646356062
残差平方和(SSE):0.182386430233526
相关系数(R): 0.990490604242565
相关系数之平方(R^2): 0.981071637092801
确定系数(DC): 0.97907392333002
F统计(F-Statist ...

真的很感谢你如此迅速的回复，多次认真的解答我的问题。
根据实际的条件，已经确定了一些参数。
[/CODE]
Parameters a1=[0,],a2=[0,],a3=[0,],b1=[0,],b2=[0,],b3=[0,],c1=[0,],c2=[0,],c3=[0,];
Variable t,x,y,z;
ODEFunction  x'=a1*x*(y/33.214-1)-a2*y+a3*z;
y'=-b1*x+b2*y*(1-y/22.854)+b3*z*(1-z/45.695);
z'=c1*x*(x/1.789-1)-c2*y-c3*z;
Data;
t=[0:1:12];
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];
求帮忙再跑一下上面的代码，还有就是是采用的什么算法啊？
赞一下 回复此楼
6楼2014-02-24 14:46:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua99dushu

铁虫 (小有名气)
有一行错了，正确的是：ODEFunction  x'=a1*x*(y/13.214-1)-a2*y+a3*z;
y/M，m=13.214.
赞一下 回复此楼
7楼2014-02-24 15:15:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hua99dushu 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭