matlab的freqs函数画幅频、相频特性曲线.docVIP

一 clear a=[1]; b=[1 0.4 0.08]; w=linspace(-2*pi,2*pi,512); h=freqs(a,b,w); subplot(2,1,2),plot(w/pi,angle(h)); subplot(2,1,1),plot(w/pi,abs(h)); 二 clear a=[2 0]; b=[1 2,10001]; w=linspace(-2*pi,2*pi,512); h=freqs(a,b,w); subplot(211),plot(w/pi,abs(h)); subplot(212),plot(w/pi,angle(h)); clear syms t; v1=(1+cos(t))*cos(100*t); V=laplace(v1); syms s; H=2*s/(s^2+2*s+10001); V2=V*H; v2=ilaplace(V2); ezplot(v2,[-0.1*pi,0.1*pi])