发布时间 : 星期三 文章数字信号处理实验一(上机)报告更新完毕开始阅读5e0c9066770bf78a652954c3
（二）序列的基本运算和时域变换 1、加法：x1(n)+x2(n)
序列的加法运算为对应位置处量值的相加，在MATLAB中可用运算符“+”实现，但要求参与运算的序列的长度必须相等。如果长度不等或者长度相等但采样位置不同，则不能直接应用该运算符，此时需要先给定参数使序列具有相同的位置向量和长度。下面给出sigadd函数实现任意两序列的加法运算。
例1.7：
function [y,n] = sigadd(x1,n1,x2,n2) % implements y(n) = x1(n)+x2(n) %
% [y,n] = sigadd(x1,n1,x2,n2)
% y = sum sequence over n, which includes n1 and n2 % x1 = first sequence over n1
% x2 = second sequence over n2 (n2 can be different from n1) %
n = min(min(n1),min(n2)):max(max(n1),max(n2)); % duration of y(n)
y1 = zeros(1,length(n)); y2 = y1; % initialization
y1(find((n>=min(n1))&(n<=max(n1))==1))=x1; % x1 with
duration of y
y2(find((n>=min(n2))&(n<=max(n2))==1))=x2; % x2 with duration of y
y = y1+y2; % sequence addition
其中x1和x2为参与加法运算的两序列，n1和n2分别为x1和x2的位置向量。
2、乘法：x1(n)·x2(n)
序列的乘法运算为对应位置处量值的相乘，在MATLAB中由数组运算符“.*”实现，也受到“+”运算符同样的限制。
3、反折：x(n)→x(-n)
序列的反折指序列的每个量值都对n=0做一个对称操作，从而得到一个新序列。在MATLAB中可由fliplr(x)函数实现，此时序列位置的反折则由-fliplr(n)实现。
4、平移：x(n)→x(n-m)
平移操作是将序列的每个量值都移动m个位置，在得到的新序列中，量值和原序列相同，只是位置向量n发生变化，当m>0时，表示序列向右平移，此时新序列的位置向量为n+m；当m<0时，表示序列向左平移，此时新序列的位置向量为n-m。
三、 实验内容及结果 （一）内容：
参考示例程序，产生一个有延迟的单位抽样序列：? (n-11)，5?n?15，绘出序列的图形。 （二）结果： （1）源程序代码：
（2）运行结果：
（一）内容：
参考示例程序，产生一个向前时移7个时刻的单位阶跃序列：u(n+7)，
?10?n?10，绘出序列的图形。
（二）结果： （1）源程序代码：