| 查看: 1118 | 回复: 5 | ||
kirkendall木虫 (小有名气)
|
[求助]
求大神用高版本1stopt跑一下复数拟合问题 已有2人参与
|
|
方程和数据如下，代码不知写得对不对，只用过1stopt1.5的版本。function中的i是虚部。按照方程的物理意义，b拟合出来应该大概在450左右。请大神们帮忙跑一下，在线等。 Parameters a,b,c; Variable x,y; Function y=((1+2.88*a^(-2.375))*8e-7*b^(1/a))/((1/c+i*x*12*8.8452e-14)^(1-1/a)); Data; 0.062831854 129739.05+128412.33i 0.082561056 107207.2+122846.48i 0.108510612 79490.176+109031.98i 0.142628309 58922.001+93550.392i 0.187364589 44202.459+80578.558i 0.246238036 31708.82+67055.782i 0.32364688 23197.226+55479.787i 0.42530882 17016.607+45832.97i 0.558889341 12532.017+37091.828i 0.734504373 8362.3418+28251.222i 0.965222941 6262.8977+23358.677i 1.268449469 4613.0619+19371.435i 1.666866255 3192.479+14077.537i 2.190569758 2343.8024+11081.69i 2.878704223 1730.1589+8715.4978i 3.783043097 1285.1065+6851.4073i 4.971444784 957.3885+5362.7526i 6.533256179 750.95267+4265.047i 8.58565869 556.60047+3345.4825i 11.28277885 434.50485+2613.6751i 14.8272494 329.96661+2047.5791i 19.48522607 253.4584+1598.6387i 25.60643095 197.95765+1248.7251i 33.65066755 150.24615+967.96302i 44.22193849 116.74476+761.28861i 58.11412424 90.84942+598.97904i 76.37048491 70.89787+460.96938i 100.3620744 57.76081+361.51781i 131.8905961 44.94704+281.36163i 173.3236799 35.85633+219.77967i 227.7728221 22.65702+170.79197i 299.3270569 25.2096+138.26227i 393.3598901 17.62742+103.39121i 516.9328053 14.75816+84.24969i 679.3259417 11.69736+63.07746i 892.7343338 8.98806+49.85415i 1173.184503 8.0246+40.88545i 1541.737309 6.16703+31.65768i 2026.069995 4.56576+25.38283i 2662.554603 3.65626+19.4245i 3498.989082 2.57597+16.18998i 4598.187254 2.36727+12.57802i 6042.695662 1.69159+10.16771i 7940.992558 1.84486+8.00604i |
回复此楼» 猜你喜欢
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
物理学I论文润色/翻译怎么收费?
已经有158人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 这样的动力学方程能够有1stopt拟合出参数来么
已经有4人回复
- 看了下1stOpt的介绍，感觉略显浮夸
已经有14人回复
- 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
已经有6人回复
- 求大神跑下1stOpt，一个sharemodel，谢谢！
已经有7人回复
- 求助--用1stopt拟合蓝宝石衬底的GaN薄膜红外反射光谱
已经有6人回复
- 求大神帮忙拟合一个非线性方程，求出模型参数
已经有15人回复
- 用1stOpt进行复数曲线拟合时，结果出错。急求教！！！
已经有5人回复
- 谁能帮我用 高版本的 1stopt拟合一下，需要涉及复数
已经有7人回复
- 求大哥大姐帮跑一下1stopt~新手啊，版本低啊
已经有9人回复
- 用1stopt拟合非线性方程结果与其他软件拟合结果差异大
已经有5人回复
- 运算过程有复数出现，求高版本1stOpt
已经有6人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- 1stOpt求解常微分方程边值问题
已经有5人回复
xxz903
金虫 (小有名气)
- 应助: 33 (小学生)
- 金币: 1045.6
- 散金: 190
- 红花: 4
- 帖子: 253
- 在线: 66.7小时
- 虫号: 1360336
- 注册: 2011-08-03
- 性别: GG
- 专业: 化工系统工程
2楼2014-12-03 22:16:44
kirkendall
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2136.9
- 散金: 4
- 帖子: 106
- 在线: 160.7小时
- 虫号: 848554
- 注册: 2009-09-15
- 专业: 等离子体物理
|
不好意思，只用过1.5版的，不知道复数拟合的格式。在本版搜索了一下，重写的格式如下，不知对否： ComplexStr = i; Parameter a,b,c; Variable x,y,y; Function y=((1+2.88*a^(-2.375))*8e-7*b^(1/a))/(1/c+i*x*2*pi*12*8.8452e-14)^(1-1/a); Data; //x real(y) imag(y) 0.01 129739.05 128412.33 0.01314 107207.2 122846.48 0.01727 79490.176 109031.98 0.0227 58922.001 93550.392 0.02982 44202.459 80578.558 0.03919 31708.82 67055.782 0.05151 23197.226 55479.787 0.06769 17016.607 45832.97 0.08895 12532.017 37091.828 0.1169 8362.3418 28251.222 0.15362 6262.8977 23358.677 0.20188 4613.0619 19371.435 0.26529 3192.479 14077.537 0.34864 2343.8024 11081.69 0.45816 1730.1589 8715.4978 0.60209 1285.1065 6851.4073 0.79123 957.3885 5362.7526 1.0398 750.95267 4265.047 1.36645 556.60047 3345.4825 1.79571 434.50485 2613.6751 2.35983 329.96661 2047.5791 3.10117 253.4584 1598.6387 4.07539 197.95765 1248.7251 5.35567 150.24615 967.96302 7.03814 116.74476 761.28861 9.24915 90.84942 598.97904 12.15474 70.89787 460.96938 15.97312 57.76081 361.51781 20.99104 44.94704 281.36163 27.58532 35.85633 219.77967 36.25117 22.65702 170.79197 47.63938 25.2096 138.26227 62.60517 17.62742 103.39121 82.27241 14.75816 84.24969 108.11808 11.69736 63.07746 142.08308 8.98806 49.85415 186.71811 8.0246 40.88545 245.37511 6.16703 31.65768 322.45905 4.56576 25.38283 423.75872 3.65626 19.4245 556.8814 2.57597 16.18998 731.82422 2.36727 12.57802 961.72487 1.69159 10.16771 1263.8482 1.84486 8.00604 |
3楼2014-12-03 22:42:03
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
感谢参与，应助指数 +1
|
参数都视为复数： 均方差(RMSE):593.045693716506 残差平方和(SSE):30949881.1455409 相关系数(R): 0.999824052026595 相关系数之平方(R^2): 0.999648135010879 确定系数(DC): 0.999627438981964 F统计(F-Statistic): 14789.3401273682 参数 最佳估算 -------------------- ------------- a.实部 6.54549367633696 a.虚部 -5.89692343813397 b.实部 -3.29031694755989E-7 b.虚部 -2.73162943837164E-20 c.实部 -14965721627473.6 c.虚部 -9689614065205.55  c135.jpg  c136.jpg |
» 本帖已获得的红花（最新10朵）
kirkendall4楼2014-12-03 22:49:33
kirkendall
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2136.9
- 散金: 4
- 帖子: 106
- 在线: 160.7小时
- 虫号: 848554
- 注册: 2009-09-15
- 专业: 等离子体物理
5楼2014-12-04 09:23:50
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
6楼2014-12-04 10:29:58
精华III:
20