信号与系统课程设计实验报告.doc

《信号与系统课程设计实验报告.doc》由会员分享，可在线阅读，更多相关《信号与系统课程设计实验报告.doc（32页珍藏版）》请在三一文库上搜索。

1、 北京信息科技大学实验报告封面实 验 报 告课程名称： 信号与系统 系 （院）： 光电学院 专 业： 光信息科学与技术 班 级： 光信0801 学生姓名： 靖鑫 学 号： 2008010750 指导教师： 开课时间： 2009-2010学年 二 学期第一次上机 实验一 信号与系统的时域分析第六题dt=0.01;t=-2:dt:2;x=sin(2*pi*t);subplot(211)plot(t,x)title(sinusoidal sign x(t)xlabel(time t(see)dt=0.2t=-2:dt:2;x=sin(2*pi*t);subplot(212)plot(t,x)titl

2、e(sinusoidal sign x(t)xlabel(time t(see)第七题n=0:0.1:4,x1=(-2).n.*u(n),subplot(211),stem(n,x1),axis(0 4 -8 16)grid on,title(-2).n.*u(n)xlabel(time index);t=0:0.01:4;x2=cos(2*pi*t).*(u(t)-u(t-3),subplot(212),plot(t,x2),axis(0 4 -1 1)grid on,title(cos(2*pi*t).*u(t)-u(t-3)xlabel(time index)第八题t=-2:0.01:2

3、0x=exp(-0.5*t).*u(t)subplot(211),plot(t,x),grid on,y=exp(-0.5*(1.5*t+3).*u(1.5*t+3)subplot(212),plot(t,y),grid on,第九题t=0:0.01:10x=(exp(-2*t)-exp(-3*t).*u(t)den=1 3 2;num=8;lsim(num,den,x,t)第十题t0=-2t1=4dt=0.01t=t0:dt:t1x=exp(-2*t).*(u(t+1)-u(t-2)h=cos(2*pi*t).*(u(t)-u(t-3)y=dt*conv(x,h);subplot(311)p

4、lot(t,x),grid on,title(signal x(t),subplot(312)plot(t,h),grid on,title(signal h(t),subplot(313)t=2*t0:dt:2*t1plot(t,y),grid on,title(The convolution of x(t) and h(t)xlabel(Time t sec) 北京信息科技大学实验报告封面实 验 报 告课程名称： 信号与系统 系 （院）： 光电学院 专 业： 光信息科学与技术 班 级： 光信0801 学生姓名： 靖鑫 学 号： 2008010750 指导教师： 开课时间： 2009-201

5、0学年 二 学期第二次上机 实验二 连续时间信号的频率分析第六题t=0:0.01:8w=0.5*pix1=cos(w*t)subplot(411)plot(t,x1)grid ontitle(cos(wt)x2=cos(3*w*t)subplot(412)plot(t,x2)grid ontitle(cos(3wt)x3=cos(5*w*t)subplot(413)plot(t,x3)grid ontitle(cos(5wt)x4=cos(w*t)-(1/3)*cos(3*w*t)+(1/5)*cos(5*w*t)subplot(414)plot(t,x4)grid ontitle(x)第七题

6、（1）T=2;dt=0.0001;t=-4:dt:4;x1=(t+1).*(u(t+1)-u(t)+(1-t).*(u(t)-u(t-1)x=0;for m=-2:2 x=x+(t-m*T+1).*(u(t-m*T+1)-u(t-m*T)+(1-t+m*T).*(u(t-m*T)-u(t-m*T-1);end w0=2*pi/T;N=10;L=2*N+1;for k=-N:1:N; ak(N+1+k)=(1/T)*x1*exp(-j*k*w0*t)*dt;endphi=angle(ak)y=0for q=1:L; y=y+ak(q)*exp(j*(-(L-1)/2+q-1)*2*pi*t/T)

7、;end;subplot(221),plot(t,x),title(The original signal x(t), axis(-4,4,-0.2,1.2),subplot(223),plot(t,y),title(The synthesis signal y(t), axis(-4,4,-0.2,1.2), xlabel(Time t),subplot(222),k=-N:N; stem(k,abs(ak),k.),title(The amplitude ak of x(t),axis(-N,N,-0.1,0.6),subplot(224)stem(k,phi,r.),title(The

8、phase phi(k) of x(t),axis(-10,10,-2,2),xlabel(Index k)第七题（2）T=2;dt=0.0001;t=-5:dt:5;x1=u(t+0.2)-u(t-0.2);x=0;for m=-2:2 x=x+u(t-m*T+0.2)-u(t-m*T-0.2);end w0=2*pi/T;N=10;L=2*N+1;for k=-N:1:N; ak(N+1+k)=(1/T)*x1*exp(-j*k*w0*t)*dt;endphi=angle(ak)y=0for q=1:L; y=y+ak(q)*exp(j*(-(L-1)/2+q-1)*2*pi*t/T);e

9、nd;subplot(221),plot(t,x),title(The original signal x(t), axis(-5,5,-0.2,1.2),subplot(223),plot(t,y),title(The synthesis signal y(t), axis(-5,5,-0.2,1.2), xlabel(Time t),subplot(222),k=-N:N; stem(k,abs(ak),k.),title(The amplitude ak of x(t),axis(-N,N,-0.1,0.6),subplot(224)stem(k,phi,r.),title(The ph

10、ase phi(k) of x(t),axis(-10,10,-2,2),xlabel(Index k)第八题T=2;dt=0.0001;t=-4:dt:4;x1=(t+2).*(u(t+2)-u(t+1)+(u(t+1)-u(t-1)+(-t+2).*(u(t-1)-u(t-2);x=0;w0=2*pi/T;N=10;L=2*N+1;for k=-N:1:N; ak(N+1+k)=(1/T)*x1*exp(-j*k*w0*t)*dt;endphi=angle(ak);y=0;for q=1:L; y=y+ak(q)*exp(j*(-(L-1)/2+q-1)*2*pi*t/T);end;sub

11、plot(311),plot(t,x1),title(The original signal x(t),axis(-2,2,0,1.1);subplot(312),k=-N:N;stem(k,abs(ak),k.);title(The amplitude |ak| of x(t);axis(-N,N,-1,1.5);subplot(313),stem(k,phi,r.);title(The phase phi(k) of x(t);axis(-N,N,-4,4),xlabel(Index k); 北京信息科技大学实验报告封面实 验 报 告课程名称： 信号与系统 系 （院）： 光电学院 专 业：

12、 光信息科学与技术 班 级： 光信0801 学生姓名： 靖鑫 学 号： 2008010750 指导教师： 开课时间： 2009-2010学年 二 学期第三次上机实验三 连续时间LIT系统的频域分析第一题 （1）b=1 0;a=1 1 25;H,w=freqs(b,a);Hm=abs(H);phai=angle(H);Hr=real(H);Hi=imag(H);subplot(221)plot(w,Hm),grid on,subplot(223)plot(w,phai),grid on,subplot(222)plot(w,Hr),grid on,subplot(224)plot(w,Hi),g

13、rid on,（2）b=1 -1;a=1 1;H,w=freqs(b,a);Hm=abs(H);phai=angle(H);Hr=real(H);Hi=imag(H);subplot(221)plot(w,Hm),grid on,subplot(223)plot(w,phai),grid on,subplot(222)plot(w,Hr),grid on,subplot(224)plot(w,Hi),grid on,（3）b=262 ;a=1 10 48 148 306 401 262;H,w=freqs(b,a);Hm=abs(H);phai=angle(H);Hr=real(H);Hi=i

14、mag(H);subplot(221)plot(w,Hm),grid on,subplot(223)plot(w,phai),grid on,subplot(222)plot(w,Hr),grid on,subplot(224)plot(w,Hi),grid on,第二题b=262;a=1 10 48 148 306 401 262;T=0.01;dw=0.1t=0:0.01:40w=-4*pi:dw:4*pih=impulse(b,a,40)x=sin(t)+sin(8*t);y=lsim(b,a,x,t)subplot(321)plot(t,x),grid on,subplot(323)p

15、lot(h),grid on,subplot(325)plot(y),grid on,X=x*exp(-j*t*w)*Tx1=abs(X);subplot(322)plot(x1),grid on,H=freqs(b,a,w)h1=abs(H)subplot(324)plot(h1),grid on,z=x1.*h1subplot(326)plot(z),grid on,第三题T=0.01;dw=0.1;t=-10:T:10;w=-4*pi:dw:4*pi;b=1 0;a=1 1 25;tao=grpdelay(b,a,w);subplot(311);plot(w,tao);b=1 -1;a=

16、1 1;tao=grpdelay(b,a,w);subplot(312);plot(w,tao);b=262;a=1 10 48 148 306 401 262;tao=grpdelay(b,a,w);subplot(313)plot(w,tao) 北京信息科技大学实验报告封面实 验 报 告课程名称： 信号与系统 系 （院）： 光电学院 专 业： 光信息科学与技术 班 级： 光信0801 学生姓名： 靖鑫 学 号： 2008010750 指导教师： 开课时间： 2009-2010学年 二 学期第四次上机实验四 抽样定理第三题t=0:0.01:10Ts=1/4n=0:Ts:10;x=cos(0.

17、5*pi*t);xn=cos(0.5*pi*n);subplot(211);plot(t,x),title(A continous-time signal x(t),xlabel(Time t);subplot(212);stem(n,xn,.),title(The sampled version xn of x(t),xlabel(Time index n)第四题 tmax=4;dt=0.01;t=0:dt:tmax;Ts=1/2;ws=2*pi/Ts;w0=20*pi;dw=0.1;w=-w0:dw:w0;n=0:1:tmax/Ts;x=exp(-4*t).*u(t);xn=exp(-4*

18、n*Ts);subplot(221);plot(t,x),title(A continuous-time signal x(t),xlabel(Time t),axis(0,tmax,0,1),grid onsubplot(223)stem(n,xn,.),title(The sampled version xn of x(t),xlabel(Time index n),axis(0,tmax/Ts,0,1),grid onXa=x*exp(-j*t*w)*dt;X=0;for k=-8:8; X=X+x*exp(-j*t*(w-k*ws)*dt;endsubplot(222)plot(w,a

19、bs(Xa)title(Magnitude spectrum of x(t),grid onaxis(-60,60,0,1.8*max(abs(Xa)subplot(224)plot(w,abs(X)title(Magnitude spectrum of xn),xlabel(Frequency in radians/s),grid onaxis(-60,60,0,1.8*max(abs(Xa) 当Ts=1/2时当Ts=1/4时当Ts=1/8时当Ts=1/10时 实验八 离散时间系统的Z域分析第一题f=sym(cos(a*k);F=ztrans(f)f1=sym(ak);F1=ztrans(f

20、1)第二题b=0.2 0.1 0.3 0.1 0.2;a=1 -1.1 1.5 -0.7 0.3;figure(1);r,p,k=tf2zp(b,a);zplane(b,a);b1=3 -5 10 0;a1=1 5 2 -10;figure(2);r,p,k=tf2zp(b1,a1);zplane(b1,a1);第三题b=1 1 0;a=2 2 1;freqz(b,a)第四题b=1 5 5 1;a=1 0 0 0;figure(1);r,p,k=tf2zp(b,a)zplane(b,a)figure(2)subplot(312)freqz(b,a)第五题clearb=1 0 0;a=1 5 6;h=impz(b,a)stem(h)第七题b=5 -5;a=4 -1;freqz(b,a)