实验1

实 验 报 告 课程名称 信号处理实验 实验项目名称 基本信号 实验类型 实验学时 2 班级 20120815 学号 2012081512 姓名 邱岳 指导教师 陈立伟 实验室名称 实验时间 11.24 14:00-15:50

实验成绩 预习部分 实验过程 表现 实验报告 部分 总成绩

教师签字 日期

哈尔滨工程大学教务处 制 实验一 基本信号 一．实验原理 本实验主要使用MATLAB内部向量程序来产生信号。用MATLAB的stem指令绘出离散时间信号。下面的MATLAB代码产生31点离散时间正弦信号。 nn=0:30; %时间标号向量 sinus=sin(nn/2+1);

依据MATLAB编址约定，标号0nn必须对应nn(1);同样，sinus（1）是正弦信号的第一个数值。Stem指令产生离散时间信号图形。使用stem指令绘出正弦波。 stem(nn,sinus);必须给定向量的第一个参数以得到正确的n轴。 1、 冲激信号 最简单的信号是单位冲激信号：

000

1[](3.1)0nnnnnn





在MATLAB中产生冲激信号，必须先确定所关注信号部分的长度。 2、正弦信号 另一个非常基本的信号是正弦信号。一般，完备的描述实正弦信号需要三个参数：振幅

A、频率0和相位。

0[]cos()(3.2)xnAn 3、指数信号 （1）衰减的指数信号是数字信号处理中的基本信号。因为它是线性常微分差分方程的解。 （2）指数序列在信号处理中常出现的另一个原因是，时移并不改变其信号特征。

（3）可以用差分方程给出的递归表示式来产生指数信号。当输入[]xn是一个冲激信号

的时候，信号[][]nynaun 是下面差分方程的解： [][1][],ynaynxn初始条件y-1=0(3.3) 由于假定差分方程以因果方式递归（即n增大），1n的初始条件是必须的。 用MATLAB的filter函数可实现差分方程。 函数filter的用法如下所示

(,,)youtfilterbaxin 假设描述离散时间系统输入与输出关系的线性常系数差分方程为 00()()nnkkkkaynkbxnk 那么01[,,]Mbbbb和01[,,]Naaaa是差分方程中的系数组成的向量，xin是输入信号向量（filter（）函数只向y返回与x中样本个数一样多的样本）。 4、复值信号 复值信号通常是非常有用的。通过把信号配对，分别作为复数的实部与虚部，再将此信号对与其他复值信号依据复数运算规则进行运算处理来实现。以这种方式使用信号对是许多信号处理系统的重要部分，尤其是那些涉及调制的系统。 在MATLAB中，real函数和image函数分别提取复数的实部和虚部。画复向量时，plot指令和stem指令的缺省设置容易使人困惑。如果z是复数，那么plot（z）指令绘出z实部与虚部的对应点；plot（n,z）指令绘出z实部与n的对应点（连线），stem（z）指令既标注实部又标注虚部。如果想同时看到绘出的实部与虚部，在stem指令之前的subplot可将屏分成几部分。 5、复指数信号 实指数信号表示式可以扩展成包含正弦和余弦信号的复指数信号形式。这些信号构成了傅里叶变换的基础。 （1）．在MATLAB中，复信号是实指数信号的自然扩展。复指数信号表示为：

00(ln)0[]()(cosjsin)(3.4)zjznnnjnnxnzerernn

此处0jnzrer。 （2）式（3.4）不足以产生所有复指数信号的通式，通式为： ()0[cos()sin()]njnjnnjnnGzAereAreArnjn

这里0njGzAeA是复指数信号的复振幅。 （3）复指数信号也可由一阶差分方程（利用filter函数）产生： 0[][1][](3.5)ynzynxn

滤波器系数0jzre 是复数，容易看出序列中连续项之比为0z;但是正确的振幅和相位必须由选择激励差分方程的冲激信号的复振幅（即[][]xnGn ）来设置。 二．实验内容 1、冲激信号： 产生并绘出下面的序列。在每种情形下，水平n轴应该只在指定的区间展开并相应标注。使用stem指令使每个序列显示成离散时间信号。 1234

[]0.9[5]120[]0.8[]1525[]1.5[333]300350[]4.5[7]100xnnnxnnnxnnnxnnn 实验程序： L=20;%信号长度 nn=1:20;%n的取值为0-20 imp=zeros(L,1);%产生一个L行1列的零矩阵 imp(5)=0.9; %n=5时，信号值为0.9 subplot(2,2,1); %将屏幕分为2行2列，指定接下来要画的图放在第一个拼块里 stem(nn,imp); %以离散形式绘制x[n] xlabel('n'); ylabel('x1[n]'); L=31; nn=-15:15; imp=zeros(L,1); imp(16)=0.8; subplot(2,2,2); stem(nn,imp); xlabel('n'); ylabel('x2[n]'); L=51; nn=300:350; imp=zeros(L,1); imp(34)=1.5; subplot(2,2,3); stem(nn,imp); xlabel('n'); ylabel('x3[n]'); L=11; nn=-10:0; imp=zeros(L,1); imp(4)=4.5; subplot(2,2,4); stem(nn,imp); xlabel('n'); ylabel('x4[n]');

实验结果：0246810121416182000.20.40.60.81nx1[n]-15-10-505101500.20.40.60.8nx2[n]

30030531031532032533033534034535000.511.5nx3[n]-10-9-8-7-6-5-4-3-2-10012345

nx4[n]

图1 四个冲激信号 2、正弦信号 产生并绘出下列每一个序列。使用MATLAB的向量功能求解此问题，将向量参数赋予余弦或正弦函数，再利用一个函数调用。在每种情形下，水平n轴应该只在指定的区间展开并相应标注。使用stem指令使每个序列显示成离散时间信号。

1234

()sin(/17)025()sin(/17)1525()sin(3/2)1010()cos(/23)050xnnnxnnnxnnnxnnn 实验程序： n=[0:25];%n的取值 y=sin(pi*n/17); subplot(2,2,1); stem(n,y);%离散形式画出x[n] xlabel('n'); ylabel('x1[n]'); n=[-15:25]; y=sin(pi*n/17); subplot(2,2,2); stem(n,y); xlabel('n'); ylabel('x2[n]'); n=[-10:10]; y=sin(3*pi*n+pi/2); subplot(2,2,3); stem(n,y); xlabel('n'); ylabel('x3[n]'); n=[0:50]; y=cos(pi/sqrt(23)*n); subplot(2,2,4); stem(n,y); xlabel('n'); ylabel('x4[n]'); 实验结果：

0102030-1-0.500.51nx1[n]-2002040-1-0.500.51

nx2[n]

-10-50510-1-0.500.51nx3[n]0204060-1-0.500.51

nx4[n]

图2 四个正弦信号 结果分析: 从上图可以看出4()xn不是周期序列，因为由周期序列的定义，当2NM 为有理数时，

序列是周期的，而这里22223/23 ，是无理数，故它不是周期序列。 3、指数信号 a. 绘制指数信号[](0.9)020nxnn b. 对(a)中产生的指数信号求和，验证 101,1(3.6)1LLnnaaaa





