| 查看: 2483 | 回復(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ò)誤,新手,不知如何解決,還請高手再次相助! ??? 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ā)表 | |
|---|---|---|---|---|
|
[考研] 306求調(diào)劑 +9 | chuanzhu川燭 2026-03-18 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 9/450 |
|
|
[考研] 298求調(diào)劑 +8 | 上岸6666@ 2026-03-20 | 8/400 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
[考研] 求調(diào)劑 +6 | 十三加油 2026-03-21 | 6/300 |
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
[考研]
|
.6lL 2026-03-18 | 8/400 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[論文投稿]
申請回稿延期一個(gè)月,編輯同意了。但系統(tǒng)上的時(shí)間沒變,給編輯又寫郵件了,沒回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|