diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /1076/CH9/EX9.12 | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '1076/CH9/EX9.12')
| -rwxr-xr-x | 1076/CH9/EX9.12/9_12.sce | 97 |
1 files changed, 97 insertions, 0 deletions
diff --git a/1076/CH9/EX9.12/9_12.sce b/1076/CH9/EX9.12/9_12.sce new file mode 100755 index 000000000..a8ef3be00 --- /dev/null +++ b/1076/CH9/EX9.12/9_12.sce @@ -0,0 +1,97 @@ +clear
+clc
+
+ey=1/(.05 + %i*.15)
+y=[
+1 2 ey
+1 5 ey
+2 5 ey
+2 3 ey
+3 4 ey
+4 5 ey
+]
+n=5
+e=6
+
+Y=zeros(n,n)
+
+for i=1:e
+ n1=real(y(i,1))
+ n2=real(y(i,2))
+ ynew=y(i,3)
+ if(n1==0)
+ Y(n2,n2)=Y(n2,n2)+ynew
+ else
+ Y(n1,n1)= Y(n1,n1)+ynew
+ Y(n1,n2)= Y(n1,n2)-ynew
+ Y(n2,n1)= Y(n2,n1)-ynew
+ Y(n2,n2)= Y(n2,n2)+ynew
+ end
+end
+// bus no | PL | QL | PG | QG | V | th | btype | Qmin | Qmax
+
+data = [
+1 1 .5 %inf %inf 1.02 0 1 %inf %inf
+2 0 0 2 %inf 1.02 0 2 .3 %inf
+3 .5 .2 0 0 1 0 3 %inf %inf
+4 .5 .2 0 0 1 0 3 %inf %inf
+5 .5 .2 0 0 1 0 3 %inf %inf
+
+]
+
+j = data(:,1); // Bus number.
+PL = data(:,2);
+QL = data(:,3);
+PG = data(:,4);
+QG = data(:,5);
+V = data(:,6);
+th = data(:,7);
+btype = data(:,8); // Type of Bus 1-Slack, 2-PV, 3-PQ.
+Qmin = data(:,9);
+Qmax = data(:,10);
+n = max(j);
+P = PG - PL;
+Q = QG - QL
+
+for i=1:n
+ V(i,1)=V(i,1) * exp (%i * th(i,1)* %pi / 180)
+end
+
+Vprev = V;
+toler = 1000;
+iteration = 1;
+while (iteration==1)
+ for i = 2:n
+ summ = 0;
+ for k = 1:n
+ if i ~= k
+ summ = summ + Y(i,k)* V(k);
+ end
+ end
+ if btype(i) == 2
+ Q(i) = -imag(conj(V(i))*(summ + Y(i,i)*V(i)));
+ if (Q(i) > Qmax(i)) | (Q(i) < Qmin(i))
+ if Q(i) < Qmin(i)
+ Q(i) = Qmin(i);
+ else
+ Q(i) = Qmax(i);
+ end
+ btype(i) = 3;
+ end
+ end
+ V(i) = (1/Y(i,i))*(((P(i)- %i*Q(i))/conj(V(i))) - summ);
+ if btype(i) == 2
+ V(i) = abs(Vprev(i))*exp(%i * atan(imag(V(i))/real(V(i))));
+ end
+ end
+ iteration = iteration + 1;
+ toler = max(abs(abs(V) - abs(Vprev)));
+ Vprev = V;
+end
+V=round(V*1e3)/1e3
+ansmat(:,1)=[1;2;3;4;5]
+ansmat(:,2)=round(abs(V(:,1))*1000)/1000
+for i=1:5
+ ansmat(i,3)=round(atand(imag(V(i))/real(V(i)))*1e3)/1e3
+end
+disp(ansmat," bus no |V| delta")
|