Page 157 - MATLAB an introduction with applications

142 ——— MATLAB: An Introduction with Applications

>> step(Ta)
>> title(‘(a)’)
>> ‘(b)’
ans =
(b)
>> numb=.045;
>> denb=[1 .025 .045];
>> Tb=tf(numb,denb)
Transfer function:
0.045
s ∧ 2 + 0.025s + 0.045
>> omeganb=sqrt(denb(3))
omeganb=0.2121
>> zetab=denb(2)/(2*omeganb)
zetab=0.0589

>> Tsb=4/ (zetab*omeganb)
Tsb = 320
>> Tpb=pi/ (omeganb*sqrt(1–zetab^2))
Tpb = 14.8354
>> Trb=(1.76*zetab^3 –.417*zetab^2 + 1.039*zetab + 1)/ omeganb
Trb = 4.9975
>> percentb=exp(–zetab*pi/ sqrt(1–zetab^2))*100
percentb = 83.0737
>> subplot(222)
>> step(Tb)
>> title(‘(b)’)
>> ‘(c)’
ans =
(c)
>> numc=10E8;
>> denc=[1 1.325*10E3 10E8];
>> Tc=tf(numc, denc)
Transfer function:
1e009
∧
s 2 + 13250s + 1e009

F:\Final Book\Sanjay\IIIrd Printout\Dt. 10-03-09