数字信号处理实验报告 MATLAB.docVIP

数字信号处理实验报告 姓名：  班级：09电信一班 学号： 实验报告（一） 姓名： 学号： 实验日期：8-12周 实验题目：离散信号的时域运算与变换 实验目的：1. 熟悉MATLAB编程特点 2. 了解离散序列的延迟、相加、相乘及平移、反折、及倒相变换 实验内容：1.设计一个实现序列移位的函数 将序列x(n)＝｛1,2,3,4,0,7｝,其中n＝0:5的每一个样本都移动3个周期,移位后的序列y(n)=x(n-3) 2.序列的奇偶分解 将序列x(n)＝｛0,1,2,3,4,3,2,1,0｝,其中n＝-3:5进行奇偶分解 用函数stem显示其奇偶序列 3.序列的加法运算 设x1(n)＝｛1,0.5,0.3,0.4｝其中n=-1:2 ;x2(n)＝｛0.2,0.3,0.4,0.5,0.8,1｝其中n=-2:3 4.实现2中序列的翻转 实验地点：4305 实验结果：1.将如下文件另存为：segshift.m文件， function [y,n] = segshift(x,n,n0) % 功能 y(n) = x(n-n0) % [y,n] = segshift(x,n,n0) n = n+n0; y = x; 在命令窗口输入：x=[1,2,3,4,0,7];n=0:5;n0=3;[y,n] = segshift(x,n,n0)按ENTER键得到如下结果， y = 1 2 3 4 0 7 n = 3 4 5 6 7 8 再在命令窗口输入：[y,n] = segshift(x,n,n0);n=0:5;n0=3;stem(n,y)按ENTER键出图如下图1。 2. 在命令窗口输入: m = -fliplr(n); m1 = min([m,n]); m2 = max([m,n]); m = m1:m2; x1 =[ zeros(1,(length(m)-length(n))),x]; xe = 0.5*(x1 + fliplr(x1)) xo = 0.5*(x1 - fliplr(x1)) subplot(1,2,1),stem(m,xe),ylabel(xe), subplot(1,2,2),stem(m,xo),ylabel(xo), 得下图二， 图二 图一 3．将如下文件另存为：sigadd.m文件 function [y,n] = sigadd(x1,n1,x2,n2) % 实现 y(n) = x1(n)+x2(n) % ----------------------------- % [y,n] = sigadd(x1,n1,x2,n2) % y = 在包含n1 和 n2 的n点上求序列和, % x1 = 在 n1上的第一序列 % x2 = 在 n2上的第二序列(n2可与 n1不等) n = min(min(n1),min(n2)):max(max(n1),max(n2)); % y(n)的长度 y1 = zeros(1,length(n)); y2 = y1; % 初始化 y1(find((n=min(n1))(n=max(n1))==1))=x1; % 具有y的长度的 x1 y2(find((n=min(n2))(n=max(n2))==1))=x2; % 具有y的长度的x2 y = y1+y2; 在命令窗口输入：x1=[1,0.5,0.3,0.4];n1=-1:2;x2=[0.2,0.3,0.4,0.5,0.8,1];n2=-2:3; [y,n] = sigadd(x1,n1,x2,n2) 得：y = 0.2000 1.3000 0.9000 0.8000 1.2000 1.0000 n = -2 -1 0 1 2 3 再在命令窗口输入：x1=[1,0.5,0.3,0.4];n1=-1:2;x2=[0.2,0.3,0.4,0.5,0.8,1];n2=-2:3;[y,n]=sigadd(x1,n1,x2,n2) ；stem(n,y)，得到下图： 4. 将如下文件另存为.m文件 n=-1:10; x=sin(0.4*pi*n); y=fliplr(x); n1=-fliplr(n); subplot(2,1,1),stem(