| 查看: 1349 | 回复: 1 | |||
[交流]
1stopt拟合求助，好难啊 已有1人参与
|
求大神帮忙拟合一下啊，本人用matlab非线性拟合出来的效果不理想，据说1stopt可以拟合带复数的公式，我尝试用1stopt拟合，结果定义的复数符号j也被拟合出来了 待拟合的方程式电力线传输函数幅频特性的一个4阶方程（4个只变gp和dp的表达式之和），所以我取了绝对值换算成对数使其单位变成dB，测试的数据放在附件的excel里，第一列是变量x，第二列是对应因变量y. ![\"1stopt拟合求助，好难啊\"](\"https://muchongimg.xmcimg.com/data/bcs/2016/0116/w428h3189795_1452916705_322.png#opennewwindow\") 就是这个函数取得对数，vp=2e8,跪求大神帮帮忙，万分感激啊！（如果拟合效果不好可增加阶数，但是待求拟合系数也会增加） Title "H_1"; Parameters gp01,a0,a1,K,dp01,gp02,dp02,gp03,dp03,gp04,dp04; ComplexStr = j; Variable x , y; Function y=20*log10(abs(gp01*exp(-(a0+a1*x^K+2*pi*j*x/2e8)*dp01)+gp02*exp(-(a0+a1*x^K+2*pi*j*x/2e8)*dp02)+gp03*exp(-(a0+a1*x^K+2*pi*j*x/2e8)*dp03)+gp04*exp(-(a0+a1*x^K+2*pi*j*x/2e8)*dp04))) ; data; 500000 -1.582418505 519500 -1.620108522 539000 -1.657877333 558500 -1.695696353 ..... "H_1" 迭代数: 0 计算用时(时:分:秒:微秒): 00:00:42:814 优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法 计算结束原因: 达到收敛判断标准 均方差(RMSE): 3.82641829064014 残差平方和(SSE): 14656.1184118803 相关系数(R): 0.89236767319423 相关系数之平方(R^2): 0.796320064162084 决定系数(DC): 0.79632005624803 卡方系数(Chi-Square): -453.261902781698 F统计(F-Statistic): 3905.75409809133 参数 最佳估算 ---------- ------------- gp01 0.0905058486354724 a0 -3.21061454592761 a1 0.0303224962664343 K 0.00164251742836774 dp01 0.24425464628118 gp02 0.162608273800758 dp02 0.244245680941401 gp03 0.000818566218494661 dp03 0.244738688504121 gp04 0.000143396963924432 dp04 2.87024769969935 j 18.9404153768973 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : code.xlsx
2016-01-16 11:44:07, 32.59 K
» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有127人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有25人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有27人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2016-01-17 21:36:05
