数字信号滤波器的重要函数

1）freqz()：通常使用格式 [h,w] = freqz(b,a,n,fs)

b和a分别为系统函数的分子和分母，n为要显示的频率等分点，默认为512，fs为要显示的频率的2倍。

例：如果要显示一个高通滤波器100hz以内系统函数的频响特性，代码：freqz(b,a,512,200);显示512个频率等分点，显示200/2以内的频

率。如下图：

 

2）tf(): 通常使用格式为sys = tf(Numerator,Denominator) 

当不知道分子分母的值，只知道传递函数时，同样可以采用该函数计算出分子分母。

例如：若已知传递函数为 H = ((1-z^(-5))/(1-z^(-1)))^2，对数字信号signal进行滤波处理，代码如下

z = tf('z');

H1 = ((1-z^(-K))/(1-z^(-1)))^M;

b = cell2mat(H1.Numerator);
a = cell2mat(H1.Denominator);

signalFiltered=filter(b,a,signal);

未完待续。。。