功率譜估計淺談(編輯修改稿)
2025-07-16 13:07 本頁面
【文章內(nèi)容簡介】 譜的主瓣能量泄漏到旁瓣中，導(dǎo)致弱信號的主瓣被強信號的旁瓣所湮沒，造成譜的模糊。 3) 平均周期圖法的收斂性較好，曲線平滑，但是功率譜主瓣較寬，分辨率低。這是由于對隨機序列的分段處理引起了長度有限所帶來的 Gibbs 現(xiàn)象而造成的。 參考文獻1. 劉志剛,李錄明,趙冬梅. 現(xiàn)代譜估計法及應(yīng)用效果[J]. 石油地球物理勘探,2009,S1:59+167+9.2. 樊劍,劉鐵,胡亮. 基于現(xiàn)代時頻分析技術(shù)的地震波時變譜估計[J]. 振動與沖擊,2007,11:7982+98+184.3. 蔡希玲,劉學(xué)偉,呂英梅,曹錫娜. 統(tǒng)計FX譜估計方法及應(yīng)用[J]. 勘探地球物理進展,2008,03:181186+163.4. 李剛. 寬帶信號空間譜估計算法研究[D].哈爾濱工程大學(xué),2007.5. 姚武川,姚天任. 經(jīng)典譜估計方法的MATLAB分析[J]. 華中理工大學(xué)學(xué)報,2000,04:4547.6. 王鳳瑛,張麗麗. 功率譜估計及其MATLAB仿真[J]. 微計算機信息,2006,31:287289.7. 王福杰,潘宏俠. MATLAB中幾種功率譜估計函數(shù)的比較分析與選擇[J]. 電子產(chǎn)品可靠性與環(huán)境試驗,2009,06:2831.附錄：譜估計的Matlab實現(xiàn)程序1：經(jīng)典法譜估計clf。Fs=1000。N=256。Nfft=256。%數(shù)據(jù)的長度和FFT所用的數(shù)據(jù)長度n=0:N1。t=n/Fs。%采用的時間序列xn=sin(2*pi*50*t)+2*sin(2*pi*120*t)+randn(1,N)。Pxx=10*log10(abs(fft(xn,Nfft).^2)/N)。%Fourier振幅譜平方的平均值，并轉(zhuǎn)化為dBf=(0:length(Pxx)1)*Fs/length(Pxx)。%給出頻率序列subplot(2,1,1),plot(f,Pxx)。%繪制功率譜曲線xlabel(39。頻率/Hz39。)。ylabel(39。功率譜/dB39。)。title(39。周期圖 N=25639。)。grid on。程序2：經(jīng)典譜加窗分析 Fs=1000。%采樣頻率n=0:1/Fs:1。%產(chǎn)生含有噪音的信號xn=sin(2*pi*50*n)+2*sin(2*pi*120*n)+randn(size(n))。plot(n,xn)。window=boxcar(length(xn))。%矩形窗nfft=512。[Pxx,f]=periodogram(xn,window,nfft,Fs)。%直接法figure(1)plot(f,10*log10(Pxx))。window=boxcar(length(xn))。%矩形窗nfft=1024。[Pxx,f]=periodogram(xn,window,nfft,Fs)。%直接法figure(2)plot(f,10*log10(Pxx))。 程序3：Bartlett平均周期圖法fs=1000。n=0:1/fs:1。xn=sin(2*pi
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1