diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /1151/CH3 | |
download | Scilab-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-x | 1151/CH3/EX3.1/test.sce | 46 | ||||
-rwxr-xr-x | 1151/CH3/EX3.2/example2.sce | 52 | ||||
-rwxr-xr-x | 1151/CH3/EX3.3/example3.sce | 9 | ||||
-rwxr-xr-x | 1151/CH3/EX3.4/example4.sce | 51 | ||||
-rwxr-xr-x | 1151/CH3/EX3.5/example5.sce | 52 | ||||
-rwxr-xr-x | 1151/CH3/EX3.6/example6.sce | 5 | ||||
-rwxr-xr-x | 1151/CH3/EX3.8/example8.sce | 52 | ||||
-rwxr-xr-x | 1151/CH3/EX3.9/example9.sce | 52 |
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
|