最小二乘法擬合圓公式推導(dǎo)及matlab實(shí)現(xiàn).doc
《最小二乘法擬合圓公式推導(dǎo)及matlab實(shí)現(xiàn).doc》由會(huì)員分享,可在線閱讀,更多相關(guān)《最小二乘法擬合圓公式推導(dǎo)及matlab實(shí)現(xiàn).doc(5頁(yè)珍藏版)》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。
2009-01-17|最小二乘法擬合圓公式推導(dǎo)及matlab實(shí)現(xiàn) 最小二乘法(least squares analysis)是一種數(shù)學(xué)優(yōu)化技術(shù),它通過(guò) 最小化誤差的平方和找到一組數(shù)據(jù)的最佳函數(shù)匹配。最小二乘法是用最 簡(jiǎn)的方法求得一些絕對(duì)不可知的真值,而令誤差平方之和為最小。 最 小二乘法通常用于曲線擬合 (least squares fitting) 。 這里有擬合圓曲線 的公式推導(dǎo)過(guò)程和vc實(shí)現(xiàn)。 matlab 實(shí)現(xiàn): function[R,A,B]=irc(x,y,N) %x,y是平面點(diǎn)的坐標(biāo),N是點(diǎn)個(gè)數(shù) %R是擬合半徑,A,B是圓心的平面坐標(biāo) x1=0; x2=0; x3=0; y1=0; y2=0; y3=0; x1y1=0; x1y2=0; x2y1=0; for i=1:N x1=x1+x(i); x2=x2+x(i)*x(i); x3=x3+x(i)*x(i)*x(i); y1=y1+y(i); y2=y2+y(i)*y(i); y3=y3+y(i)*y(i)*y(i); x1y1=x1y1+x(i)*y(i); x1y2=x1y2+x(i)*y(i)*y(i); x2y1=x2y1+x(i)*x(i)*y(i); end C=N*x2-x1*x1; D=N*x1y1-x1*y1; E=N*x3+N*x1y2-(x2+y2)*x1; G=N*y2-y1*y1; H=N*x2y1+N*y3-(x2+y2)*y1; a=(H*D-E*G)/(C*G-D*D); b=(H*C-E*D)/(D*D-G*C); c=-(a*x1+b*y1+x2+y2)/N; A=a/(-2); B=b/(-2); R=sqrt(a*a+b*b-4*c)/2; VC void CViewActionImageTool::LeastSquaresFitting() { if (m_nNum<3) { return; } int i=0; double X1=0; double Y1=0; double X2=0; double Y2=0; double X3=0; double Y3=0; double X1Y1=0; double X1Y2=0; double X2Y1=0; for (i=0;i- 1.請(qǐng)仔細(xì)閱讀文檔,確保文檔完整性,對(duì)于不預(yù)覽、不比對(duì)內(nèi)容而直接下載帶來(lái)的問(wèn)題本站不予受理。
- 2.下載的文檔,不會(huì)出現(xiàn)我們的網(wǎng)址水印。
- 3、該文檔所得收入(下載+內(nèi)容+預(yù)覽)歸上傳者、原創(chuàng)作者;如果您是本文檔原作者,請(qǐng)點(diǎn)此認(rèn)領(lǐng)!既往收益都?xì)w您。
下載文檔到電腦,查找使用更方便
9.9 積分
下載 |
- 配套講稿:
如PPT文件的首頁(yè)顯示word圖標(biāo),表示該P(yáng)PT已包含配套word講稿。雙擊word圖標(biāo)可打開(kāi)word文檔。
- 特殊限制:
部分文檔作品中含有的國(guó)旗、國(guó)徽等圖片,僅作為作品整體效果示例展示,禁止商用。設(shè)計(jì)者僅對(duì)作品中獨(dú)創(chuàng)性部分享有著作權(quán)。
- 關(guān) 鍵 詞:
- 最小二乘法 擬合 公式 推導(dǎo) matlab 實(shí)現(xiàn)
鏈接地址:http://m.italysoccerbets.com/p-6646561.html