《最小二乘法與曲線(xiàn)擬合.ppt》由會(huì)員分享,可在線(xiàn)閱讀,更多相關(guān)《最小二乘法與曲線(xiàn)擬合.ppt(38頁(yè)珍藏版)》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。
1、如果已知函數(shù)f(x)在若干點(diǎn)xi(i=1,2,,n)處 的值yi,便可根據(jù)插值原理來(lái)建立插值多項(xiàng)式作為f(x)的近似。但在科學(xué)實(shí)驗(yàn)和生產(chǎn)實(shí)踐中,往往會(huì)遇到這樣一種情況,即節(jié)點(diǎn)上的函數(shù)值并不是很精確的,這些函數(shù)值是由實(shí)驗(yàn)或觀測(cè)得到的數(shù)據(jù),不可避免地帶有測(cè)量誤差,如果要求所得的近似函數(shù)曲線(xiàn)精確無(wú)誤地通過(guò)所有的點(diǎn)(xi,yi),就會(huì)使曲線(xiàn)保留著一切測(cè)試誤差。,,最小二乘法與曲線(xiàn)擬合,為此,我們希望從給定的數(shù)據(jù)(xi,yi)出發(fā),構(gòu)造一個(gè)近似函數(shù) ,不要求函數(shù) 完全通過(guò)所有的數(shù)據(jù)點(diǎn),只要求所得的近似曲線(xiàn)能反映數(shù)據(jù)的基本趨勢(shì), 如圖5-7所示。,,,圖5-7 曲線(xiàn)擬合示意圖,也就是說(shuō)擬合函數(shù) 在xi處的
2、偏差(亦稱(chēng)殘差) 不都嚴(yán)格地等于零。即為矛盾方程組。,,,,,,,曲線(xiàn)擬合函數(shù) 不要求嚴(yán)格地通過(guò)所有數(shù)據(jù)點(diǎn),但是,為了使近似曲線(xiàn)能盡量反映所給數(shù)據(jù)點(diǎn)的 變化趨勢(shì),要求 按某種度量標(biāo)準(zhǔn)最小。若記向量,即要求向量 的某種范數(shù) 最小,如 的1-范數(shù) 或-范數(shù),即,為最小。這種要求誤差(偏差)平方和最小的擬合稱(chēng)為曲線(xiàn)擬合的最小二乘法。,為了便于計(jì)算、分析與應(yīng)用,通常要求 的2-范數(shù),實(shí)質(zhì)仍然是求矛盾方程組的最小二乘解。,作擬合直線(xiàn),,,,,,,,,(1)直線(xiàn)擬合,該直線(xiàn)不是通過(guò)所有的數(shù)據(jù)點(diǎn) , 而是使偏差平方和,設(shè)已知數(shù)據(jù)點(diǎn) , 分布大致為一條直線(xiàn)。,為最小,,其中每組數(shù)據(jù)與擬合曲線(xiàn)
3、的偏差為,根據(jù)最小二乘原理,應(yīng)取 和 使 有極小值,,故 和 應(yīng)滿(mǎn)足下列條件:,解法一:,即得如下正規(guī)方程組,求解該方程組,解得,代人 即得擬合曲線(xiàn)。,也可將條件帶入構(gòu)成矛盾方程組,其中,利用,解法二:,即得如下正規(guī)方程組,求解該方程組,解得,代人 即得擬合曲線(xiàn)。,例:某種合成纖維的強(qiáng)度與其拉伸倍數(shù)有直接關(guān)系,下表是實(shí)際測(cè)定的24個(gè)纖維樣品的強(qiáng)度與相應(yīng)拉伸倍數(shù)的記錄。試確定這種關(guān)系。,,(提示:將拉伸倍數(shù)作為x, 強(qiáng)度作為y, 在座標(biāo)紙上標(biāo)出各點(diǎn),可以發(fā)現(xiàn)什么?),,解:設(shè) y=a+bx,從上圖中可以看出強(qiáng)度與拉伸倍數(shù)大致成線(xiàn)形關(guān)系, 可用一條直線(xiàn)來(lái)表示兩者之間的關(guān)
4、系。,則:,解得: a=0.15 , b=0.859 直線(xiàn)方程為:y=0.15+0.859x,計(jì)算出它的正規(guī)方程得,,1 2 3 4 1.36 1.37 1.95 2.28 14.094 16.844 18.475 20.963,,,,,,,,,用最小二乘法求以上數(shù)據(jù)的擬合函數(shù),例 設(shè)有某實(shí)驗(yàn)數(shù)據(jù)如下:,解:把表中所給數(shù)據(jù)畫(huà)在坐標(biāo)紙上,將會(huì)看到數(shù)據(jù)點(diǎn)的分布可以用一條直線(xiàn)來(lái)近似地描述,,設(shè)所求的擬合直線(xiàn)為,,,,,,,則正規(guī)方程組為,解得,即得擬合直線(xiàn),將以上數(shù)據(jù)代入上式正規(guī)方程組,得,其中,,,(2)多項(xiàng)式擬合 有時(shí)所給數(shù)據(jù)點(diǎn)的分布并不一定近似地呈一條直線(xiàn),這時(shí)仍用直線(xiàn)擬合
5、顯然是不合適的,可用多項(xiàng)式擬合。,,,,對(duì)于給定的一組數(shù)據(jù), 尋求次數(shù)不超過(guò)m (m<
6、拋物線(xiàn),因此設(shè)所求的多項(xiàng)式為,,由法方程組(5.46), n=6, 經(jīng)計(jì)算得,,其法方程組為,,,解之得,所求的多項(xiàng)式為,例1 設(shè)函數(shù)y=f(x)的離散數(shù)據(jù)如下表所示,試用二次多項(xiàng)式擬和上述數(shù)據(jù),解:設(shè),則,由 可得,例:試用最小二乘法求形如 的多項(xiàng)式,使之與下列數(shù)據(jù)擬合。,解:,由題目可知:,由 可得,(3)可化為線(xiàn)性擬合的非線(xiàn)性擬合,1 2 3 4 5 6 0 0.5 1 1.5 2 2.5 2.0 1.0 0.9 0.6 0.4 0.3,,,,,,,,,用最小二乘法求擬合曲線(xiàn),例 設(shè)某實(shí)驗(yàn)數(shù)據(jù)如下:,解:將已給數(shù)據(jù)點(diǎn)描在坐標(biāo)系中下圖所示,,,,,可以看出這些點(diǎn)接近指數(shù)
7、曲線(xiàn),因而可取指數(shù)函數(shù) 作為擬合函數(shù):,,,,對(duì)函數(shù) 兩邊取對(duì)數(shù)得.,令,則就得到線(xiàn)性模型,得,則正規(guī)方程組為,,其中,,,,,,,,,,,將以上數(shù)據(jù)代入上式正規(guī)方程組,得,解得,由 得,于是得到擬合指數(shù)函數(shù)為,由 得,有些非線(xiàn)性擬合曲線(xiàn)可以通過(guò)適當(dāng)?shù)淖兞刻鎿Q轉(zhuǎn)化為線(xiàn)性曲線(xiàn),從而用線(xiàn)性擬合進(jìn)行處理,對(duì)于一個(gè)實(shí)際的曲線(xiàn)擬合問(wèn)題,一般先按觀測(cè)值在直角坐標(biāo)平面上描出散點(diǎn)圖,看一看散點(diǎn)的分布同哪類(lèi)曲線(xiàn)圖形接近,然后選用相接近的曲線(xiàn)擬合方程。再通過(guò)適當(dāng)?shù)淖兞刻鎿Q轉(zhuǎn)化為線(xiàn)性擬合問(wèn)題,按線(xiàn)性擬合解出后再還原為原變量所表示的曲線(xiàn)擬合方程。,下表列舉了幾類(lèi)經(jīng)適當(dāng)變換后化為線(xiàn)性擬合求解的 曲線(xiàn)擬合方程及變換關(guān)系,曲線(xiàn)擬合方程 變換關(guān)系 變換后線(xiàn)性擬合方程,,,,,,,,,,,,,,,,,,,,,,,,,,,