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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 关于1STOPT拟合的小问题
151/2返回列表12下一页
查看: 986  |  回复: 14
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

heisekaishui

银虫 (小有名气)

[求助] 关于1STOPT拟合的小问题已有2人参与

各位前辈，我最近用1stopt拟合了一组数据，但是发现拟合的数据使用不同阶的龙格库塔方程结果不同，但是精度都差不多，不知道是什么原因，希望大家不吝赐教。
代码如下：
Parameters k1,k2,k3;
Variable t,a,b,c;
ODEFunction a'=-k1*a-k2*a;
b'=k1*a;
c'=k2*a;
Data;
t=0,2,4,6,8,10;
a=1.161575096,1.098850041,1.022186084,0.952491579,0.885120223,0.825879893;
b=0,0.060153328,0.132977117,0.197793007,0.254614938,0.306825415;
c=0,0.002571727,0.006411895,0.01129051,0.021839935,0.028869787;

同时，我在设定参数的范围时，并没有起到效果，而且1stopt的结果纵坐标范围是0-4而不是0-10不知道是什么原因
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
Don'tPanic!
1楼2015-07-17 19:41:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ybkooo

至尊木虫 (著名写手)

^^

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
heisekaishui: 金币+10, ★★★很有帮助, 非常感谢，我还想问一下如果我要设置参数的取值范围应该怎么设置，比如我想让所有的 2015-07-18 13:32:29
heisekaishui: 金币+10, ★★★很有帮助, 非常感谢，我还想问一下如果我要设置参数的取值范围应该怎么设置，比如我想让所有的 2015-07-18 13:32:48
看看这样能不能解决你的问题
Parameters k1,k2,k3;
Variable t,a,b,c;
SharedModel;
ODEFunction a'=-k1*a-k2*a;
b'=k1*a;
c'=k2*a;
RowData;
0,2,4,6,8,10;
1.161575096,1.098850041,1.022186084,0.952491579,0.885120223,0.825879893;
0,0.060153328,0.132977117,0.197793007,0.254614938,0.306825415;
0,0.002571727,0.006411895,0.01129051,0.021839935,0.028869787;
赞一下 回复此楼
abcd
2楼2015-07-17 21:15:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ybkooo

至尊木虫 (著名写手)

^^
Parameters k1,k2,k3,k4;
修改方程，这个可以得到稳定结果
Variable t,a,b,c;
ODEFunction a'=-k1*a-k2*a;
b'=k3*a;
c'=k4*a;
RowData;
0,2,4,6,8,10;
1.161575096,1.098850041,1.022186084,0.952491579,0.885120223,0.825879893;
0,0.060153328,0.132977117,0.197793007,0.254614938,0.306825415;
0,0.002571727,0.006411895,0.01129051,0.021839935,0.028869787;
赞一下 回复此楼
abcd
3楼2015-07-17 21:41:22
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ybkooo

至尊木虫 (著名写手)

^^
坐标显示0-4可能是他的作图功能不好，或不是正版软件，只看他的计算值就行了（计算值不是有，每时刻的计算值吗）
2        1.098850041        .060153328        .002571727        1.08633992640023        6.97307002743706E-02        5.60446897778746E-03        For a:
4        1.022186084        .132977117        .006411895        1.01597773545576        .134944946131689        1.08459367441951E-02        DC: 0.99487234139559
6        .952491579        .197793007        .01129051        .950172900633618        .195935268051574        1.57479148657442E-02        R: 0.999959508427279
8        .885120223        .254614938        .021839935        .888630242171105        .25297524940825        2.03323920723293E-02        RMSE: 0.00692414160103415
10        .825879893        .306825415        .028869787        .831073699086233        .306320753592402        2.46199328847424E-02        SSE: 0.000239718684555859
赞一下 回复此楼
abcd
4楼2015-07-17 21:48:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ybkooo

至尊木虫 (著名写手)

^^
11
关于1STOPT拟合的小问题
1111.png
赞一下 回复此楼
abcd
5楼2015-07-17 21:51:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
稳定唯一：

均方差(RMSE):0.00523400858904149
残差平方和(SSE):0.000410922688652401
相关系数(R): 0.99229618847388
相关系数之平方(R^2): 0.98465172565979
确定系数(DC): 0.949184328781669
F统计(F-Statistic): 585.146344863956

参数                  最佳估算
--------------------        -------------
k1        0.0310189175554483
k2        0.0024794201096689
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

heisekaishui
6楼2015-07-17 23:23:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

heisekaishui

银虫 (小有名气)
引用回帖:
2楼: Originally posted by ybkooo at 2015-07-17 21:15:36
看看这样能不能解决你的问题
Parameters k1,k2,k3;
Variable t,a,b,c;
SharedModel;
ODEFunction a'=-k1*a-k2*a;
b'=k1*a;
c'=k2*a;
RowData;
0,2,4,6,8,10;
1.161575096,1.098850041,1.022186084,0.95249 ...

我还想问一下如果我要设置参数的取值范围应该怎么设置，比如我想让所有的k大于0，应该怎样设置。
赞一下 回复此楼
Don'tPanic!
7楼2015-07-18 13:33:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

heisekaishui

银虫 (小有名气)
送红花一朵
引用回帖:
6楼: Originally posted by dingd at 2015-07-17 23:23:35
稳定唯一：

均方差(RMSE):0.00523400858904149
残差平方和(SSE):0.000410922688652401
相关系数(R): 0.99229618847388
相关系数之平方(R^2): 0.98465172565979
确定系数(DC): 0.949184328781669
F统计(F-St ...

非常感谢
回复此楼
Don'tPanic!
8楼2015-07-18 13:37:22
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

heisekaishui

银虫 (小有名气)
引用回帖:
3楼: Originally posted by ybkooo at 2015-07-17 21:41:22
Parameters k1,k2,k3,k4;
修改方程，这个可以得到稳定结果
Variable t,a,b,c;
ODEFunction a'=-k1*a-k2*a;
b'=k3*a;
c'=k4*a;
RowData;
0,2,4,6,8,10;
1.161575096,1.098850041,1.022186084,0.952491579,0. ...

我发现这样也是不稳定，采用龙格库塔一阶拟合的和版主的就一致，但是阶数越高差别越大，可能跟版本有关吧。
赞一下 回复此楼
Don'tPanic!
9楼2015-07-18 13:39:07
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

heisekaishui

银虫 (小有名气)
引用回帖:
6楼: Originally posted by dingd at 2015-07-17 23:23:35
稳定唯一：

均方差(RMSE):0.00523400858904149
残差平方和(SSE):0.000410922688652401
相关系数(R): 0.99229618847388
相关系数之平方(R^2): 0.98465172565979
确定系数(DC): 0.949184328781669
F统计(F-St ...

我采用龙格库塔一阶拟合的和版主的就一致，但是阶数越高差别越大，四阶结果就是k1=1.3，K2=0.1了，不知道是不是跟版本有关。
赞一下 回复此楼
Don'tPanic!
10楼2015-07-18 13:41:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 heisekaishui 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭