机械工程测试与信号处理-第七章.ppt

机械工程测试?信息?信号分析 主讲：闫静 南京航空航天大学 第七章 数字信号分析- Z变换 第七章 数字信号分析- Z变换 第七章 数字信号分析- Z变换 第七章 数字信号分析- Z变换 第七章 数字信号分析- Z变换 有限长序列x(n) Z变换的ROC 第七章 数字信号分析-数字滤波 右边序列因果序列x(n) Z变换的ROC 右边序列非因果序列x(n) Z变换的ROC 双边序列x(n) Z变换的ROC 第七章 数字信号分析-数字滤波 第七章 数字信号分析-数字滤波 第七章 数字信号分析-数字滤波 第七章 数字信号分析-数字滤波 逆Z变换-----幂级数法 逆Z变换-----幂级数法 逆Z变换-----部分分式法 逆Z变换-----部分分式法 LSI系统的转移函数 系统稳定性判据1： 一个LSI系统是稳定的充要条件是 由此，FIR系统的单位抽样响应h(n)为有限长，所以这一类系统是稳定的。IIR系统的h(n)为无限长，可能是稳定的，也可能是不稳定的。 即h(0)=b (0), h(1)=b (1),…, h(M)=b (M);对nM,h(n)=0。 所以该系统为FIR系统。 FIR系统由于其h(n)为有限长，在输入端不包含输出对输入的反馈，因此总是稳定的。 若a(k),k=1,2,…,N不全为零，那么输入端包含输出对输入的反馈，因此h(n)为无限长，故该系统为IIR系统，存在稳定性问题。 对于一LSI系统，给定其转移函数，我们可以求出差分方程； 反之，给定差分方程，也可求出其转移函数。 差分方程是时域的表示，转移函数是复频域的表示。 离散系统的极零分析 由极零图估计离散系统的频率响应 离散系统的极零分析 在计算机上计算相频特性时，要用到反正切函数ATAN2(HI，HR), HI， HR分别是H(ejw)的虚部和实部。ATAN2规定，在一、二象限的角度为0~∏，而在三、四象限的角度为0~-∏。由此，若一个角度从0变到2∏，但实际得到的结果是0~∏，再由-∏~0，在 处出现了跳变，跳变的幅度为2 ∏，这种现象称为相位的卷绕（wrapping）。图c的两处跳变就是 由相位的卷绕所引起的。 为了得到连续的相频曲线，可在发生2 ∏跳变的以后各处都加上（或减去） 2 ∏，这种做法称为相位的解卷绕（unwrapping）。 图b即是图c解卷绕的结果。 IIR系统的信号流图 用Z变换求解差分方程 第七章 数字信号分析-数字滤波 冲激响应不变法设计IIR数字低通滤波器 第七章 数字信号分析-数字滤波 第七章 数字信号分析-数字滤波 第七章 数字信号分析-数字滤波 对于低通和带通滤波器，当 足够小时，冲激响应不变法可给出较为满意的结果 由于高通、带阻滤波器不是限带的，因此不能用冲激响应不变法实现G(s)到H(z)的转换。 第七章 数字信号分析-数字滤波 第七章 数字信号分析-数字滤波 试设计一低通数字滤波器，要求在通带 内衰减不大于3dB，在阻带 内衰减不小于20dB。给定 第七章 数字信号分析-数字滤波 第七章 数字信号分析-数字滤波 第七章 数字信号分析-数字滤波 用双线性z变换法设计IIR数字低通滤波器 克服冲激响应不变法的缺点，寻找由s平面到z平面的另外的映射关系，这种关系应保证： ①s平面的整个 轴只映射为z平面的单位圆一周； ②若G(s)是稳定的，由G(s)映射得到的H(z)也应该是稳定的； ③这种映射是可逆的，既能由G(s)得到H(z)，也能由H(z)得到G(s)； ④如果G(j0)＝1，那么 也应等于1 用双线性z变换法设计IIR数字低通滤波器 用双线性z变换法设计IIR数字低通滤波器 用双线性z变换法设计IIR数字低通滤波器 用双线性z变换法设计IIR数字低通滤波器 FI = fft2(I); figure;imshow(log(abs(FI)),[]);title(Lena log(abs(FI))); figure;mesh(log(abs(FI)));title(FFT P); figure;mesh(log(abs(fftshift(FI))));title(fftshift P); figure;imshow(log(abs(fftshift(FI))),[]);title(fftshift P); mFilter = ones