内蒙古科技大学信号与系统实验2,3.doc

实验2作业 clear;close all; t0=-2;t1=4;dt=0.01; t=t0:dt:t1; u0=(t=0); u1=((t-1)=0); u2=((t-2)=0); u3=((t-3)=0); x=u0-u2; h=u1-u3; y=dt*conv(x,h); subplot(221) plot(t,x),grid on,title(Signal x(t)),axis([t0,t1,-0.2,1.2]) subplot(222) plot(t,h),grid on,title(Signal h(t)),axis([t0,t1,-0.2,1.2]) subplot(212) t=2*t0:dt:2*t1; plot(t,y),grid on,title(The convolution of x(t) and h(t)), axis([2*t0,2*t1,-0.1,3]), xlabel(Time t sec) 实验2例题 clear ;close all; t0=-2;t1=4;dt=0.01; t=t0:dt:t1; x=u(t)-u(t-1); h=t.*(u(t)-u(t-1)); y=dt*conv(x,h); subplot(221) plot(t,x),grid on,title(Signal x(t)),axis([t0,t1,-0.2,1.2]) subplot(222) plot(t,h),grid on,title(Signal h(t)),axis([t0,t1,-0.2,1.2]) subplot(212) t=2*t0:dt:2*t1; plot(t,y),grid on,title(The convolution of x(t) and h(t)),axis([2*t0,2*t1,-0.1,0.6]), xlabel(Time t sec) 实验3: 系统1带通滤波器 %program3_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 ,title(Magnitude response), subplot(223) plot(w,phai),grid on ,title(Phase response), xlabel(Frequency in rad/sec) subplot(222) plot(w,Hr),grid on ,title(Real part of frequency response), xlabel(Frequency in rad/sec) subplot(224) plot(w,Hi),grid on ,title(Imaginary part of frequency response), xlabel(Frequency in rad/sec) 系统2 全通滤波器 %program3_1 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 ,title(Magnitude response), subplot(223) plot(w,phai),grid on ,title(Phase response), xlabel(Frequency in rad/sec) subplot(222) plot(w,Hr),grid on ,title(Real part of frequency response), xlabel(Frequency in rad/sec) subplot(224) plot(w,Hi),grid on ,title(Imaginary part of frequency response), xlabel(Frequency in rad/sec) 系统3 低通滤波器 %program3_1 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=imag(H); subplot(221) plot(w,Hm),grid on ,title(Magnitude response), subplot(223) plot(w,phai),grid on ,title(Phase response), xlabel(Freq