| 查看: 578 | 回复: 1 | ||
[求助]
求助用1stopt代跑数据 已有1人参与
|
|
求助大神用高版本1stopt代跑一下代码： 十分感谢！！！ Parameter a1<1000; Parameter 0.002<a2<0.0045; Parameter a3>0.0000001; Parameter a4<0.1 Parameter a5<0.01 Parameter a6>100 Parameter a7>100 Parameter a8<0.0001 Parameter a9>100 Parameter a10>100 Parameter a11<0.001 Parameter a12>100 ComplexStr = i; Variable x, e[realPart], e[imagPart]; Function e=a1+(0.0045-a2)/(1/a3+i*2*pi*x*a4)+1/(i*2*pi*x*a5+1/(a6+a7+i*2*pi*x*a8+a9*sqrt((a10+1/(i*2*pi*x*a11+1/a12))/(2*a9))*(asinh(a2/sqrt((a10+1/(i*2*pi*x*a11+1/a12))/(2*a9)))+coth(a2/sqrt((a10+1/(i*2*pi*x*a11+1/a12))/(2*a9)))))); Data; //f e1 e2 10010 51.2394 -20.38674 8252 53.8227 -23.53173 6812 56.8326 -26.3307 5625 60.1506 -29.0799 4643 63.7056 -32.0898 3833 67.4739 -35.4789 3164 71.5977 -39.2709 2612 76.1955 -43.2762 2153 81.2673 -47.4474 1777 86.6946 -51.7371 1470 92.5248 -56.1216 1211 98.8527 -60.7668 1001 105.702 -65.6253 825.2 113.1201 -70.7682 683.6 121.1544 -76.3614 561.5 129.6627 -82.5234 463.9 138.8346 -89.1831 380.9 148.599 -96.3405 317.4 159.0033 -104.0667 258.8 170.1897 -112.3854 214.8 182.3004 -121.4862 175.8 195.1695 -131.3691 146.5 208.086 -141.8682 122.1 222.9459 -153.9078 97.66 239.844 -166.6821 82.54 257.145 -180.5466 68.13 276.105 -196.71 56.23 296.013 -214.8168 46.42 317.58 -235.578 38.31 341.043 -258.567 31.62 367.587 -283.452 26.1 396.975 -310.707 21.54 428.97 -341.28 17.78 464.994 -375.645 14.68 504.099 -413.091 12.12 546.759 -457.884 10 593.922 -505.995 8.254 645.825 -555.765 6.813 703.416 -611.697 5.623 765.747 -674.739 4.642 837.084 -745.839 3.831 915.768 -823.812 3.162 1005.117 -913.872 2.61 1103.709 -1014.36 2.154 1221.972 -1137.6 1.778 1357.299 -1280.511 1.468 1493.574 -1425.792 1.212 1635.3 -1570.836 1 1791.483 -1717.539 0.8254 1959.753 -1873.248 0.6813 2134.659 -2037.963 0.5623 2322.126 -2218.32 0.4642 2528.79 -2422.14 0.3831 2763.42 -2649.66 0.3162 3019.38 -2898.51 0.261 3308.52 -3166.32 0.2154 3642.69 -3467.31 0.1778 4007.67 -3782.52 0.1468 4405.83 -4109.58 0.1212 4853.76 -4467.45 0.1 5387.01 -4856.13 0.08254 5986.62 -5254.29 0.06813 6624.15 -5659.56 0.05623 7311.45 -6074.31 0.04642 8053.26 -6491.43 0.03831 8868.54 -6913.29 0.03162 9754.92 -7325.67 0.0261 10714.77 -7721.46 0.02154 11745.72 -8103.03 0.01778 12835.92 -8444.31 0.01468 13985.37 -8738.19 0.01212 15186.96 -8965.71 0.01 16428.84 -9115.02@月只蓝 |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有245人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
clhclh: 金币+10, ★★★很有帮助, 十分感谢！ 2019-01-09 22:15:27
独孤神宇: 金币+5, 鼓励交流 2019-01-10 14:53:25
感谢参与，应助指数 +1
clhclh: 金币+10, ★★★很有帮助, 十分感谢！ 2019-01-09 22:15:27
独孤神宇: 金币+5, 鼓励交流 2019-01-10 14:53:25
|
参考下，似乎过拟合，多解： 1： 均方差(RMSE): 413.848207127873 残差平方和(SSR): 25005469.4272714 相关系数(R): 0.993262258305013 相关系数之平方(R^2): 0.986569913773174 修正R平方(Adj. R^2): 0.985986128777314 确定系数(DC): 0.985121768798174 F统计(F-Statistic): 392.096545775466 参数 最佳估算 -------------------- ------------- a1 462.198090627634 a2 0.00408963492348397 a3 13435297.7758474 a4 4.34611028583419E-8 a5 0.000629953572178909 a6 2811.22829859182 a7 127.837745740069 a8 6.34021810351535E-5 a9 100.002191476285 a10 115.082841542099 a11 -603273406.375728 a12 100.001250047267 2： 均方差(RMSE): 413.848207055404 残差平方和(SSR): 25005469.418514 相关系数(R): 0.993262258472133 相关系数之平方(R^2): 0.986569914105161 修正R平方(Adj. R^2): 0.985986129126611 确定系数(DC): 0.985121772183895 F统计(F-Statistic): 392.096349528075 参数 最佳估算 -------------------- ------------- a1 462.198852162116 a2 0.00343061299542595 a3 5155649.13583405 a4 1.13257201861152E-7 a5 0.000629954327168894 a6 378.856209184058 a7 100.000000110878 a8 0.0001 a9 100.000000000001 a10 113.418539877976 a11 -4620210365.48644 a12 100 |
2楼2019-01-09 21:21:02