信号系统实验.docVIP

实验一： (1.1): 方波，T=4,A=1 t=(0:0.02:10); y=square(0.5*pi*t); plot(t,y) (1.2):时移，反转，尺度变换 syms t; y=sym(square(0.5*pi*t)); y1=subs(f,t,t-1); y2=subs(f,t,-t); y3=subs(f,t,0.5*t-1); subplot(2,2,1);ezplot(y,[-8,8]);grid on; subplot(2,2,2);ezplot(y1,[-8,8]);grid on; subplot(2,2,3);ezplot(y2,[-8,8]);grid on; subplot(2,2,4);ezplot(y3,[-8,8]);grid on; (2.1) f1=exp(-0.5*t),f2=cos(2*pi*t)相加、相乘 t=0:0.01:2; f1=exp(-0.5*t); f2=cos(2*pi*t); f3=f1+f2; f4=f1.*f2; subplot(2,2,1);plot(t,f1);title(f1(t)); subplot(2,2,2);plot(t,f2);title(f2(t)); subplot(2,2,3);plot(t,f3);title(f1+f2); subplot(2,2,4);plot(t,f4);title(f1*f2); x1=2*exp(-0.5*n); x2=[1,2,3,4,0]; x1=2*exp(-0.5*n); x2=[1,2,3,4,0]; x1=2*exp(-0.5*n); x2=[1,2,3,4,0]; (2.2) x1=2*exp(-0.5*n)，x2=[1,2,3,4,0]相加、相乘 n=0:4; x1=2*exp(-0.5*n); x2=[1,2,3,4,0]; m=1:length(x2); subplot(2,2,1);stem(n,x1); subplot(2,2,2);stem(m,x2); subplot(2,2,3);stem(n,x1+x2); subplot(2,2,4);stem(n,x1.*x2); 实验二 （1） y1=[ones(1,10),zeros(1,10)]; y2=[ones(1,24),zeros(1,24)]; y=conv(y1,y2); n1=1:length(y1); n2=1:length(y2); L=length(y) subplot(3,1,1);plot(n1,y1);axis([1,L,0,2]); subplot(3,1,2);plot(n2,y2);axis([1,L,0,2]); n=1:L; subplot(3,1,3);plot(n,y);axis([1,L,0,20]); （3） y1=[ones(1,10),zeros(1,10)]; t=0:0.1:5; y2=exp(-1*t); y=conv(y1,y2); n1=1:length(y1); n2=1:length(y2); L=length(y) subplot(3,1,1);plot(n1,y1);axis([1,L,0,2]); subplot(3,1,2);plot(t,y2); n=1:L; subplot(3,1,3);plot(n,y);axis([1,L,0,20]); 实验三 已知频率特性函数为：，求其幅频特性和相频特性。 w=linspace(0,5,400); b=[2 0 1 4]; a=[1 3 2 5 2]; H=freqs(b,a,w); subplot(2,1,1); plot(w,abs(H)); set(gca,xtick,[0 1 2 3 4]); set(gca,ytick,[0 0.4 0.707 1]); xlabel(\omega); ylabel(|H(j\omega)); title(幅频特性曲线,FontSize,15); subplot(2,1,2); plot(w,angle(H)); set(gca,xtick,[0 1 2 3 4 5]); xlabel(\omega); ylabel(\phi(\omega)); title(相频特性曲线,FontSize,15); 结果： 实验四 设设（1）p1=-2,p2=-30;(2)p1=-2,p3=3 a：针对极点参数（1）（2），画出系统零、极点分布图，判断该系统稳定性。 （1） %程序部分,当p1=-2,p2=-30时 num=[1,0]; den=[1,32,60];%这些可以运用MATLAB中simple对分母多项式进行简化求得 [r,pole