summaryrefslogtreecommitdiff
path: root/1055/CH21/EX21.1/ch21.sce
blob: d9f058496299f75792231d09b2f9e1c720ea98cd (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
//Determine the load bus voltage
clear 
clc;
load1=10+%i*15;//load per phase(MVA)
SCC=250/3;
V=11/sqrt(3);
P=30;
Q=45;
Z=(11/sqrt(3))^2/(250/3);//Equivalent short circuit impedence
dsc=atand(5);
R=.0949;
X=.4746;
//Using equation: V^2= (Vcosd+PR/V)^2 + (Vsind+QX/V)^2, we get 
y=poly([51.7 0 -27.5 0 1],'V','c');
disp(y,"we get equation :");
X=roots(y);
disp(X,"Roots of above equation are ");
V=5.046;
mprintf("V=%.3f\n",V);
dV=6.35-V;
Ssc=250;
//using expression ,a=dV/v=1(Pcos(dsc)+Qsin(dsc))/Ssc +j(Psin(dsc)-Qcos(dsc))/Ssc
a=(P*cosd(dsc)+Q*sind(dsc))/Ssc +%i*(P*sind(dsc)-Q*cosd(dsc))/Ssc;
disp(abs(a),"dV/V= ");