blob: dda24d54dd1a22dd58be9ed9e0c10705f90336ea (
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
|
clear
clc
clear
clc
Z=[.2 .2 .05]' * %i
S=30
V=11
I=round(S*1e8/(sqrt(3)*V*1e3) )/1e2
E=1
Ia1=E/(Z(1)+(Z(2)*Z(3)/(Z(2)+Z(3))))
Ia1=round(Ia1*1000)/1000
Ia2=-Ia1 * Z(3)/(Z(2)+Z(3))
Ia0=-Ia1 * Z(2)/(Z(2)+Z(3))
a=exp(%i * 2 * %pi/3)
A=[1 1 1
1 a^2 a
1 a a^2
]
Ia=A*[ Ia0 Ia1 Ia2]' * I
mprintf("Line currents: ")
mprintf("\nIa= %.2f ang(0) A",abs(Ia(1)))
mprintf("\nIb= %.2f ang(%.2f) A",abs(Ia(2)),atand(imag(Ia(2))/real(Ia(2)))+180)
mprintf("\nIc= %.2f ang(%.2f) A",abs(Ia(3)),atand(imag(Ia(3))/real(Ia(3))))
If=Ia(2)+Ia(3)
mprintf("\nFault current= %.0f ang(%.2f) A",abs(If),atand(imag(If)/real(If)))
Va1=1-(Ia1*Z(1))
Va2=Va1
Va0=Va1
Va=A*[Va0 Va1 Va2]' * V/sqrt(3)
mprintf("\nLine to neutral voltages: ")
mprintf("\nVa= %.3f ang(%.2f) kV",abs(Va(1)),atand(imag(Va(1))/real(Va(1))))
mprintf("\nVb= %.3f KV",abs(Va(2)))
mprintf("\nVc= %.3f KV",abs(Va(3)))
|
