ex6_40.m

ex5_1 I-型线性相位FIR滤波器 ex5_2 II-型线性相位FIR滤波器 ex5_3 III-型线性相位FIR滤波器 ex5_4 IV-型线性相位FIR滤波器 ex5_5 矩

Fs = 1000;     % 采样频率
NFFT = 256;
t = 0:1/Fs:1;  　% 采样时间
x = sin(2*pi*100*t) + sin(2*pi*200*t) + sin(2*pi*400*t) + randn(size(t));
w = hanning(NFFT)'; % 256个点的海宁窗，注意取转置
Pxx = (abs(fft(w.*x(1:NFFT))).^2 + ...
       abs(fft(w.*x(NFFT*1/2 + 1:NFFT*3/2))).^2 + ...
       abs(fft(w.*x(NFFT*2/2 + 1:NFFT*4/2))).^2 + ...
       abs(fft(w.*x(NFFT*3/2 + 1:NFFT*5/2))).^2 + ...
       abs(fft(w.*x(NFFT*4/2 + 1:NFFT*6/2))).^2 + ...
       abs(fft(w.*x(NFFT*5/2 + 1:NFFT*7/2))).^2) / (norm(w)^2*6);
f = (0:(NFFT-1))./NFFT*Fs;       % 横坐标变换
PXX = 10 * log10(Pxx);           % 纵坐标变换
figure
plot(f, PXX);
xlabel('频率（Hz）');
ylabel('功率谱（dB）');
grid;