发布时间 : 星期日 文章数字信号处理实验-熟悉MATLAB环境更新完毕开始阅读
西安交通大学实验报告
成绩
第 页 共 页
课 程 数字信号处理
系 别 实 验 日 期 年 月 日
专业班级 组别 交报告 日 期 年 月 日 姓 名 学号 报 告 退 发 (订正、重做) 同 组 者 教师审批签字
实验名称 熟悉MATLAB环境
一、 实验目的
（1） 熟悉MATLAB的主要操作命令； （2） 学会简单的矩阵输入和数据读写； （3） 掌握简单的绘图命令； （4） 学习用MATLAB编程。
二、 实验内容
（1）已知位置向量同为n=0:6的以下三个序列：
x=[2,-1,4,7,0,-3,5];y=[-2,8,0,-3,6,2,-4];w=[7,0,-9,-3,4,2,-1]; 请用MATLAB计算 (a)u=x+y (b)v=x.*w (c)z=x-y.*w
（2）用MATLAB计算下列序列的线性卷积并作图：（可用conv函数求卷积，用stem函数作图）
(a)x(n)=δ(n)，h(n)=R5(n) (b)x(n)= R3(n)，h(n)=R4(n) (c)x(n)=δ(n-2)，h(n)=0.5R3(n)
（3）设滤波器差分方程为y(n)=x(n)+x(n-1)+1/3y(n-1)+1/4y(n-2), 请用
n
MATLAB计算并作图。
(a)系统的频率响应（数值与图形0～π共512点）。（可用freqz函数） (b)设抽样频率为10kHz，输入正弦波幅度为5，频率为1kHz，记录长度5ms，求系统输出。（可用filter函数）
（4）一个特定的线性时不变系统，描述它的差分方程如下： y(n)+0.1y(n-1)-0.06y(n-2)=x(n)-2x(n-1)
(a) 在0≤n≤10之间求得并画出系统得单位脉冲响应，从脉冲响应确定系统的稳定性。
(b)如果此系统的输入为x(n)=[5+3cos(0.2πn)+4sin(0.6πn)]u(n) 0≤n≤20，求响应y(n)并将x(n)与y(n)作图。
三、实验内容程序及结果
（1）已知位置向量同为n=0:6的以下三个序列：
x=[2,-1,4,7,0,-3,5] y=[-2,8,0,-3,6,2,-4] w=[7,0,-9,-3,4,2,-1];
用MATLAB计算(a)u=x+y (b)v=x.*w (c)z=x-y.*w 运行程序为：
x=[2,-1,4,7,0,-3,5]; y=[-2,8,0,-3,6,2,-4]; w=[7,0,-9,-3,4,2,-1]; n=0:6; u=x+y v=x.*w z=x-y.*w
subplot(3,1,1);stem(n,u) subplot(3,1,2);stem(n,v) subplot(3,1,3);stem(n,z)
运行结果如下图：
u = 0 7 4 4 6 -1 1 v = 14 0 -36 -21 0 -6 -5 z = 16 -1 4 -2 -24 -7 1
实验结果分析：由图所示结果可知u,v,z的结果，计算正确，并有图示序列表示。 （2）用MATLAB计算序列的线性卷积并作图。 (a) x(n)=δ(n)，h(n)=R5(n) 运行程序为：
x=[1];h=[1,1,1,1,1]; y=conv(x,h);
n1=0;n2=0:(length(y)-1);n3=0:4; subplot(3,1,1);stem(n1,x)
subplot(3,1,2);stem(n2,h) subplot(3,1,3);stem(n3,y) 运行结果如下图：
(b) x(n)= R3(n)，h(n)=R4(n) 运行程序为：
x=[1,1,1]; h=[1,1,1,1]; y=conv(x,h) n1=0:2;n2=0:3; n3=0:(length(y)-1); subplot(3,1,1);stem(n1,x) subplot(3,1,2);stem(n2,h) subplot(3,1,3);stem(n3,y) 运行结果如下图：
(c) x(n)=δ(n-2)，h(n)=0.5R3(n) 运行程序为：
x=[1]; h1=[1,1,1]; n1=2; n2=0:2; h=h1.*(0.5.^n2); y=conv(x,h); n3=0:(length(y)-1); subplot(3,1,1);stem(n1,x) subplot(3,1,2);stem(n2,h) subplot(3,1,3);stem(n3,y) 运行程序如下图：
n