24小时热门版块排行榜    

查看: 2399  |  回复: 15
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yujianmiao

新虫 (初入文坛)

[求助] 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,请问为什么会出现这种情况,如何解决?最好提供一下程序,谢谢!
回复此楼

» 猜你喜欢

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

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

xqaiqz

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by 月只蓝 at 2016-01-11 16:58:41
加不加点,结果应该也是2楼那样。MATLAB里涉及数组计算要加点,别的软件没有这样的规则。...

想问这道题的Matlab程序代码
matlab 非线性最小二乘法问题



发自小木虫IOS客户端
9楼2016-01-12 01:48:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +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
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2016-01-11 16:20:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yujianmiao

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 月只蓝 at 2016-01-11 16:20:09
如果追求的是效果而不是代码,就用1stopt吧:

相关系数(R): 0.99999998677724
相关系数之平方(R^2): 0.999999973554479
决定系数(DC): 0.999999972874056
卡方系数(Chi-Square): 6.68692837741245E-6
F统计( ...

谢谢!
3楼2016-01-11 16:37:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yujianmiao

新虫 (初入文坛)

自己不细心,少写了个点哈,程序应该是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)
正在重新计算中,还是希望有大神帮忙看下matlab该如何编这个程序
4楼2016-01-11 16:38:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 11:49 by jCd0dEvKHShX
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:40 by jCd0dEvKHShX
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 11:18 by jCd0dEvKHShX
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 8/400 2026-08-29 11:17 by jCd0dEvKHShX
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +14 晴天加油 2026-08-26 15/750 2026-08-29 11:06 by biily
[基金申请] 中青基了要发朋友圈吗? +3 349506619 2026-08-28 3/150 2026-08-29 10:49 by 木风9012
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 10:18 by winsaint
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 05:45 by gy1nBQXYQJqL
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见