乐吧的博客

（注意：正频率对应系数1/2(a-bj)与前文中ak=a+bj不一致，因为前边2（ac-bd）=2(a cos(kw0t)+b sin(kw0t))这一步符号错了，应该为减号。谢谢网友指正。）Matlab代码：clear;fs=32000;t=[0:1/fs:0.1];f0 = 1000;f0_mirror = 1000;a =1;b ...

我们已经知道，IIR 滤波器可以表示为 多个 1/(1-c*z^(-1)) 相乘（这里假定只有 x(n) 参与滤波）, 每个单元我们称之为子滤波器，多个以上公式的相乘，表示在时域里面，依次经过这些子滤波器过滤，因此，只要任何一个子滤波器不稳定，则导致整个滤波器不稳定，换句话说，要保证系统稳定，则需要保证每个子滤波器稳定。下面，我们以一个子滤波器为例来分析：Y(z) = 1 / ( 1 - c*z^(

在af_aneuqalizer.c这个模块中用到这个宏定义FILTER_ORDER，按照字面理解，可能会被认为是滤波器的阶数，实际上是指用了几个双一阶IIR滤波器来级联。

代码如下：%一个最简单的IIR低通滤波器 % H(z) = (1-a)/(1-b*z^(-1))，其中b=1-a %时域表达式 % y(n) = (1-a)*x(n) + b*y(n-1) %下面的代码在[-pi,pi]之间来计算幅度 %通过两种方式计算 % 1）线下计算出来幅度响应为a./((1-2*b*cos(w)+b^2).^0.5) % 2）直接用e^jw来计算，然后用abs求

针对butterworth类型的双二阶滤波器，给出了求解系数的计算公式，并给出了matlab参考代码。

matlab的函数firceqrip，可以实现FIR高通或者低通滤波器。 需要给定的参数包括：以下四个是滤波器本身的参数： N: 滤波器的阶数 dev: the maximum deviation ，指留下的频段的部分的纹波 rip: ripple 指滤除掉的部分的纹波 fo: 截止频率其他辅助参数： ‘high’ : 指明是高通，如果不给这个参数，缺省是低通 ‘passed

前言双二阶滤波器说明补充说明和一些建议三种滤波器类型的计算过程系数如何使用Direct Form I 的计算过程Direct Form II 的计算过程Transposed direct forms附录英文原文第一个链接来源的英文原文第二个链接来源的英文原文前言本文翻译自这个博客： http://www.earlevel.com/main/2003/02/28/biquads

在学习《信号与系统》时，对四个公式的关系，以前一直比较困惑，尤其是在matlab中实现得到的值与真实世界中要表达的频率是什么关系。四个公式为：公式1）连续且周期的时间函数===》傅里叶级数（非周期的）公式2）离散且周期的时间函数===》傅里叶级数（周期的）公式3）连续且非周期的时间函数===》傅里叶变换（非周期的）公式4）离散且非周期的时间函数===》傅里叶变换（周

简单的解释奈奎斯特采样定理中的镜像频谱问题