MATLAB 畫訊號頻譜的子函式

2022-09-08 13:51:13 字數 521 閱讀 6836

輸入訊號序列和取樣率,該子函式可以畫出該訊號的頻譜圖。

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 首先定義...