信号处理-基础知识 fft、滤波器

fft 频谱分析实例

fs=600;
N=1024;
n = 0:N-1;
x1=sin(2*pi*5*n/fs);
x2=sin(2*pi*10*n/fs);
x = x1+x2;
X = fft(x,N);
f=fs/N* (0:N/2-1);
figure;
subplot(221);plot(n,x);grid on;title('原信号x(n)');
subplot(222);plot(f,abs(X(1:N/2)));grid on;title('x(n)的幅频谱'); 

滤波器频谱响应 freqz()

 f0 =50; 
 fs=500;
 r =0.9;N=1024;
 n=0:N-1;
 w0=2*pi*f0/fs;
 b=[1 -2*cos(w0) 1];
a=[1 -2*r*cos(w0) r*r]; 
[H,w] = freqz(b,a);
subplot(221);plot(w,abs(H));grid on;title('陷波器的幅频响应');
subplot(222);plot(w,angle(H));grid on;title('陷波器的相频响应');
subplot(223);zplane(b,a);grid on;title('陷波器的零极点图');

陷波器设计
转载：https://blog.csdn.net/limindaihong/article/details/9205211