数字信号实验报告.docVIP

数字信号处理实验报告 专 业： 姓 名： 学 号： 指导教师： 实验一 序列、频谱、卷积 实验目的 掌握序列的输入方法； 熟悉不同序列的特征； 了解确定性信号谱分析的方法； 验证卷积的计算过程； 实验要求 利用matlab程序，生成几种常用的序列，如矩形序列，单位脉冲序列； 绘制图形，观察序列特征； 研究其频率特性，绘制图形，观察频率响应特征； 利用matlab程序，验证卷积的过程； 实验步骤 矩形序列 （1）生成长度为N的矩形序列， 观察并记录生成的图形； n=1:50 x=sign(sign(10-n)+1); close all; subplot(3,1,1);stem(x);title(单位矩形信号序列); （2）研究其频率特性，，分别研究其幅频特性和相频特性，观察并记录生成的图形； k=-25:25; X=x*(exp(-j*pi/25)).^(n*k); magX=abs(X); %绘制x(n)的幅度谱 subplot(3,1,2);stem(magX);title(单位矩形信号的幅度谱); angX=angle(X); %绘制x(n)的相位谱 subplot(3,1,3);stem(angX) ; title (单位矩形信号的相位谱) 答：根据以上两个源程序得到的图形如下图1.1： 图1.1 分析：幅度谱可以发现其幅度呈偶对称，而相位谱与幅度普均为离散信号。其频率特性就是通过傅里叶变换得来的，是用频率信息表示时域信息，与连续信号的相同的矩形函数的频率特性对应，只不过一个是离散的一个是连续的。 单位脉冲序列 生成单位脉冲序列，观察并记录生成的图形； n=1:50; %定义序列的长度是50 x=zeros(1,50); %注意：MATLAB中数组下标从1开始 x(1)=1; close all; subplot(3,1,1);stem(x);title(单位冲击信号序列); 研究其频率特性，，分别研究其幅频特性和相频特性，观察并记录生成的图形； k=-25:25; X=x*(exp(-j*pi/12.5)).^(n*k); magX=abs(X); %绘制x(n)的幅度谱 subplot(3,1,2);stem(magX);title(单位冲击信号的幅度谱); angX=angle(X); %绘制x(n)的相位谱 subplot(3,1,3);stem(angX) ; title (单位冲击信号的相位谱) 答：根据以上两个源程序得到的图形如下图1.2： 图1.2 分析：由上图观察单位冲激信号的频率响应可知，其幅度谱为矩形函数，相位谱呈奇对称。与连续的冲激响应对应，只是这里是有限序列。 卷积过程 ， n=1:50; %定义序列的长度是50 hb=zeros(1,50); %注意：MATLAB中数组下标从1开始 hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1; close all; subplot(3,1,1);stem(hb);title(系统hb[n]); m=1:50; %定义序列的长度是50 A=444.128; %设置信号有关的参数 a=50*sqrt(2.0)*pi; T=0.001; %采样率 w0=50*sqrt(2.0)*pi; x=A*exp(-a*m*T).*sin(w0*m*T); subplot(3,1,2);stem(x);title(输入信号x[n]); y=conv(x,hb); subplot(3,1,3);stem(y);title(输出信号y[n]); 答：根据以上两个源程序得到的图形如下图1.3： 图1.3 分析：离散信号的卷积为相对移动然后累计相加，所以卷积的结果最高可达到1000. 思考： 如何生成实指数序列？写出对应的matlab程序； 答：（1）程序如下所示： a1=2 n=1:50 x1=(a1.^n) subplot(1,1,1) stem(x1); title(实指数序列) 所得图形如下图1.4所示： 编写程序验证卷积定律。 答：编写的源程序如下： x=[1 2 0 3 1];%原序列 y=[1 2 0 3 1]; %直接得到结果 z=conv(x,y); figure(1),subplo