24小时热门版块排行榜    

查看: 874  |  回复: 4

东景辉

新虫 (正式写手)

[求助] 1stopt微分方程组参数拟合(各位大哥求帮忙,没有高版本的软件) 已有1人参与

初始金币 5
Parameter  a11=0.01[0,10000];
Parameter  a12=0.01[0,10000];
Parameter  a13=0.01[0,10000];
Parameter  b1=0.01[0,10000];
Parameter  a21=0.01[0,10000];
Parameter  a22=0.01[0,10000];
Parameter  a23=0.01[0,10000];
Parameter  b2=0.01[0,10000];
Parameter  a31=0.01[0,10000];
Parameter  a32=0.01[0,10000];
Parameter  a33=0.01[0,10000];
Parameter  b3=0.01[0,10000];
Variable t,x,y;
ODEFunction x'=a11* x+a12*y+a13*z-b1*x*x;
         y'= a21* x+a22*y+a23*z-b2*y*y;
z'= a31* x+a32*y+a33*z-b2*y*y;
Data;
1     0  0      0
33    15 12     6
43    39 30     30
72   84   60     63
92    174  102    119
102   334  194    188
109   414  258    284
123  608  452    404 @beefly
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
内容已删除
2楼2017-05-11 13:48:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

东景辉

新虫 (正式写手)

非常感谢!!参数a11为0有点问题啊。 @beefly

发自小木虫Android客户端
3楼2017-05-11 15:16:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

东景辉

新虫 (正式写手)

能否再算一下Parameter??a11=0.01[0,10000];
Parameter??a12=0.01[0,10000];
Parameter??b1=0.01[0,10000];
Parameter??a21=0.01[0,10000];
Parameter??a22=0.01[0,10000];
Parameter??b2=0.01[0,10000];
Variable t,x,y;
ODEFunction x'=a11* x+a12*y-b1*x*x);
? ?? ?? ?y'= a21* x+a22*y-b2*y*y);
Data;
1? ???0? ?? ?0
33? ? 12? ???6
43? ? 30? ???30
72? ? 60? ???63
92? ? 102? ? 119
102? ?194? ? 188
109? ?258? ? 284
123? ?452? ? 404 @beefly

发自小木虫Android客户端
4楼2017-05-11 15:18:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

东景辉

新虫 (正式写手)

引用回帖:
2楼: Originally posted by dingd at 2017-05-11 13:48:49
参考下:
均方差(RMSE): 22.4162459568179
残差平方和(SSE): 7034.83315915174
相关系数(R): 0.99124454515698
相关系数之平方(R^2): 0.982565748303469
修正R平方(Adj. R^2): 0.965161343617583
确定系数(DC ...

非常感谢楼主,不过参数a11理论上不能是0,能否再重新算一下,谢谢
Parameter a11=0.01[0,10000];
Parameter a12=0.01[0,10000];
Parameter b1=0.01[0,10000];
Parameter a21=0.01[0,10000];
Parameter a22=0.01[0,10000];
Parameter b2=0.01[0,10000];
Variable t,x,y;
ODEFunction x'=a11* x+a12*y-b1*x*x);
y'= a21* x+a22*y-b2*y*y);
Data;
1 0 0
33 12 6
43 30 30
72 60 63
92 102  119
102 194 188
109 258  284
123 452 404 @beefly
5楼2017-05-11 15:32:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 东景辉 的主题更新
信息提示
请填处理意见