| 查看: 2410 | 回复: 15 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
matlab 非线性最小二乘法问题 已有3人参与
|
||||
|
现有两组散点x、y，想拟合为y=C1.*x.^2/(1+(1-(K+1).*C1.^2.*x.^2).^(1/2))+a*(x.^2)+b*(x.^4)+c*(x.^6)+d*(x.^8)形式， x=[0.2:0.2:10]; y=[0.0004 0.0018 0.0041 0.0074 0.0116 0.0166 0.0226 0.0294 0.0371 0.0456 0.0549 0.0649 0.0757 0.0872 0.0993 0.1120 0.1254 0.1393 0.1537 0.1685 0.1838 0.1995 0.2155 0.2318 0.2484 0.2653 0.2823 0.2994 0.3167 0.3340 0.3514 0.3689 0.3863 0.4037 0.4210 0.4383 0.4555 0.4726 0.4895 0.5063 0.5230 0.5396 0.5560 0.5722 0.5883 0.6042 0.6199 0.6355 0.6508 0.6661]; 运行程序如下： syms C1 K a b c d x=[0.2:0.2:10]; fi=C1.*x.^2/(1+(1-(K+1).*C1.^2.*x.^2).^(1/2))+a*(x.^2)+b*(x.^4)+c*(x.^6)+d*(x.^8) 出现警告： Warning: System is inconsistent. Solution does not exist. fi = [ Inf + a/25 + b/625 + c/15625 + (1511157274518287*d)/590295810358705651712, Inf + (4*a)/25 + (16*b)/625 + (64*c)/15625 + (1511157274518287*d)/2305843009213693952, Inf + (9*a)/25 + (81*b)/625 + (729*c)/15625 + (1210290878676085*d)/72057594037927936, ...... 后面过多，在此不复制了，总之就是多项式的第一项都是Inf，请问为什么会出现这种情况，如何解决？最好提供一下程序，谢谢！ |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有265人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
13楼2016-01-13 10:21:19
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★
感谢参与，应助指数 +1
yujianmiao: 金币+1, ★★★很有帮助, 谢谢帮助，新人金币少，见谅！ 2016-01-11 20:34:26
感谢参与，应助指数 +1
yujianmiao: 金币+1, ★★★很有帮助, 谢谢帮助，新人金币少，见谅！ 2016-01-11 20:34:26
|
如果追求的是效果而不是代码，就用1stopt吧： 相关系数(R): 0.99999998677724 相关系数之平方(R^2): 0.999999973554479 决定系数(DC): 0.999999972874056 卡方系数(Chi-Square): 6.68692837741245E-6 F统计(F-Statistic): 1815052131.43651 参数 最佳估算 ---------- ------------- c1 -0.000413575526215245 k -17351671.8653803 a 0.0118084009534888 b -8.27043394152969E-5 c 4.06674236596116E-7 d -9.21591370383961E-10 |
2楼2016-01-11 16:20:09
3楼2016-01-11 16:37:05
4楼2016-01-11 16:38:47
30