自相關估計

estimation of autocorrelation

?? 隨機信號()的相關函數是在時間域內描述隨機過程的重要特征。自相關函數是隨機信號在不同時刻的值之間的依賴性的量度，是一個很有用的統計平均量，其定義為自相關函數

????????[1003-01]?????? (1)?? 式中[·]表示數學期望，表示共軛值,為時間滯后數。

?? 在隨機信號處理中，自相關函數可以用來檢測淹沒在隨機噪聲干擾中的信號，隨機信號的自功率譜等于它的自相關函數的傅里葉變換。因此，通過自相關估計可求得信號的功率譜。

?? 利用計算機計算自相關估值有兩種方法。一種是直接方法，先計算出隨機信號和它的滯后序列的乘積，再取其平均值即得相關函數的估計值。另一種是間接方法，先用快速變換算法計算隨機序列的功率譜密度，再作反變換計算出相關函數。

?? 直接算法??設離散隨機信號序列()是平穩的，其長度為，自相關函數的估值記作()，定義為

?????? [1003-05]??????

???????????????? [1003-10]??????????????????(3)它的均值[1003-06],因而是無偏的估計。顯然，設為有限值，當→∞,則從式(2)可以得到漸近無偏估計計算個滯后數時的自相關估計約需次實數乘加運算。

?? 間接算法??間接方法是利用快速傅里葉變換的方法計算出功率譜密度函數的估值，然后再計算它的傅里葉反變換，即得自相關函數估值。由于采用了快速傅里葉變換算法，計算速度較快。如當=2時，間接算法所需要的運算量約為8[kg2]次實數乘加運算[kg2]因此,兩種方法的速度比是

????????????????[1003-07]如[kg2]＝13，＝0.1＝819，則[1003-08]，即間接算法比直接算法約快8倍。在用間接算法計算相關函數時,需要把隨機信號序列的長度補零擴大到2-1之后再計算其相關函數。

?? 參考書目

何振亞：《數字信號處理的理論與應用》下冊，人民郵電出版社，北京，1983。

J.S.Bendat et al. ,Random Data: Analysis??and Measurement??Procedures，Wiley-Interscience，New York,1971.