summaryrefslogtreecommitdiff
path: root/1151/CH3
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /1151/CH3
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '1151/CH3')
-rwxr-xr-x1151/CH3/EX3.1/test.sce46
-rwxr-xr-x1151/CH3/EX3.2/example2.sce52
-rwxr-xr-x1151/CH3/EX3.3/example3.sce9
-rwxr-xr-x1151/CH3/EX3.4/example4.sce51
-rwxr-xr-x1151/CH3/EX3.5/example5.sce52
-rwxr-xr-x1151/CH3/EX3.6/example6.sce5
-rwxr-xr-x1151/CH3/EX3.8/example8.sce52
-rwxr-xr-x1151/CH3/EX3.9/example9.sce52
8 files changed, 319 insertions, 0 deletions
diff --git a/1151/CH3/EX3.1/test.sce b/1151/CH3/EX3.1/test.sce
new file mode 100755
index 000000000..a3b21be32
--- /dev/null
+++ b/1151/CH3/EX3.1/test.sce
@@ -0,0 +1,46 @@
+s=%s
+p= poly([0 0 50],'s','coeff');
+q= poly([0.1 2 10],'s','coeff');
+M=p/q
+ disp (M,"M( s )=")
+ H =1;
+ R =1;
+ b= coeff (p)
+a= coeff (q)
+// s t e p input
+if (a(1 ,1) ==b(1 ,1)) then
+ printf (" f o r u n i t s t e p input Ess=0 \nn" )
+ else
+ Ess =1/ H*(1 -(b(1 ,1)*H/a(1 ,1)))*R
+disp (Ess ," f o r u n i t s t e p input Ess=")
+ end
+// ramp input
+ c=0
+ for i =1:2
+if(a(1,i)-b(1,i)*H ==0) then
+c=c+1
+ end
+end
+if(c ==2)
+printf (" f o r u n i t ramp input Es s=0 \nn")
+else if(c ==1) then
+Ess =(a(1 ,2) -b(1 ,2)*H)/a(1 ,1)*H
+disp (Ess ," for unit ramp input Ess=")
+ else printf (" for unit ramp input Es s=infinite\n")
+ end
+ end
+// p a r a b o l i c input
+c=0
+for i =1:3
+if(a(1,i)-b(1,i)*H ==0) then
+c=c+1
+end
+end
+if(c ==3)
+printf (" f o r u n i t p a r a b o l i c input Es s=0 \n")
+ else if(c ==2) then
+Ess =(a(1 ,3) -b(1 ,3)*H)/a(1 ,1)*H
+ disp (Ess ," for unit parabolic input Ess=")
+ else printf (" for unit parabolic input Ess=infinite \n")
+ end
+ end
diff --git a/1151/CH3/EX3.2/example2.sce b/1151/CH3/EX3.2/example2.sce
new file mode 100755
index 000000000..df1d8eced
--- /dev/null
+++ b/1151/CH3/EX3.2/example2.sce
@@ -0,0 +1,52 @@
+s=%s
+p= poly([5 10 500],'s','coeff' );
+q= poly([1.8 25 10 0 0],'s','coeff' );
+M=p/q
+ disp (M,"M( s )=")
+ H =1;
+ R =1;
+ b= coeff (p)
+a= coeff (q)
+// s t e p input
+if (a(1 ,1) ==b(1 ,1)) then
+ printf (" f o r u n i t s t e p input Ess=0 \nn" )
+ else
+ Ess =1/ H*(1 -(b(1 ,1)*H/a(1 ,1)))*R;
+ kp=(1-Ess)/Ess;
+disp (Ess ," f o r u n i t s t e p input Ess=")
+disp(kp,"Kp=")
+ end
+// ramp input
+ c=0
+ for i =1:2
+if(a(1,i)-b(1,i)*H ==0) then
+c=c+1
+ end
+end
+if(c ==2)
+printf (" f o r u n i t ramp input Es s=0 \nn")
+else if(c ==1) then
+Ess =(a(1 ,2) -b(1 ,2)*H)/a(1 ,1)*H;
+kv=1/Ess;
+disp (Ess ," for unit ramp input Ess=")
+disp(kv,"Kv=")
+ else printf (" for unit ramp input Ess=infiite \n")
+ end
+ end
+// p a r a b o l i c input
+c=0
+for i =1:3
+if(a(1,i)-b(1,i)*H ==0) then
+c=c+1
+end
+end
+if(c ==3)
+printf (" f o r u n i t p a r a b o l i c input Ess=0 \nn")
+ else if(c ==2) then
+Ess =(a(1 ,3) -b(1 ,3)*H)/a(1 ,1)*H;
+ka=1/Ess;
+ disp (Ess ," for unit parabolic input Ess=")
+ diisp(ka,"Ka=")
+ else printf (" for unit parabolic input Ess=infinite \nn")
+ end
+ end
diff --git a/1151/CH3/EX3.3/example3.sce b/1151/CH3/EX3.3/example3.sce
new file mode 100755
index 000000000..6f63bc25b
--- /dev/null
+++ b/1151/CH3/EX3.3/example3.sce
@@ -0,0 +1,9 @@
+ // g i v e n G(s)H(s) =K/ s (s+20) H(s)=1 R(s)=1/s^2
+G =1+400/(s*(s+20));
+h=1/G;
+r=1/s^2;
+e= 0.05;
+disp (e," Steady state error=")
+//if ess=0.02 then k=?
+k=20/0.02;
+disp(k,"K=");
diff --git a/1151/CH3/EX3.4/example4.sce b/1151/CH3/EX3.4/example4.sce
new file mode 100755
index 000000000..09c3505fb
--- /dev/null
+++ b/1151/CH3/EX3.4/example4.sce
@@ -0,0 +1,51 @@
+s=%s
+p= poly([0 0 20],'s','coeff');
+q= poly([1 2 22 4 40 0], 's','coeff');
+ disp (M,"M( s )=")
+ H =1;
+ R =1;
+ b= coeff (p)
+a= coeff (q)
+// s t e p input
+if (a(1 ,1) ==b(1 ,1)) then
+ printf (" f o r u n i t s t e p input Ess=0 \nn" )
+ else
+ Ess =1/ H*(1 -(b(1 ,1)*H/a(1 ,1)))*R;
+ kp=(1-Ess)/Ess;
+disp (Ess ," f o r u n i t s t e p input Ess=")
+disp(kp,"Kp=")
+ end
+// ramp input
+ c=0
+ for i =1:2
+if(a(1,i)-b(1,i)*H ==0) then
+c=c+1
+ end
+end
+if(c ==2)
+printf (" f o r u n i t ramp input Ess=0 \nn")
+else if(c ==1) then
+Ess =(a(1 ,2) -b(1 ,2)*H)/a(1 ,1)*H;
+kv=1/Ess;
+disp (Ess ," for unit ramp input Ess=")
+disp(kv,"Kv=")
+ else printf (" for unit ramp input Ess=infinite \nn")
+ end
+ end
+// p a r a b o l i c input
+c=0
+for i =1:3
+if(a(1,i)-b(1,i)*H ==0) then
+c=c+1
+end
+end
+if(c ==3)
+printf (" for unit parabolic input Ess=0 \nn")
+ else if(c ==2) then
+Ess =(a(1 ,3) -b(1 ,3)*H)/a(1 ,1)*H;
+ka=1/Ess;
+ disp (Ess ," for unit parabolic input Ess=")
+ diisp(ka,"Ka=")
+ else printf (" for unit parabolic input Ess=infinite\nn")
+ end
+ end
diff --git a/1151/CH3/EX3.5/example5.sce b/1151/CH3/EX3.5/example5.sce
new file mode 100755
index 000000000..e3af87f96
--- /dev/null
+++ b/1151/CH3/EX3.5/example5.sce
@@ -0,0 +1,52 @@
+s=%s
+p= poly([0 0 108],'s','coeff');
+q= poly([1 7 24 48 0 0],'s','coeff');
+M=p/q
+ disp (M,"M( s )=")
+ H =1;
+ R =1;
+ b= coeff (p)
+a= coeff (q)
+// s t e p input
+if (a(1 ,1) ==b(1 ,1)) then
+ printf (" f o r u n i t s t e p input Ess=0 \nn" )
+ else
+ Ess =1/ H*(1 -(b(1 ,1)*H/a(1 ,1)))*R;
+ kp=(1-Ess)/Ess;
+disp (Ess ," f o r u n i t s t e p input Ess=")
+disp(kp,"Kp=")
+ end
+// ramp input
+ c=0
+ for i =1:2
+if(a(1,i)-b(1,i)*H ==0) then
+c=c+1
+ end
+end
+if(c ==2)
+printf (" f o r u n i t ramp input Es s=0 \nn")
+else if(c ==1) then
+Ess =(a(1 ,2) -b(1 ,2)*H)/a(1 ,1)*H;
+kv=1/Ess;
+disp (Ess ," for unitramp input Ess=")
+disp(kv,"Kv=")
+ else printf (" for unit ramp input Ess=infinite \nn")
+ end
+ end
+// p a r a b o l i c input
+c=0
+for i =1:3
+if(a(1,i)-b(1,i)*H ==0) then
+c=c+1
+end
+end
+if(c ==3)
+printf (" f o r u n i t p a r a b o l i c input Es s=0 \nn")
+ else if(c ==2) then
+Ess =(a(1 ,3) -b(1 ,3)*H)/a(1 ,1)*H;
+ka=1/Ess;
+ disp (Ess ," for unit parabolic input Ess=")
+ diisp(ka,"Ka=")
+ else printf (" for unit parabolic input Ess=infinite\nn")
+ end
+ end
diff --git a/1151/CH3/EX3.6/example6.sce b/1151/CH3/EX3.6/example6.sce
new file mode 100755
index 000000000..49a69ac1b
--- /dev/null
+++ b/1151/CH3/EX3.6/example6.sce
@@ -0,0 +1,5 @@
+//determine the type of following feedback system
+printf("a) G(s)=K/((s+10)*(s+5) \n solution : Since in denominator the power of s is zero i.e.m=0,henece its type zero system");
+printf("\n b) G(s)=20/(s*(0.05*s+1)*(0.1*s+1) \n solution : Since in denominator the power of s is one i.e.m=1,henece its type one system\n");
+printf("c) G(s)=20*(s+2)/(s^3*(s+2)*(s^2+5*s+5) \n solution : Since in denominator the power of s is three i.e.m=3,henece its type three system\n");
+printf("d) G(s)=K/(s*(s+1)*(s+2) \n solution : Since in denominator the power of s is one i.e.m=1,henece its type one system\n");
diff --git a/1151/CH3/EX3.8/example8.sce b/1151/CH3/EX3.8/example8.sce
new file mode 100755
index 000000000..8572fbdca
--- /dev/null
+++ b/1151/CH3/EX3.8/example8.sce
@@ -0,0 +1,52 @@
+s=%s
+p= poly([0 0 1],'s','coeff');
+q= poly([1 2 3 2],'s','coeff');
+M=p/q
+ disp (M,"M( s )=")
+ H =1;
+ R =1;
+ b= coeff (p)
+a= coeff (q)
+// s t e p input
+if (a(1 ,1) ==b(1 ,1)) then
+ printf (" f o r u n i t s t e p input Ess=0 \nn" )
+ else
+ Ess =1/ H*(1 -(b(1 ,1)*H/a(1 ,1)))*R;
+ kp=(1-Ess)/Ess;
+disp (Ess ," f o r u n i t s t e p input Ess=")
+disp(kp,"Kp=")
+ end
+// ramp input
+ c=0
+ for i =1:2
+if(a(1,i)-b(1,i)*H ==0) then
+c=c+1
+ end
+end
+if(c ==2)
+printf (" f o r u n i t ramp input Es s=0 \nn")
+else if(c ==1) then
+Ess =(a(1 ,2) -b(1 ,2)*H)/a(1 ,1)*H;
+kv=1/Ess;
+disp (Ess ," f o r u n i t ramp input Es s=")
+disp(kv,"Kv=")
+ else printf (" f o r u n i t ramp input Es s=i n f \nn")
+ end
+ end
+// p a r a b o l i c input
+c=0
+for i =1:3
+if(a(1,i)-b(1,i)*H ==0) then
+c=c+1
+end
+end
+if(c ==3)
+printf (" f o r u n i t p a r a b o l i c input Es s=0 \nn")
+ else if(c ==2) then
+Ess =(a(1 ,3) -b(1 ,3)*H)/a(1 ,1)*H;
+ka=1/Ess;
+ disp (Ess ," f o r u n i t p a r a b o l i c input Es s=")
+ diisp(ka,"Ka=")
+ else printf (" f o r u n i t p a r a b o l i c input Es s=i n f \nn")
+ end
+ end
diff --git a/1151/CH3/EX3.9/example9.sce b/1151/CH3/EX3.9/example9.sce
new file mode 100755
index 000000000..75f107ad9
--- /dev/null
+++ b/1151/CH3/EX3.9/example9.sce
@@ -0,0 +1,52 @@
+s=%s
+p= poly([0 0 20],'s','coeff');
+q= poly([1 2 22 4 4 0],'s','coeff');
+M=p/q
+ disp (M,"M( s )=")
+ H =1;
+ R =1;
+ b= coeff (p)
+a= coeff (q)
+// s t e p input
+if (a(1 ,1) ==b(1 ,1)) then
+ printf (" f o r u n i t s t e p input Ess=0 \nn" )
+ else
+ Ess =1/ H*(1 -(b(1 ,1)*H/a(1 ,1)))*R;
+ kp=(1-Ess)/Ess;
+disp (Ess ," f o r u n i t s t e p input Ess=")
+disp(kp,"Kp=")
+ end
+// ramp input
+ c=0
+ for i =1:2
+if(a(1,i)-b(1,i)*H ==0) then
+c=c+1
+ end
+end
+if(c ==2)
+printf (" f o r u n i t ramp input Es s=0 \nn")
+else if(c ==1) then
+Ess =(a(1 ,2) -b(1 ,2)*H)/a(1 ,1)*H;
+kv=1/Ess;
+disp (Ess ," for unitramp input Es s=")
+disp(kv,"Kv=")
+ else printf (" for unit ramp input Ess=infinite \nn")
+ end
+ end
+// p a r a b o l i c input
+c=0
+for i =1:3
+if(a(1,i)-b(1,i)*H ==0) then
+c=c+1
+end
+end
+if(c ==3)
+printf (" for unit parabolic input Ess=0 \nn")
+ else if(c ==2) then
+Ess =(a(1 ,3) -b(1 ,3)*H)/a(1 ,1)*H;
+ka=1/Ess;
+ disp (Ess ," for unit parabolic input Ess=")
+ diisp(ka,"Ka=")
+ else printf (" for unit parabolic input Ess=infinite \nn")
+ end
+ end