| 查看: 2479 | 回復(fù): 28 | ||||
liangyufei0木蟲 (著名寫手)
|
[求助]
求助Matlab曲線擬合參數(shù)求解 已有3人參與
|
|||
|
求各位大俠幫助。。∈指屑ぃ。 問題描述如下: 已知實(shí)驗(yàn)數(shù)據(jù)x=[0.0019 0.0031 0.0049 0.0093 0.0168 0.0330]; y=[0.698 0.609 0.520 0.350 0.223 0.130]; 要用函數(shù)http://muchongimg.xmcimg.com/data/b ... _1456729728_303.png進(jìn)行擬合,求擬合參數(shù)a、b、c。 圖片1.png |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗(yàn): +52 |

木蟲 (小有名氣)
|
x=[0.0019 0.0031 0.0049 0.0093 0.0168 0.0330]; y=[0.698 0.609 0.520 0.350 0.223 0.130]; f=@(c,x) 1-((c(3)*x)^c(2)/(1+(c(3)*x)*c(2)))^c(1); c0=[1 1 1]; [c,resnorm]=lsqcurvefit(f,c0,x,y); plot(x,y,'.-',x,f(c,x),'r:x') legend('原始數(shù)據(jù)','擬合數(shù)據(jù)') |
木蟲 (著名寫手)
|
首先感謝hellosu應(yīng)助! 但程序運(yùn)行出現(xiàn)如下錯(cuò)誤,新手,不知如何解決,還請(qǐng)高手再次相助! ??? Error using ==> mpower Matrix must be square. Error in ==> @(c,x)1-((c(3)*x)^c(2)/(1+(c(3)*x)^c(2)))^c(1) Error in ==> lsqncommon at 152 fuser = funfcn{3}(xargin{:},varargin{:}); Error in ==> lsqcurvefit at 186 [x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ... Caused by: Failure in initial user-supplied objective function evaluation. LSQCURVEFIT cannot continue. |
木蟲 (小有名氣)
|
clear clc x=[0.0019 0.0031 0.0049 0.0093 0.0168 0.0330]; y=[0.698 0.609 0.520 0.350 0.223 0.130]; f=@(c,x) (1-((c(3).*x).^c(2)/(1+(c(3).*x).*c(2))).^c(1)); c0=[1 1 1]; [c,resnorm]=lsqcurvefit(f,c0,x,y); plot(x,y,'.-',x,f(c,x),'r:x') legend('原始數(shù)據(jù)','擬合數(shù)據(jù)') %%%這樣運(yùn)行后還是出錯(cuò),可能你的數(shù)據(jù)點(diǎn)太少了 |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗(yàn): +52 |
|
二樓的方法正確,但是擬合是數(shù)值運(yùn)算,所以在構(gòu)造函數(shù)f時(shí)應(yīng)該采用點(diǎn)除和點(diǎn)乘這類元素運(yùn)算,而不是矩陣運(yùn)算。樓主直接運(yùn)行下面的函數(shù)不會(huì)出錯(cuò) clear clc x=[0.0019 0.0031 0.0049 0.0093 0.0168 0.0330]; y=[0.698 0.609 0.520 0.350 0.223 0.130]; f=@(c,x) (1-((c(3).*x).^c(2)./(1+(c(3).*x).*c(2))).^c(1)); c0=[1 1 1]; [c,resnorm]=lsqcurvefit(f,c0,x,y); plot(x,y,'.-',x,f(c,x),'r:x') legend('原始數(shù)據(jù)','擬合數(shù)據(jù)') |

專家顧問 (正式寫手)
![]() |
專家經(jīng)驗(yàn): +52 |
|
matlab擬合程度的選擇與初始值選擇有很大的關(guān)系。我改下初始值,這個(gè)擬合程度非常好。 a=0.770311398382774 b=1.02743657844741 c=142.278261696316 clear,clc x=[0.0019 0.0031 0.0049 0.0093 0.0168 0.0330]; y=[0.698 0.609 0.520 0.350 0.223 0.130]; f=@(c,x) (1-((c(3).*x).^c(2)./(1+(c(3).*x).*c(2))).^c(1)); c0=[1 1 100]; c=lsqcurvefit(f,c0,x,y); plot(x,y,'.-',x,f(c,x),'r:x') legend('原始數(shù)據(jù)','擬合數(shù)據(jù)') |

版主 (知名作家)

銀蟲 (小有名氣)

銀蟲 (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 5/250 |
|
|
[考研] 085601求調(diào)劑 +3 | Du.11 2026-03-16 | 3/150 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 311求調(diào)劑 +5 | 26研0 2026-03-15 | 5/250 |
|
|
[考研] 312求調(diào)劑 +3 | 陌宸希 2026-03-16 | 4/200 |
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂 2026-03-14 | 5/250 |
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧, 2026-03-12 | 4/200 |
|
|
[基金申請(qǐng)]
有必要更換申報(bào)口嗎
20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
[考研] 材料工程,326分,求調(diào)劑 +6 | KRSLSR 2026-03-10 | 6/300 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|