实验一 离散时间信号时域表示

《实验一 离散时间信号时域表示》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、数字信号处理实验报告实验一离散时间信号的时域表示一、实验目的1、熟悉Matlab命令，掌握离散时间信号－序列的时域表示方法。2、掌握用Matlab描绘二维图像的方法。3、掌握用Matlab对序列进行基本的运算和时域变换的方法。二、实验内容（1）参考示例程序，产生一个有延迟的单位抽样序列：d(n-11)，，绘出序列的图形。>>n=[5:15];x=(n>=11)-(n>=12);stem(n,x);>>（2）参考示例程序，产生一个向前时移7个时刻的单位阶跃序列：u(n+7)，，绘出序列的图形。>>n=[-10:

2、10];x=(n>=-7);stem(n,x);（3）产生一个指数为[-0.1+(pi/6)*i]n的复指数序列，并绘出序列的实部、虚部、幅度和相位的波形。>>n=[-10:10];x=exp((-0.1+(pi/6)*i).*n);subplot(221);stem(n,real(x));subplot(222);stem(n,imag(x));subplot(223);stem(n,abs(x));subplot(224);stem(n,angle(x));（4）已知x(n)={1,2,3,4,5,6,7

3、,6,5,4,3,2,1}，，参考示例程序，绘出下列序列的波形。a.x1(n)=2x(n-5)-3x(n+4)b.x2(n)=x(3-n)+x(n)x(n-2)A．>>n=-2:10;n01=5;n02=-4;n03=3;n04=2;x=[1234567654321];n1=n+n01;n2=n+n02;n3=n+n03;n4=n+n04;n5=-fliplr(n3);%x(3-n)ÐòÁÐλÖõĕ´ÕÛnn1=min([n1,n2]):max([n1,n2]);N1=length(nn1);y=zer

4、os(1,N1);y1=y;y2=y;y(find((nn1>=min(n))&(nn1<=max(n))))=x;y1(find((nn1>=min(n1))&(nn1<=max(n1))))=x;%y1=x(n-5)y2(find((nn1>=min(n2))&(nn1<=max(n2))))=x;%y2=x(n+4)x1=2*y1-3*y2;figure(1);subplot(2,2,1),stem(nn1,y);title('x(n)');axis([min(nn1),max(nn1),1.2*min

5、(y),1.2*max(y)]);subplot(2,2,2),stem(nn1,y1);title('x(n-5)');axis([min(nn1),max(nn1),1.2*min(y1),1.2*max(y1)]);subplot(2,2,3),stem(nn1,y2);title('x(n+4)');axis([min(nn1),max(nn1),1.2*min(y2),1.2*max(y2)]);subplot(2,2,4),stem(nn1,x1);title('2x(n-5)-3x(n+4)')

6、;axis([min(nn1),max(nn1),1.3*min(x1),1.3*max(x1)]);B．n=-2:10;n01=5;n02=-4;n03=3;n04=2;x=[1234567654321];n1=n+n01;n2=n+n02;n3=n-n03;n4=n+n04;n5=-fliplr(n3);%x(3-n)序列位置的反折nn2=min([n4,n5]):max([n4,n5]);N2=length(nn2);z=zeros(1,N2);z1=z;z2=z;z(find((nn2>=min(n)

7、)&(nn2<=max(n))))=x;z1(find((nn2>=min(n4))&(nn2<=max(n4))))=x;%z1=x(n-2)z2(find((nn2>=min(n5))&(nn2<=max(n5))))=x;%z2=x(3-n)x2=z2+z.*z1;figure(2);subplot(2,2,1),stem(nn2,z);title('x(n)');axis([min(nn2),max(nn2),1.2*min(z),1.2*max(z)]);subplot(2,2,2),stem(nn

8、2,z1);title('x(n-2)');axis([min(nn2),max(nn2),1.2*min(z1),1.2*max(z1)]);subplot(2,2,3),stem(nn2,z2);title('x(3-n)');axis([min(nn2),max(nn2),1.2*min(z2),1.2*max(z2)]);subplot(2,2,4),stem(nn2,x1);title('