輸入訊號序列和取樣率,該子函式可以畫出該訊號的頻譜圖。
function [f,spec,nfft]=spec_fft_plot(sample,l,fs)% 輸入資料說明:
% sample:訊號序列;
% l:訊號序列的長度;
% fs:該訊號的取樣頻率。
% 輸出資料說明:
% f:頻率;
% spec:頻譜圖縱座標;
nfft = 2^nextpow2(l);
% nfft = l;
spec=abs(fft(sample,nfft)./l);
flag = 1;
if flag ==1
f = fs/2*linspace(0,1,nfft/2+1);
spec=spec(1:length(f));
else
f = fs/2*linspace(-1,1,nfft);
spec = fftshift(spec);
end
訊號的頻譜,頻譜密度,能力譜區別
訊號的頻譜,頻譜密度,能力譜區別 訊號可以分成能量訊號與功率訊號,非週期能量訊號具有能量譜密度,是傅利葉變換的平方,功率訊號具有功率譜密度,其與自相關函式是一對傅利葉變換對,等於傅利葉變換的平方 區間長度。不能混淆。能量訊號是沒有功率譜的。隨機訊號在時間上是無限的,在樣本上也是無窮多,因此隨機訊號的...
matlab窗函式的頻譜 如何理解頻譜洩露?
這只是我個人的理解,錯誤之處還請指正。如果原始訊號中有頻率成分處於兩個基本函式的頻率之間,會怎麼樣呢?fig.1 a 解釋了答案,原始訊號包含兩路不同頻率的正弦波,一路頻率與基本函式匹配,一路頻率不匹配。前者用乙個點就能表示峰值,而後者則會出現乙個峰值伴隨著兩個尾巴的頻譜,即發生了譜洩露 拖尾 可以...
數字訊號處理 觀察訊號的頻譜
今天學習將時域訊號通過fft轉換為頻域訊號之後,將其各個頻率分量的幅值繪製成圖,可以很直觀地觀察訊號的頻譜。重點理解fft變換的過程。程式來自參考書 python科學計算 import numpy as np import pylab as pl from pylab import mpl 首先定義...