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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求帮忙运行一下下面的1stopt的程序，微分方程-参数估计
51/1返回列表
查看: 783  |  回复: 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

铁虫 (小有名气)
引用回帖:
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的回帖
查看全部 7 个回答

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的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭