数字信号处理实验报告MATLAB

1文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.

2文档来源为:从网络收集整理.word 版本可编辑.欢迎下载支持. 00.2

0.4

0.6

0.8

1

1.2

1.4

再在命令窗口输入：[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

再在命令窗口输入：