From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 1904/CH11/EX11.1/11_1.sce | 10 ++++++++++ 1904/CH11/EX11.2/11_2.sce | 33 +++++++++++++++++++++++++++++++ 1904/CH11/EX11.3/11_3.sce | 47 ++++++++++++++++++++++++++++++++++++++++++++ 1904/CH11/EX11.4/11_4.sce | 42 +++++++++++++++++++++++++++++++++++++++ 1904/CH11/EX11.5/11_5.sce | 49 ++++++++++++++++++++++++++++++++++++++++++++++ 1904/CH11/EX11.6/11_6.sce | 38 +++++++++++++++++++++++++++++++++++ 1904/CH11/EX11.7/11_7.jpg | Bin 0 -> 20974 bytes 1904/CH11/EX11.7/11_7.sce | 21 ++++++++++++++++++++ 1904/CH11/EX11.8/11_8.jpg | Bin 0 -> 43698 bytes 1904/CH11/EX11.8/11_8.sce | 32 ++++++++++++++++++++++++++++++ 1904/CH11/EX11.9/11_9.sce | 22 +++++++++++++++++++++ 11 files changed, 294 insertions(+) create mode 100755 1904/CH11/EX11.1/11_1.sce create mode 100755 1904/CH11/EX11.2/11_2.sce create mode 100755 1904/CH11/EX11.3/11_3.sce create mode 100755 1904/CH11/EX11.4/11_4.sce create mode 100755 1904/CH11/EX11.5/11_5.sce create mode 100755 1904/CH11/EX11.6/11_6.sce create mode 100755 1904/CH11/EX11.7/11_7.jpg create mode 100755 1904/CH11/EX11.7/11_7.sce create mode 100755 1904/CH11/EX11.8/11_8.jpg create mode 100755 1904/CH11/EX11.8/11_8.sce create mode 100755 1904/CH11/EX11.9/11_9.sce (limited to '1904/CH11') diff --git a/1904/CH11/EX11.1/11_1.sce b/1904/CH11/EX11.1/11_1.sce new file mode 100755 index 000000000..a702b9037 --- /dev/null +++ b/1904/CH11/EX11.1/11_1.sce @@ -0,0 +1,10 @@ +//To Determine the Approximate value of the component reliability +//Page 598 +clc; +clear; +Rsys=0.99 //Minimum Acceptable System Reliabilty +n=15; //Number of identical Components +q=(1-Rsys)/n; //Probability of component failure +Ri=1-q; //Approximate value of the component reliability + +printf('The Approximate Value of The component reliability is %g\n',Ri) diff --git a/1904/CH11/EX11.2/11_2.sce b/1904/CH11/EX11.2/11_2.sce new file mode 100755 index 000000000..8e47035de --- /dev/null +++ b/1904/CH11/EX11.2/11_2.sce @@ -0,0 +1,33 @@ +//To Determine the fault components of the system +//Page 606 +clc; +clear; +L=4; //Total Length of the cable +Lov=3;//Length of Overhead Cable +Lu=L-Lov; //Length of Underground Cable +Nct=2; //Number of circuit terminations +T=10; //No of years for which the record is shown + +Fov=2; // Faults Per Mile of the Over Head Cable +Fu=1; //Faults Per Mile of The Underground cable + +Ct=0.3/100// Cable Termination Fault Rate + +//Repair Time +Tov=3; //Over Head +Tu=28; //Underground +Tct=3; //Cable Termination + +lamdaFDR= (Lov*Fov/T)+(Lu*Fu/T)+(2*Ct); //Total Annual Fault Rate + +rFDR=((Tov*Lov*Fov/T)+(Tu*Lu*Fu/T)+(2*Ct*Tct))/lamdaFDR; //Annual Fault Restoration Time + +mFDR=8760-rFDR; //Annual Mean Time of Failure + +UFDR=rFDR*100/(rFDR+mFDR); //Unavailability of Feeder +AFDR=100-UFDR; //Availability of Feeder + +printf('a) The Total Annual Fault Rate is %g faults per year\n',lamdaFDR) +printf('b) The Annual Fault Restoration Time is %g hours per fault per year\n',rFDR) +printf('c) Unavailability of the feeder is %g percent\n',UFDR) +printf('d) Availability of the feeder is %g percent\n',AFDR) diff --git a/1904/CH11/EX11.3/11_3.sce b/1904/CH11/EX11.3/11_3.sce new file mode 100755 index 000000000..2902cd7a2 --- /dev/null +++ b/1904/CH11/EX11.3/11_3.sce @@ -0,0 +1,47 @@ +//To Determine the Annual Fault properties for A B C Customers +//Page 608 +clc; +clear; + +//Annual average Fault rates +Fm=0.08; +Fl=0.2; + + +//Average Repair Times +Rm=3.5; //Main +Rl=1.5; //Lateral +Rs=0.75;//Manual Sections + +// Distances of the Lateral Feeders of A,B, and C respectively +Lla=2; +Llb=1.5; +Llc=1.5; + +// Distances of the Main Feeders of A,B, and C respectively +Lma=1; +Lmb=1; +Lmc=1; + +TFm=(Lma*Fm)+(Lmc*Fm)+(Lmb*Fm); //Annual Fault of the Main Sections + +deff('x=SusInt(y)','x=TFm+(Fl*y)') //Function to find the Total Annual Sustained Interruption rates + +//Sustained Interruption Rates for A,B and C +IrA=SusInt(Lla); +IrB=SusInt(Llb); +IrC=SusInt(Llc); + +//Annual Repair time for A,B and C +rA=((Lma*Fm*Rm)+(Lmb*Fm*Rs)+(Lmc*Fm*Rs)+(Lla*Fl*Rl))/IrA; +rB=((Lma*Fm*Rm)+(Lmb*Fm*Rm)+(Lmc*Fm*Rs)+(Llb*Fl*Rl))/IrB; +rC=((Lma*Fm*Rm)+(Lmb*Fm*Rm)+(Lmc*Fm*Rm)+(Llc*Fl*Rl))/IrC; + +printf('\ni) The Annual Sustained Interruption Rates for:\n') +printf('Customer A : %g faults per year\n',IrA) +printf('Customer B : %g faults per year\n',IrB) +printf('Customer C : %g faults per year\n',IrC) +printf('\nii) The Average Annual Repair Time (Restoration Time) for:\n') +printf('Customer A : %g hours per fault per year\n',rA) +printf('Customer A : %g hours per fault per year\n',rB) +printf('Customer A : %g hours per fault per year\n',rC) diff --git a/1904/CH11/EX11.4/11_4.sce b/1904/CH11/EX11.4/11_4.sce new file mode 100755 index 000000000..4a33576fa --- /dev/null +++ b/1904/CH11/EX11.4/11_4.sce @@ -0,0 +1,42 @@ +//To Determine the Equivalent System Reliability of Each configuration +//Page 612 +clc; +clear; + +Ri=0.85; + +deff('x=relp(y,z)','x=1-((1-(Ri^y))^z)')//Equal Parallel Combination + +deff('x=rels(y,z)','x=(1-((1-Ri)^y))^z') //Equal Series Combination + +//Case 1: 4 elements in series + +Req1= rels(1,4); + +//Case 2: Two Comination of 4 elements in series, parallel to each other + +Req2=relp(4,2); + +//Case 3 : ((two elements in series)//(two elements in series))in series with ((two elements in series)//(two elements in series)) + +//Two Segments +R1=relp(2,2); +R2=relp(2,2); +Req3=R1*R2; + +//Case 4 : (two elements in parallel)in series with ((three elements in series)//(three elements in series)) + +//Two Segments +R1=relp(1,2); +R2=relp(3,2); +Req4=R1*R2; + +//Case 5, 4 groups of (2 elements in parallel) connected in series to each other +Req5=rels(2,4); + +printf('The Equivalent System reliability for:\n') +printf('a) Configuration A : %g\n',Req1) +printf('b) Configuration B : %g\n',Req2) +printf('c) Configuration C : %g\n',Req3) +printf('d) Configuration D : %g\n',Req4) +printf('e) Configuration E : %g\n',Req5) diff --git a/1904/CH11/EX11.5/11_5.sce b/1904/CH11/EX11.5/11_5.sce new file mode 100755 index 000000000..dd26b14ab --- /dev/null +++ b/1904/CH11/EX11.5/11_5.sce @@ -0,0 +1,49 @@ +//To Design the system to meet the given Equivalent System Reliability +//Page 614 +clc; +clear; + +//Individual System Reliabilities +Ra=0.8; +Rb=0.95; +Rc=0.99; +Rd=0.90; +Re=0.65; + +//When All Are Connected in Series + +Req=Ra*Rb*Rc*Rd*Re; //Equivalent System Reliability + +Rr=0.8; //Required + +Rae=Rr/(Rb*Rc*Rd); + +//Since Connecting the elements in parallel will increase their reliability +deff('x=rel(Ri,y,)','x=(1-((1-Ri)^y))') //Equal Only Parallel Combination + +//Since Connecting the elements in parallel will increase their reliability +//Conditions to Find The Number of Elements to be used +for i= 1:10 + L=i; //Number of Time Element A is used + R1=rel(Ra,i); + X=R1-Rae; + if(abs(X)+X==0) + continue; + else + break; + end +end + +for i= 1:10 + M=i;//Number of Time Element E is used + R2=rel(Re,i); + X=R2-Rae; + if(abs(X)+X==0) + continue; + else + break; + end +end + +printf('a) The Equivalent system Reliability is %g\n',Req) +printf('b) One Each of B,C and D all connected in series are connected in series\nwith the series combination of X(Comination of %g elements of A, All Connected in Parallel)\nand Y(Comination of %g elements of E, All Connected in Parallel) to achieve \n%g Equivalent System Realibility\n',L,M,Rr) diff --git a/1904/CH11/EX11.6/11_6.sce b/1904/CH11/EX11.6/11_6.sce new file mode 100755 index 000000000..e1c726c17 --- /dev/null +++ b/1904/CH11/EX11.6/11_6.sce @@ -0,0 +1,38 @@ +//To Find The Probability on the reliability of transformers +//Page 614 +clc; +clear; + +//Reliabilities of The Three Transformers +Pa=0.9; +Pb=0.95; +Pc=0.99; + +//Faliures of Three Transformers +Qa=1-Pa; +Qb=1-Pb; +Qc=1-Pc; + +//Probability of NO Transformer Failing +Pnf=Pa*Pb*Pc; + +PfA=Qa*Pb*Pc//Probability of Transformer A Failing +PfB=Pa*Qb*Pc//Probability of Transformer B Failing +PfC=Pa*Pb*Qc//Probability of Transformer C Failing + +PfAB=Qa*Qb*Pc//Probability of Transformer A and B Failing +PfBC=Pa*Qb*Qc//Probability of Transformer B and C Failing +PfCA=Qa*Pb*Qc//Probability of Transformer C and A Failing + +Pf=Qa*Qb*Qc; //Probability of All Transformers failing + +printf('\na) Probability of No Transformer Failing is %g\n',Pnf) +printf('\nb)\n') +printf('Probability of Transformer A Failing is %g\n',PfA) +printf('Probability of Transformer B Failing is %g\n',PfB) +printf('Probability of Transformer C Failing is %g\n',PfC) +printf('\nc)\n') +printf('Probability of Transformers A and B Failing is %g\n',PfAB) +printf('Probability of Transformers B and C Failing is %g\n',PfBC) +printf('Probability of Transformers C and A Failing is %g\n',PfCA) +printf('\nd) Probability of All Three Transformers Failing is %g\n',Pf) diff --git a/1904/CH11/EX11.7/11_7.jpg b/1904/CH11/EX11.7/11_7.jpg new file mode 100755 index 000000000..0d7127da8 Binary files /dev/null and b/1904/CH11/EX11.7/11_7.jpg differ diff --git a/1904/CH11/EX11.7/11_7.sce b/1904/CH11/EX11.7/11_7.sce new file mode 100755 index 000000000..d1bb40baf --- /dev/null +++ b/1904/CH11/EX11.7/11_7.sce @@ -0,0 +1,21 @@ +//To Determine Probabilities Using Markovian Principle +//Page 619 +clc; +clear; + +//Conditional Probabilites Present Future +Pdd=2/100; //Down Down +Pud=5/100; //Up Down +Pdu=1-Pdd; //Down up +Puu=1-Pud; //Up Up + +P=[Pdd,Pdu;Pud,Puu]; //Transition Matrix + +printf('\na) The Conditional Probabilites for\n') +printf('Transformers Down in Present and Down in Future is %g\n',Pdd) +printf('Transformers Down in Present and Up in Future is %g\n',Pdd) +printf('Transformers Up in Present and Down in Future is %g\n',Pdd) +printf('Transformers Up in Present and Up in Future is %g\n',Pdd) +printf('\nb) The Transition Matrix is\n') +disp(P) +printf('\nc) The Transition Diagram can be viewed with the result file attached to this code\n') diff --git a/1904/CH11/EX11.8/11_8.jpg b/1904/CH11/EX11.8/11_8.jpg new file mode 100755 index 000000000..f2dec90fa Binary files /dev/null and b/1904/CH11/EX11.8/11_8.jpg differ diff --git a/1904/CH11/EX11.8/11_8.sce b/1904/CH11/EX11.8/11_8.sce new file mode 100755 index 000000000..1e5219a79 --- /dev/null +++ b/1904/CH11/EX11.8/11_8.sce @@ -0,0 +1,32 @@ +//To Determine the Conditional Outage Probabilites +//Page 620 +clc; +clear; + +//Conditional Outage Probabilites From The Table Given +P11=40/100; +P12=30/100; +P13=30/100; +P21=20/100; +P22=50/100; +P23=30/100; +P31=25/100; +P32=25/100; +P33=50/100; + +//Transition Matrix +P=[P11,P12,P13;P21,P22,P23;P31,P32,P33]; + +printf("\na) The Conditional Outage Probabilites for:\n") +printf("Presently Outaged Feeder is 1, Next Outaged Feeder is 1 is %g\n",P11) +printf("Presently Outaged Feeder is 1, Next Outaged Feeder is 2 is %g\n",P12) +printf("Presently Outaged Feeder is 1, Next Outaged Feeder is 3 is %g\n",P13) +printf("Presently Outaged Feeder is 2, Next Outaged Feeder is 1 is %g\n",P21) +printf("Presently Outaged Feeder is 2, Next Outaged Feeder is 2 is %g\n",P22) +printf("Presently Outaged Feeder is 2, Next Outaged Feeder is 3 is %g\n",P23) +printf("Presently Outaged Feeder is 3, Next Outaged Feeder is 1 is %g\n",P31) +printf("Presently Outaged Feeder is 3, Next Outaged Feeder is 2 is %g\n",P32) +printf("Presently Outaged Feeder is 3, Next Outaged Feeder is 3 is %g\n",P33) +printf("\nb) Transition Matrix is\n") +disp(P) +printf("\nc) The Transition figure is displayed in the result file attached to this code\n") diff --git a/1904/CH11/EX11.9/11_9.sce b/1904/CH11/EX11.9/11_9.sce new file mode 100755 index 000000000..b77aa9a50 --- /dev/null +++ b/1904/CH11/EX11.9/11_9.sce @@ -0,0 +1,22 @@ +//To Determine the vector of state probabilities at a specific time +//Page 624 +clc; +clear; + +P=[0.6,0.4;0.3,0.7]; //One Step Transition Matrix + +Po=[0.8,0.2]; //Initial State Probability Vector + +//Funtion to determine the Vector of State Probability +deff('x=VSP(y)','x=(Po*(P^y))') + +P1=VSP(1); //Vector of State Probability at Time t1 +P4=VSP(4); //Vector of State Probability at Time t4 +P8=VSP(8); //Vector of State Probability at Time t8 + +printf('\na) The Vector of State Probability at time t1 is\n') +disp(P1) +printf('\na) The Vector of State Probability at time t4 is\n') +disp(P4) +printf('\na) The Vector of State Probability at time t8 is\n') +disp(P8) -- cgit