DSP实验报告离散时间信号与系统时频域表示离散傅立叶变换和z变换数字滤波器频域分析和实现数字滤波器设计.doc

南京邮电大学 实 验 报 告 实验名称：离散时间信号与系统的时、频域表示 离散傅立叶变换和z变换 数字滤波器的频域分析和实现 数字滤波器的设计 课程名称 数字信号处理A(双语) 班级学号 姓 名 陈志豪 开课时间 　2015/2016学年，第1学期 实验名称：离散时间信号与系统的时、频域表示 实验目的和任务： 熟悉Matlab基本命令，理解和掌握离散时间信号与系统的时、频域表示及简单应用。在Matlab环境中，按照要求产生序列，对序列进行基本运算；对简单离散时间系统进行仿真，计算线性时不变（LTI）系统的冲激响应和卷积输出；计算和观察序列的离散时间傅立叶变换（DTFT）幅度谱和相位谱。 实验内容： 基本序列产生和运算： Q1.1～1.3，Q1.23，Q1.30～1.33 离散时间系统仿真： Q2.1～2.3 LTI系统：Q2.19，Q2.21，Q2.28 DTFT：Q3.1，Q3.2，Q3.4 实验过程与结果分析： Q1.1运行程序P1.1，以产生单位样本序列u[n]并显示它。 clf; n = -10:20; u = [zeros(1,10) 1 zeros(1,20)]; stem(n,u); xlabel(Time index n); ylabel(Amplitude); title(Unit Sample Sequence); axis([-10 20 0 1.2]); Q1.2 命令clf，axis，title，xlabel和ylabel 命令的作用是什么？ 答：clf命令的作用：清除图形窗口上的图形； axis命令的作用：设置坐标轴的范围和显示方式； title命令的作用：给当前图片命名； xlabel命令的作用：添加x坐标标注； ylabel c命令的作用：添加y坐标标注； Q1.3修改程序P1.1，以产生带有延时11个样本的延迟单位样本序列ud[n]。运行修改的程序并显示产生的序列。 clf; n = -10:20; u = [zeros(1,21) 1 zeros(1,9)]; stem(n,u); xlabel(Time index n); ylabel(Amplitude); title(Unit Sample Sequence); axis([-10 20 0 1.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([0 50 -3 3]); grid; title(Sinusoidal Sequence); xlabel(Time index n); ylabel(Amplitude); axis; 答：周期为：T====22.5。 Q1.30未污染的信号s[n]是什么样的形式？加性噪声d[n]是什么样的形式？ 答：未污染的信号：s[n]=×。 加性噪声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(Desired length of the filter = ); 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(