24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1461  |  回复: 23

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朵)

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朵)

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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 307求调剂 +3 所念及所望 2026-04-06 3/150 2026-04-06 17:30 by 土木硕士招生
[考研] 305求调剂 +3 77Qi 2026-04-06 3/150 2026-04-06 15:36 by going home
[考研] 材料专硕322 +11 哈哈哈吼吼吼哈 2026-04-05 11/550 2026-04-06 14:07 by lqwchd
[考研] 285求调剂 +8 AZMK 2026-04-04 11/550 2026-04-06 13:56 by BruceLiu320
[考研] 302分 085601求调剂推荐 +11 zyx上岸! 2026-04-05 11/550 2026-04-05 22:13 by dongzh2009
[考研] 材料工程302分求调剂 +9 zyx上岸! 2026-04-04 9/450 2026-04-05 22:08 by 醉翁wl
[考研] 315求调剂 +5 &123456789 2026-04-05 5/250 2026-04-05 19:55 by nepu_uu
[考研] 328分调剂 +6 门men 2026-04-04 6/300 2026-04-05 13:40 by imissbao
[考研] 求调剂 +3 小沢 2026-04-03 3/150 2026-04-05 09:10 by sihailian3
[考研] 考研调剂 +6 15615482637 2026-04-04 6/300 2026-04-04 22:43 by yu221
[考研] 求调剂 +3 ffyyu 2026-04-02 3/150 2026-04-04 19:03 by 蓝云思雨
[考研] 292求调剂 +11 2022080213 2026-04-04 13/650 2026-04-04 18:38 by macy2011
[考研] 一志愿东北大学085901土木专硕345求调剂 +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by 土木硕士招生
[考研] 400分求调剂 +3 尴尬且挠头 2026-04-04 3/150 2026-04-04 08:41 by jp9609
[考研] 一志愿武汉理工0856,初试334 +3 26考研材料 2026-04-02 3/150 2026-04-02 21:22 by dongzh2009
[考研] 348求调剂 +11 zzzzyk123 2026-04-01 11/550 2026-04-02 16:52 by Wang200018
[考研] 初试301,代码085701环境工程,本硕一致,四六级已过,有二区一作,共发表5篇论文 +6 axibli 2026-04-01 6/300 2026-04-02 13:42 by Ecowxq666!
[考研] 一志愿北交大材料工程总分358 +3 cs0106 2026-04-02 5/250 2026-04-02 11:37 by olim
[考研] 材料专硕322分 +11 哈哈哈吼吼吼哈 2026-04-01 11/550 2026-04-02 10:52 by lnilvy
[考研] 物理学调剂 +4 小羊36 2026-03-30 4/200 2026-03-31 16:16 by lishahe
信息提示
请填处理意见