资源描述：
《信号与系统上机实验报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、姓名：学号：学院：Q1:0<=n<=31x1(n)=sin(pi*n/4)*cos(pi*n/4);x2(n)=cos(pi*n/4)*cos(pi*n/4)x3(n)=sin(pi*n/4)*cos(pi*n/8)分别画出图形,求出其周期。代码：n=0:31;x1=sin(pi*n/4).*cos(pi*n/4);x2=cos(pi*n/4).*cos(pi*n/4);x3=sin(pi*n/4).*cos(pi*n/8);stem(n,x1)stem(n,x2)stem(n,x3)结果：由图形可知周期T1=4由图形可知周期T2=4由图形可知周期T3=16Q2:当0<=n<=5时，h
2、(n)=n;其他h(n)=0;x(n)=h(n);求y(n)=x(n)*h(n);用stem函数画出y(n).代码：n=0:5;h=n;x=h;y=conv(x,h);stem(y)结果：Q3：（a）.定义用向量a1和b1描述差分方程y(n)-0.8y(n-1)=2x(n)-x(n-2)表征的因果LTI系统，(b).用在（a）中的系数向量，利用freqz定义H1是在0和pi之间4个等份频率上频率响应的值，omega1是这些频率值。(c).用在（a）中的系数向量，利用freqz定义H2是在0和2*pi之间4个等份频率上频率响应的值，omega2是这些频率值。代码：n=4;a1=[5,0,
3、-4]b1=[10,0,-5][H1,W1]=freqz(b1,a1,n)[H2,W2]=freqz(b1,a1,n,'whole')结果：a1=50-4b1=100-5H1=5.00001.7073-0.3659i1.66671.7073+0.3659iW1=00.78541.57082.3562H2=5.00001.66675.00001.6667W2=01.57083.14164.7124Q4:X1(n)=u(n)-u(n-8);其周期N1=8,X2(n)=u(n)-u(n-8);其周期N2=16,X3(n)=u(n)-u(n-8);其周期N3=32,(1)画出这些周期信号在0<
4、=n<=63的图形(2)求其对应的付氏级数，(分别为a1,a2,a3)并画图。(3)综合x3a.x31(n)=a3(1)*exp(j*2*pi*n/32)+a3(2)*exp(j*2*2*pi*n/32)+a3(30)*exp(j*2*30*pi*n/32)+a3(31)*exp(j*2*31*pi*n/32)+a3(32)*exp(j*2*32*pi*n/32)b.x32(n)=a3(1)*exp(j*2*pi*n/32)+……+a3(8)*exp(j*2*8*pi*n/32)+a3(24)*exp(j*2*24*pi*n/32)+……+a3(32)*exp(j*2*32*pi*n/
5、32)c..x33(n)=a3(1)*exp(j*2*pi*n/32)+……+a3(12)*exp(j*2*12*pi*n/32)+a3(20)*exp(j*2*20*pi*n/32)+……+a3(32)*exp(j*2*32*pi*n/32)c..x34(n)=a3(1)*exp(j*2*pi*n/32)+……+a3(32)*exp(j*2*32*pi*n/32)yi=abs(x3i),用stem对yi作图，试比较其于x3的区别。(1)(2)代码：N=64;n=0:63;x1=zeros(1,64);x2=zeros(1,64);x3=zeros(1,64);flag1=0;flag
6、2=0;fori=1:64x1(i)=1;endfori=1:64if((i>=1&&i<=8)
7、
8、(i>=17&&i<=24)
9、
10、(i>=33&&i<=40)
11、
12、(i>=49&&i<=56))x2(i)=1;elsex2(i)=0;endendfori=1:64if((i>=1&&i<=8)
13、
14、(i>=33&&i<=40))x3(i)=1;elsex3(i)=0;endendstem(n,x1);figure;stem(n,x2);figure;stem(n,x3);a1=fft(x1)figure;stem(a1)figure;a2=fft(x2)stem(a2)figure;a3
15、=fft(x3)stem(a3)结果：(3)代码：x31=a3(1)*exp(1i*2*pi*n/32)+a3(2)*exp(1i*2*2*pi*n/32)+a3(30)*exp(1i*2*30*pi*n/32)+a3(31)*exp(1i*2*31*pi*n/32)+a3(32)*exp(1i*2*32*pi*n/32)stem(abs(x31))figure;x32=sum(a3(1:8))*exp(1i*2*8*pi*n/32)+sum(
