initial commit / add all books

10 files changed, 172 insertions, 0 deletions

diff --git a/416/CH17/EX17.1/exp17_1PP.sce b/416/CH17/EX17.1/exp17_1PP.sce
new file mode 100755
index 000000000..c8940faf5
--- /dev/null
+++ b/416/CH17/EX17.1/exp17_1PP.sce
@@ -0,0 +1,15 @@
+clc

+clear

+disp("example 17.1")

+//given

+n=2 //number of generating station

+f=0.03 //F.O.R

+a=1-f

+p=40 //generation station power 

+function [y]=comb(m,r)

+y=factorial(m)/(factorial(m-r)*factorial(r))

+endfunction

+for i=0:n

+    pg(i+1)=comb(n,i)*((f)^i)*((a)^(n-i))

+    printf("\nnumber of units out %d ,capacity out  %dMW ,capacity available %dMW ,probability %4f ",i,p*i,p*(n-i),pg(i+1))

+end

diff --git a/416/CH17/EX17.1/exp17_1PP.txt b/416/CH17/EX17.1/exp17_1PP.txt
new file mode 100755
index 000000000..2340f124f
--- /dev/null
+++ b/416/CH17/EX17.1/exp17_1PP.txt
@@ -0,0 +1,6 @@
+ 

+ example 17.1   

+

+number of units out 0 ,capacity out  0MW ,capacity available 80MW ,probability 0.940900 

+number of units out 1 ,capacity out  40MW ,capacity available 40MW ,probability 0.058200 

+number of units out 2 ,capacity out  80MW ,capacity available 0MW ,probability 0.000900  
\ No newline at end of file
diff --git a/416/CH17/EX17.2/exp17_2PP.sce b/416/CH17/EX17.2/exp17_2PP.sce
new file mode 100755
index 000000000..d07980936
--- /dev/null
+++ b/416/CH17/EX17.2/exp17_2PP.sce
@@ -0,0 +1,38 @@
+clc

+clear

+disp("example 17 2")

+//given

+n1=2 //number of generating station

+f1=0.03 //F.O.R

+a1=1-f1

+p1=40 //genetaion station power

+n2=1  //number of genreting station 

+f2=0.03 //F.O.R for second set

+a2=1-f2

+p2=30 //generating station power in second set 

+function [y]=comb(m,r)

+y=factorial(m)/(factorial(m-r)*factorial(r))

+endfunction

+for i=0:n2

+    pg2(i+1)=comb(n2,i)*((f2)^i)*((a2)^(n2-i))

+   co2(i+1)=p2*i;ca2(i+1)=p2*(n2-i)

+    printf("\nnumber of units out %d ,capacity out  %dMW ,capacity available  %dMW ,probability %4f ",i,co2(i+1),ca2(i+1),pg2(i+1))

+end

+printf("\nfor exp 17 1 ")

+for i=0:n1

+    pg1(i+1)=comb(n1,i)*((f1)^i)*((a1)^(n1-i))

+    co1(i+1)=p1*i;ca1(i+1)=p1*(n1-i)

+    printf("\nnumber of units out %d ,capacity out  %dMW ,capacity available %dMW ,probability %4f ",i,co1(i+1),ca1(i+1),pg1(i+1))

+end

+printf("\ncombination of 2 set of stations")

+tp=1

+pocg=0

+for i=0:n1

+    for j=0:n2

+        og=co1(i+1)+co2(j+1) //now total system capacity out

+        cg=ca1(i+1)+ca2(j+1)  //now total system capacity available

+        tp=tp-pocg

+        pocg=pg1(i+1)*pg2(j+1) //individual stste probability

+        printf("\ncapacity out %dMW ,capacity available %dMW ,individual state probability %.6f ,cumulative probability %.6f",og,cg,pocg,tp)

+    end

+end
\ No newline at end of file
diff --git a/416/CH17/EX17.2/exp17_2PP.txt b/416/CH17/EX17.2/exp17_2PP.txt
new file mode 100755
index 000000000..003c6878b
--- /dev/null
+++ b/416/CH17/EX17.2/exp17_2PP.txt
@@ -0,0 +1,16 @@
+

+ example 17 2   

+

+number of units out 0 ,capacity out  0MW ,capacity available  30MW ,probability 0.970000 

+number of units out 1 ,capacity out  30MW ,capacity available  0MW ,probability 0.030000 

+for exp 17 1 

+number of units out 0 ,capacity out  0MW ,capacity available 80MW ,probability 0.940900 

+number of units out 1 ,capacity out  40MW ,capacity available 40MW ,probability 0.058200 

+number of units out 2 ,capacity out  80MW ,capacity available 0MW ,probability 0.000900 

+combination of 2 set of stations

+capacity out 0MW ,capacity available 110MW ,individual state probability 0.912673 ,cumulative probability 1.000000

+capacity out 30MW ,capacity available 80MW ,individual state probability 0.028227 ,cumulative probability 0.087327

+capacity out 40MW ,capacity available 70MW ,individual state probability 0.056454 ,cumulative probability 0.059100

+capacity out 70MW ,capacity available 40MW ,individual state probability 0.001746 ,cumulative probability 0.002646

+capacity out 80MW ,capacity available 30MW ,individual state probability 0.000873 ,cumulative probability 0.000900

+capacity out 110MW ,capacity available 0MW ,individual state probability 0.000027 ,cumulative probability 0.000027 
\ No newline at end of file
diff --git a/416/CH17/EX17.3/exp17_3.jpg b/416/CH17/EX17.3/exp17_3.jpg
new file mode 100755
index 000000000..32e9b017a
--- /dev/null
+++ b/416/CH17/EX17.3/exp17_3.jpg
diff --git a/416/CH17/EX17.3/exp17_3PP.sce b/416/CH17/EX17.3/exp17_3PP.sce
new file mode 100755
index 000000000..ab917f249
--- /dev/null
+++ b/416/CH17/EX17.3/exp17_3PP.sce
@@ -0,0 +1,31 @@
+clc

+clear

+disp("example 17 3")

+//given

+n=4 //number of generating station

+f=0.05 //F.O.R

+a=1-f

+p=50 //generation station power 

+mp=150 //maximum alowable  power

+lf=50 //load factor in persentage

+function [y]=comb(m,r)

+y=factorial(m)/(factorial(m-r)*factorial(r))

+endfunction

+for i=0:n

+    pg(i+1)=comb(n,i)*((f)^i)*((a)^(n-i))

+    co(i+1)=p*i;ca(i+1)=p*(n-i)

+    printf("\nnumber of units out %d ,capacity out %dMW ,capacity available %dMW ,probability %4f ",i,co(i+1),ca(i+1),pg(i+1))

+end

+ld=mp:-lf:0

+[m n]=size(ld)

+plot(ld)

+tg(n-1)=round(10000/(n-1))/100

+tg(n)=tg(n-1)*2

+tg(n+1)=100

+tg(2)=0;tg(1)=0 //maximum load limit

+for i=0:n

+    el(i+1)=pg(i+1)*tg(i+1)

+    printf("\nnumber of units out %d ,capacity out %dMW ,capacity available %dMW ,probability %4f ,tg in persentage %.2f ,expected load %.6fMW",i,co(i+1),ca(i+1),pg(i+1),tg(i+1),el(i+1))

+end

+lt=sum(el)

+printf("\n\nexpected loss of load is %.6fMW percent of time. assuming 365 days in a year, then expected loss of load is %.3fMW days per year",lt,lt*365/100)
\ No newline at end of file
diff --git a/416/CH17/EX17.3/exp17_3PP.txt b/416/CH17/EX17.3/exp17_3PP.txt
new file mode 100755
index 000000000..9708c03c9
--- /dev/null
+++ b/416/CH17/EX17.3/exp17_3PP.txt
@@ -0,0 +1,15 @@
+  

+ example 17 3   

+

+number of units out 0 ,capacity out 0MW ,capacity available 200MW ,probability 0.814506 

+number of units out 1 ,capacity out 50MW ,capacity available 150MW ,probability 0.171475 

+number of units out 2 ,capacity out 100MW ,capacity available 100MW ,probability 0.013538 

+number of units out 3 ,capacity out 150MW ,capacity available 50MW ,probability 0.000475 

+number of units out 4 ,capacity out 200MW ,capacity available 0MW ,probability 0.000006 

+number of units out 0 ,capacity out 0MW ,capacity available 200MW ,probability 0.814506 ,tg in persentage 0.00 ,expected load 0.000000MW

+number of units out 1 ,capacity out 50MW ,capacity available 150MW ,probability 0.171475 ,tg in persentage 0.00 ,expected load 0.000000MW

+number of units out 2 ,capacity out 100MW ,capacity available 100MW ,probability 0.013538 ,tg in persentage 33.33 ,expected load 0.451205MW

+number of units out 3 ,capacity out 150MW ,capacity available 50MW ,probability 0.000475 ,tg in persentage 66.66 ,expected load 0.031664MW

+number of units out 4 ,capacity out 200MW ,capacity available 0MW ,probability 0.000006 ,tg in persentage 100.00 ,expected load 0.000625MW

+

+expected loss of load is 0.483493MW percent of time. assuming 365 days in a year, then expected loss of load is 1.765MW days per year 
\ No newline at end of file
diff --git a/416/CH17/EX17.4/exp17_4.jpg b/416/CH17/EX17.4/exp17_4.jpg
new file mode 100755
index 000000000..98afe1019
--- /dev/null
+++ b/416/CH17/EX17.4/exp17_4.jpg
diff --git a/416/CH17/EX17.4/exp17_4PP.sce b/416/CH17/EX17.4/exp17_4PP.sce
new file mode 100755
index 000000000..3f5416516
--- /dev/null
+++ b/416/CH17/EX17.4/exp17_4PP.sce
@@ -0,0 +1,34 @@
+clc 

+clear

+disp("example 17 4")

+//given

+n=4 //number of generating station

+f=0.02 //F.O.R

+a=1-f

+p=50 //generation station power 

+mp=150 //maximum alowable  power

+minp=30 //minimum power 

+lf=60 //load factor in persentage

+function [y]=comb(m,r)

+y=factorial(m)/(factorial(m-r)*factorial(r))

+endfunction

+for i=0:n

+    pg(i+1)=comb(n,i)*((f)^i)*((a)^(n-i))

+    co(i+1)=p*i;ca(i+1)=p*(n-i)

+    printf("\nnumber of units out %d ,capacity out %dMW ,capacity available %dMW ,probability %.7f ",i,co(i+1),ca(i+1),pg(i+1))

+end

+ld=mp:-lf:minp

+[m n1]=size(ld)

+[mm m]=max(co)

+plot(ld)

+tg(1)=0

+for i=2:n+1

+    tg(i)=(mp-ca(i))*100/(2*lf) //percentage time 

+end

+disp("")

+for i=1:n+1

+    el(i)=pg(i)*tg(i)

+    printf("\nnumber of units out %d ,capacity out %dMW ,capacity available %dMW ,probability %4f ,tg in persentage %.2f ,expected load %.6fMW",i-1,co(i),ca(i),pg(i),tg(i),el(i))

+end

+lt=sum(el)

+printf("\n\nexpected loss of load is %.6fMW percent of time. assuming 365 days in a year, then expected loss of load is %.3fMW days per year ,some times the loss of load is also expressed as reciprocal of this figure and then the units are years per day this result is %.4fMW years per day.",lt,lt*365/100,100/(lt*365))
\ No newline at end of file
diff --git a/416/CH17/EX17.4/exp17_4PP.txt b/416/CH17/EX17.4/exp17_4PP.txt
new file mode 100755
index 000000000..6121c721b
--- /dev/null
+++ b/416/CH17/EX17.4/exp17_4PP.txt
@@ -0,0 +1,17 @@
+   

+ example 17 4   

+

+number of units out 0 ,capacity out 0MW ,capacity available 200MW ,probability 0.9223682 

+number of units out 1 ,capacity out 50MW ,capacity available 150MW ,probability 0.0752954 

+number of units out 2 ,capacity out 100MW ,capacity available 100MW ,probability 0.0023050 

+number of units out 3 ,capacity out 150MW ,capacity available 50MW ,probability 0.0000314 

+number of units out 4 ,capacity out 200MW ,capacity available 0MW ,probability 0.0000002  

+    

+

+number of units out 0 ,capacity out 0MW ,capacity available 200MW ,probability 0.922368 ,tg in persentage 0.00 ,expected load 0.000000MW

+number of units out 1 ,capacity out 50MW ,capacity available 150MW ,probability 0.075295 ,tg in persentage 0.00 ,expected load 0.000000MW

+number of units out 2 ,capacity out 100MW ,capacity available 100MW ,probability 0.002305 ,tg in persentage 41.67 ,expected load 0.096040MW

+number of units out 3 ,capacity out 150MW ,capacity available 50MW ,probability 0.000031 ,tg in persentage 83.33 ,expected load 0.002613MW

+number of units out 4 ,capacity out 200MW ,capacity available 0MW ,probability 0.000000 ,tg in persentage 125.00 ,expected load 0.000020MW

+

+expected loss of load is 0.098673MW percent of time. assuming 365 days in a year, then expected loss of load is 0.360MW days per year ,some times the loss of load is also expressed as reciprocal of this figure and then the units are years per day this result is 2.7766MW years per day. 
\ No newline at end of file