信号与系统上机实验 电子科技大学.docx

信号与系统上机实验Q10=n=31 x1(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) 分别画出图形,求出其周期。（1）x1(n)=sin(pi*n/4)*cos(pi*n/4);程序如下：n=0:31x1=sin(pi*n/4).*cos(pi*n/4)stem(n,x1)MATLAB画出图形如下：由上图图可知周期T=4（2）x2(n)=cos(pi*n/4)*cos(pi*n/4)程序如下：n=0:31x2=cos(pi*n/4).*cos(pi*n/4) stem(n,x2)MATLAB画出图形如下：由上图可知周期T=4（3）x3(n)=sin(pi*n/4)*cos(pi*n/8) 程序如下：n=0:31x3=sin(pi*n/4).*cos(pi*n/8) stem(n,x3)MATLAB画出图形如下：由上图可知周期T=16Q2当0=n=5时，h(n)=n; 其他h(n)=0; x(n)=h(n);求y(n)=x(n)*h(n); 用stem函数画出y(n).程序如下：n=0:5; y=n.^2; stem(y)MATLAB画出图形如下：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=4a1=[5,0,-4]b1=[10,0,-5][H1,W1]=freqz(b1,a1,n)[H2,W2]=freqz(b1,a1,n,whole)输出结果如下：n = 4a1 = 5 0 -4b1 = 10 0 -5H1 = 5.0000 + 0.0000i 1.7073 - 0.3659i 1.6667 + 0.0000i 1.7073 + 0.3659iW1 = 0 0.7854 1.5708 2.3562H2 = 5.0000 1.6667 5.0000 1.6667W2 = 0 1.5708 3.1416 4.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=n=63的图形 程序如下：N=64;n=0:63;x1=zeros(1,64);x2=zeros(1,64);x3=zeros(1,64);flag1=0;flag2=0;for i=1:64x1(i)=1;endfor i=1:64if((i=1i=8)||(i=17i=24)||(i=33i=40)||(i=49i=56))x2(i)=1;elsex2(i)=0;endendfor i=1:64if((i=1i=8)||(i=33i=40))x3(i)=1;elsex3(i)=0;endendstem(n,x1);figure;stem(n,x2);figure;stem(n,x3);用MATLAB画出图形分别如下：(2) 求其对应的付氏级数，(分别为a1,a2,a3)并画图。 程序如下：N=64;n=0:63;x1=zeros(1,64);x2=zeros(1,64);x3=zeros(1,64);flag1=0;flag2=0;for i=1:64x1(i)=1;endfor i=1:64if((i=1i=8)||(i=17i=24)||(i=33i=40)||(i=49i=56))x2(i)=1;elsex2(i)=0;endendfor i=1:64if((i=1i=8)||(i=33i=40))x3(i)=1;elsex3(i)=0;endenda1=fft(x1)figure;stem(a1)figure;a2=fft(x2)stem(a2)figure;a3=fft(x3)stem(a3)MATLAB画出波形如下：(3) 综合x3 a. 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