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