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
|
//Chapter 10
//Example 10_10
//Page 241
clear;clc;
r=0.25;
xl=0.8;
y=14*1e-6;
vr=66000;
pd=15000;
pfr=0.8;
ll=100;
tr=ll*r;
txl=ll*xl;
ty=y*ll;
z=tr+%i*txl;
l=pd/vr/pfr;
i=pd*1000/vr/pfr;
vr_phasor=vr+%i*0;
ir_phasor=i*(pfr-%i*sin(acos(pfr)));
ic=%i*ty*vr;
is_phasor=ir_phasor+ic;
mag_is=abs(is_phasor);
vd=is_phasor*z;
vs_phasor=vr_phasor+is_phasor*z;
mag_vs=abs(vs_phasor);
reg=(mag_vs-vr)/vr*100;
theta1=atan(-imag(is_phasor)/real(is_phasor));
theta2=atan(imag(vs_phasor)/real(vs_phasor));
thetas=abs(theta1)+theta2;
pfs=cos(thetas);
printf("Total resistance = %.0f ohm \n\n", tr);
printf("Total reactance = %.0f ohm \n\n", txl);
printf("Total susceptance = %.0f ohm \n\n", ty);
printf("Recieving end voltage = %.0f V \n\n", vr);
printf("Load current = %.0f A \n\n", i);
printf("Vr phasor = %.2f+j%.2f \n\n", real(vr_phasor), imag(vr_phasor));
printf("Load current phasor = %.2f+j%.2f \n\n", real(ir_phasor), imag(ir_phasor));
printf("Capacitive current = j%.2f \n\n", imag(ic));
printf("(i) Sending end current = %.2f+j%.2f \n", real(is_phasor), imag(is_phasor));
printf(" Magnitude = %.0f V \n\n", mag_is);
printf("(ii) Voltage drop = %.2f+j%.2f \n", real(vd), imag(vd))
printf(" Sending end voltage = %.2f+j%.2f \n", real(vs_phasor), imag(vs_phasor));
printf(" Magnitude of Vs = %.0f V \n\n", mag_vs);
printf("(iii) %% Regulation = %.2f %% \n\n", reg);
printf("(iv) Phase angle between Vr and Ir = %.2f degrees \n\n", theta1*180/%pi);
printf(" Phase angle between Vr and Vs = %.2f degrees \n\n", theta2*180/%pi);
printf(" Supply power factor angle = %.2f degrees \n\n", thetas*180/%pi);
printf(" Supply power factor = %.2f lag \n\n", pfs);
|
