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

1、电子科技大学信号与系统上机实验报告 3信号与系统实验报告姓名：学号：学院：指导教师：实验时间： 1 第一题1.1 实验目的画出离散时间正弦信号并确定基波周期1.2 实验程序n=0:31;x1=sin(pi*n/4).*cos(pi*n/4);x2=(cos(pi*n/4).2;x3=sin(pi*n/4).*cos(pi*n/8);subplot(3,1,1);stem(n,x1)title(figures of signals)xlabel(n)ylabel(x1n)subplot(3,1,2);stem(n,x2)xlabel(n)ylabel(x2n)subplot(3,1,3);ste

2、m(n,x3)xlabel(n)ylabel(x3n)1.3 实验图像1.4 实验结论可得信号x1n与x2n的周期均为4，x3n的周期为16.2 第二题2.1 实验目的确定离散时间系统的性质证明 不满足线性性质证明不是因果的2.2 实验程序2.2.1 系统1x1=1 zeros(1,10);x2=2 zeros(1,10);y1=sin(pi/2).*x1);y2=sin(pi/2).*x2);x3=2*x1+6*x2;y3=sin(pi/2).*x3);y4=2*y1+6*y2;subplot(2,1,1);stem(n,y3)title(figure of y3)xlabel(n)ylab

3、el(y3)subplot(2,1,2);stem(n,y4)title(figure of y4)xlabel(n)ylabel(y4)2.2.2 系统2nx=-5:9;x1=zeros(1,5) 1 ones(1,9);x2=zeros(1,4) 1 ones(1,10);y=x1+x2;subplot(3,1,1);stem(nx,x1);xlabel(n);ylabel(x1n);subplot(3,1,2);stem(nx,x2);xlabel(n);ylabel(x2n);subplot(3,1,3) ; stem(nx,y);xlabel(n);ylabel(yn);2.3 实验

4、图像2.3.1 系统1 2.3.2 系统22.4 实验结论系统不是线性的。系统不是因果的。3 第三题3.1 实验目的计算卷积3.2 实验程序x,n1=stepseq(0,0,5);h=0 1 2 3 4 5;y=conv(x,h);3.3 实验图像3.4 实验结论计算得xn和hn的卷积如图所示。4 第四题4.1 实验目的使用filter函数，并与上题比较结果是否一致用freqz函数求差分方程的频率响应用freqs函数画出微分方程的频率响应4.2 实验程序4.2.1 程序1h=0 1 2 3 4 5;x2=1 1 1 1 1 1 0 0 0 0 0;y2=filter(h,1,x2);stem(

5、y2)4.2.2 程序2b=2 0 -1;a=1 -0.8;H omega=freqz(b,a,4)4.2.3 程序3a=1 3;b=3;figure;freqs(b,a);w=linspace(0,3*pi);h=freqs(b,a,w);figure;plot(w,abs(h);4.3 实验图像或输出4.3.1 程序14.3.2 程序2 H = 5.0000 2.8200 - 1.3705i 1.8293 - 1.4634i 0.9258 - 0.9732iomega = 0 0.7854 1.57082.35624.3.3 程序34.4 实验结论用filter函数所得结果，与上题一致用f

6、reqz函数求差分方程的频率响应如上用freqs函数画出微分方程的频率响应图如上5 第五题5.1 实验目的计算离散时间傅里叶级数。5.2 实验程序n=0:31;N=32;x=1 1 zeros(1,30);a=(1/N)*fft(x);subplot(2,1,1);stem(n,real(a)title(Real part of a)xlabel(n)ylabel(real(a)subplot(2,1,2);stem(n,imag(a)title(Image part of a)xlabel(n)ylabel(image(a)5.3 实验图像5.4 实验结论函数实部与虚部的图像如图所示。6 第

7、六题6.1 实验目的验证LTI系统的特征函数6.2 实验程序n=-20:100;x1=exp(i*0.25*pi*n);x2=sin(pi*n/8+pi/16);x3=power(0.9,n);x4=n+1;figuresubplot(2,1,1)stem(real(x1),.)subplot(2,1,2) stem(imag(x1),.)figurestem(x2,.)figurestem(x3,.)figuren=0:100;x1=exp(i*0.25*pi*n);x2=sin(pi*n/8+pi/16);x3=power(0.9,n);x4=n+1;a=1,-0.25;b=1,0.9;y

8、1=filter(b,a,x1);y2=filter(b,a,x2);y3=filter(b,a,x3);y4=filter(b,a,x4);figuresubplot(2,1,1)stem(real(y1),.)subplot(2,1,2)stem(imag(y1),.)figurestem(y2,.)figurestem(y3,.)figurestem(y4,.)stem(x4,.) h1=y1./x1;h2=y2./x2;h3=y3./x3;h4=y4./x4;figure subplot(2,1,1)stem(real(h1),.)title(the real part of h1)x

9、label(n)ylabel(h1)subplot(2,1,2)stem(imag(h1),.)title(the imag part of h1)xlabel(n)ylabel(h1)figurestem(h2,.)xlabel(n)ylabel(h2)figurestem(h3,.)xlabel(n)ylabel(h3)figurestem(h4,.)xlabel(n)ylabel(h4)6.3 实验图像x1x2x3x4y1y2y3y4h1h2h3h46.4 实验结论x1,x2,x3,x4,y1,y2,y3,y4,h1,h2,h3,h4的图像分别如图所示。7 第七题7.1 实验目的用离散时

10、间傅里叶级数综合信号7.2 实验程序x1=ones(1,8) ones(1,8) ones(1,8) ones(1,8) ones(1,8) ones(1,8) ones(1,8) ones(1,8);x2=ones(1,8) zeros(1,8) ones(1,8) zeros(1,8) ones(1,8) zeros(1,8) ones(1,8) zeros(1,8) ones(1,8) zeros(1,8);x3=ones(1,8) zeros(1,24) ones(1,8) zeros(1,24);stem(x1)title(singal x1)xlabel(n);ylabel(x1)

11、stem(x2,.)title(singal x2)xlabel(n);ylabel(x2)stem(x3,.)title(singal x3)xlabel(n);ylabel(x3)x1=ones(1,8);x2=ones(1,8) zeros(1,8);x3=ones(1,8) zeros(1,24);a1=1/8*fft(x1);a2=1/16*fft(x2);a3=1/64*fft(x3);stem(abs(a1)title(a1 Amplitude)xlabel(k);ylabel(abs(a1);stem(abs(a2)title(a2 Amplitude)xlabel(2);ylabel(abs(a2);stem(abs(a3)title(a3 Amplitude)xlabel(k);ylabel(abs(a3);7.3 实验图像