G(S)=100/S(0.04S+1),Increase the Closed-loop Phase margin:PM=28.0243 deg at Wcp=46 rad/sec Consider the Lead Compensator Gc(S)=(1+0.0262S)/(1+0.0106S),Draw the Bode Plot?
>> G=tf(100,[0.04,1,0])
Transfer function:
100
------------
0.04 s^2 + s
>> [Gm,Pm,Wcg,Wcp]=margin(G);
>> Gc1=tf([0.0262,1],[0.0106,1]);
>> GGc1=G*Gc1;
>> [Gm1,Pm1,Wcg1,Wcp1]=margin(GGc1);
>> figure(1),
>> bode(G,GGc1);
>> grid off;
>> bode(G,GGc1);
MATLAB CODE FOR DESIGN OF LEAD COMPENSATOR:
>> G=tf(100,[0.04,1,0])
Transfer function:
100
------------
0.04 s^2 + s
>> [Gm,Pm,Wcg,Wcp]=margin(G);
>> Gc1=tf([0.0262,1],[0.0106,1]);
>> GGc1=G*Gc1;
>> [Gm1,Pm1,Wcg1,Wcp1]=margin(GGc1);
>> figure(1),
>> bode(G,GGc1);
>> grid off;
>> bode(G,GGc1);
