sy3_fft_hp.m

MATLAB实现什么什么滤波

%滤除cos(2*pi/3*n)
omegap=3*pi/5;Apmax=3;
omegar=3*pi/8;Armin=50;
T=5*pi*10^(-6);
wp=[omegap/pi];wr=[omegar/pi];
[N,wc]=buttord(wp,wr,Apmax,Armin,'s');
[b,a]=butter(N,wc,'high');
sysl=tf(b,a)
w=linspace(0,2*pi,500);
h=freqz(b,a,w);
length(h)
mag=fftshift(20*log10(abs(h)));
n=0:23;
x=1+cos(pi/4*n)+cos(2*pi/3*n);
X=fft(x,500);
X1=fftshift(X);
omega=2*pi/500*([0:499]);
Y=fftshift(X.*h);


title('原信号');
ylabel('幅度');xlabel('频率');
subplot(3,1,1);plot(omega,abs(X1));axis([0,2*pi,0,35]);

title('滤波器幅频响应');
ylabel('Gain(dB)');xlabel('Frequency(rad)');
subplot(3,1,2);plot(w,mag);axis([0,2*pi,-100,5]);

subplot(3,1,3);
plot(omega,abs(Y));axis([0,2*pi,0,25]);