数字信号处理a(双语)dspa实验报告资料.docVIP

南京邮电大学 实 验 报 告 实验名称：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.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); title(单位样本序列); axis([-10 20 0 1.2]); 结果： Q1.2 clf 清除波形 axis 设置坐标轴范围，可读比例等 title 给图形加标题 xlabel给x加标注 Ylabel给y加标注 Q1.3 clf; n= -10:20; u = [zeros(1,10) 1 zeros(1,20)]; stem(n+11,u); xlabel(时间序列 n);ylabel(振幅); title(单位样本序列); axis([0 32 0 1.2]); 结果： Q1.23 n = 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([0 50 -3 3]); grid; title(正弦序列); xlabel(时间序号n); ylabel(振幅); axis; 结果： Q1.30 s[n]是线性增加伴随着实指数缓慢衰减的图像加性噪声d[n]是均匀分布在-0.4和+0.4之间的随机序列 Q1.31 不能，因为d是列向量，s是行向量 Q1.32 x1是x的延时，x2和x相等，x3超前于x Q1.33 产生图例说明 Q2.1 clf; 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.2 n = 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); plo