ex6_36.m

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

Fs=2000;
NFFT=1024;
n=0:1/Fs:1;
x=sin(2*pi*100*n)+4*sin(2*pi*500*n)+randn(size(n));  % 噪声序列
Cx=xcorr(x,’unbiased’);        % 计算序列的自相关函数
Cxk=fft(Cx,NFFT);            % 计算FFT
pxx=abs(Cxk);                % 求解PSD
t=0:round(NFFT/2-1);
k=t*Fs/NFFT;
P=10*log(pxx(t+1));           % 纵坐标的单位为dB
plot(k,P);