首页 > 代码库 > 变换域中的离散时间系统

变换域中的离散时间系统

某些实际信号不存在傅里叶变换。就像拉普拉斯变换的引入一样,加了个衰减因子就满足条件了。而从拉普拉斯到z变换,可以理解为连续到离散的映射。

z变换是一个无穷级数,无穷级数就有收敛域的问题。可以理解为收敛域就是傅里叶变换存在的区域。


z变换求反变换的部分分式法有函数可以计算:[r,p,C] = residuez(b,a)

其中b和a为按z-1升幂序列排列的多项式的分子和坟墓的系数向量;

r为各个根的留数向量;p为极点向量。C先不管。

也可以用h = impz(b,a,N)。这个之前有介绍过,就是已知多项式分子分母求h(n)的。也就是说可以来求反变换。


至于求解差分方程,之前介绍过filter(b,a,x,xic),xic是初始条件输入序列。

其中初始条件计算:xic = filtic(b,a,Y,X)

b和a是分子分母系数数组。Y和X是初始条件数组,Y=[y(-1),y(-2),...],X=[x(-1),x(-2)...]。


接下来讲讲z平面上的谱分析。

之前学过DTFT的几何画法。可以发现,如果极点靠单位圆很近,频率特性在靠近极点附近会出现大的谐振峰,分母迅速减小。

由于稳定性要求,极点要在单位圆内,这样阐释的都是负相移。当零点也在单位圆内,系统的负相移最小(零点可产生正相移抵消),称最小相位系统。


非单位圆周上的频谱分析。

例如语音信号处理中,常常需要知道极点所对应的频率。如果极点里单位圆较远,则单位圆上的频谱就很平滑。

如果使采样点轨迹沿一条接近这些极点的弧线或圆周进行,则采样结果会在极点对应的频率上出现明显的尖峰。


关于理想滤波器,其脉冲响应是sa函数。为了因果,只能截取n>=0部分。考虑到线性相位要求,截取的序列必须对称。

为了使更接近于理想情况,应该尽可能增加延迟时间,加大截取长度(阶数)。

截取的序列越短,幅频特性与理想情况差别越大。

截取的序列若是对称的,则相频为线性。若不对称,相频特性则非线性。


用零极点分析滤波器。

规律是:离零点越近的频率,幅度越小。离极点越近的频率,幅度越大。

由z = eiw,z=-1离低频最远。因此取零点z=-1可以得到更高的低频幅度。

z=-1后,对一阶低通滤波器,通带宽度与极点a的关系近似是wp = 1-a。注意wp是数字频率。

二阶则更加灵活。为了滤波或者陷波,可以直接把零点配置在这个角频率的单位圆上ejw0。

同理,梳状滤波器就是把零点均匀分布在单位圆上。极点位置很靠近零点位置,能将陷波特性做的很窄。

不过陷波器相频特性不好,一般要级联全通滤波器进行校正。

变换域中的离散时间系统