diff options
Diffstat (limited to '1151/CH1')
32 files changed, 150 insertions, 0 deletions
diff --git a/1151/CH1/EX1.1/example1.sce b/1151/CH1/EX1.1/example1.sce new file mode 100755 index 000000000..e30b57367 --- /dev/null +++ b/1151/CH1/EX1.1/example1.sce @@ -0,0 +1,5 @@ +printf("\n trnsfer function of the given network ")
+printf("\n Vi(s)=R*I(s)+s*L*I(s)")
+printf("\n V0=s*L*I")
+printf("\n V0(s)/Vi(s)=s*L*I/(R*I(s)+s*L*I(s))")
+printf("V0(s)/Vi(s)=s*L/(R+s*L) is the required transfer function")
diff --git a/1151/CH1/EX1.10/example10.sce b/1151/CH1/EX1.10/example10.sce new file mode 100755 index 000000000..7dc11966d --- /dev/null +++ b/1151/CH1/EX1.10/example10.sce @@ -0,0 +1,4 @@ +//find out transfer function of the system
+
+printf("syms K1 K2 M1 M2 B12 fc1 fc2");printf(" the differential equations in laplace form are \n F(s)=X1(s)[s^2*M1+s*B12+s*fc1+K1]-B12*s*X2(s)")
+printf("\n X1(s)=X2(s)*(s^2*M2+B12*s+s*fc2+K2)/(B12*s)\nf=B12*s/((s^2*M2+B12*s+s*fc2+K2)*(s^2*M1+B12*s+s*fc1+K1)-s^2*B12);\n Transfer function X2(s)/F(s)=f")
diff --git a/1151/CH1/EX1.18/example18.sce b/1151/CH1/EX1.18/example18.sce new file mode 100755 index 000000000..4eb583a53 --- /dev/null +++ b/1151/CH1/EX1.18/example18.sce @@ -0,0 +1,5 @@ +//write down systems of equations and find F(s)/X2(s)
+
+
+printf(" the differential equatios are \n F(s)=s^2*M1**X1(s)+(K1+K2)*X2(s)-K2*X2(s)")
+printf("\n X1(s)=X2(s)*(s^2*M2+K3+K2)/K2\n f=((s^2*M1+K1+K2)*(s^2*M2+K2+K3)-K2^2)/K2\n Transfer function F(s)/X2(s)=f")
diff --git a/1151/CH1/EX1.19/example18.sce b/1151/CH1/EX1.19/example18.sce new file mode 100755 index 000000000..4eb583a53 --- /dev/null +++ b/1151/CH1/EX1.19/example18.sce @@ -0,0 +1,5 @@ +//write down systems of equations and find F(s)/X2(s)
+
+
+printf(" the differential equatios are \n F(s)=s^2*M1**X1(s)+(K1+K2)*X2(s)-K2*X2(s)")
+printf("\n X1(s)=X2(s)*(s^2*M2+K3+K2)/K2\n f=((s^2*M1+K1+K2)*(s^2*M2+K2+K3)-K2^2)/K2\n Transfer function F(s)/X2(s)=f")
diff --git a/1151/CH1/EX1.2/example2.sce b/1151/CH1/EX1.2/example2.sce new file mode 100755 index 000000000..a6b7aa52a --- /dev/null +++ b/1151/CH1/EX1.2/example2.sce @@ -0,0 +1,4 @@ +printf("\n trnsfer function of the given network ")
+printf("\n E0(s)=(1/(C*s))*I(s)")
+printf("\n E0(s/Ei(s)=(I(s)/C*s)*((C*s)/(I(s)*(s^2*LC+s*R*C+1)))")
+printf(" E0(s/Ei(s)= 1/(s^2*LC+s*R*C+1)is the required transfer function")
diff --git a/1151/CH1/EX1.20/example20.sce b/1151/CH1/EX1.20/example20.sce new file mode 100755 index 000000000..2eeb445a1 --- /dev/null +++ b/1151/CH1/EX1.20/example20.sce @@ -0,0 +1,29 @@ +printf(" given \n N1=300\n N2=50 \n w1=3 rad/sec\n Torque for wheel 1=10Nm \n displacement for wheel 1=3 radian(clockwise)\n angular accerleration of wheel 1 is 4 rad/sec^2 ")
+printf("find a) r1:r2 and r1:r3\n b)w3,w4 \n c)torque for wheel 2,wheel 3 \n displacement for wheel 2,3and 4\n e) angular displacement of wheel3 ")
+n1=300;
+n2=50;
+n3=150;
+n4=30;
+r=n1/n2;
+r1=n1/n3;
+disp(r,"r1:r2=");
+disp(r1,"r1:r3=")
+w1=15;
+w3=r1*w1; w4=n1*w1/n2;
+disp(w3,"w3=")
+disp(w4,"w4=")
+t1=10;
+t2=t1/r;
+t3=t1/r1;
+disp(t2,"torque for wheel 2 (in Nm) is")
+disp(t3,"torque for wheel 3(in Nm) is")
+disp(d2,"displacement for wheel 2 (in radian) is")
+
+d1=3;
+d2=r*d1;
+d3=d2*n2/n3;
+d4=d3*n3/n1;
+disp(d3,"displacement for wheel 3 (in radian) is")
+disp(d4,"displacement for wheel 4 (in radian) is")
+a=r1*4;
+disp(a,"angular displacement for wheel 3(in rad/sec^2) is")
diff --git a/1151/CH1/EX1.23/example23.sce b/1151/CH1/EX1.23/example23.sce new file mode 100755 index 000000000..1f9497c38 --- /dev/null +++ b/1151/CH1/EX1.23/example23.sce @@ -0,0 +1,2 @@ +//determine the transfer function
+printf("z1=(1+R*C*s)/(C*s)\n z2=1/(C*s) \n g=z2/z2 \nTransfer function Eo(s)/Ei(s)=g")
diff --git a/1151/CH1/EX1.24/example24.sce b/1151/CH1/EX1.24/example24.sce new file mode 100755 index 000000000..99df02246 --- /dev/null +++ b/1151/CH1/EX1.24/example24.sce @@ -0,0 +1,5 @@ +//find out transfer function of the system
+
+printf("z1=(1+R*C2*s)*C1/(C1+C2+C1*C2*R*s)\n f=K1/((s^2*M2+K1+K2)*(K1+s^2*M1)-K1^2 \n Transfer function of electrical network Vo(s)/Vi(s)=f");
+printf("g=(s*B+K2)/(K1+K2+s*B\n");
+printf("Transfer function ofmechanical network analogous to electrical network X2(s)/X1(s)=g");
diff --git a/1151/CH1/EX1.26/example26.sce b/1151/CH1/EX1.26/example26.sce new file mode 100755 index 000000000..e633cfc7b --- /dev/null +++ b/1151/CH1/EX1.26/example26.sce @@ -0,0 +1,3 @@ +//determine the transfer function
+printf("syms R1 R2 C1 C2\nz1=R2/(R2*C2*s+1);\n z2=(R1*C1*s+1)/(C1*s);\n g=z2/(z1+z2);Transfer function Eo(s)/Ei(s)=g")
+
diff --git a/1151/CH1/EX1.3/example3.sce b/1151/CH1/EX1.3/example3.sce new file mode 100755 index 000000000..8976872ec --- /dev/null +++ b/1151/CH1/EX1.3/example3.sce @@ -0,0 +1,3 @@ +//find transfer function of V2(s)/V1(s
+printf("syms R1 R2 C")
+printf("G=(R2+R1*R2*C*s)/(R1+R2+R1*R2*C*s)\n transfer function V2(s)/V1(s)=G")
diff --git a/1151/CH1/EX1.30/example30.sce b/1151/CH1/EX1.30/example30.sce new file mode 100755 index 000000000..9fa30d593 --- /dev/null +++ b/1151/CH1/EX1.30/example30.sce @@ -0,0 +1,4 @@ +//derive transfer function using block diagram reduction
+
+printf("syms G1 H1 G2 H2 H3 \n l1=G1/(1+G1*H1);//feedback configuration\n l2=G2/(1+G2*H2);//feedback configuration\n l3=l1*l2;//cascading combination\n g=l3/(1+l3*H3");
+printf("\nC(s)/R(s)=g")
diff --git a/1151/CH1/EX1.31/example31.sce b/1151/CH1/EX1.31/example31.sce new file mode 100755 index 000000000..ad0a4d933 --- /dev/null +++ b/1151/CH1/EX1.31/example31.sce @@ -0,0 +1,2 @@ +//find the transfer function using block diagram reduction
+printf("syms R1 R2 C1 C2 \n l1=1/(1+R2*C2*s);//unity feedback\nl2=1/(s*C1); \n l3=1/R1;\n l4=s*C2;\n l5=1+R2*C2*s;\n g1=l1*l2;//cascading of blocks\ng2=g1/(1+g1*l4);//feedback configuration\n g3=l3*g2;//cascading of blocks\n g4=g3/(1+g3*l5);\nTransfer function Vo(s)/Vi(s)=g4")
diff --git a/1151/CH1/EX1.32/example32.sce b/1151/CH1/EX1.32/example32.sce new file mode 100755 index 000000000..93c711f27 --- /dev/null +++ b/1151/CH1/EX1.32/example32.sce @@ -0,0 +1,2 @@ +//derive transfer function using block diagram reduction
+printf("syms G1 H1 G2 H2 \nG3l1=G2*G3;//cascading of blocks\nl2=l1/(1+l1*H2);//feedback configuration\nl3=l2*G1;//cascading of blocksl4=H1/G3;l5=l3/(1+l3*l2)\n C(s)/R(s)=l5")
diff --git a/1151/CH1/EX1.34/example34.sce b/1151/CH1/EX1.34/example34.sce new file mode 100755 index 000000000..f96cbfc61 --- /dev/null +++ b/1151/CH1/EX1.34/example34.sce @@ -0,0 +1,2 @@ +//derive transfer function using block diagram reduction
+printf("syms G1 H1 G2 H2 G3 G4\nl1=(G1-G1*G2*G4)/(1-G1*G2*G3*G4);//X1/R\nl2=G2*(1-G1*G3)/(1-G1*G2*G3*G4);//X2/R\n g=l1+l2;//X1/R+X2/R\n C(s)/R(s)=g");
diff --git a/1151/CH1/EX1.35/example35.sce b/1151/CH1/EX1.35/example35.sce new file mode 100755 index 000000000..447cdc54c --- /dev/null +++ b/1151/CH1/EX1.35/example35.sce @@ -0,0 +1,2 @@ +//determine transfer function using bloack diagram reduction techniques
+printf("syms H1 H2 G\n l1=H1/G;\n l2=1+l1;//parallel blocks\nl3=G/(1+G*H2);\nl4=l3*l2;\nC(s)/R(s)=l3");
diff --git a/1151/CH1/EX1.37/example37.sce b/1151/CH1/EX1.37/example37.sce new file mode 100755 index 000000000..91810aa14 --- /dev/null +++ b/1151/CH1/EX1.37/example37.sce @@ -0,0 +1,2 @@ +//derive transfer function using block diagram reduction
+printf("syms G1 H1 G2 H2 G3\n l1=G3/(1+G3*H1);//feedback blocks\n l2=G2*l1;//cascading of blocks\n l3=H2/G3;\n l4=l2/(1+l2*l3);//-ve feedback\n l5=l4*G1;//cascading blocks\n l6=l5/(1+l5);C(s)/R(s)=l6")
diff --git a/1151/CH1/EX1.38/example38.sce b/1151/CH1/EX1.38/example38.sce new file mode 100755 index 000000000..bc34dc297 --- /dev/null +++ b/1151/CH1/EX1.38/example38.sce @@ -0,0 +1,2 @@ +//derive transfer function using block diagram reduction
+printf("\nl1=3/(s+1);\n l2=-6/(s+2);\n l3=3/(s+3);\n g=l1+l2+l3;\nC(s)/R(s)=g")
diff --git a/1151/CH1/EX1.39/example39.sce b/1151/CH1/EX1.39/example39.sce new file mode 100755 index 000000000..0a83e270b --- /dev/null +++ b/1151/CH1/EX1.39/example39.sce @@ -0,0 +1,2 @@ +//derive transfer function using block diagram reduction
+printf("syms G1 H1 G2 H2\nL1=G1*G2/(1+G2*H2+G1*G2*H1) //WHEN R3=0,R2=0,R4=0\nL2=G2/(1+G2*H2+G1*G2*H1) //WHEN R3=0,R1=0,R4=0\nL3=-G2/(1+G2*H2+G1*G2*H1) //WHEN R2=0,R1=0,R4=0\nL4=-G2*G1*H1/(1+G2*H2+G1*G2*H1) //WHEN R2=0,R1=0,R3=0\nl=L1+L2+L3+L4;\nC(s)/R(s)=l")
diff --git a/1151/CH1/EX1.4/example4.sce b/1151/CH1/EX1.4/example4.sce new file mode 100755 index 000000000..9d0c9a78d --- /dev/null +++ b/1151/CH1/EX1.4/example4.sce @@ -0,0 +1,2 @@ +//find transfer function of lag network
+printf("syms R1 R2 C\n ei=(R1+R2+(1/C)*s);\n eo=(R2+(1/C)*s);\n G=eo/ei;\n transfer function Eo(s)/Ei(s)is:G")
diff --git a/1151/CH1/EX1.41/example41.sce b/1151/CH1/EX1.41/example41.sce new file mode 100755 index 000000000..0a34370b1 --- /dev/null +++ b/1151/CH1/EX1.41/example41.sce @@ -0,0 +1,5 @@ +//to find transfer function using mason gain formula
+
+printf("syms a12 a23 a32 a34 a42 a44 a35 a52 a45\n");
+printf("there are two forward paths:\n 1) x1 tox2 to x3 to x4 to x5\n 2) x1 to x2 to x3 to x5")
+printf("\nP1=a12*a23*a34*a45;//forward path2 gain\nP2=a12*a23*a35;//forward path1 gain\n//gain of individual loops\nL1=a23*a32;\nL2=a23*a34*a42;\nL3=a23*a34*a45*a52;\nL4=a23*a35*a52;\nL5=a44;\n//gain of two non touching loops\ng1=L1*L5;\n g2=L4*L5;\n//since all the loops touches the forward path1 and loop 5 do not touch forward path 2\nd1=1;\nd2=1-L5;\n d=1-(L1+L2+L3+L4+L5)+(g1+g2);\nG=(P1*d1+P2*d2)/d;\ntransfer function x5/x1=G")
diff --git a/1151/CH1/EX1.42/example42.sce b/1151/CH1/EX1.42/example42.sce new file mode 100755 index 000000000..c8dde05aa --- /dev/null +++ b/1151/CH1/EX1.42/example42.sce @@ -0,0 +1,3 @@ +//to find transfer function using mason gain formula
+
+printf("syms G1 G2 G3 G4 H1 H2 H3 H4\n//gains of forward paths\nP1=G1*G2*G3*G4;//forward path1 gain\nP2=G1*G5;//forward path2 gain\n//gain of individual loops\nL1=-G1*H1;\nL2=-G3*H2;\nL3=-G1*G2*G3*H3;\nL4=-H4;\n//gain of two non touching loops\ng1=G1*G3*H1*H2;\ng2=G1*H1*H4;\ng3=G3*H2*H4;\ng4=G1*G2*G3*H2*H4;\n//gain of three non touching loops\nl1=-G1*G3*H1*H2*H4;\n//since all the loops touches the forward path1 and loop 2 do not touch forward path 2\nd1=1;\nd2=1-L2;\nd=1-(L1+L2+L3+L4)+(g1+g2+g3+g4)-l1;\nG=(P1*d1+P2*d2)/d;\ntransfer function C/R=G")
diff --git a/1151/CH1/EX1.43/example43.sce b/1151/CH1/EX1.43/example43.sce new file mode 100755 index 000000000..4250148da --- /dev/null +++ b/1151/CH1/EX1.43/example43.sce @@ -0,0 +1,3 @@ +//to find transfer function using mason gain formula
+
+printf("syms G1 G2 G3 G4 G5 G6 H1 H2 \n //gains of forward paths\n P1=G1*G2*G3*G6;//forward path1 gain\n P2=G1*G2*G4*G6;//forward path2 gain\nP3=G1*G2*G5*G6;//forward path3 gain\n //gain of individual loops\n L1=-G1*H1;\n L2=-G1*G2*G3*H2;\n L3=-G1*G2*G4*H2;\n L4=-G1*G2*G5*H2;\n //NO TWO NON TOUCHING LOOPS ARE THERE\n d1=1;\n d2=1;\n d3=1;\n d=1-(L1+L2+L3+L4);\n G=(P1*d1+P2*d2+P3*d3)/d;\ntransfer function C/R=G")
diff --git a/1151/CH1/EX1.44/example44.sce b/1151/CH1/EX1.44/example44.sce new file mode 100755 index 000000000..966fe9aeb --- /dev/null +++ b/1151/CH1/EX1.44/example44.sce @@ -0,0 +1 @@ +printf("//to find transfer function using mason gain formula\nsyms G1 G2 G3 G4 \n//gains of forward paths\nP1=G1*G2*G3;//forward path1 gain\n P2=G4;//forward path2 gain\n//gain of individual loops\nL1=-G1*G2;\nL2=G2*G3;\n//NO TWO NON TOUCHING LOOPS ARE THERE\n//since all the loops touches the forward path1 and loop1 and 2 do not touch forward path 2\\nd1=1;\nd2=1-(L1+L2);\nG=(P1*d1+P2*d2)/d2;\ntransfer function C/R=G")
diff --git a/1151/CH1/EX1.45/example45.sce b/1151/CH1/EX1.45/example45.sce new file mode 100755 index 000000000..ae68be815 --- /dev/null +++ b/1151/CH1/EX1.45/example45.sce @@ -0,0 +1 @@ +printf("//to find transfer function using mason gain formula\n syms G1 G2 G3 G4 H1 H2 \n //gains of forward paths\n P1=G1*G2*G3;//forward path1 gain\n P2=-G1*G4;//forward path2 gain\n //gain of individual loops\nL1=-G1*G2*H1;\nL2=-G2*G3*H2;\nL3=-G1*G2*G3;\nL4=G1*G4;\nL5=G4*H2;\n//NO TWO NON TOUCHING LOOPS ARE THERE\nd1=1;\nd2=1;\nd=1-(L1+L2+L3+L4+L5);\nG=(P1*d1+P2*d2)/d;\ntransfer function C/R=G")
diff --git a/1151/CH1/EX1.46/example46.sce b/1151/CH1/EX1.46/example46.sce new file mode 100755 index 000000000..db0d6b4c7 --- /dev/null +++ b/1151/CH1/EX1.46/example46.sce @@ -0,0 +1,3 @@ +//to find transfer function using mason gain formula
+
+printf("syms G1 G2 G3 G4 H1 H2 \n //gains of forward paths\n P1=G1*G2*G3;//forward path1 gain\n P2=G4;//forward path2 gain\n //gain of individual loops\n L1=-G2*H1;\n L2=-G1*G2*H1;\n L3=-G2*G3*H2;\n //there are no two non touching loops\n //since all the loops touches the forward path1 and loop 1,2,3 do not touch forward path 2\n d1=1;\n d2=1-(L1+L3+L2);\n G=(P1*d1+P2*d2)/d2;\n transfer function C/R=G")
diff --git a/1151/CH1/EX1.48/example48.sce b/1151/CH1/EX1.48/example48.sce new file mode 100755 index 000000000..11950933f --- /dev/null +++ b/1151/CH1/EX1.48/example48.sce @@ -0,0 +1,27 @@ +//to find transfer function using mason gain formula
+
+syms G1 G2 G3 G4 G5 G6 G7 G8 H1 H2
+//gains of forward paths
+P1=G2*G6*G4;//forward path1 gain
+P2=G3*G5*G7;//forward path2 gain
+P3=G1*G2*G7;//forward path3 gain
+P4=G3*G8*G6;//forward path4 gain
+P5=-G1*G2*G8*G6*H2;//forward path5 gain
+P6=-G3*G8*G1*G7*H1;//forward path6 gain
+//gain of individual loops
+L1=-G4*H1;
+L2=-G5*H2;
+L3=H1*H2*G1*G8;
+//gain of two non touching loops
+g1=G5*G4*H1*H2;
+
+//since all the loops touches the forward path3,4,5,6 and loop 2 do not touch forward path 1 and loop1 donot touch forward path2
+d1=1-L2;
+d2=1-L1;
+d3=1;
+d4=1;
+d5=1;
+d6=1;
+d=1-(L1+L2+L3)+g1;
+G=(P1*d1+P2*d2+P4*d4+P3*d3+P5*d5+P6*d6)/d;
+disp(G,"transfer function C/R=")
diff --git a/1151/CH1/EX1.49/example49.sce b/1151/CH1/EX1.49/example49.sce new file mode 100755 index 000000000..8f31d3149 --- /dev/null +++ b/1151/CH1/EX1.49/example49.sce @@ -0,0 +1,3 @@ +//to find transfer function using mason gain formula
+
+printf("syms R1 R2 C1 C2 \n //gains of forward path\n P1=1/(R1*R2*C1*C2*s^2);//forward path1 gain\n //gain of individual loops\n L1=-1/(R1*C1*s);\n L2=-1/(R2*C1*s);\n L3=-1/(R2*C2*s);\n //gain of two non touching loops\n g1=1/(s^2*R1*R2*C1*C2);\n //since all the loops touches the forward path1 so\n d1=1\n d=1-(L1+L2+L3)+g1;\n G=(P1*d1)/d;\n transfer function C/R=G")
diff --git a/1151/CH1/EX1.5/example5.sce b/1151/CH1/EX1.5/example5.sce new file mode 100755 index 000000000..a1518e17b --- /dev/null +++ b/1151/CH1/EX1.5/example5.sce @@ -0,0 +1,2 @@ +//determine the transfer function
+printf("syms R1 R2 C1 C2 \n z1=R1/(R1*C1*s+1);\n z2=(R2*C2*s+1)/(C2*s)\n g=z2/(z1+z2)\n Transfer function Eo(s)/Ei(s)= g")
diff --git a/1151/CH1/EX1.51/example51.sce b/1151/CH1/EX1.51/example51.sce new file mode 100755 index 000000000..ec240bdfa --- /dev/null +++ b/1151/CH1/EX1.51/example51.sce @@ -0,0 +1,6 @@ +printf("\n A thermometer has time constant =15.33sec\n its quickly taken from 0 (degree C) to water bath having a temperatur 100(degree C)\n find the temperature indicated after 60 sec ");
+printf("temperature indicated by thermometer can be given as =100*(1-e^-(t/T))");
+T=15.33;
+t=60;
+h=100*(1-%e^-(t/T));
+disp(h,"temperature indicated by thermometer(in degree C)) after 60 seconds=")
diff --git a/1151/CH1/EX1.6/example6.sce b/1151/CH1/EX1.6/example6.sce new file mode 100755 index 000000000..f13e07979 --- /dev/null +++ b/1151/CH1/EX1.6/example6.sce @@ -0,0 +1,2 @@ +//determine the transfer function of the transformer coupled network
+printf("syms R1 R2 C2 C1 L1 L2 M L3 \n G=s^3*R2*C1*C2*M/(((s^2*C2*(L3+L2))+1+s*R2*C2)*(S^2*L1*C1+s*C1*R1+1)-M^2*s^4*C1*C2)\n Transfer function=G")
diff --git a/1151/CH1/EX1.8/example8.sce b/1151/CH1/EX1.8/example8.sce new file mode 100755 index 000000000..519df72f4 --- /dev/null +++ b/1151/CH1/EX1.8/example8.sce @@ -0,0 +1,5 @@ +//find out transfer function of the system
+
+printf("syms K1 K2 M1 M2")
+printf(" the differential equatios in laplace form are \n F(s)=M1*S^2*X1(s)+K1*X1(s)-K2*X2(s)")
+printf("\n X1(s)=X2(s)*(s^2*M2+K1+K2)/K1 \n f=K1/((s^2*M2+K1+K2)*(K1+s^2*M1)-K1^2)\nTransfer function X2(s)/F(s)=f")
diff --git a/1151/CH1/EX1.9/example9.sce b/1151/CH1/EX1.9/example9.sce new file mode 100755 index 000000000..b42fac7fe --- /dev/null +++ b/1151/CH1/EX1.9/example9.sce @@ -0,0 +1,4 @@ +//write down equations for block diagram
+printf("F(s)=M*s^2*X(s)+B*S*X(s)+K*X(s) \n given\n M=10 kg \n B=30N/m/sec \n K=20N/m")
+printf("putting the values in the above equation")
+printf("F(s)=10*s^2*X(s)+30*S*X(s)+20*X(s")
|