应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求大神用高版本1stopt跑一下复数拟合问题
61/1返回列表
查看: 971  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-12-03 20:04:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxz903

金虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
你的编写代码有问题  建议照着手册重新写一下！
赞一下 回复此楼
2楼2014-12-03 22:16:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kirkendall

木虫 (小有名气)
引用回帖:
2楼: Originally posted by xxz903 at 2014-12-03 22:16:44
你的编写代码有问题  建议照着手册重新写一下！

不好意思，只用过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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +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
求大神用高版本1stopt跑一下复数拟合问题
c135.jpg


求大神用高版本1stopt跑一下复数拟合问题-1
c136.jpg
赞一下(1人) 回复此楼

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

kirkendall
4楼2014-12-03 22:49:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kirkendall

木虫 (小有名气)
送红花一朵
引用回帖:
4楼: Originally posted by dingd at 2014-12-03 22:49:33
参数都视为复数：

均方差(RMSE):593.045693716506
残差平方和(SSE):30949881.1455409
相关系数(R): 0.999824052026595
相关系数之平方(R^2): 0.999648135010879
确定系数(DC): 0.999627438981964
F统计(F-S ...

感谢大神回复！汗，忘了说参数都是实数，可否麻烦您再帮忙跑一下？
赞一下 回复此楼
5楼2014-12-04 09:23:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kirkendall: 金币+20, ★★★很有帮助, 好的，谢谢了。 2014-12-04 10:46:52
都是实数的话，拟合效果非常差。
赞一下 回复此楼
6楼2014-12-04 10:29:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kirkendall 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭