资源描述：
《武汉理工大学信号与系统实验答案》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、1-1-1clearclosealldt=0.01;(dt=0.2;)t=-2:dt:2;x=sin(2*pi*t);plot(t,x)title('sinusoidalsingnalx(t)')xlabel('Timet(sec)')1-2t=-5:0.001:5;x=exp(-0.5*t).*u(t);y=exp(-0.5*(1.5*t+3)).*u(1.5*t+3);subplot(211)plot(t,x)title('x(t)')subplot(212)plot(t,y)title('y(t)=x(1,5*t+3)')1-3clear;closeall;t0=0
2、;t1=3.5;dt=0.001;t=t0:dt:t1;x1=u(t)-u(t-2);x2=u(t-1)-u(t-3);y=dt*conv(x1,x2);subplot(221)plot(t,x1)gridontitle('Signalx1(t)')axis([t0,t1,0,1.5])subplot(222)plot(t,x2)gridontitle('signalx2(t)')axis([t0,t1,0,1.5])subplot(212)t=t0:dt:2*t1;plot(t,y)gridontitle('Theconvolutionofx1(t)andx2(t)')
3、axis([2*t0,2*t1,0,2.5])xlabel('Timet(sec)')+2-1clearcloseallnum=input('Typeintherightcoeffcientvectorofdifferentialequation:');den=input('Typeintheleftcoefficientvectorofdifferentialequation:');t=0:0.01:10;x=input('Typeintheexpressionoftheinputsignalx(t):');lsim(num,den,x,t)运行过程：Typeinthe
4、rightcoeffcientvectorofdifferentialequation:8Typeintheleftcoefficientvectorofdifferentialequation:[1,3,2]Typeintheexpressionoftheinputsignalx(t):(exp(-2*t)-exp(-3*t)).*u(t)3-1clearcloseallT=2;dt=0.00001;t=-2:dt:2;x1=(u(t+1)-u(t-1-dt)).*abs(t);x=0;form=-1:1x=x+u(t-m*T)-u(t-1-m*T-dt);endw0=
5、2*pi/T;N=10;L=2*N+1;fork=-N:N;ak(N+1+k)=(1/T)*x1*exp(-j*k*w0*t')*dt;endphi=angle(ak);运行结果：ak=Columns1through60.0000-0.0000i-0.0025-0.0000i0.0000+0.0000i-0.0041+0.0000i0.0000+0.0000i-0.0081-0.0000iColumns7through120.0000-0.0000i-0.0225+0.0000i0.0000-0.0000i-0.2026+0.0000i0.5000-0.2026-0.00
6、00iColumns13through180.0000+0.0000i-0.0225-0.0000i0.0000+0.0000i-0.0081+0.0000i0.0000-0.0000i-0.0041-0.0000iColumns19through210.0000-0.0000i-0.0025+0.0000i0.0000+0.0000i3-2b=input('b=');a=input('a=');[H,w]=freqs(b,a);Hm=abs(H);phai=angle(H);hr=real(H);Hi=imag(H);subplot(221)plot(w,Hm)grid
7、ontitle('Magnituderesponse');xlabel('Frequencyinrad/sec');subplot(223)plot(w,phai)gridontitle('Phaseresponse')xlabel('Frequencyinrad/sec');subplot(222)plot(w,hr),gridontitle('Realpartoffrequencyresponse')xlabel('Frequencyinrad/sec');subplot(224)plot(w,Hi),gridontitl
