数字信号处理A(双语)DSPA实验报告.doc

南京邮电大学实验报告实验名称：1、离散时间信号与系统的时、频域表示2、离散傅立叶变换和z变换3、数字滤波器的频域分析和实现4、数字滤波器的设计课程名称数字信号处理A(双语)班级学号姓名指导老师开课时间2016/2017学年第一学期实验一：离散时间信号与系统的时、频域表示一、实验目的和任务：主要内容：序列的生成、序列DTFT谱计算、简单滤波器的仿真，对应实验指导书Page1－34。实验要求：（1）按照要求产生基本序列和复杂序列，对序列进行运算，求序列的DTFT幅度谱和相位谱。（2）计算简单的离散时间系统的输出，包括简单滤波器的仿真、计算单个系统和互连系统的输出。二、实验内容：具体包括基本序列产生和运算：Q1.1－Q1.3，Q1.23，Q1.30－Q1.33离散时间系统仿真：Q2.1－Q2.3LTI系统：Q2.19、Q2.21，Q2.28DTFT：Q3.1，Q3.2，Q3.4三、实验过程与结果分析：Q1.1运行P1_1产生单位样本序列u[n]的程序与显示的波形如下:clf;n=-10:20;u=[zeros(1,10)1zeros(1,20)];stem(n,u);xlabel('时间序列n');ylabel('振幅');title('单位样本序列');axis([-102001.2]);结果：Q1.2clf清除波形axis设置坐标轴范围，可读比例等title给图形加标题xlabel给x加标注Ylabel给y加标注Q1.3clf;n=-10:20;u=[zeros(1,10)1zeros(1,20)];stem(n+11,u);xlabel('时间序列n');ylabel('振幅');title('单位样本序列');axis([03201.2]);结果：Q1.23n=0:50;f=0.08;phase=pi/2;A=2.5;arg=2*pi*f*n-phase;x=A*cos(arg);clf;stem(n,x);axis([050-33]);grid;title('正弦序列');xlabel('时间序号n');ylabel('振幅');axis;结果：Q1.30s[n]是线性增加伴随着实指数缓慢衰减的图像加性噪声d[n]是均匀分布在-0.4和+0.4之间的随机序列Q1.31不能，因为d是列向量，s是行向量Q1.32x1是x的延时，x2和x相等，x3超前于xQ1.33产生图例说明Q2.1clf;n=0:100;s1=cos(2*pi*0.05*n);s2=cos(2*pi*0.47*n);x=s1+s2;M=input('滤波器所需的长度=');num=ones(1,M);y=filter(num,1,x)/M;subplot(2,2,1);plot(n,s1);axis([0,100,-2,2]);xlabel('时间序列n');ylabel('振幅');title('信号#1');subplot(2,2,2);plot(n,s2);axis([0,100,-2,2]);xlabel('时间序列n');ylabel('振幅');title('信号#2');subplot(2,2,3);plot(n,x);axis([0,100,-2,2]);xlabel('时间序列n');ylabel('振幅');title('输入信号');subplot(2,2,4);plot(n,y);axis([0,100,-2,2]);xlabel('时间序列n');ylabel('振幅');title('输出信号');axis;结果：S[n]被离散时间系统抑制Q2.2n=0:100;s1=cos(2*pi*0.05*n);s2=cos(2*pi*0.47*n);x=s1+s2;M=input('滤波器所需长度=');num=(-1).^[0:M-1];y=filter(num,1,x)/M;clf;subplot(2,2,1);plot(n,s1);axis([0,100,-2,2]);xlabel('时间序号n');ylabel('振幅');title('信号#1');subplot(2,2,2);plot(n,s2);axis([0,100,-2,2]);xlabel('时间序号n');ylabel('振幅');title('信号#2');subplot(2,2,3);plot(n,x);axis([0,100,-2,2]);xlabel('时间序号n');ylabel('振幅');title('输入信号');subplot(2,2,4);plot(n,y);axis([0,100,-2,2]);xlabel('时间序号n');ylabel('振幅');title('输出信号');axis;显示的波形如下:改变LTI系统对输入的影响是,系统现在是