数字信号处理尹为民16、17-51频率响应系统函数+52全通最小相位课件教学.pptVIP

第5章 离散时间系统分析 § 5.1 离散时间系统的频率响应和系统函数 系统的频域分析引入 一、频率响应 一、频率响应 一、频率响应 一、频率响应 一、频率响应 物理可实现系统的条件 实际低通滤波器特性 二、系统函数 二、系统函数 二、系统函数 二、系统函数 三、频率响应的几何确定法 三、频率响应的几何确定法 三、频率响应的几何确定法 三、频率响应的几何确定法 三、频率响应的几何确定法 四、有理系统函数的单位冲激响应 ——IIR系统和FIR系统 四、有理系统函数的单位冲激响应 ——IIR系统和FIR系统 五、 MATLAB实现 五、 MATLAB实现 § 5.2 全通系统与最小相位系统 一、全通系统 一、全通系统 一、全通系统 一、全通系统 一、全通系统 一、全通系统 一、全通系统 一、全通系统 二、最小相位系统 二、最小相位系统 二、最小相位系统 作业: 5-7; 5-10 馋死 PPT研究院 POWERPOINT ACADEMY 例： 给定系统 求： 频率响应 极－零图 极－零图 频率响应 2、从系统函数看： 若为无限长，称为IIR系统； 若为有限长，则称为FIR系统 存在极点（IIR) 无极点（FIR) 1、从离散时域看： 3、从差分方程看： 递归 非递归 利用MATLAB画出z平面的零极点分布和频率响应 zplane(b,a) freqz(b,a) y = filter(b, a, x); b，a ：H(z)分子和分母多项式系数矩阵 %zeros and poles of the transfer function %frequency response b =[1 2 1]; a =[1 -0.5 -0.005 0.3]; figure(1);zplane(b,a); figure(2);freqz(b,a) [例] 已知一离散因果的LTI系统的系统函数为 画出z平面的零极点分布图及频率响应。 [例] 已知一离散因果的LTI系统的系统函数为 画出z平面的零极点分布图及频率响应。 2、设计简单FIR高通滤波器 [例2]已知输入信号是角频率分别为0.1和0.4的离散正弦信号。 设计简单FIR高通滤波器, 滤除低频分量，保留高频分量。 [解] 假设该DF是一个具有如下形式的长度为3的FIR系统 h[0]=h[2]=a, h[1]=b 系统的频率响应为 解上述方程组得 满足要求的FIR DF的差分方程为 y (n) =x (n) *h (n) =x (n) *(a d (n) +b d (n-1) +a d (n-2)) = a x (n) + b x (n-1)+ a x (n-2) [例2]已知输入信号是角频率分别为0.1和0.4的离散正弦信号。 设计简单FIR高通滤波器, 滤除低频分量，保留高频分量。 为了滤除低频分量，保留高频分量则要求 %信号通过简单FIR DF的响应 % Set up the filter coefficients a=1; b=[-6.76195 13.456335 -6.76195]; % Generate the two sinusoidal sequences N=100; k = 0:N-1; x1 = cos(0.1*k); x2 = cos(0.4*k); %Generate the filter output sequence y = filter(b, a, x1+x2); plot(k,y,r,k,x2,b--,k,x1+x2,k:); axis([0 100 -2 4]); ylabel(Amplitude); xlabel(Time index k); legend(y[k], x2[k], x1[k]+x2[k]); [例2]已知输入信号是角频率分别为0.1和0.4的离散正弦信号。 设计简单FIR高通滤波器, 滤除低频分量，保留高频分量。 0 10 20 30 40 50 60 70 80 90 100 -2 -1 0 1 2 3 4 Amplitude Time index k y[k] x2[k] x1[k]+x2[k] [例2]已知输入信号是角频率分别为0.1和0.4的离散正弦信号。 设计简单FIR高通滤波器, 滤除低频分量，保留高频分量。 [例2]已知输入信号是角频率分别为0.