第09章-离散时间滤波器的实现--数字信号处理[刘兴钊][电子教案]幻灯片.pptVIP

第9章 离散时间滤波器的实现 9.1 信号流图表示 9.2 IIR滤波器的基本实现结构1. 直接I型 2. 直接II型（典范型） 3. 级联型 4. 并联型 9.3 FIR滤波器的基本实现结构1.直接型 2. 级联型 3. 线性相位型 4. DFT实现 1.重叠相加法 2. 重叠保留法 有限字长效应简介 FIR和IIR比较 FIR数字滤波器 IIR数字滤波器 单位取样响应有限长 单位取样响应无限长 总是稳定(只有原点处的极点) 要考虑稳定性 采用卷积或递归实现 不能用卷积实现，采用递归 极易实现线性相位 不能实现真正的线性相位 要达到同样的性能指标 由于递归的采用，能以较低阶数比IIR数字滤波器高很多 阶数逼近指标 设计过程需反复迭代 不需迭代 可以采用FFT快速实现 无快速算法 注意，不是互为逆系统 例如典范型比直接I型存储空间少， 线性相位型比直接型乘法少 级联型偏差小，易调整 * 9.1 信号流图表示 9.2 IIR滤波器的基本实现结构 9.3 IIR滤波器的基本实现结构 信号流图是由连接节点的有向支路组成的网络。 带箭头的线段表示一条支路； 支路的两个端点称为节点，一个节点代表一个信号变量，用圆点或圆圈表示； 支路上的箭头表示信号流动的方向； 支路上的增益表示输入节点到输出节点的线性变换，其中常数增益表示乘以常数，增益表示延迟k个样本，增益没有标示则表示传输比为1或恒等变换。 每个节点的变量值等于进入该节点的所有支路的输出变量之和。 没有流进支路的节点又称源节点； 仅有流进支路的节点又称汇节点。 LTI系统三种基本运算的流图表示 举例 流图转置定理：对于单输入/输出系统，将流图中所有支路的方向颠倒，但保持支路增益不变，并将源节点和汇节点互换，则所得流图与原流图具有相同的系统函数。 先实现零点，后实现极点 优点：简单；缺点：延迟多；对字长敏感；调整零极点不便 double x[4],y[3]; while(!eof(in_file)) { for(k=3;k0;k--) //M=3 x[k]=x[k-1]; x[0]=getc(in_file)-128; //8bit wav文件是单极性的 for(k=2;k0;k--) y[k]=y[k-1]; for(k=0,y[0]=0;k=3;k++) y[0]+=b[k]*x[k]; for(k=1;k=2;k++) y[0]+=a[k]*y[k] putc(out_file,y[0]+128); //写文件 } 直接形式I的C语言实现 先实现极点，后实现零点 M=N 优点：延迟少一半为MAX{M，N} double w[4],y; while(!eof(in_file)) { for(k=3;k0;k--) //M=3 w[k]=w[k-1]; w[0]=getc(in_file)-128; //8bit wav文件是单极性的 for(k=1;k=2;k++) w[0]+=a[k]*w[k]; for(k=0,y=0;k=3;k++) y+=b[k]*w[k] putc(out_file,y+128); } 直接形式II的C语言实现 举例 互为共轭的复数零点两两配对，或实数零点任意两两配对, 互为共轭的复数极点两两配对，或实数极点任意两两配对,形成二阶IIR子系统的系统函数 4乘2阶子系统 5乘2阶子系统 优点：对系数量化的灵敏度比直接型低， 可寻找因有限字长效应造成的误差最小的 实数零点极点配对和级联顺序； 易于调整零点和极点； 可用一个二阶环时分复用。 优点：对系数量化的灵敏度比直接型低， 由有限字长效应造成的总误差较小； 易于调整极点； 硬件实现快。 缺点：不能调整零点，不能用于对零点位置精度要求 高的滤波器，例如点阻和窄带带阻。 又称横向结构 double x[4],y; while(!eof(in_file)) { for(k=3;k0;k--) //M=3 x[k]=x[k-1]; x[0]=getc(in_file)-128; for(k=0,y=0;k=3;k++) y+=h[k]*x[k]; putc(out_file,y+128); } 直接形式的C语言实现 优点：易于调整零点； 可用一个二阶环时分复用