数字信号处理实验1离散序列的基本运算

1、实验1离散序列的基本运算学号： 姓名：一、实验目的学会运用matlab进行离散序列的运算，并掌握程序的编写方法。二、实验涉及的matlab子函数1、find功能：寻找非零元素的索引号。调用格式：find ( (n>=min(nl) &(n<=max(n 1)；在符合关系运算条件的范围内寻找非零 元素的索引号。2、fliplr功能：对矩阵行元素进行左右翻转。调用格式：xl= fliplr (x)；将x的行元素左右翻传，赋给变量xl。三、实验原理离散序列的时域运算包括信号的相加、相乘，信号的时域变换包括信号的移 位、反折等。在matlab中，离散序列的相加、相乘等运算是两个向量

2、之间的运算，因 此参加运算的两个序列向量必须具冇相同的维数，否则应进行相应的处理。下面介绍齐种离散序列的时域运算和时域变换的性质。(n) =x (nm)%牛成离散信号的时间序列%生成离散信号xo(n)%生成离散信号x 1 (n)%生成离散信号x2(n)1、序列的移位x例、x0=u (n),(显示-10<n<10) xi=u (n+6)(显示-10<n<10), x2= u (n-4)(显示-10<n<10)matlab源程序为： k0=0;kl=-6;k2=4; n=-10:10;xo=n>=ko;xl=(n-kl)>=0;x2=(n-k2)&g

3、t;=0;subplot(3,l, 1 ),stem(n, xo,'filled','k'); axis(n 1 ,n2,l .1 *min(xo),l .1 *max(xo); ylabelcucn)*)；subplot(3,1,2),stem(n, x 1,'filled',v); axis(n 1 ,n2,1.1 *min(x 1 ),1.1 *max(x 1); ylabel('u(n+6)t);subplot(3,l ,3),stem(n, x2,filled',k,); axis(n 1 ,n2,l.l *min(x2

4、),l.l *max(x2); ylabel('u(n4)');例、已知x (n) =2sm,求x (n-2)和x (n+2)在-210区间的波形。10matlab源程序为：n=-2:10 ;n0=2; n 1=-2;x=2*sin(2*pi*n/10);%建立原信号 x(n)xl=2*sin(2*pi*(n-n0)/10); %建立 x(n2)信号 x2=2*sin(2*pi*(n-n 1)/10); %建立 x(n+2)信号 subplot(3,l, 1 ),stem(n,x,filled7k,); ylabek'x(ny);subplot(3,l ,2),stem

5、(n,x 1,'filled1,v); ylabel('x(n2);subplot(3,l ,3)?stem(n,x2,filled,k,); ylabel('x(n+2y);2、序列的相加情况1两序列具有相同的维数例、求 x(n)=(7(n - 2)+ cr(n - 4)(0<n< 10)matlab源程序为：n 1 =0;n2= 10;n01 =2;n02=4;n=nl:n2;xl=(n-n01)=0;x2=(n-n02)=0;x3=xl+x2;subplot(3,l ,1 );stem(n,xl，filled','k'); ax

6、is(n 1 ,n2,1.1 *min(x 1),1.1 *max(x 1); ylabel(,5(n-2)t);subplot(3,l ,2);stem(n,x2/filled7k'); axis(nl ,n2,1.1 *min(x2),l .1 *max(x2); ylabel(3(n4y);subplot(3,l ,3)；stem(n,x3/filled7k'); axis(n 1 ,n2,1.1 *min(x3),l.l *max(x3); ylabel(,5(n-2)+8(n-4)');,5 a euio1367810,5 a (truio1236810-.5

7、910情况2两序列具有不同的维数例、xi (n) =r4n+2-2<=n<=5x2(n)=r5 n-11 <=n<=5 求 x (n) =xn+x2nmatlab源程序为：nl=-4:6;x 1=(n 1 +2)>=0)&(n 1 -5)<=0);%建立 xl 信号n2=-l:5;x2=(n2-1 )>=0)&(n2-5)<=0);%建立 x2 信号n=min(n 1 ,n2):max(n 1 ,n2); %为 x 信号建立时间序列 n n=length(n);%求时间序列n的点数nyl=zeros(l,n);y2=zeros(l

8、,n);%新建一维n列的yl全0数组%新建一维n列的y2全。数组y 1 (find(n>=min(n 1)&(n<=max(n 1 )=x 1; %为 yl 赋值 y2(find(n>=min(n2)&(n<=max(n2)=x2; %为 y2 赋值 x=yl+y2;subplot(3,l ,1 ),stem(n 1 ,x 1 /filled'/k1); ylabelfxl(n)');axis(min(n),max(n), 1.1 *min(x), 1.1 *max(x); subplot(3,l ,2),stem(n2,x2,fille

9、d7k,);ylabel('x2(n)')；axis(min(n),max(n), 1*min(x), 1.1 *max(x);3、序列相乘例、xl (n) =3e'°'25n ri5n+4(-4<=n<= 10)x2(n)=r9(n+2)(-2<=n<=6)求 x (n) = xl (n) *x2(n)matlab源程序为：nl=-4:10;xl=3*exp(-0.25*n);%建立 xl 信号n2=-2:6;n02=-2;x2=(n2-n02)>=0;%建立 x2 信号n=min(nl,n2):max(nl,n2);

10、%为 x 信号建立时间序列 n n=length(n);%求时间序列n的点数nyl=zeros(l,n);%新建一维n列的yl全0数组y2=zeros(l,n);%新建一-维n列的y2全0数组y 1 (find(n>=min(n 1)&(n<=max(n 1 )=x 1; %为 yl 赋值 y2(find(n>=min(n2)&(n<=max(n2)=x2; %为 y2 赋值 x=yl.*y2;subplot(3,l ,l),stem(n,x 1 /filled'/k1); ylabelcxhn)*); axis(min(n),max(n), 1

11、*min(x 1),1.1 *max(x 1); subplot(3,l,2),stem(n,y2,'filled',显示 x (n) = cr(n + 3)+ 2cr(n -4) 显示 x (n) =u (n-2) +u (n+2) 已知 xl=exp (- n /16) rs(n+2), x2=5sin (2*pi*n/10) r7(n-2), yn=xl 乘以x2«);ylabelcx2(n)*);axis(min(n),max(n), 1.1 *m in(y2),l.l * max(y2); subplot(3,l ,3),stem(n,x/filled7k');ylabelcx(ny);axis(min(n),max(n), 1.1 *min(x), 1.1 *max(x);4、序列反折例、x (n) =e0'3n(-4<n<4),求 x (-n)omatlab源程序为: n=-4:4;x=exp(-0.3*n); xl=fliplr(x);subplot( 1,2,1 ),stem(n,x,'filled7