blob: 51822c2305bd78a4b78f5a6806f294b717d56b5d (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
s=%s
syms K J f
K=60; // given
J=10; // given
p=K/J
q=K/J+(f/J)*s+s^2
G=p/q;
disp(G,"Qo(s)/Qi(s)=")
zeta=0.3; // given
cof1=coeffs(q,'s',0)
// on comparing the coefficients
Wn=sqrt(cof1)
cof2=coeffs(q,'s',1)
// 2*zeta*Wn=cof2
f/J=2*zeta*Wn
r=s^2+f/J
s=s^2+f/J+K/J
H=r/s;
disp(H,"Qe(s)/Qi(s)=")
|
