24小时热门版块排行榜    

查看: 495  |  回复: 2

yejianling

新虫 (正式写手)

[求助] 求助高版本的1stOpt代跑个程序 已有1人参与

代码如下:
Parameter p1,p2,p3,p4;
Variable t,x1,x2,x3,x4;
ODEFunction
x2'=(x3/270)*(x1-x2)-p1*exp((-p2+p3*96485*(x4-p4))/(8.314*363))*x2^2;
data;
t,x1,x2,x3,x4;
0.1        165        69        100.00         361.00
0.3        115        71        110.00         356.00
0.5        103        62        112.50         358.00
0.7        150        69        118.75         380.00
0.9        128        95        112.50         380.00
1.1        140        63        100.00         383.00
1.3        110        36        100.00         383.00
1.5        103        33        118.75         370.00
1.7        96        31        137.50         367.00
1.9        70        25        100.00         356.00
2.1        83        31        125.00         373.00
2.3        90        81        137.50         363.00
2.5        100        88        125.00         353.00
2.7        112        60        137.50         346.00
2.9        119        35        137.50         353.00
3.1        66        30        112.50         353.00
3.3        70        35        125.00         373.00
3.5        160        63        112.50         370.00
3.7        147        74        112.50         370.00
3.9        140        86        112.50         373.00
4.1        108        76        107.50         373.00
4.3        109        79        75.00         363.00
4.5        76        50        125.00         353.00
4.7        79        47        131.00         363.00
4.9        82        43        112.00         362.00
5.1        136        61        135.00         373.00
5.3        145        69        125.00         363.00
5.5        143        52        110.00         366.00
5.7        159        63        125.00         370.00
5.9        178        76        100.00         370.00
6.1        112        56        100.00         360.00
6.3        121        70        81.25         364.00
6.5        152        47        112.50         363.00
6.7        157        60        100.00         373.00
6.9        162        82        93.75         383.00
7.1        170        63        118.75         363.00
7.3        62        49        112.50         291.00
7.5        83        78        112.50         356.00
7.7        227        95        137.50         359.00
7.9        226        75        125.00         360.00
8.1        223        47        125.00         355.00
8.3        116        44        143.75         366.00
8.5        130        70        143.75         370.00
回复此楼
学习
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
楼主的公式明显有问题:
p1*exp((-p2+p3*96485*(x4-p4))/(8.314*363))
完全等同于:
exp((-p2+p3*96485*(x4))/(8.314*363))

p1和p4多余,去掉后效果是一样的,否则过拟合,参数组值不唯一。

计算用时(时:分:秒:微秒): 00:01:01:312
均方差(RMSE): 17.9229682919401
残差平方和(SSE): 13491.7772805435
相关系数(R): 0.290562575469603
相关系数之平方(R^2): 0.0844266102635286
修正R平方(Adj. R^2): 0.0374741287385814
确定系数(DC): 0.0816922624187681
F统计(F-Statistic): 2.44181044608464

参数                  最佳估算
--------------------        -------------
p2        20003.4712902898
p3        0.000151432517221469

» 本帖已获得的红花(最新10朵)

2楼2016-11-07 21:01:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yejianling

新虫 (正式写手)

送红花一朵
引用回帖:
2楼: Originally posted by dingd at 2016-11-07 07:01:17
楼主的公式明显有问题:
p1*exp((-p2+p3*96485*(x4-p4))/(8.314*363))
完全等同于:
exp((-p2+p3*96485*(x4))/(8.314*363))

p1和p4多余,去掉后效果是一样的,否则过拟合,参数组值不唯一。

计算用时(时:分 ...

谢谢了,不过残差平方和有点大了

发自小木虫Android客户端
学习
3楼2016-11-08 11:35:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yejianling 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +3 匹克i 2026-03-23 3/150 2026-03-23 23:18 by peike
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +4 石能伟 2026-03-17 4/200 2026-03-23 17:51 by 17862566385
[考研] 333求调剂 +6 87639 2026-03-21 10/500 2026-03-23 10:41 by Iveryant
[考研] 276求调剂 +3 YNRYG 2026-03-21 4/200 2026-03-23 08:31 by 醉在风里
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 318求调剂 +4 plum李子 2026-03-21 7/350 2026-03-22 14:17 by ColorlessPI
[考研] 286分人工智能专业请求调剂愿意跨考! +4 lemonzzn 2026-03-17 8/400 2026-03-21 22:49 by lemonzzn
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
信息提示
请填处理意见