diff options
Diffstat (limited to '257/CH11/EX11.23/example_11_23.sce')
-rw-r--r-- | 257/CH11/EX11.23/example_11_23.sce | 20 |
1 files changed, 20 insertions, 0 deletions
diff --git a/257/CH11/EX11.23/example_11_23.sce b/257/CH11/EX11.23/example_11_23.sce new file mode 100644 index 000000000..1f8344570 --- /dev/null +++ b/257/CH11/EX11.23/example_11_23.sce @@ -0,0 +1,20 @@ +s=poly(0,'s');
+F=syslin('c',[1/((1+s*0.02)*s*(0.05*s+1))]) //without k
+
+fmin=0.1; //Min freq in Hz
+fmax=20; //Max freq in Hz
+
+scf(1);clf;
+bode(F,fmin,fmax); //Plots frequency response of open-loop system in Bode diagram
+
+show_margins(F) //display gain and phase margin and associated crossover frequencies
+
+// TO FIND VALUE OF K
+
+// for GM=10 dB, the point on the plot without k is 26dB away from 0dB line.
+k1=10^(26/20) // 20*log(k1)=26
+disp(k1,"k for GM=30 is ")
+F1=syslin('c',[(k1)/((1+s*0.02)*s*(0.05*s+1))])
+
+[PhaseMargin,freqPM]=p_margin(F1)
+disp(PhaseMargin,"PM=")
|