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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 哪位有1stopt软件把我拟合一下数据好吗？我用matlab非线性曲线拟合lsqnonli误差太大
191/2返回列表12下一页
查看: 2904  |  回复: 18
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者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）。
谢谢啦！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-12-12 23:54:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

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

feixiaolin

荣誉版主 (文坛精英)

优秀版主
感觉模型不对。应该具有类似于1/f^m的形式。
赞一下(1人) 回复此楼
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
赞一下(1人) 回复此楼
4楼2013-12-13 10:21:09
 已阅   回复此楼   关注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
赞一下(1人) 回复此楼
6楼2013-12-13 15:47:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

不知道为什么不一开始就将条件写清楚？

均方差(RMSE):13022.2141607379
残差平方和(SSE):16957806164.8122
相关系数(R): 0.999999236591807
相关系数之平方(R^2): 0.999998473184197
确定系数(DC): 0.999998472851297
F统计(F-Statistic): 921693.387349975

参数                  最佳估算
--------------------        -------------
par8        -73034.4843701816
par2        1102649.0406055
par3        -555.829042669526
par5        -57900.3471518271
par6        106.28243222226
par9        25147768.6160925
par10        -70.3275713005673
par12        75428.7724722937
par13        -9.99023890475422
par14        -22.1155742342423
par1        -0.199654791037603
par4        0.905553964621441
par7        0.999995295119092
par11        0.999562964175576
赞一下(1人) 回复此楼
10楼2013-12-14 23:08:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
参数名        参数值
par8        1960.22460122061
par2        2034520.31550587
par3        -418.577772243193
par5        942776.811502915
par6        6.07601523064739
par9        28642022.2437356
par10        -69.6342108036788
par12        -5072891.66359184
par13        -66.3096238997002
par1        -0.200265686103831
par4        0.845402442610382
par7        0.355788573157042
par11        0.990825917868463
par14        0.799847352320103
哪位有1stopt软件把我拟合一下数据好吗？我用matlab非线性曲线拟合lsqnonli误差太大-2
c30.jpg


哪位有1stopt软件把我拟合一下数据好吗？我用matlab非线性曲线拟合lsqnonli误差太大-3
c31.jpg
赞一下(1人) 回复此楼
13楼2013-12-15 21:53:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

lilytian_11

新虫 (小有名气)
呵呵！做点说的仓促！f是频率，e1和e2是我要你和的复介电常数的实部与虚部，(即 e=e1-j*e2），上面的e就是复介电常数，只不过没画成实部与虚部   而已，我要拟合的就是这个公式里的par1-par14个参数。麻烦给位帮帮忙。我在matlab中用非线性拟合效果十分不好！下载了个1stopt不能用。
赞一下 回复此楼
5楼2013-12-13 12:08:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

这组结果似乎更好（多花点时间效果会更好）：
均方差(RMSE):7127.64005487103
残差平方和(SSE):5080325275.18019
相关系数(R): 0.999999748138628
相关系数之平方(R^2): 0.99999949627732
确定系数(DC): 0.999999494594063
F统计(F-Statistic): 2777695.33037448

参数                  最佳估算
--------------------        -------------
par1        -0.199968813555673
par2        1339288.89927689
par3        -495.020793280607
par4        0.895286908990819
par5        -6384.18174038482
par6        -61.7284581663964
par7        -1.99044328205707
par8        8400.54248223923
par9        224823.16878659
par10        16.1681045888452
par11        21.2506313470173
par12        24549005.7016819
par13        -70.9911667494408
par14        1.00929177593573
赞一下 回复此楼
7楼2013-12-13 16:17:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lilytian_11

新虫 (小有名气)
感谢dingd，我要拟合得这14个参数中par(4)、par(7)、par(11)应该在（0——1）之间，能不能再麻烦您给我再算一下！谢谢
赞一下 回复此楼
8楼2013-12-14 19:54:58
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lilytian_11

新虫 (小有名气)
再次谢dingd，上边我说的少了一个参数
我要拟合得这14个参数中par(4)、par(7)、par(11)、par(14)的值应该在（0——1）之间，能不能再麻烦您给我再算一下！谢谢
赞一下 回复此楼
9楼2013-12-14 20:11:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lilytian_11 的主题更新
191/2返回列表12下一页
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 每次骚扰女学生的都是院系领导，而不是普通教师，小编们要注意措辞正确 +9 zju2000 2024-06-15 11/550 2024-06-16 14:49 by appleapple2
[基金申请] 希望今年自己国自然面上项目和老婆青年项目能中！ +7 恐龙爸爸 2024-06-14 7/350 2024-06-16 14:48 by redfish105
[教师之家] 请问事业编制和年薪制冲突吗？ +7 ZHONGWU_U 2024-06-14 7/350 2024-06-16 14:27 by 听风Ming
[博后之家] 山东大学(青岛)“天然药物生物智造”课题组 招聘“博士后”(年薪20.4-55.6万元) +5 第二种态度 2024-06-11 8/400 2024-06-16 10:14 by 午睡未进行
[考博] 34岁读博士晚吗 +26 emitdne 2024-06-13 26/1300 2024-06-16 07:16 by liyeqik
[找工作] 江西双非一本和四川双一流高校如何选择？ 5+7 寒山敲钟 2024-06-12 23/1150 2024-06-16 06:52 by puterde
[基金申请] 面青地会评时间 +5 tanjydd 2024-06-15 5/250 2024-06-15 23:08 by andywei1028
[基金申请] 有没有机械的前辈分享一下评上海优都是什么成果啊 +7 wulala800 2024-06-10 7/350 2024-06-15 09:33 by 晓目崇
[论文投稿] 投了一篇4区的SCI，审稿人一个拒稿，一个小修，编辑给了大修。 +9 安稳22123 2024-06-13 10/500 2024-06-14 23:45 by jurkat.1640
[基金申请] 面上基金有一个（两个）C是不是就没戏了.... 5+3 zzzm116 2024-06-13 17/850 2024-06-14 22:26 by lzt8076
[基金申请] 博士后基金需要结题吗？ +8 zhouchuck 2024-06-13 8/400 2024-06-14 17:27 by liuyupu132
[考博] 申博找导师 +4 疏影横斜水清浅3 2024-06-13 6/300 2024-06-14 14:31 by zxl_1105
[论文投稿] 最近写了一篇控制优化领域的文章，可以投哪里啊？有没有水一些的期刊推荐 +7 香瓜木香 2024-06-12 13/650 2024-06-14 07:05 by 香瓜木香
[论文投稿] 文章proof要求使用机构的邮箱 5+3 不可不信缘 2024-06-11 11/550 2024-06-14 07:00 by 3001160025
[基金申请] 工材E10口函评结束了吗 10+3 我1的飞翔 2024-06-13 5/250 2024-06-14 06:35 by nono2009
[基金申请] 博士后面上项目状态还是专家评审吗 10+9 Thatcheremu 2024-06-13 55/2750 2024-06-13 21:23 by 乌合麒麟
[基金申请] 连续两年医学口青年项目初审体会 +11 进击的荣耀 2024-06-09 18/900 2024-06-13 17:27 by 进击的荣耀
[硕博家园] 科研求助 +5 杲www 2024-06-12 6/300 2024-06-13 16:16 by 姓李名明
[论文投稿] with editor日期变更 +3 慎独的小花卷 2024-06-12 8/400 2024-06-13 11:00 by 慎独的小花卷
[论文投稿] 摩擦磨损论文投稿 +3 jmysan 2024-06-12 3/150 2024-06-13 08:36 by 莱茵润色
信息提示
关闭
请填处理意见
关闭