yujianmiao

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 12
帖子: 14
在線: 4.7小時
蟲號: 1749967
注冊: 2012-04-11
專業(yè): 光學(xué)

[求助] matlab 非線性最小二乘法問題已有3人參與

現(xiàn)有兩組散點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)

出現(xiàn)警告：
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, ......

后面過多，在此不復(fù)制了，總之就是多項式的第一項都是Inf，請問為什么會出現(xiàn)這種情況，如何解決？最好提供一下程序，謝謝！

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有221人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有72人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

求matlab非線性規(guī)劃程序已經(jīng)有4人回復(fù)
matlab非線性回歸問題已經(jīng)有8人回復(fù)
matlab非線性擬合模型參數(shù) 已經(jīng)有3人回復(fù)
matlab 非線性擬合數(shù)據(jù)擬合已經(jīng)有3人回復(fù)
MATLAB非線性擬合已經(jīng)有0人回復(fù)
matlab非線性擬合求助已經(jīng)有10人回復(fù)
matlab 非線性擬合已經(jīng)有16人回復(fù)
Matlab非線性規(guī)劃編程求助已經(jīng)有2人回復(fù)
matlab非線性最小二乘法求解已經(jīng)有12人回復(fù)
matlab 非線性微分方程求解已經(jīng)有3人回復(fù)
matlab非線性擬合求助已經(jīng)有7人回復(fù)
matlab 非線性模擬已經(jīng)有1人回復(fù)
matlab非線性擬合的問題已經(jīng)有2人回復(fù)
matlab非線性擬合計算求助已經(jīng)有10人回復(fù)
matlab非線性擬合已經(jīng)有3人回復(fù)
matlab 非線性參數(shù)擬合問題已經(jīng)有1人回復(fù)
matlab 非線性擬合已經(jīng)有2人回復(fù)
MATLAB 非線性擬合 fminsearch中包含非線性方程的求解問題已經(jīng)有2人回復(fù)
MATLAB 非線性擬合已經(jīng)有2人回復(fù)
跪求matlab 非線性回歸程序已經(jīng)有5人回復(fù)
matlab非線性參數(shù)擬合問題已經(jīng)有7人回復(fù)
matlab 非線性方程求解函數(shù) 已經(jīng)有2人回復(fù)
求解 matlab非線性約束編程已經(jīng)有2人回復(fù)

1樓 2016-01-11 15:49:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
yujianmiao: 金幣+1, ★★★很有幫助, 謝謝幫助，新人金幣少，見諒！ 2016-01-11 20:34:26

如果追求的是效果而不是代碼，就用1stopt吧：

相關(guān)系數(shù)(R): 0.99999998677724
相關(guān)系數(shù)之平方(R^2): 0.999999973554479
決定系數(shù)(DC): 0.999999972874056
卡方系數(shù)(Chi-Square): 6.68692837741245E-6
F統(tǒng)計(F-Statistic): 1815052131.43651

參數(shù) 最佳估算
---------- -------------
c1          -0.000413575526215245
k       -17351671.8653803
a       0.0118084009534888
b       -8.27043394152969E-5
c       4.06674236596116E-7
d       -9.21591370383961E-10

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2016-01-11 16:20:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yujianmiao

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 12
帖子: 14
在線: 4.7小時
蟲號: 1749967
注冊: 2012-04-11
專業(yè): 光學(xué)

引用回帖:

2樓: Originally posted by 月只藍 at 2016-01-11 16:20:09
如果追求的是效果而不是代碼，就用1stopt吧：

相關(guān)系數(shù)(R): 0.99999998677724
相關(guān)系數(shù)之平方(R^2): 0.999999973554479
決定系數(shù)(DC): 0.999999972874056
卡方系數(shù)(Chi-Square): 6.68692837741245E-6
F統(tǒng)計( ...

謝謝！

回復(fù)此樓

3樓2016-01-11 16:37:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yujianmiao

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 12
帖子: 14
在線: 4.7小時
蟲號: 1749967
注冊: 2012-04-11
專業(yè): 光學(xué)

自己不細(xì)心，少寫了個點哈，程序應(yīng)該是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該如何編這個程序

贊一下

回復(fù)此樓

4樓2016-01-11 16:38:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

引用回帖:

4樓: Originally posted by yujianmiao at 2016-01-11 16:38:47
自己不細(xì)心，少寫了個點哈，程序應(yīng)該是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該如何編這個程序

加不加點，結(jié)果應(yīng)該也是2樓那樣。MATLAB里涉及數(shù)組計算要加點，別的軟件沒有這樣的規(guī)則。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

5樓2016-01-11 16:58:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yujianmiao

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 12
帖子: 14
在線: 4.7小時
蟲號: 1749967
注冊: 2012-04-11
專業(yè): 光學(xué)

引用回帖:

5樓: Originally posted by 月只藍 at 2016-01-11 16:58:41
加不加點，結(jié)果應(yīng)該也是2樓那樣。MATLAB里涉及數(shù)組計算要加點，別的軟件沒有這樣的規(guī)則。...

是的，但是因為那個軟件自己沒用過，還是想弄個程序以后自己改數(shù)算也方便
那可否幫我計算一下幫擬合函數(shù)改成
fi=C1.*x.^2/(1+(1-(K+1).*C1.^2.*x.^2).^(1/2))+b*(x.^4)+c*(x.^6)+d*(x.^8)后C1 K b c d 的值？

贊一下

回復(fù)此樓

6樓2016-01-11 19:01:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yujianmiao

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 12
帖子: 14
在線: 4.7小時
蟲號: 1749967
注冊: 2012-04-11
專業(yè): 光學(xué)

引用回帖:

下載了那個軟件但是不太會用，能幫忙看下代碼哪兒有問題嗎，無法執(zhí)行，報錯說變量不能超過4個？
Title "NLReg Demo - 1";
Parameter c1 k a b c d;
Variable x, y;
Function 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);
Data;

贊一下

回復(fù)此樓

7樓2016-01-11 20:21:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xqaiqz

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 186.5
帖子: 9
在線: 7.2小時
蟲號: 4103032
注冊: 2015-09-25
專業(yè): 地理學(xué)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

大神可不可以幫忙寫下這道題的程序

發(fā)自小木蟲IOS客戶端

贊一下

回復(fù)此樓

8樓2016-01-12 01:39:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xqaiqz

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 186.5
帖子: 9
在線: 7.2小時
蟲號: 4103032
注冊: 2015-09-25
專業(yè): 地理學(xué)

引用回帖:

想問這道題的Matlab程序代碼

發(fā)自小木蟲IOS客戶端

贊一下

回復(fù)此樓

9樓2016-01-12 01:48:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

【答案】應(yīng)助回帖

引用回帖:

9樓: Originally posted by xqaiqz at 2016-01-12 01:48:17
想問這道題的Matlab程序代碼

...

函數(shù)表達式: (66505-65486*x)^2+(66505-66847*(1-x))^2
目標(biāo)函數(shù)值(最小): 2234258262.01286
x: 0.499947390458923

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

10樓2016-01-12 09:27:23

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 280分求調(diào)劑一志愿085802 +4	PUMPT 2026-03-22	7/350	2026-03-22 22:13 by 星空星月
[考研] 311求調(diào)劑 +3	26研0 2026-03-20	3/150	2026-03-22 14:46 by ColorlessPI
[考研] 291 求調(diào)劑 +3	化工2026屆畢業(yè)?/a> 2026-03-21	3/150	2026-03-22 14:26 by ColorlessPI
[考研] 285求調(diào)劑 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +4	lemonzzn 2026-03-17	8/400	2026-03-21 22:49 by lemonzzn
[考研] 一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑 +3	Grand777 2026-03-21	3/150	2026-03-21 19:23 by 簡之-
[考研] 求調(diào)劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學(xué)員8dgXkO
[考研] 材料學(xué)碩333求調(diào)劑 +3	北道巷 2026-03-18	3/150	2026-03-21 18:17 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +6	Mqqqqqq 2026-03-19	6/300	2026-03-21 08:04 by JourneyLucky
[考研] 22408 344分求調(diào)劑一志愿華電計算機技術(shù) +4	solanXXX 2026-03-20	4/200	2026-03-20 23:49 by alg094825
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3	niko- 2026-03-19	6/300	2026-03-20 21:58 by luoyongfeng
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11	adaie 2026-03-19	11/550	2026-03-20 21:30 by laoshidan
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考研] 08工學(xué)調(diào)劑 +5	用戶573181 2026-03-20	5/250	2026-03-20 15:47 by xia_2003
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨