Fz2=(1/2)fp1; Fp2=fs; Fp3=fs;
[marg_G0,phase_G0]=bode(G0,fg*2*pi); Marg_G=1/marg_G0; AV1=fz2/fg*marg_G; AV2=fp2/fg*marg_G; R2=10*10^3; R3=R2/AV2;
C1=1/(2*pi*fz1*R2); C3=1/(2*pi*fzp2*R3); C2=1/(2*pi*fp3*R2); R1=1/(2*pi*C3*fz1);
Num=conv([C1*R2 1],[(R2+R3)*C3 1]); Den1=conv([(C1+C2)*R1 0],[R3*C3 1]); Den=conv(den1,[R2*C1*C2/(C1+C2) 1]); Gc=tf(num,den); Figure(2); Bode(Gc);
G=series(Gc,G0); Figure(3) Margin(G)
11