【答案】应助回帖

» 猜你喜欢

求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有169人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有22人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复

1楼2019-03-27 20:34:14

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31014.8
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 855.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

感谢参与，应助指数 +1

1stOpt 的国外试用版 5.5版 4参数限制
http://mob.muchong.com/t-8050245-1

[ 发自手机版 http://muchong.com/3g ]

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

嬛嬛爱豆豆

数值计算

2楼2019-03-27 21:13:25

嬛嬛爱豆豆

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 4192
散金: 602
红花: 1
帖子: 605
在线: 125.4小时
虫号: 3514884
注册: 2014-11-02
专业: 无机化学

送红花一朵

引用回帖:

2楼: Originally posted by 独孤神宇 at 2019-03-27 21:13:25
1stOpt 的国外试用版 5.5版 4参数限制
http://mob.muchong.com/t-8050245-1

感谢回复，有6个参数，想试下1.5版本的

3楼2019-03-28 09:43:39

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31014.8
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 855.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

引用回帖:

3楼: originally posted by 嬛嬛爱豆豆 at 2019-03-28 09:43:39
感谢回复，有6个参数，想试下1.5版本的...

试用版，算不了6个参数的。你可以将代码贴上来

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

嬛嬛爱豆豆

数值计算

4楼2019-03-28 09:53:33

嬛嬛爱豆豆

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 4192
散金: 602
红花: 1
帖子: 605
在线: 125.4小时
虫号: 3514884
注册: 2014-11-02
专业: 无机化学

送红花一朵

引用回帖:

4楼: Originally posted by 独孤神宇 at 2019-03-28 09:53:33
试用版，算不了6个参数的。你可以将代码贴上来...

代码奉上，十分感谢！请大侠帮忙代跑一次

Title " Model";
Parameters b,c,d,e,f,h;
Variable T,x1,y;
Constant R=8.314;
Constant a=0.3;
ConstStr x2=0.1*(1-x1);
ConstStr x3=0.9*(1-x1);
ConstStr t12=b/(R*T);
ConstStr t21=c/(R*T);
ConstStr t31=d/(R*T);
ConstStr t13=e/(R*T);
ConstStr t32=f/(R*T);
ConstStr t23=h/(R*T);
Function y=exp((exp(-a*t21)*x2+exp(-a*t31)*x3)*(t21*exp(-a*t21)*x2+t31*exp(-a*t31)*x3)/(x1+x2*exp(-a*t21)+x3*exp(-a*t31))^2
+(t13*exp(-a*t13)*x3*x3+exp(-a*t13)*exp(-a*t23)*x2*x3*(t13-t23))/(x3+x1*exp(-a*t13)+x2*exp(-a*t23))^2+(t12*exp(-a*t12)*x2*x2+exp(-a*t12)*exp(-a*t32)*x2*x3*(t12-t32))/(x2+x1*exp(-a*t12)+x3*exp(-a*t32))^2);
Data;
T,x1,y=
293.15 0.01625 4.16012
298.15 0.01864 4.16046
303.15 0.02192 4.04024
308.15 0.02591 3.88661

5楼2019-03-28 16:07:51

嬛嬛爱豆豆

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 4192
散金: 602
红花: 1
帖子: 605
在线: 125.4小时
虫号: 3514884
注册: 2014-11-02
专业: 无机化学

引用回帖:

4楼: Originally posted by 独孤神宇 at 2019-03-28 09:53:33
试用版，算不了6个参数的。你可以将代码贴上来...

5楼数据没有贴完，完整数据如下：
T,x1,y=
293.15 0.01625 4.16012
298.15 0.01864 4.16046
303.15 0.02192 4.04024
308.15 0.02591 3.88661
313.15 0.03076 3.70732
318.15 0.03629 3.54462
323.15 0.04234 3.41427
328.15 0.04979 3.25127
333.15 0.05738 3.14852

6楼2019-03-28 16:09:39

独孤神宇

版主 (知名作家)

应助: 490 (硕士)
贵宾: 0.008
金币: 31014.8
散金: 802
红花: 122
沙发: 1
帖子: 5600
在线: 855.5小时
虫号: 3522474
注册: 2014-11-06
性别: GG
专业: 机械动力学
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
嬛嬛爱豆豆: 金币+20, ★★★★★最佳答案, 十分感谢！ 2019-03-28 20:15:24

引用回帖:

6楼: Originally posted by 嬛嬛爱豆豆 at 2019-03-28 16:09:39
5楼数据没有贴完，完整数据如下：
T,x1,y=
293.15 0.01625 4.16012
298.15 0.01864 4.16046
303.15 0.02192 4.04024
308.15 0.02591 3.88661
313.15 0.03076 3.70732
318.15 0.03629 3.54462
323.15 0.042 ...

优化算法: 通用全局优化算法(UGO1)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.0015578099942271
残差平方和(SSR): 3.15480357154798E-5
相关系数(R): 0.999990115350376
相关系数之平方(R^2): 0.999980230798458
修正R平方(Adj. R^2): 0.999973641064611
确定系数(DC): 0.99998023064809
卡方系数(Chi-Square): 4.72247457185271E-6
F统计(F-Statistic): 70813.5551328973

参数最佳估算
---------- -------------
b -5788.70837047638
c 11752.3000158702
d 33274.8006880531
e 25847.2491575168
f 51644.7184858163
h 792529371.588477

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

No. 实测值 y 计算值 y
1 4.16012 4.15989392027371
2 4.16046 4.16072052222905
3 4.04024 4.04075052758997
4 3.88661 3.88571920524604
5 4.16012 4.15989392027371
6 4.16046 4.16072052222905
7 4.04024 4.04075052758997
8 3.88661 3.88571920524604
9 3.70732 3.70779757506783
10 3.54462 3.54561945681809
11 3.41427 3.41131070742377
12 3.25127 3.25526126939215
13 3.14852 3.14670695704689

数值计算

7楼2019-03-28 17:08:08