陈数字信号实验嘉俊.docVIP

广东海洋大学学生实验报告书 实验名称 离散时间信号和离散时间系统 课程名称 数字信号处理 成绩 学院（系） 信息学院 专业 电气工程及其自动化 班级 电气1093 学生姓名 陈嘉俊 学号 200911632203 实验地点 微机2室 实验日期 一、实验目的 1.掌握计算线性时不变系统的冲激响应的方法。 2.理解时域采样的概念及方法。 3.掌握离散时间信号的z变换和z逆变换分析 4.了解离散时间傅里叶变换(DTFT) 二、实验原理 （一）信号采样 采样就是利用周期性抽样脉冲序列pT(t),从连续信号xa(t)中抽取一系列的离散值，得到抽样信号（或称抽样数据信号）即离散时间信号。 （二）线性时不变离散时间系统 线性系统：满足线性叠加原理的系统。若y1(n)和y2(n)分别是输入序列x1(n)和x2(n)的响应，则输入x(n)=ax1(n)+bx2(n)的输出响应为y(n)=ay1(n)+by2(n)。 时不变系统：即系统参数不随时间变化的系统，亦即系统对于输入信号的响应与信号加于系统的时间无关。即满足 ：若y(n)是x(n)的响应,则y(n-m)是输入x(n-m)的响应，其中m是任意整数。 数字滤波器对单位样本序列的响应称为冲激响应，用h(n)表示。线性时不变离散系统对输入信号x(n)的响应y(n)可用h(n)来表示：。 （三）z变换和逆z变换 序列的z变换定义为： 其中，z是复变量。相应地，单边z变换定义为： MATLAB提供了计算离散时间信号单边z变换的函数ztrans和z反变换函数iztrans： Z=ztrans(x)，x=iztrans(z)。上式中的x和Z分别为时域表达式和z域表达式的符号表示，可通过sym函数来定义。 如果信号的z域表示式是有理函数，进行z反变换的另一个方法是对进行部分分式展开，然后求各简单分式的z反变换。设的有理分式表示为 （4-3） MATLAB信号处理工具箱提供了一个对进行部分分式展开的函数residuez，其语句格式为 [R,P,K]=residuez(B,A) 其中，B，A分别表示X(z)的分子与分母多项式的系数向量；R为部分分式的系数向量；P为极点向量；K为多项式的系数。若X(z)为有理真分式，则K为零。 离散时间傅里叶变换(DTFT) 1.序列x(n)的离散时间傅里叶变换定义为： 是变量的连续函数。 并可写为实部和虚部相加的形式： 也可以表示为：。其中，。称为幅度函数，称为相位函数，又分别称为幅度谱和相位谱，都是的实函数。 2.的离散时间傅里叶逆变换为： 3.由于是连续函数，而在MATLAB中数据只能以向量的形式存在，所以只能在一个给定L个离散频率点的离散频率集合中计算，需要尽可能大地选取L的值以表示连续函数 三、实验内容 （一）线性时不变系统的冲激响应的计算 设系统为y(n)-0.5y(n-1)+0.75y(n-2)=2.5x(n)+2.5x(n-1)+2x(n-2)，计算上述系统的冲激响应。 参考程序如下： N=40; num=[2.5 2.5 2]; den=[1 -0.5 0.75]; y=impz(num,den,N); %画出冲激响应 stem(y); xlabel(时间序号n); ylabel(振幅); title(冲激响应); grid; （二）时域采样 对连续正弦时间信号x(t)=cos(2πft)进行采样，其中f=13。 t=0:0.0005:1; f=13; xa=cos(2*pi*f*t); subplot(2,1,1) plot(t,xa);grid xlabel(时间，msec); ylabel(振幅); title(连续时间信号); axis([0 1 -1.2 1.2]) subplot(2,1,2); T=0.1; n=0:T:1; xs=cos(2*pi*f*n); k=0:length(n)-1; stem(k,xs);grid xlabel(时间序号n); ylabel(振幅); title(离散时间信号); axis([0 length(n)-1 -1.2 1.2]) （三）z变换和z反变换 1.用ztrans函数求函数的z变换。 MATLAB参考程序如下： x=sym(a^n*cos(pi*n)); Z1=ztrans(x); Z=simplify(Z1) Z =z/(z+a) 2.用iztrans函数求函数的z反变换。 MATLAB参考程序如下： Z=sym(z*(2*z^2-11*z+12)/(z-1)/(z-2)^3); x=iztrans(Z); simplify(x) X（n）=-3+3*2^n-1/4*2^n*n-1/4*2^n