11 files changed, 294 insertions, 0 deletions
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
--- /dev/null
+++ b/1904/CH11/EX11.7/11_7.jpg
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
--- /dev/null
+++ b/1904/CH11/EX11.8/11_8.jpg
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)