matlab 万能,matlab 万能实用的线性曲线拟合方法
在科學(xué)計(jì)算和工程應(yīng)用中,經(jīng)常會(huì)遇到需要擬合一系列的離散數(shù)據(jù),最近找了很多相關(guān)的文章方法,在這里進(jìn)行總結(jié)一下其中最完整、幾乎能解決所有離散參數(shù)線性擬合的方法
第一步:得到散點(diǎn)數(shù)據(jù)
根據(jù)你的實(shí)際問題得到一系列的散點(diǎn)
例如:
x=[3.2,3.6,3.8,4,4.2,4.8,5,5.4,6.2,6.4,6.6,6.9,7.1]';%加上一撇表示對(duì)矩陣的轉(zhuǎn)置
y=[0.38,0.66,1,0.77,0.5,0.66,0.83,1,0.71,0.71,1,0.87,0.83]';
第二步:確定函數(shù)模型
根據(jù)上述的實(shí)際散點(diǎn)確定應(yīng)該使用什么樣的曲線,或者說是想要模擬的曲線
t=[3.2,3.6,3.8,4,4.2,4.8,5,5.4,6.2,6.4,6.6,6.9,7.1]';
tt=[0.38,0.66,1,0.77,0.5,0.66,0.83,1,0.71,0.71,1,0.87,0.83]';
plot(t,tt,'.');%得到散點(diǎn)圖
matlab中現(xiàn)成的函數(shù):ployfit()
定義:polyfit函數(shù)是matlab中用于進(jìn)行曲線擬合的一個(gè)函數(shù)。其數(shù)學(xué)基礎(chǔ)是最小二乘法曲線擬合原理。曲線擬合:已知離散點(diǎn)上的數(shù)據(jù)集,即已知在點(diǎn)集上的函數(shù)值,構(gòu)造一個(gè)解析函數(shù)(其圖形為一曲線)使在原離散點(diǎn)上盡可能接近給定的值。
調(diào)用方法:polyfit(x,y,n)。用多項(xiàng)式求過已知點(diǎn)的表達(dá)式,其中x為源數(shù)據(jù)點(diǎn)對(duì)應(yīng)的橫坐標(biāo),可為行向量、矩陣,y為源數(shù)據(jù)點(diǎn)對(duì)應(yīng)的縱坐標(biāo),可為行向量、矩陣,n為你要擬合的階數(shù),一階直線擬合,二階拋物線擬合,并非階次越高越好,看擬合情況而定。
matlab polyfit 做出來的值從左到右表示從高次到低次的多項(xiàng)式系數(shù)。
例子:
x = (0: 0.1: 2.5)';
y = erf(x);
p = polyfit(x,y,6)
p =
0.0084 -0.0983 0.4217 -0.7435 0.1471 1.1064 0.0004
則y=0.0084x^6-0.0983x^5+0.4217x^4-0.7435x^3+0.1471x^2+1.1064x+0.0004
完整例子程序:
%二項(xiàng)式擬合
t=[3.2,3.6,3.8,4,4.2,4.5,4.8,5,5.3,5.4,5.6,5.8,6,6.2,6.4,6.6,6.9,7.1];
y=[0.38,0.66,1,0.77,0.5,0.33,0.66,0.83,0.33,1,0.33,0.5,0.33,0.71,0.71,1,0.87,0.83];
plot(t,y,'.');
hold on
p=polyfit(t,y,2)
y1=polyval(p,t);
axis([0,8,0,1]);
plot(t,y1)
%plot是畫圖函數(shù)
%polyval是求值函數(shù)
%polyfit是曲線擬合函數(shù)
%polyfit用于多項(xiàng)式曲線擬合
%p=polyfit(x,y,m)
%其中, x, y為已知數(shù)據(jù)點(diǎn)向量, 分別表示橫,縱坐標(biāo), m為擬合多項(xiàng)式的次數(shù), 結(jié)果返回m次擬合多項(xiàng)式系數(shù), 從高次到低次存放在向量p中.
%y0=polyval(p,x0)
%可求得多項(xiàng)式在x0處的值y0
結(jié)果:另外可以自行加上對(duì)應(yīng)的橫縱坐標(biāo)內(nèi)容,這里就不多說了。
總結(jié)一下matlab非線性擬合散點(diǎn)圖的過程:得到散點(diǎn)數(shù)據(jù)=>確定線性函數(shù)模型=>求解函數(shù)模型的待定系數(shù)=>得到擬合函數(shù)的具體形式=>畫出擬合圖像
總結(jié)
以上是生活随笔為你收集整理的matlab 万能,matlab 万能实用的线性曲线拟合方法的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: linux网卡固件名,修改CentOS7
- 下一篇: 上海城市坐标