24小时热门版块排行榜

返回列表

misaka86

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 627.3
帖子: 60
在线: 28.6小时
虫号: 2339266
注册: 2013-03-12
性别: GG
专业: 光学

[求助] 1stopt模拟求助 6参数复数模拟有程序代码已有1人参与

Title "Debye water";
ComplexStr = i;
Parameter a, b, c, d, e, f;
Variable x, y[realPart],y[imagPart];
Function y=a+(b-c)/(1+i*2*pi*x*e)+(c-d)/(1+i*2*pi*x*f);
Data;
2.92969E11 5.0493       4.7848
3.51563E11 5.1434       4.3787
4.10157E11 5.13917       3.89761
4.6875E11 5.03814       3.70969
5.27344E11 4.99643       3.50592
5.85938E11 4.8553       3.32553
6.44532E11 4.73974       3.2601
7.03126E11 4.57753       3.09131
7.6172E11 4.40119       2.93438
8.20313E11 4.31773       2.76128
8.78907E11 4.22822       2.56528
9.37501E11 4.21345       2.42731
9.96095E11 4.21551       2.28549
1.05469E12 4.21018       2.20991
1.11328E12 4.23369       2.14796
1.17188E12 4.21764       2.09597
1.23047E12 4.20956       2.08097
1.28906E12 4.19519       2.05444
1.34766E12 4.15119       2.0478
1.40625E12 4.11325       2.04438
1.46485E12 4.0412       2.01188
1.52344E12 3.99479       1.96784
1.58203E12 3.9394       1.92388
1.64063E12 3.91918       1.87774
1.69922E12 3.89392       1.76959
1.75781E12 3.89626       1.75711
1.81641E12 3.8712       1.64121
1.875E12          3.9422       1.78256
1.9336E12 3.94556       1.57435
1.99219E12 3.85584       1.73505

回复此楼

» 猜你喜欢

求助中医院临床中医中标的青年标书学习，谢谢已经有6人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有259人回复
物理学I论文润色/翻译怎么收费? 已经有263人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有47人回复
25届硕士申请，材料计算方向已经有2人回复
课题组招1-2名硕士已经有1人回复
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上已经有0人回复
斯坦福大学与Elsevier世界前2%学者新排名查询已经有1人回复
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有19人回复
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计已经有1人回复
apl投稿状态已经有0人回复

Нипуханипера

1楼 2016-08-04 09:59:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
misaka86: 金币+20, ★★★★★最佳答案, 感谢您的帮助 2016-08-05 13:17:57

参数组值似乎不唯一：
1：
均方差(RMSE): 0.1151868529408
残差平方和(SSE): 0.796080665424332
相关系数(R): 0.979421705966726
相关系数之平方(R^2): 0.959266878118772
修正R平方(Adj. R^2): 0.954754164383851
确定系数(DC): 0.952126547248776
F统计(F-Statistic): 71.670774807762

参数       最佳估算
-------------------- -------------
a 3.45529988534589
b -0.297845276953659
c -2.89632987059642
d 18.9910742182084
e -2.09276637084744E-13
f 2.87573889727082E-12

2：
均方差(RMSE): 0.1151868529408
残差平方和(SSE): 0.796080665424332
相关系数(R): 0.979421705967326
相关系数之平方(R^2): 0.959266878119948
修正R平方(Adj. R^2): 0.954754164385121
确定系数(DC): 0.952126547252169
F统计(F-Statistic): 71.6707748132263

参数       最佳估算
-------------------- -------------
a 3.45529988541862
b -22.3048300048426
c -0.417425925914286
d -3.01591051979003
e 2.875738896093E-12
f -2.09276637110207E-13

3：
均方差(RMSE): 0.1151868529408
残差平方和(SSE): 0.796080665424331
相关系数(R): 0.979421705966883
相关系数之平方(R^2): 0.959266878119079
修正R平方(Adj. R^2): 0.954754164384169
确定系数(DC): 0.952126547251118
F统计(F-Statistic): 71.6707748109358

参数       最佳估算
-------------------- -------------
a 3.45529988540051
b -21.8938777471039
c -0.0064736660315961
d -2.60495825985104
e 2.8757388963447E-12
f -2.09276637103644E-13

公式有过拟合现象吧：(b-c)可以用b表示，(c-d)用c表示，d去掉，公式变为：
y=a+b/(1+i*2*pi*x*e)+c/(1+i*2*pi*x*f);

结果稳定唯一：
均方差(RMSE): 0.1151868529408
残差平方和(SSE): 0.796080665424331
相关系数(R): 0.979421705967872
相关系数之平方(R^2): 0.959266878121017
修正R平方(Adj. R^2): 0.954754164386264
确定系数(DC): 0.95212654725604
F统计(F-Statistic): 93.3213213792955

参数       最佳估算
-------------------- -------------
a 3.45529988550685
b 2.59848459417536
e -2.09276637141881E-13
c -21.8874040662521
f 2.8757388945819E-12

赞一下