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

回复此楼

» 猜你喜欢

散金币祈福已经有186人回复
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有274人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
深圳大学应用技术学院招聘凝聚态物理博士后已经有113人回复
深圳大学2027级光电信息工程/物理学研究生推免已经有69人回复
求助中医院临床中医中标的青年标书学习，谢谢已经有4人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有257人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有39人回复

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

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

赞一下