热工试验技术与数据处理课件

1、热工试验技术与数据处理课件热工试验技术与数据处理课件第三节 模拟滤波模拟滤波器是对模拟信号实行线性滤波的一种线性非时变系统H(j)xa(t)X(j)ha(t)ya(t)Y(j)X(j)Y(j)H(j) 动态特性的描述（1）单位冲激响应（2）传递函数（3）频率特性9/11/20222热能工程系第三节 模拟滤波模拟滤波器是对模拟信号实行线性滤波的一种 一、滤波器的传递函数理想幅频特性模拟滤波器无源滤波 R, C, L有源滤波 R, C, 运算放大器H(j)H0 2H(j)H0 0H(j)H0 0H(j)H0 2 1 1二阶传输函数9/11/20223热能工程系 一、滤波器的传递函数模拟滤波器无源滤

2、波 R, C, L有低通带通高通带阻其中： H0：任意增益因子0：特征频率， 对低通，高通， 0是截止频率 带通，带阻， 0是截止频率q： 选择性因子近似方法巴特沃思近似切比雪夫近似贝塞尔近似只讲归一化低通9/11/20224热能工程系低通带通高通带阻其中： H0：任意增益因子近似方法巴特沃思 通过保角变换得到高通滤波器得到带通滤波器得到带阻滤波器例：低通归一化二阶滤波器9/11/20225热能工程系 通过保角变换得到高通滤波器得到带通滤波器得到带阻滤波器例：1. 巴特沃思近似幅频函数传递函数例： n = 1 s1= -1 Bn(s)=s+1 H(s)=1/(s+1)n=1,2,.n=2巴特沃

3、思多项式9/11/20226热能工程系1. 巴特沃思近似幅频函数传递函数例： n = 1 2. 切比雪夫近似Tn() n阶切比雪夫多项式Vn() 切比雪夫滤波多项式0.5db波纹n=1S+2.863n=2S2+1.425S+1.5161db波纹n=1S+1.965n=2S2+1.098S+1.103切比雪夫多项式Vn() 9/11/20227热能工程系2. 切比雪夫近似Tn() n阶切比雪夫多项式0.53. 贝塞尔近似n=1S+1n=2S2+3S+3En()贝塞尔滤波多项式9/11/20228热能工程系3. 贝塞尔近似n=1S+1En()贝塞尔滤波多项式9/4. 三种近似方法的特点巴特沃思：通

4、带内幅频曲线的幅度平坦，最平幅度逼近，相移与频率的关系不是很线性的，阶跃响应有过冲。切比雪夫：下降最陡，但通带之间幅频曲线有波纹。贝塞尔： 相移和频率之间有良好的线性关系，阶跃响应过冲小，但幅频曲线的下降陡度较差。巴特沃思贝塞尔切比雪夫H()/09/11/20229热能工程系4. 三种近似方法的特点巴特沃思：通带内幅频曲线的幅度平坦二、RC有源滤波器无限增益多路反馈电路vov2Y3viv1Y1Y2Y4Y5设A0 , Ib 根据电流守恒定律，则由此解得9/11/202210热能工程系二、RC有源滤波器无限增益多路反馈电路vov2Y3viv1Y1. 低通滤波器将Y1,Y3,Y4用电阻1/R，Y2,

5、Y5用电容CS，带入上式高频下：C2,C5相当于短路，传输函数为0低频下： C2,C5相当于开路，传输函数-R4/R1viv1R1C2R4C5v2voR3与比较得截止角频率增益因子选择性因子9/11/202211热能工程系1. 低通滤波器将Y1,Y3,Y4用电阻1/R，Y2,Y5用2. 高通滤波器将Y1,Y3,Y4电容，Y2,Y5用电阻，带入上式高频下：C1,C3,C4相当于短路，传输函数为1低频下： C1,C3,C4相当于开路，传输函数= 0viv1C1R2C4R5C3vo与比较得截止角频率增益因子选择性因子9/11/202212热能工程系2. 高通滤波器将Y1,Y3,Y4电容，Y2,Y5用

6、电阻，带三、设计方法1. 选类型，阶数2. 去归一化（所有手册给出的都是归一化传输函数3. 多路反馈实现，求元件值归一化H()0H0H()1H0去归一化推广到复数9/11/202213热能工程系三、设计方法1. 选类型，阶数归一化H()0H0H(例1：二阶低通巴特沃思滤波器，截止频率f0=10Hz去归一化用多路反馈电路来实现若令R=R1=R3=R4=10K则9/11/202214热能工程系例1：二阶低通巴特沃思滤波器，截止频率f0=10Hz去归一化例2：二阶高通切比雪夫滤波器，截止频率f0=10Hz通带波纹3db解：(1) 查表（二阶低通）(2) 保角变换（二阶高通）(3) 去归一化9/11/

7、202215热能工程系例2：二阶高通切比雪夫滤波器，截止频率f0=10Hz通带波 (4) 用多路反馈电路实现若令C=C1=C3=C4=1F四、高阶滤波器的设计奇次阶：一个单极点RC与若干个二阶级联偶次阶：若干个二阶级联9/11/202216热能工程系 (4) 用多路反馈电路实现若令C=C1=C3=C4=1F 练习：设计一个二阶高通滤波器，采用巴特沃思近似，滤波器截止频率为50Hz,采用无限增益多路反馈电路，C1=C3=C4=1F, 求R2,R59/11/202217热能工程系 练习：设计一个二阶高通滤波器，采用巴特沃思近似，滤波器截止第四节 数字滤波器一、数字滤波器的描述数字滤波器：离散时间系

8、统，它接收一个输入序列，对该序列进行某种修正后，作为输出序列送出。H(z)x(n)X(z)h(n)y(n)Y(z)二、数字滤波器的分类（一）根据冲激响应函数的时间特性分无限冲激响应 IIR有限冲激响应 FIR（二）根据滤波的实现方法和型式分递归型非递归型FFT9/11/202218热能工程系第四节 数字滤波器一、数字滤波器的描述H(z)x(n)X(z三、数字滤波器的设计和实现设计：由给定的规定求找滤波器的转移函数实现：按照转移函数获得实际的数字滤波系统设计模拟H(s)数字H(z)差分方程H(z)实现软件硬件软硬结合9/11/202219热能工程系三、数字滤波器的设计和实现设计：由给定的规定求找

9、滤波器的转移四、差分方程和 z 变换（一）差分方程离散系统只能用差分方程描述，对于一个因果系统，用常系数线性差分方程来描述。用途直接得到系统的结构求解系统的瞬态响应9/11/202220热能工程系四、差分方程和 z 变换（一）差分方程用途直接得到系统的结构 例：y(n)=a0 x(n)+a1x(n-1)-b1y(n-1)z-1a0a1-b1y(n)x(n)设：a0=0, a1=1,b1=1/2 y(n)=x(n)+y(n-1)/2 x(n)=(n)= 1 n=0 0 n0设：y(n) = 0 n2 y(-1)=1则：y(0) = x(0)+y(-1)/2 = 1+1/2 = 1.5 y(1)

10、= x(1)+y(0)/2 = 0+1.5/2 = 0.75 y(2) = x(2)+y(1)/2 = 0+0.75/2 = 0.375.n0-1 0 1 2y(n)n9/11/202221热能工程系 例：y(n)=a0 x(n)+a1x(n-1)-b1y(n-(二) Z 变换1. Z 变换的定义双边Z变换Z是一个复变量单边Z变换例1：x(n)=u(n)= 1n0 0n0收敛域ReImZ平面9/11/202222热能工程系(二) Z 变换1. Z 变换的定义双边Z变换Z是一个复变量 对两个不同的序列，表示Z变换的求和表达式可能相同，而收敛域不同。一个序列的Z变换要用它的闭合形式和收敛域来共同描

11、述。如果Z变换的收敛域包括复数平面上的单位圆，则当Z沿单位圆取值时，Z变换就成为傅立叶变换。当实变量T在(-,)间连续变化时，ejT沿Z平面上的单位圆周期性变化，X(ejT)也是周期性变化。从这里可以直观地看出离散信号的频谱X(ejT)随频率T变化的周期性。9/11/202223热能工程系 对两个不同的序列，表示Z变换的求和表达式可能相同，而收敛域2. Z反变换Z反变换的目的在于使我们恢复离散时刻的原始时间函数值，也就是说，能恢复X(nT)若则一般求解反变换有三种常用办法长除法用留数定律求解部分分式展开法9/11/202224热能工程系2. Z反变换Z反变换的目的在于使我们恢复离散时刻的原始时

12、(1)长除法（通过长除法把x(n)表示成Z-1的幂级数根据收敛域判断x(n)为正时间序列，把分子分母按z的降幂排列例：然后进行长除结果由此得 x(0)=1, x(1)=-1, x(2)=4,.9/11/202225热能工程系(1)长除法（通过长除法把x(n)表示成Z-1的幂级数根据收(2)利用留数定律求解若X(z)Zn-1在围线C以内的所有极点集合为Zk，则根据留数定律其中：ResX(z)Zn-1, Zk表示X(z)Zn-1在极点Zk上的留数值(3)部分分式法对有理Z变换9/11/202226热能工程系(2)利用留数定律求解若X(z)Zn-1在围线C以内的所有极例：根据收敛域可以判断x(n)是一个正时间序列，则根据Z变换的线性特性，可得9/11/202227热能工程系例：根据收敛域可以判断x(n)是一个正时间序列，则根据Z变换（三）拉氏变换与Z变换对连续信号进行理想采样这个变换的一个重要性质是它的周期性，周期为2/T9/11/202228热能工程系（三）拉氏变换与Z变换对连续信号进行理想采样这个变换的一个重 对照采样序列x(n)=x(nT)的Z变换当将S表达为直角坐标形式，Z表达成极坐标形式Z的模对应S的实部Z的幅角对应S的虚部9/11/202229热能工程系 对照采样序列x(n)=x(nT)的Z变换当将S表达为直