blob: 82784ffddff2f7ca7ca7b0514322752fc4ab154e (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
s=%s;
syms Kp Ka Kt J f
// given
J=0.4;
Kp=0.6;
Kt=2;
f=2;
p=Kp*Ka*Kt
q=s^2+f/J+(Kp*Ka*Kt)/J
G=p/q;
disp(G,"Qm(s)/Qr(s)=")
cof_1=coeffs(q,'s',0)
// on comparing the coefficients
// Wn=sqrt(cof_1)
Wn=10;
Ka=(Wn)^2*J/(Kp*f)
disp(Ka,"Amplifier Constant=")
|