实验二 Z 变换和傅立叶变换
一、实验目的
1、掌握Z 变换和反变换和傅立叶变换在MATLAB 中的处理过程
2、应用MATLAB 信号处理工具箱提供的求连续和离散系统频率响应的函数。
二、实验内容
(1)有限长序列的Z 变换和逆Z 变换
原理: 根据已知求出1x 和2x 的z 变换X 1(z)和X 2(z)。根据z 变换的时域卷积定理,只要求出X(z)= X 1(z) X 2(z)即可。这是两个多项式相乘,可用conv 函数来求得序列x(n)。conv 函数只能给出x 数组,n 的起点ns=ns1+ns2=-3,终点nf=nf1+nf2=2。由x,n 即可的出X(z)。
离散序列的卷积计算:两序列的卷积:
数学形式:)()()(11n h n x n y ?=
MATLAB 表述:y=conv(x1,x2)
(2)求Z 多项式分式的逆变换
内容:设系统函数为2115.02.223)(---+--=z z z z W ,输入为(1)中的2x 信号,用z 变换计算输出)(n y ,)]([)(z Y IZT n y =
原理:用留数法求z 的反变换函数residuez ,
MATLAB 描述:[r,p,k]=residuez(B,A)
向量B 和A 分别为系统函数H (z )的分子和分母多项式系数。
...)1()2()()1()()()(...)()2()2()()1()1()(+-+++++=n k n k n u N p N r n u p r n u p r n y n n n δδ k 是当N M ≥时的直接项,也即有限序列,而其余的则是无限序列。有(1)可知)(2z X ,调用[r,p,k]=residuez(B,A),可由B ,A 求出r,p,k,进而求z 的反变换。
(3)离散时间傅立叶变换
内容:离散时间傅立叶变换
取一个周期的正弦信号,作8点采样,求它的连续频谱。然后对该信号进行N 个周期延拓,再求它的连续频谱。把N 无限增大,比较分析结果。 原理:利用了离散时间傅立叶变换子程序dtft.m
function X=dtft(x,w)
X=x*exp(-j*[1:length(x)]’*w);
(4)由离散序列恢复模拟信号
内容:已知两序列],3,2,1[1=x ]1:1[1-=n 及]5,3,4,2[2x ,]1:2[2-=n ,求出1x 与2x 及卷积x 的z 变换。