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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stopt 语句修改
242/3返回列表上一页123下一页
查看: 1326  |  回复: 23
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xxz903

金虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
你这个可以直接用微分方程拟合啊！不过1.5版本的应该不能运行！代码如下：
Parameters k;
Variables a,y;
odeFunction y'=k*y;
RowData;
0, 1, 2, 3;
1, 3, 4, 5;
赞一下 回复此楼
11楼2014-12-08 09:15:10
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chengxl

木虫 (小有名气)
引用回帖:
11楼: Originally posted by xxz903 at 2014-12-08 09:15:10
你这个可以直接用微分方程拟合啊！不过1.5版本的应该不能运行！代码如下：
Parameters k;
Variables a,y;
odeFunction y'=k*y;
RowData;
0, 1, 2, 3;
1, 3, 4, 5;

您好。我的是高版本的，不能运行，请问a值如何体现？不用体现？
赞一下 回复此楼
12楼2014-12-08 10:50:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chengxl

木虫 (小有名气)
引用回帖:
11楼: Originally posted by xxz903 at 2014-12-08 09:15:10
你这个可以直接用微分方程拟合啊！不过1.5版本的应该不能运行！代码如下：
Parameters k;
Variables a,y;
odeFunction y'=k*y;
RowData;
0, 1, 2, 3;
1, 3, 4, 5;

请指导，我运行结果都是0.   a是初值，y是某一点（t=2）函数值，咋么写程序嫩狗狗体现出来这些。请高手指点
赞一下 回复此楼
13楼2014-12-08 10:56:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
chengxl: 金币+2 2014-12-09 08:33:43
引用回帖:
11楼: Originally posted by xxz903 at 2014-12-08 09:15:10
你这个可以直接用微分方程拟合啊！不过1.5版本的应该不能运行！代码如下：
Parameters k;
Variables a,y;
odeFunction y'=k*y;
RowData;
0, 1, 2, 3;
1, 3, 4, 5;

Parameters c, k;
Variable t, y;
Function  y=c*(exp(k*t)-1)+1;
data;
0  1
1  3
2  4
3  5
// ----------------------------------------
"Type your title here"

迭代数: 66
计算用时(时:分:秒:微秒): 00:00:00:94
优化算法: 粒子群算法PSO
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.0895640117970001
残差平方和(SSE): 0.0320868488366927
相关系数(R): 0.998194643831721
相关系数之平方(R^2): 0.996392546974336
决定系数(DC): 0.996332931561521
卡方系数(Chi-Square): 0.00435149314601671
F统计(F-Statistic): 552.407773509935

参数        最佳估算
----------        -------------
c                  -5.45183726787717
k                  -0.428004961121383
====== 结果输出 =====
No        实测值y        计算值y
1        1        1.0000000
2        3        2.8982851
3        4        4.1356029
4        5        4.9420968
赞一下 回复此楼

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

chengxl
14楼2014-12-08 15:08:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chengxl

木虫 (小有名气)
送红花一朵
引用回帖:
14楼: Originally posted by feixiaolin at 2014-12-08 15:08:18
Parameters c, k;
Variable t, y;
Function  y=c*(exp(k*t)-1)+1;
data;
0  1
1  3
2  4
3  5
// ----------------------------------------
"Type your title here"

迭代数: 66
计算用时( ...

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;
赞一下 回复此楼
15楼2014-12-08 19:40:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
15楼: Originally posted by chengxl at 2014-12-08 19:40:52
y'=k*y
已知不同初始条件下（t=0）初始值及其对应的t=10时函数值，如何用1stopt 语句拟合k值。下列语句运行不了，请高手指导，感激不尽。假定函数不能积分得到解析解。
Parameters k;
Variables a,y;
Function ...

我的版本也很低，跑不了。
赞一下 回复此楼
16楼2014-12-08 20:35:10
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chengxl

木虫 (小有名气)
引用回帖:
16楼: Originally posted by feixiaolin at 2014-12-08 20:35:10
我的版本也很低，跑不了。...

谢谢了，麻烦您了
回复此楼
17楼2014-12-09 08:33:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雕刻你的拾光

新虫 (初入文坛)
SharedModel;
Parameters X0, Xm, um,a,b,P0,Yx,Yp,m,S0;
Variables t, X, P,S;
Function
           X= X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))
           P=P0+a*(X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))-X0)+(b*Xm/um)*ln((Xm-X0+X0*exp(um*t))/Xm)
           S=S0-(1/Yx+a/Yp)*(X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))-X0)-Xm/um*(b/Yp+m)*ln((Xm-X0+X0*exp(um*t))/Xm);
Data;
t    X          P        S
0   0.05        0        79
4   0.16        0        78
8   0.6515      0        76
12  1.5426      1.5      74
16  2.4047      4.9      69
20  3.1066      8.7      60
24  3.6693      12.4     50
28  3.8165      16.1     43
32  3.8874      17.2     40
36  3.9618      17.9     38
三个函数表达式分开可以算出来，可是合在一起就说表达式有误，分开算出的数据是错的，应该怎么弄，求大神指点
赞一下 回复此楼
18楼2014-12-16 22:30:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

CODE:
Parameters k;
Variable a,y;
Function y=int(k*y,y=0,10)+a;
RowData;
0,1,2,3;
1,3,4,5;

用微分方程的形式拟合应该也可以实现。
1stopt 语句修改
附图1.png
赞一下 回复此楼

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

chengxl
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
19楼2014-12-17 09:14:43
 已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 chengxl 的主题更新
242/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭