| 查看: 2411 | 回复: 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论文润色/翻译怎么收费?
已经有189人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求matlab非线性规划程序
已经有4人回复
- matlab非线性回归问题
已经有8人回复
- matlab非线性拟合模型参数
已经有3人回复
- matlab 非线性拟合 数据拟合
已经有3人回复
- MATLAB非线性拟合
已经有0人回复
- matlab非线性拟合求助
已经有10人回复
- matlab 非线性拟合
已经有16人回复
- Matlab非线性规划编程求助
已经有2人回复
- matlab非线性最小二乘法求解
已经有12人回复
- matlab 非线性微分方程求解
已经有3人回复
- matlab非线性拟合求助
已经有7人回复
- matlab 非线性模拟
已经有1人回复
- matlab非线性拟合的问题
已经有2人回复
- matlab非线性拟合计算求助
已经有10人回复
- matlab非线性拟合
已经有3人回复
- matlab 非线性参数拟合问题
已经有1人回复
- matlab 非线性拟合
已经有2人回复
- MATLAB 非线性拟合 fminsearch中包含非线性方程的求解问题
已经有2人回复
- MATLAB 非线性拟合
已经有2人回复
- 跪求matlab 非线性回归程序
已经有5人回复
- matlab非线性参数拟合问题
已经有7人回复
- matlab 非线性方程求解函数
已经有2人回复
- 求解 matlab非线性约束编程
已经有2人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
10楼2016-01-12 09:27:23
月只蓝
主管区长 (职业作家)
-
专家经验: +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
20