summaryrefslogtreecommitdiff
path: root/1319/CH5/EX5.5/5_5.sce
blob: e641c98e8dd927d3acbad2c31839b31b9cb2692f (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
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
//Voltage regulation at 0.8 pf lagging

clc;
clear;

E1=1100;
E2=110;

P=5*(10^3);

I=P/E1; // Primary full load current

I2=P/E2;// Secondary full load current

V=33;

pf=0.8; // Power Factor lagging, so the angle is positive

theta=acosd(pf);// Power factor angle

Pc=85;

R=Pc/(I^2);

Z=V/I;

V1=E1;

X=sqrt((Z^2)-(R^2));

// Using equation 5.22 to determine V2

V2=poly([0 1],'V2','c');

x=(V2^2)+(2*V2*I*R*pf)+(2*V2*I*X*sind(theta))+((I^2)*((R^2)+(X^2)))-(V1^2);

r=roots(x);

a1=sqrt(r(1));
a2=sqrt(r(2));

if(imag(a1))
    V2=r(2);
else
    if(imag(a2))
        V2=r(1);
    else
        disp('Error')
    end
end

reg=(V1-V2)/V2; // Voltage regulation

regper=reg*100;// Voltage regulation percent

disp(x,'The characteristic equation to find V2 equated to zero is')

disp(regper,'The percentage voltage regulation for a load at 0.8 pf lagging is')