南京邮电大学实验报告实验名称:离散时间信号与系统的时、频域表示离散傅立叶变换和z变换数字滤波器的频域分析和实现数字滤波器的设计课程名称数字信号处理A(双语)班级学号B13011025姓名陈志豪开课时间2015/2016学年,第1学期实验名称:离散时间信号与系统的时、频域表示实验目的和任务:熟悉Matlab基本命令,理解和掌握离散时间信号与系统的时、频域表示及简单应用。在Matlab环境中,按照要求产生序列,对序列进行基本运算;对简单离散时间系统进行仿真,计算线性时不变(LTI)系统的冲激响应和卷积输出;计算和观察序列的离散时间傅立叶变换(DTFT)幅度谱和相位谱。实验内容:基本序列产生和运算:Q1.1~1.3,Q1.23,Q1.30~1.33离散时间系统仿真:Q2.1~2.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('Timeindexn');ylabel('Amplitude');title('UnitSampleSequence');axis([-102001.2]);Q1.2命令clf,axis,title,xlabel和ylabel命令的作用是什么?答:clf命令的作用:清除图形窗口上的图形;axis命令的作用:设置坐标轴的范围和显示方式;title命令的作用:给当前图片命名;xlabel命令的作用:添加x坐标标注;ylabelc命令的作用:添加y坐标标注;Q1.3修改程序P1.1,以产生带有延时11个样本的延迟单位样本序列ud[n]。运行修改的程序并显示产生的序列。clf;n=-10:20;u=[zeros(1,21)1zeros(1,9)];stem(n,u);xlabel('Timeindexn');ylabel('Amplitude');title('UnitSampleSequence');axis([-102001.2]);Q1.23修改上述程序,以产生长度为50、频率为0.08、振幅为2.5、相移为90度的一个正弦序列并显示它。该序列的周期是多少?n=0:50;f=0.08;phase=90;A=2.5;arg=2*pi*f*n-phase;x=A*cos(arg);clf;stem(n,x);axis([050-33]);grid;title('SinusoidalSequence');xlabel('Timeindexn');ylabel('Amplitude');axis;答:周期为:T=2=f1=0.081=22.5。Q1.30未污染的信号s[n]是什么样的形式?加性噪声d[n]是什么样的形式?答:未污染的信号:s[n]=n2×n0.9。加性噪声d[n]是均匀分布在-04到+0.4之间的随机序列。Q1.31使用语句x=s+d能产生被噪声污染的信号吗?若不能,为什么?答:不能,因为d是列向量,s是行向量。Q1.32信号x1,x2和x3与信号x之间的关系是什么?答:X1是x的延时一个单位,x2和x相等,x3超前于x一个单位。Q1.33legend命令的作用是什么?答:产生图例说明。Q2.1对M=2,运行上述程序,生成输入x[n]=s1[n]+s2[n]的输出信号。输出x[n]的哪个分量被该离散时间系统抑制?答:输入x[n]被该离散时间系统抑制的分量为Signal2的高频分量。Q2.2若线性时不变系统由y[n]=0.5(x[n]+x[n–1])变成y[n]=0.5(x[n]–x[n–1]),对输入x[n]=s1[n]+s2[n]的影响是什么?n=0:100;s1=cos(2*pi*0.05*n);s2=cos(2*pi*0.47*n);x=s1+s2;M=input('Desiredlengthofthefilter=');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('Timeindexn');ylabel('Amplitude');title('Signal#1');subplot(2,2,2);plot(n,s2);axis([0,100,-2,2]);xlabel('Timeindexn');ylabel('Amplitude');title('Signal#2');subplot(2,2,3);plot(n,x);axis([0,100,-2,2]);xlabel('Timeindexn');ylabel('Amplitude');title('InputSignal');subplot(2,2,4);plot(n,y);axis([0,100,-2,2]);xlabel('Timeindexn');ylabel('Amplitude');title('OutputSignal');axis;答:对于输入的影响是-该系统是一个高通滤波器,它通过高频率的输入分量S2,而不是低频的输入分量S1。Q2.3对滤波器长度M和正弦信号s1[n]和s2[n]的频率取其他值,运行程序P2.1,算出结果。M=3,f1=0.1,f2=0.2M=8,f1=0.25,f2=0.5Q2.19运行P2_5,生成式(2.15)所给离散时间系统的冲激响应。Q2.21利用filter命令编写一个MATLAB程序,生成式(2.17)给出的因果线性时不变系统的冲激响应,计算并画出前40个样本。把你的结果和习题Q2.20中得到的结果相比较。clf;N=40;num=[0.9-0.450.350....