查看: 1749  |  回复: 23
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

chengxl

木虫 (小有名气)

[求助] 1stopt 语句修改 已有4人参与

y'=k*y
已知不同初始条件下(t=0)初始值及其对应的t=10时函数值,如何用1stopt 语句拟合k值。下列语句运行不了,请高手指导,感激不尽。假定函数不能积分得到解析解。
Parameters k;
Variables a,y;
Function y=int(k*y,t=0,10)+a;
RowData;
0 1 2 3;
1 3 4 5;
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

6.0运行结果(稳定唯一):

均方差(RMSE):0.483374374706425
残差平方和(SSE):7.00952358368483
相关系数(R): 0.981724799896839
相关系数之平方(R^2): 0.963783582732488
确定系数(DC): 0.937297250322943
F统计(F-Statistic): 0

参数                  最佳估算
--------------------        -------------
x0        0.0334076309586811
xm        4.64024768139639
um        0.241057508591346
a        4.1046342230345
b        -0.00415265747668047
p0        -0.287687677048546
yx        0.127392818019361
yp        5.42816890517109
m        0.0260506659407243
s0        78.5148748503121

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

文件: 数据文件-1
No        目标 x        计算 x        目标 p        计算 p        目标 s        计算 s
1        0.05        0.0334076309586811        0        -0.287687677048546        79        78.5148748503121
2        0.16        0.0866085836080747        0        -0.0702457278129587        78        78.0513786216106
3        0.6515        0.220476302701316        0        0.476847080403129        76        76.8848015882738
4        1.5426        0.536863715228011        1.5        1.76956429951415        74        74.1258474300948
5        2.4047        1.18549279397836        4.9        4.41819545471953        69        68.4600563772828
6        3.1066        2.19800926890448        8.7        8.54648023793028        60        59.5776135163638
7        3.6693        3.25942090148948        12.4        12.8576041989121        50        50.165648205155
8        3.8165        3.99496251671989        16.1        15.8159219548877        43        43.4653609261927
9        3.8874        4.3710537088627        17.2        17.2897532001951        40        39.8032194191773
10        3.9618        4.53378908930007        17.9        17.8835669783185        38        37.9511990644952
20楼2014-12-17 09:58:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
chengxl: 金币+5 2014-12-06 09:33:20
直接进行微分方程拟合不行吗?
2楼2014-12-05 22:43:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
如果用的是1.5版本的1stOpt软件,是不支持 int 函数的。
这种拟合用1stOpt软件的话,得是高版本的。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
3楼2014-12-06 08:58:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chengxl

木虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2014-12-05 22:43:26
直接进行微分方程拟合不行吗?

我在拟合方面知道的少,请多指导。如何直接进行微分方程拟合,请指导
4楼2014-12-06 09:02:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见