選帶傅里葉變換的原理大家可以看書。大致的步驟為

移頻 （將選帶的中心頻率移動到零頻）

數字低通濾波器? （防止頻率混疊）

重新采樣? （將采樣的數據再次間隔采樣，間隔的數據取決于分析的帶寬，就是放大倍數）

復FFT （由于經過了移頻，所以數據不是實數了）

頻率調整 （將負半軸的頻率成分移到正半軸）

程序如下：

function [f, y] = zfft(x, fi, fa, fs) % x為采集的數據 % fi為分析的起始頻率 % fa為分析的截止頻率 % fs為采集數據的采樣頻率 % f為輸出的頻率序列 % y為輸出的幅值序列(實數）f0 = (fi + fa) / 2; %中心頻率 N = length(x); %數據長度r = 0:N-1; b = 2*pi*f0.*r ./ fs; x1 = x .* exp(-1j .* b); %移頻bw = fa - fi; B = fir1(32, bw / fs); %濾波 截止頻率為0.5bw x2 = filter(B, 1, x1); c = x2(1:floor(fs/bw):N); %重新采樣 N1 = length(c); f = linspace(fi, fa, N1); y = abs(fft(c)) ./ N1 * 2; y = circshift(y, [0, floor(N1/2)]); %將負半軸的幅值移過來 end

應用實例：

fs = 2048; T = 100; t = 0:1/fs:T; x = 30 * cos(2*pi*110.*t) + 30 * cos(2*pi*111.45.*t) + 25*cos(2*pi*112.3*t) + 48*cos(2*pi*113.8.*t)+50*cos(2*pi*114.5.*t); [f, y] = zfft(x, 109, 115, fs); plot(f, y);
效果：

《新程序員》：云原生和全面數字化實踐50位技術專家共同創作，文字、視頻、音頻交互閱讀

總結

以上是生活随笔為你收集整理的选带傅里叶变换(zoom-fft)的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。