相关代码 |
|
w1s=2;w1p=1;N=5; w2s=0.5;w2p=1; for i=1:3 if m==1 for i=0:1:30 w1=0:0.1:i*0.1; if i==0 w2=3.2; elseif i==1 w2(2)=3.1 elseif i==2 w2(3)=3.0 elseif i==3 w2(4)=2.8 else w2(i+1)=1/w1(i+1); end G=1./((j.*w1).^5+3.236.*(j.*w1).^4+5.2361.*(j.*w1).^3+5.2361.*(j.*w1).^2+3.2361.*(j.*w1)+1); G1=abs(G); subplot(2,1,1) plot(w1,G1); axis([0,3.5,0,1.2]) subplot(2,1,2) plot(w2,G1) axis([0,3.5,0,1.2]) pause(1) end