matlab h(z)变换成h(t),自动控制实验二

实验二 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 变换。