| 5 | 1/1 | 返回列表 |
| 查看: 2158 | 回復(fù): 15 | ||||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
yujianmiao新蟲(chóng) (初入文壇)
|
[求助]
matlab 非線性最小二乘法問(wèn)題 已有3人參與
|
|||
|
現(xiàn)有兩組散點(diǎ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]; 運(yùn)行程序如下: 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, ...... 后面過(guò)多,在此不復(fù)制了,總之就是多項(xiàng)式的第一項(xiàng)都是Inf,請(qǐng)問(wèn)為什么會(huì)出現(xiàn)這種情況,如何解決?最好提供一下程序,謝謝! |
新蟲(chóng) (初入文壇)
主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專(zhuān)家經(jīng)驗(yàn): +1059 |
|
如果追求的是效果而不是代碼,就用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)計(jì)(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 |

新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 083000學(xué)碩274求調(diào)劑 +9 | Li李魚(yú) 2026-03-26 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 279求調(diào)劑 +12 | j的立方 2026-03-29 | 12/600 |
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +5 | 朵朵- 2026-03-26 | 7/350 |
|
|
[考研]
|
鐘llll 2026-03-26 | 5/250 |
|
|
[考研] 一志愿南開(kāi)大學(xué)0710生物學(xué)359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 7/350 |
|
|
[考研] 一志愿廈門(mén)大學(xué)材料工程專(zhuān)碩354找調(diào)劑!! +5 | 貝唄鋇鋇 2026-03-30 | 5/250 |
|
|
[考研] 08工科,295,接受跨專(zhuān)業(yè)調(diào)劑 +3 | lmnlzy 2026-03-30 | 3/150 |
|
|
[考研] 材料專(zhuān)碩調(diào)劑 +11 | 椰椰。 2026-03-29 | 11/550 |
|
|
[考研] 085701求調(diào)劑初試286分 +5 | secret0328 2026-03-28 | 5/250 |
|
|
[考研] 283求調(diào)劑(080500) +14 | A child 2026-03-27 | 14/700 |
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
[考研] 317求調(diào)劑 +10 | 蛋黃咸肉粽 2026-03-26 | 10/500 |
|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見(jiàn)α 2026-03-25 | 3/150 |
|
|
[材料工程] 一志愿C9材料與化工專(zhuān)業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
[考研] 266求調(diào)劑 +11 | 陽(yáng)陽(yáng)哇塞 2026-03-27 | 12/600 |
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺(jué)_16 2026-03-24 | 4/200 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|