基于比時(shí)法的晶振頻率測(cè)量建模與分析

　　用最小二乘法們來(lái)估計(jì)參數(shù)β0，β。設(shè)b0，b分別是參數(shù)β0，β的最小二乘估計(jì)，于是得到一元線性回歸的回歸方程：

　　式中：b0，b是回歸方程的回歸系數(shù)，分別表示晶振相對(duì)于標(biāo)準(zhǔn)頻率的初始誤差和累積誤差。應(yīng)用最小二乘法可求得回歸系數(shù)b，b0為：

　　3 數(shù)據(jù)分析與處理

　　3.1 回歸系數(shù)估計(jì)

　　為了定量分析數(shù)據(jù)，從而確定晶振頻率誤差的組成，首先對(duì)上述測(cè)得的數(shù)據(jù)進(jìn)行歸一化處理。實(shí)際測(cè)量中得到的是晶振脈沖的計(jì)數(shù)個(gè)數(shù)，設(shè)測(cè)量系統(tǒng)所用晶振頻率為10 MHz，可將計(jì)數(shù)數(shù)據(jù)轉(zhuǎn)化為晶振相對(duì)于標(biāo)準(zhǔn)時(shí)間每秒的時(shí)間之差。例如在x=30 s時(shí)，y=2 349，表示在第30 s時(shí)，晶振頻率相對(duì)于標(biāo)準(zhǔn)頻率的計(jì)數(shù)值為2 349，若晶振頻率f=10 MHz，則可得到在第30 s時(shí)晶振相對(duì)于GPS時(shí)間的誤差為t=y／f=234.9μs。

　　用Matlab對(duì)歸一化數(shù)據(jù)進(jìn)行處理，依照最小二乘原理，得到計(jì)數(shù)時(shí)間x與時(shí)間差值y的均值，以及x的自相關(guān)、x和y的互相關(guān)、y的自相關(guān)及回歸方程如表2所示。

　　以上求得了回歸方程，但是該方程是否基本上符合y與x之間的客觀規(guī)律，是否符合晶振頻率誤差變化的實(shí)際特點(diǎn)，還需要對(duì)回歸方程做進(jìn)一步的分析。在回歸分析法中，通常采用方差分析法們對(duì)回歸方程的顯著性進(jìn)行檢驗(yàn)，其實(shí)質(zhì)是將N個(gè)測(cè)量值的影響從數(shù)量上區(qū)分開(kāi)，然后用F檢驗(yàn)法對(duì)所求回歸方程進(jìn)行顯著性檢驗(yàn)。