24小时热门版块排行榜    

查看: 2908  |  回复: 18
【悬赏金币】回答本帖问题,作者lilytian_11将赠送您 30 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lilytian_11

新虫 (小有名气)

[求助] 哪位有1stopt软件把我拟合一下数据好吗?我用matlab非线性曲线拟合lsqnonli误差太大已有3人参与

大家好!哪位有1stopt软件把我拟合一下数据好吗?我用matlab非线性曲线拟合lsqnonli误差太大,对初值太敏感。
f        e1        e2
10        2.57E+07        3.63E+08
15.849        2.50E+07        2.32E+08
25.119        2.33E+07        1.51E+08
39.811        2.01E+07        1.01E+08
63.096        1.49E+07        6.96E+07
100        9.33E+06        4.80E+07
158.49        5.11E+06        3.23E+07
251.19        2.70E+06        2.12E+07
398.11        1.46E+06        1.38E+07
630.96        8.04E+05        8.94E+06
1000        4.35E+05        5.78E+06
1584.9        2.31E+05        3.71E+06
2511.9        1.24E+05        2.37E+06
3981.1        68552        1.51E+06
6309.6        40496        9.64E+05
10000        25909        6.13E+05
15849        18055        3.90E+05
25119        13634        2.48E+05
39811        11009        1.58E+05
63096        9324.4        1.01E+05
100000        8089.2        65074
158490        6973.4        42501
251190        5749.1        28345
398110        4354.1        19309
630960        2962.4        13267
1000000        1836.4        9040
1584900        1082.7        6060.7
2511900        636.43        4005.7
3981100        386.74        2622.8
6309600        248.48        1708.2
10000000        170.73        1109.3
15849000        125.8        719.37
25119000        99.06        466.26
39811000        82.708        302.21
63096000        72.482        195.98
100000000        65.972        127.25
158490000        61.765        82.84
251190000        59.002        54.235
398110000        57.144        35.949
630960000        55.835        24.472
1000000000        54.811        17.591
1584900000        53.842        13.949
2511900000        52.654        12.752
3981100000        50.845        13.54
6309600000        47.801        15.909
10000000000        42.764        19.109
15849000000        35.36        21.688
25119000000        26.512        21.952
39811000000        18.313        19.44
63096000000        12.322        15.363
公式是:
e=par(8)+par(2)./(1+(j*X./par(3)).^par(4))+par(5)./(1+(j*X./par(6)).^par(7))+par(9)./(1+(j*X./par(10)).^par(11))+par(12)./(1+(j*X./par(13)).^par(14))+par(1)./(j*2*3.1416.*X*(8.85e-12));
要拟合的参数是par(1)、par(2)、par(3)、par(4)、par(5)、par(6)、par(7)、par(8)、par(9)、par(10)、par(11)、par(12)、par(13)、par(14)。
谢谢啦!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

公式中的j是虚数符号吧,1stOpt求解代码及结果:
CODE:
Parameter par(14);
ComplexStr = j;
Variable x, e[realPart], e[imagPart];
Function
e=par8+par2/(1+(j*X/par3)^par4)+par5/(1+(j*X/par6)^par7)+par9/(1+(j*X/par10)^par11)+par12/(1+(j*X/par13)^par14)+par1/(j*2*3.1416*X*(8.85e-12));
Data;
//f        e1        e2
10        2.57E+07        3.63E+08
15.849        2.50E+07        2.32E+08
25.119        2.33E+07        1.51E+08
39.811        2.01E+07        1.01E+08
63.096        1.49E+07        6.96E+07
100        9.33E+06        4.80E+07
....

均方差(RMSE):11353.4863245101
残差平方和(SSE):12890165172.0838
相关系数(R): 0.999999282992593
相关系数之平方(R^2): 0.9999985659857
确定系数(DC): 0.999998541615234
F统计(F-Statistic): 959707.257043923

参数                  最佳估算
--------------------        -------------
par1        -0.19997550069211
par2        69323.2691995449
par3        99.2809228245236
par4        -132.634485633095
par5        -1716075.69507839
par6        -399.918912293264
par7        -0.87287541809133
par8        1648129.75621163
par9        279271.025103101
par10        -16.3273493394721
par11        18.9137493447292
par12        24164131.0587702
par13        -70.3023756542921
par14        1.01534559728677
哪位有1stopt软件把我拟合一下数据好吗?我用matlab非线性曲线拟合lsqnonli误差太大
c1.jpg


哪位有1stopt软件把我拟合一下数据好吗?我用matlab非线性曲线拟合lsqnonli误差太大-1
c2.jpg

6楼2013-12-13 15:47:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
公式和数据的对应关系说清楚点啊,数据的f、e1、e2对应公式里各是谁?(e、x、j)
2楼2013-12-13 09:00:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

感觉模型不对。应该具有类似于1/f^m的形式。
3楼2013-12-13 10:17:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
2楼: Originally posted by dingd at 2013-12-13 09:00:55
公式和数据的对应关系说清楚点啊,数据的f、e1、e2对应公式里各是谁?(e、x、j)

f=frenqucy
e1=real part
e2=image part
4楼2013-12-13 10:21:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 对象没有,还非常想读博,难以抉择 +18 pvrw0224 2024-06-23 40/2000 2024-06-24 19:22 by tjushede
[硕博家园] 回家两天,不想呆了 +4 368ghnf 2024-06-22 8/400 2024-06-24 19:03 by 我是王小帅
[基金申请] 评审专家会不会很在意申请人的单位啊 +7 lancet0903 2024-06-24 7/350 2024-06-24 17:23 by LNP@mRNA
[基金申请] 这样的说辞是上会了吗 +9 学员d3zYCz 2024-06-24 12/600 2024-06-24 16:46 by yingxz
[硕博家园] 申博雅思求助 +4 刘嘉研究究 2024-06-24 4/200 2024-06-24 16:20 by 仙女科研
[基金申请] 刚刚收到科研之友邮件 +28 olivermiaoer 2024-06-19 42/2100 2024-06-24 15:37 by jiann_li84
[金属] 寻找钛合金热压缩代做 +4 liuyang358 2024-06-23 4/200 2024-06-24 13:50 by 搬砖狗不放弃
[基金申请] 说博后基金7月出的真打电话了吗? +12 antonysole 2024-06-24 14/700 2024-06-24 13:39 by sizhouyi
[基金申请] 青年和面上,哪个上会难度更大 +12 今晚推荐22 2024-06-21 18/900 2024-06-24 11:08 by 半简体
[催化] 求助原位红外 +3 jfdhj 2024-06-21 7/350 2024-06-24 10:52 by qvhm2609
[硕博家园] 数据不好 +3 Hetai 2024-06-23 5/250 2024-06-23 20:36 by 82年拉菲
[基金申请] 演员杨幂在核心期刊发表论文,中国广播电视学期刊 +8 babu2015 2024-06-22 11/550 2024-06-23 11:49 by jurkat.1640
[硕博家园] 夏至,要不要硕博联谊 +3 我是王小帅 2024-06-21 5/250 2024-06-22 22:51 by 我是王小帅
[公派出国] CSC博士联培对将来就业有用吗 +3 也就这样 2024-06-22 3/150 2024-06-22 22:04 by 326lhpqk
[基金申请] 江南大学到瑞士招聘,称取消非升即走,改预聘+长聘 +24 babu2015 2024-06-18 29/1450 2024-06-22 16:29 by 风今25
[论文投稿] ACS 编辑的意见 10+3 哈哈妞1993 2024-06-20 3/150 2024-06-21 17:06 by 投个论文
[考博] 2025年博士申请,可先做科研助理 +5 limit888 2024-06-18 8/400 2024-06-21 12:45 by 半简体
[基金申请] 青年基金会评专家到底是怎么会评的呀?主审专家是不是一般不会改动系统按函评给的顺序 5+4 他山攻玉之石 2024-06-18 18/900 2024-06-20 16:33 by 他山攻玉之石
[论文投稿] 投稿求助 +4 平凡的日子 2024-06-19 5/250 2024-06-20 16:24 by yueyueyue@
[基金申请] 基金得中 +4 woaini0218 2024-06-18 4/200 2024-06-19 17:27 by mengzl
信息提示
请填处理意见