blob: f4c63a56c8dce49d3009e97546b7ddba41dcf7b7 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
//Caption:overall_transfer_function_and_root_locus
//example 7.24.10
//page 314
syms Ka Ke Kf Rf Lf eq N1 N2 N3 N4 N5
//where Ka=amplifier_gain; Ke=error_detector_gain; Kf=motor_torque_const; Rf=field_resistance; Lf=field_inductance Jeq=moment_of_inertia; feq=coeff_of_viscous_friction;
s=%s;
d=N1/N2;
e=N4/N3;
f=N4/N5;
n=N3/N5
Ke=0.05; Kf=2; Rf=10; Lf=2; Jeq=0.5*10^-4; feq=10^-4; d=0.1; e=5; f=0.5;n=0.2;
a=(1*Ke); //in series
b=Kf/(Rf+s*Lf);
c=1/(s*(Jeq*s+feq));
g= (b*c) //in series
h= (g*a) //in series
j= (h*0.02) //in series
k=j/(1+j*0.5);
a1=1/20;
cl=a1*k
disp(cl,"C(s)/R(s)=");
clf();
G=syslin('c',cl)
evans(G,200)
xgrid(2)
|
