数字陷波器的设计和仿真(Matlab+C)

目录

一、数字陷波器的模型

二、Matlab仿真

1. 示例1

2. 示例2

三、C语言仿真

1. 由系统函数计算差分方程

2. 示例代码

一、数字陷波器的模型

二、Matlab仿真

1. 示例1

clear 
clc
f0=100;%滤掉的100Hz
fs=1000;%大于两倍的信号最高频率
r=0.9; 
w0=2*pi*f0/fs;%转换到数字域b=[1 -2*cos(w0) 1];
a=[1 -2*r*cos(w0) r*r];N=1024;
[H,w]=freqz(b,a,N);
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('陷波器的零极点图');%输入信