計(jì)量學(xué)基礎(chǔ)教學(xué)線性參數(shù)最小二乘處理

1、 最小二乘法原理是一種在多學(xué)科領(lǐng)域中獲得廣泛應(yīng)用的數(shù)據(jù)處理方法應(yīng)了解最小二乘法的基本思路和基本原理，以及在等精度或不等精度測(cè)量中線性、非線性參數(shù)的最小二乘估計(jì)方法，并科學(xué)給出估計(jì)精度。一、引入一、引入待測(cè)量（難以直接測(cè)量）：txxx,21直接測(cè)量量：nyyy,21),(),(),(212122221111tnnnttxxxfylxxxfylxxxfyl問(wèn)題：如何根據(jù)和測(cè)量方程解得待測(cè) 量的估計(jì)值？nlll,21txxx,21: tn 直接求得。txxx,21: tn 有利于減小隨機(jī)誤差，方程組有冗余，采用最小二乘原理求 。txxx,21討論：最小二乘原理：最可信賴值應(yīng)使殘余誤差平方和最小。二、

2、最小二乘原理二、最小二乘原理 設(shè)直接測(cè)量量 的估計(jì)值為 ，則有nyyy,21nyyy,21),(),(),(2121222111tnnttxxxfyxxxfyxxxfy由此得測(cè)量數(shù)據(jù) 的殘余誤差nlll,21),(),(),(212122221111tnnnttxxxflvxxxflvxxxflv殘差方程式 若 不存在系統(tǒng)誤差，相互獨(dú)立并服從正態(tài)分布，標(biāo)準(zhǔn)差分別為 ，則 出現(xiàn)在相應(yīng)真值附近 區(qū)域內(nèi)的概率為nlll,21n,21nddd,21), 2 , 1(21)2(22nidepiiiii由概率論可知，各測(cè)量數(shù)據(jù)同時(shí)出現(xiàn)在相應(yīng)區(qū)域的概率為nnniidddeppniii21)2(2111222

3、1nlll,21測(cè)量值 已經(jīng)出現(xiàn)，有理由認(rèn)為這n個(gè)測(cè)量值出現(xiàn)于相應(yīng)區(qū)間的概率p為最大。要使p最大，應(yīng)有nlll,212222222121nn最小由于結(jié)果只是接近真值的估計(jì)值，因此上述條件應(yīng)表示為2222222121nnvvv最小等精度測(cè)量的最小二乘原理：niinvvvv1222221最小 不等精度測(cè)量的最小二乘原理：niiinnvpvpvpvp122222211最小最小二乘原理（其他分布也適用）最小二乘原理（其他分布也適用） 測(cè)量結(jié)果的最可信賴值應(yīng)使殘余誤差平方和（或加權(quán)殘余誤差平方和）最小。三、等精度測(cè)量的線性參數(shù)最小二乘原理三、等精度測(cè)量的線性參數(shù)最小二乘原理線性參數(shù)的測(cè)量方程和相應(yīng)的估計(jì)

4、量為：tntnnnttttxaxaxayxaxaxayxaxaxay22112222121212121111tntnnnttttxaxaxayxaxaxayxaxaxay22112222121212121111殘差方程為)()()(2211222212122121211111tntnnnnttttxaxaxalvxaxaxalvxaxaxalv令ntnnttnnnaaaaaaaaaavvvvxxxxllll212222111211212121則殘差方程的矩陣表達(dá)式為xalv等精度測(cè)量最小二乘原理的矩陣形式：最?。ǎㄗ钚alxalvvtt不等精度測(cè)量最小二乘原理的矩陣形式：22222212

5、21000000000000nnnnpppp權(quán)矩陣最小）（）（最小xalpxalpvvtt四、不等精度測(cè)量的線性參數(shù)最小二乘原理四、不等精度測(cè)量的線性參數(shù)最小二乘原理正規(guī)方程：誤差方程按最小二乘法原理轉(zhuǎn)化得到的有確定解的代數(shù)方程組。一、等精度測(cè)量線性參數(shù)最小二乘處理的正規(guī)方程一、等精度測(cè)量線性參數(shù)最小二乘處理的正規(guī)方程tntnnnnttttxaxaxalvxaxaxalvxaxaxalv2211222212122121211111最小22221nvvv0)(0)(12112nniiniixvxv正規(guī)方程：正規(guī)方程：tniititniiitniiitiniittniitiniiiniiiinii

6、tniitiniiiniiiiniixaaxaaxaalaxaaxaaxaalaxaaxaaxaala12121111122122111212112121111111特點(diǎn)：特點(diǎn)：主對(duì)角線分布著平方項(xiàng)系數(shù)，正數(shù)相對(duì)于主對(duì)角線對(duì)稱分布的各系數(shù)兩兩相等看正規(guī)方程組中第r個(gè)方程：012121111tniitirniiirniiiriniirxaaxaaxaala02211nnrrrvavava則正規(guī)方程可寫成000221122221121221111nntttnnnnvavavavavavavavava0vat即正規(guī)方程的矩陣形式正規(guī)方程的矩陣形式將代入到中，得xalv0vat0xaalattlaxa

7、attaactlaxctlacxt1（待測(cè)量的無(wú)偏估計(jì)）例1 已知銅棒的長(zhǎng)度和溫度之間具有線性關(guān)系：，為。為獲得時(shí)銅棒的長(zhǎng)度和銅的線膨脹系數(shù)，現(xiàn)測(cè)得不同溫度下銅棒的長(zhǎng)度，如下表，求，的最可信賴值。)1 (0tyyt0y0ycti0/mmli/解： 1）列出誤差方程)(00iiitayylv令 為兩個(gè)待估參量，則誤差方程為daycy00,)(dtclviii按照最小二乘的矩陣形式計(jì)算60150140130120110160.200148.200107.200180.200072.200036.2000adcxl則有：0012. 0034. 0034. 013. 11c03654. 097.199

8、91dclacxt那么：cydmmcy000/0000183. 0/97.1999二、不等精度測(cè)量線性參數(shù)最小二乘處理的正規(guī)二、不等精度測(cè)量線性參數(shù)最小二乘處理的正規(guī) 方程方程最小niiivp120)(0)(12112nniiiniiixvpxvp由此可得不等精度測(cè)量線性參數(shù)最小二乘處理的正規(guī)方程：tniititiniiitiniiitiiniititniitiiniiiiniiiiiniiitniitiiniiiiniiiiiniiixaapxaapxaaplapxaapxaapxaaplapxaapxaapxaaplap12121111122122111212112121111111整理得

9、：000222111222221121122121111nntnttnnnnnnvapvapvapvapvapvapvapvapvap即0pvat不等精度的正規(guī)方程不等精度的正規(guī)方程將代入上式，得xalv0xpaaplattplaxpaattpaactplaxctplacxt1（待測(cè)量的無(wú)偏估計(jì)）例2 某測(cè)量過(guò)程有誤差方程式及相應(yīng)的標(biāo)準(zhǔn)差： 08. 0)5(27.1508. 0)4(22.1308. 0)3(81.1006. 0)2(60. 806. 0)(44. 652154214321322121211xxvxxvxxvxxvxxv試求 的最可信賴值。21,xx解：首先確定各式的權(quán)9:9:

10、9:16:161:1:1:1:1:252423222154321ppppp令61514131211127.1522.1381.1060. 844. 621axxxl900000900000900000160000016nnp227. 2186. 4)(121plapaaxxxtt三、非線性參數(shù)最小二乘處理的正規(guī)方程三、非線性參數(shù)最小二乘處理的正規(guī)方程針對(duì)非線性函數(shù)), 2 , 1(),(21nixxxfytii其測(cè)量誤差方程為 ),(),(),(212122221111tnnnttxxxflvxxxflvxxxflv02010,txxxtiiiititixfxfxfxxxfxxxf020210

11、10201021)()()(),(),(令 ，現(xiàn)將函數(shù)在 處展開，則有tttxxxxxx022021101,將上述展開式代入誤差方程，令002201102010)(,)(,)(),(tiitiiiitiiixfaxfaxfaxxxfll則誤差方程轉(zhuǎn)化為線性方程組)()()(2211222212122121211111tntnnnnttttaaalvaaalvaaalv于是可解得 ，進(jìn)而可得 。), 2 , 1(trr), 2 , 1(trxr近似值近似值為獲得函數(shù)的展開式，必須首先確定 02010,txxx1）直接測(cè)量2）通過(guò)部分方程式進(jìn)行計(jì)算：從誤差方程中選取 最簡(jiǎn)單的t個(gè)方程式，如令 ，由

12、此可解得 。0iv02010,txxx四、最小二乘原理與算術(shù)平均值原理的關(guān)系四、最小二乘原理與算術(shù)平均值原理的關(guān)系 為確定一個(gè)被測(cè)量x的估計(jì)值x，對(duì)它進(jìn)行n次直接測(cè)量，得n個(gè)數(shù)據(jù) ，相應(yīng)的權(quán)分別為nlll,21nppp,21，則測(cè)量的誤差方程為xlvxlvxlvnn2211按照最小二乘原理可求得niiniiiplpx11結(jié)論：最小二乘原理與算術(shù)平均值原理是一致的，結(jié)論：最小二乘原理與算術(shù)平均值原理是一致的， 算術(shù)平均值原理是最小二乘原理的特例。算術(shù)平均值原理是最小二乘原理的特例。組合測(cè)量：通過(guò)直接測(cè)量待測(cè)參數(shù)的組合量（一般是 等精度），然后對(duì)這些測(cè)量數(shù)據(jù)進(jìn)行處理， 從而求得待測(cè)參數(shù)的估計(jì)量，求其精度估計(jì)。以檢定三段刻線間距為例，要求檢定刻線a、b、c、d間的距離 。321,xxxabcd1x3x2xabcd1l3l2l4l6l5l直接測(cè)量各組合量，得mmlmmlmmlmmlmmlmml032. 3981. 1016. 2020. 1985. 0015. 1654321首先列出誤差方程)(