summaryrefslogtreecommitdiff
path: root/3432/CH6/EX6.7
diff options
context:
space:
mode:
authorprashantsinalkar2017-10-10 12:27:19 +0530
committerprashantsinalkar2017-10-10 12:27:19 +0530
commit7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch)
treedbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3432/CH6/EX6.7
parentb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff)
downloadScilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2
Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip
initial commit / add all books
Diffstat (limited to '3432/CH6/EX6.7')
-rw-r--r--3432/CH6/EX6.7/EX6_7_f0.pdfbin0 -> 34768 bytes
-rw-r--r--3432/CH6/EX6.7/Ex6_7.sce44
2 files changed, 44 insertions, 0 deletions
diff --git a/3432/CH6/EX6.7/EX6_7_f0.pdf b/3432/CH6/EX6.7/EX6_7_f0.pdf
new file mode 100644
index 000000000..9e46a121d
--- /dev/null
+++ b/3432/CH6/EX6.7/EX6_7_f0.pdf
Binary files differ
diff --git a/3432/CH6/EX6.7/Ex6_7.sce b/3432/CH6/EX6.7/Ex6_7.sce
new file mode 100644
index 000000000..0ea6d5b7f
--- /dev/null
+++ b/3432/CH6/EX6.7/Ex6_7.sce
@@ -0,0 +1,44 @@
+//Example 6.7
+//Computation of velocity error constant Kv from Bode plot
+
+xdel(winsid())//close all graphics Windows
+clear;
+clc;
+//------------------------------------------------------------------
+//System transfer function and its bode plot
+K=10;
+s=poly(0,'s');
+Gs=syslin('c',(K)/(s*(s+1)));
+//The bode plot of the system
+
+fmin=0.01/2/%pi; //mininmum frq. in Hz for response (0.1 rad/sec)
+fmax=10/2/%pi; //maximum frq. in Hz for response (100 read/sec)
+//------------------------------------------------------------------
+//Bode plot for frequency in Hz (scilab ver. 5.4.1)
+//bode(g,fmin,fmax);
+//OR
+//Bode plot for frequency in rad/sec (scilab ver. 5.5.1)
+bode(Gs,fmin,fmax,0.01,"rad")
+title(['Determination of Kv from the Bode plot for the system',...
+'$10/[s(s+1)]$'],'fontsize',3)
+//choose frequency (rad) and magnitude from bode plot and calculate Kv
+//Here at w=0.01, magngitude in db is M=60
+//i.e actual magnitude of the reponse is |A|=10^(M/20)
+w=0.01; // in rad
+M=60 // in db
+A=10^(M/20) //actual gain
+
+//Velocity error constant Kv=w*|A(w)|
+Kv=w*A;
+disp(Kv,"The Velocity error Constant from bode plot is: ")
+//------------------------------------------------------------------
+// Computation of the Kv
+[frq repf]=repfreq(Gs,fmin,fmax);
+//frq in Hz, repf is freq. response in rectangular form.
+//From bode plot, Kv=w*|A(w)|
+//i.e Kv=2*pi*f*|A(2*pi*f)|
+
+idx=1;//selecting the frequency and response at that frequency from arrays
+Kv=2*%pi*frq(idx)*abs(repf(idx))
+disp(Kv,"The Velocity error Constant is computed at 0.0015915 Hz (0.01 rad/sec) : ")
+//------------------------------------------------------------------