基于Matlab的信號(hào)平穩(wěn)性檢驗(yàn)系統(tǒng)設(shè)計(jì)
由式(6)產(chǎn)生的替代數(shù)據(jù)s(t)的時(shí)域波形和頻域波形如圖4所示。
由圖3和圖4不難看出,替代數(shù)據(jù)與原數(shù)據(jù)的傅里葉變換幅值相同,但替代數(shù)據(jù)傅里葉變換的相位是隨機(jī)的。
圖5顯示了由式(7)計(jì)算得到的原數(shù)據(jù)和替代數(shù)據(jù)的時(shí)頻分布圖。由圖5中可見(jiàn),原數(shù)據(jù)的時(shí)頻分布圖有明顯的結(jié)構(gòu)性。它表明是非平穩(wěn)的,而替代數(shù)據(jù)的時(shí)頻分布圖的結(jié)構(gòu)性較原數(shù)據(jù)有明顯減弱,表明替代數(shù)據(jù)的平穩(wěn)性增加。
由式(11)計(jì)算得到的0的概率密度函數(shù)f(0)如圖6所示,其中一共計(jì)算了1 000次替代數(shù)據(jù)。
由圖6可見(jiàn),替代數(shù)據(jù)的平穩(wěn)度主要分布在0.02附近。數(shù)據(jù)的平穩(wěn)度落在0~0.04之間可以認(rèn)為是平穩(wěn)的,而在此之外可以認(rèn)為是非平穩(wěn)的。
由式(11)計(jì)算得到的1=0.046。位于上述區(qū)間之外,因此判為非平穩(wěn)。這一結(jié)果也與圖5所示的結(jié)果相吻合。
4結(jié)語(yǔ)
利用替代數(shù)據(jù)法和時(shí)頻透視法,并采用Matlab的GUI開(kāi)發(fā)環(huán)境,設(shè)計(jì)了一個(gè)信號(hào)平穩(wěn)性檢驗(yàn)系統(tǒng)。該系統(tǒng)的用戶界面友好。利用該系統(tǒng)可以觀察信號(hào)及其替代數(shù)據(jù)的頻域和時(shí)頻域波形,檢驗(yàn)信號(hào)的平穩(wěn)性。
本文引用地址:http://butianyuan.cn/article/149834.htm
評(píng)論