西南交大信号与系统实验报告.docVIP

西南交通大学 信号与系统实验报告 班级 ____________ 学号 ____________ 姓名 ___________ 实验时间 ___ 年___ 月___ 日 实验名称__________________________________________________________ 预 习： 操 作： 报 告： 总 评： 指导老师： 一、实验预习内容 1、实验目的： 2、实验内容及实验函数说明： 电气工程基础实验中心 二、实验内容 1.2 t=-3:0.01:3; y=rectpuls(t,2); plot(t,y); axis([-5,5,-0.5,1.5]); 1.3 t=-1:0.0001:1; y1=5*rectpuls(t,0.2); y2=10*rectpuls(t,0.1); plot(t,y1,t,y2) axis([-1,1,-1,12]); 1.4（1） t=-3:0.01:3; y=5*cos(12*pi*t)+6*cos(18*pi*t); plot(t,y); axis([0,1,-9,11]); 1.4（2） syms t; y=cos(2*pi*t)*cos(2*pi*t); ezplot(y,[0,1]); axis([0,1,0,1]); 4.2 t=0:0.001:6; a1=[1 5 6]; a2=[1 2 5]; a3=[1 2 1]; b=[1 5]; y1=step(b,a1,t); y2=step(b,a2,t); y3=step(b,a3,t); plot(t,y1,t,y2,+,t,y3,*); 4.5 a=[1 5 6]; b=[2 8]; [A,B,C,D]=tf2ss(b,a); sys=ss(A,B,C,D); t=0:0.01:5; f=zeros(1,length(t)); x0=[0,-3]; subplot(1,3,1);plot(t,y1); y1=lsim(sys,f,t,x0); f=exp(-t); subplot(1,3,2);plot(t,y2); y2=lsim(sys,f,t); subplot(1,3,3);plot(t,y3); y3=lsim(sys,f,t,x0); 5.3 Ts=0.001; t=-1:Ts:5; f=(t0)-(t1); h=0.5*((t1)-(t3)); y=conv(f,h)*Ts; L=length(y); k=(0:L-1)*Ts+1; subplot(2,1,2); plot(k,y);grid on; axis([0,10,0,0.7]); %6.1 syms t T n; y1=4/T*t+1; y2=-4/T*t+1; a0=1/T*(int(y1,-T/2,0)+int(y2,0,T/2)); y3=y1*cos(n*2*pi/T*t); y4=y2*cos(n*2*pi/T*t); an=2/T*(int(y3,-T/2,0)+int(y4,0,T/2)); bn=0; an=zeros(50,1); cn=zeros(50,1); T=2*pi; an(1)=0; cn(1)=0; for(i=1:1:50) an(i+1)=-2*(-2+2*cos(i*pi)+i*pi*sin(i*pi))/i^2/pi^2; cn(i+1)=abs(an(i+1)); end syms k; k=0:1:50; stem(k,cn); subplot(2,1,1); plot(k,cn); t=-5:0.01:5; f=sawtooth((t-pi),0.5); subplot(2,1,2); plot(t,f);1 %7.3 a=[0.08 0.4 1]; b=[1]; w=linspace(0,10,300); [h,w]=freqs(b,a,w); subplot(2,1,1); plot(w,abs(h)); grid on; subplot(2,1,2); plot(w,angle(h)); grid on; %7.1 syms t; f=exp(-2*abs(t)); subplot(3,1,1) ezplot(f); y1=fourier(f); y=abs(y1); subplot(3,1,2); ezplot(y); im=imag(y1); re=real(y1); phase=atan(im/re); subplot(3,1,3); ezplot(phase); %7.2 syms tao w t; tao=1; F=(tao*sin(w*tao/2/pi))/(w*tao/2); subplot(2,2,1);ezplot(F);