diff options
Diffstat (limited to '1067/CH46')
-rwxr-xr-x | 1067/CH46/EX46.01/46_01.sce | 39 | ||||
-rwxr-xr-x | 1067/CH46/EX46.02/46_02.sce | 7 | ||||
-rwxr-xr-x | 1067/CH46/EX46.03/46_03.sce | 11 |
3 files changed, 57 insertions, 0 deletions
diff --git a/1067/CH46/EX46.01/46_01.sce b/1067/CH46/EX46.01/46_01.sce new file mode 100755 index 000000000..e10a27758 --- /dev/null +++ b/1067/CH46/EX46.01/46_01.sce @@ -0,0 +1,39 @@ +clear;
+clc;
+//for low loads
+p1(1)=20;
+p2(1)=30;
+t1(1)=.1*p1(1)+20;
+t2(1)=.12*p2(1)+16;
+//when load is further increased
+t2(4)=22;
+p2(4)=(t2(4)-16)/.12;
+t1(4)=t2(4);
+//upper limit 125MW
+p2(5)=125;
+t1(5)=1.12*p2(5)+16;
+p1(5)=(t1(5)-20)/.1;
+n=7;
+t2(1)=19.6;
+t2(2)=20;
+t2(3)=21;
+t2(4)=22;
+t2(5)=31;
+t2(6)=32;
+t2(7)=32.5;
+p1(5)=110;
+p1(6)=120;
+p1(7)=125;
+for j=1:4
+ p1(j)=20;
+end;
+mprintf("incremental cost(rs./MWhr)\tloading of unit 1(MW) \t loading of unit 2(MW)\ttotal generating power(MW)");
+for i=1:n
+ p2(i)=(-16+t2(i))/.12;
+ if(t2(i)>=31)
+ p2(i)=125;
+ end;
+ pt(i)=p1(i)+p2(i);
+ mprintf("\n%f\t\t\t%f\t\t\t%f\t\t%f",t2(i),p1(i),p2(i),pt(i));
+
+end;
diff --git a/1067/CH46/EX46.02/46_02.sce b/1067/CH46/EX46.02/46_02.sce new file mode 100755 index 000000000..0d3546cd4 --- /dev/null +++ b/1067/CH46/EX46.02/46_02.sce @@ -0,0 +1,7 @@ +clear;
+clc;
+p=180;
+p2=(20-16+(180*.1))/(.1+.12);
+p1=p-p2;
+t=.1*p1+20;
+mprintf("loading of unit 1 P1=%dMW\nthe loading of unit 2 P2=%dMW\nincremental operating cost =%dRs/MWhr",p1,p2,t);
diff --git a/1067/CH46/EX46.03/46_03.sce b/1067/CH46/EX46.03/46_03.sce new file mode 100755 index 000000000..ef8bec90e --- /dev/null +++ b/1067/CH46/EX46.03/46_03.sce @@ -0,0 +1,11 @@ +clear;
+clc;
+p11=80;
+p12=90;
+p21=100;
+p22=90;
+x=integrate('.1*x+20','x',p11,p12);
+y=integrate('.2*x+6','x',p21,p22);
+p=x+y;
+as=p*8760;
+mprintf("economic loading for unit 1=%dRs/hr\neconomic loading for unit 2=%dRs/hr\nannual savings=%dRs",x,y,as);
|