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
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
|
clc
// Fundamental of Electric Circuit
// Charles K. Alexander and Matthew N.O Sadiku
// Mc Graw Hill of New York
// 5th Edition
// Part 2 : AC Circuits
// Chapter 12 : Three Phase Circuit
// Example 12 - 10
clear; clc; close; // Clear the work space and console.
//
// Given data
V1_mag = 120.0000 * sqrt(3);
V1_angle = 30.0000;
V2_mag = 120.0000 * sqrt(3);
V2_angle = -90.000;
Van_mag = 120.0000;
Van_angle = 0.0000;
Vbn_mag = Van_mag;
Vbn_angle = -120.0000;
Vcn_mag = Van_mag;
Vcn_angle = 120.0000;
Zan = complex(0,5);
Zbn = complex(10,0);
Zcn = complex(0,-10);
V1 = complex(V1_mag*cosd(V1_angle),V1_mag*sind(V1_angle));
V2 = complex(V2_mag*cosd(V2_angle),V2_mag*sind(V2_angle));
Z11 = complex(10,5)
Z12 = complex(-10,0);
Z21 = Z12;
Z22 = complex(10,-10);
// Calculations Determinants
delta = det([Z11 Z12; Z21 Z22]);
delta_mag = norm(delta);
delta_real = real(delta);
delta_imag = imag(delta);
delta_angle = atand(delta_imag,delta_real);
// Calculations Determinants For Source Voltage V1
delta1 = det([V1 Z12; V2 Z22]);
delta1_mag = norm(delta1);
delta1_real = real(delta1);
delta1_imag = imag(delta1);
delta1_angle = atand(delta1_imag,delta1_real);
// Calculations Determinants For Source Voltage V2
delta2 = det([Z11 V1; Z12 V2]);
delta2_mag = norm(delta2);
delta2_real = real(delta2);
delta2_imag = imag(delta2);
delta2_angle = atand(delta2_imag,delta2_real);
// Calculations Mesh Current I1
I1_mag = delta1_mag/delta_mag;
I1_angle = delta1_angle - delta_angle
// Calculations Mesh Current I2
I2_mag = delta2_mag/delta_mag;
I2_angle = delta2_angle - delta_angle
// Calculations Line Current a
I_a_mag = I1_mag;
I_a_angle = I1_angle
// Calculations Line Current b
I1 = complex(I1_mag*cosd(I1_angle),I1_mag*sind(I1_angle))
I2 = complex(I2_mag*cosd(I2_angle),I2_mag*sind(I2_angle))
I_b_mag = norm(I2 - I1);
I_b_angle = atand(imag(I2 - I1),real(I2 - I1))
// Calculations Line Current c
I_c_mag = I2_mag
I_c_angle = -180.0000 + I2_angle;
// Calculations Power Absorbed by the Load for Phase A
Sal =(I_a_mag)^2*Zan;
// Calculations Power Absorbed by the Load for Phase B
Sbl =(I_b_mag)^2*Zbn;
// Calculations Power Absorbed by the Load for Phase C
Scl =(I_c_mag)^2*Zcn;
// Calculations Total Complex Power Absorbed by the Load
Stl = Sal + Sbl + Scl;
Stl_real = real(Stl);
Stl_imag = imag(Stl);
// Calculations Power Absorbed by the Source for Phase A
Sas =(complex(Van_mag*cosd(Van_angle),Van_mag*sind(Van_angle)))*conj(complex(I_a_mag*cosd(I_a_angle),I_a_mag*sind(I_a_angle)));
// Calculations Power Absorbed by the Load for Phase B
Sbs =(complex(Vbn_mag*cosd(Vbn_angle),Vbn_mag*sind(Vbn_angle)))*conj(complex(I_b_mag*cosd(I_b_angle),I_b_mag*sind(I_b_angle)));
// Calculations Power Absorbed by the Load for Phase C
Scs =(complex(Vcn_mag*cosd(Vcn_angle),Vcn_mag*sind(Vcn_angle)))*conj(complex(I_c_mag*cosd(I_c_angle),I_c_mag*sind(I_c_angle)));
// Calculations Total Complex Power Absorbed by The Source
Sts = Sas + Sbs + Scs;
Sts_real = -real(Sts);
Sts_imag = -imag(Sts);
//
disp("Example 12-10 Solution : ");
disp("a. The Line Currents: ");
printf(" \n I_a_mag = Magnitude of Line Currents a = %.3f A",I_a_mag)
printf(" \n I_a_angle = Angle of Line Currents a = %.3f degree",I_a_angle)
printf(" \n I_b_mag = Magnitude of Line Currents b = %.3f A",I_b_mag)
printf(" \n I_b_angle = Angle of Line Currents b = %.3f degree",I_b_angle)
printf(" \n I_c_mag = Magnitude of Line Currents c = %.3f A",I_c_mag)
printf(" \n I_c_angle = Angle of Line Currents c = %.3f degree",I_c_angle)
disp("")
disp("b Total Power Complex Absorbed By The Load");
printf(" \n Stl_real = Real Part of Power Complex = %.3f Watt",Stl_real)
printf(" \n Stl_imag = Imaginary Part of Power Complex = %.3f Var",Stl_imag)
disp("")
disp("c. Total Power Complex Absorbed By The Source");
printf(" \n Sts_real = Real Part of Power Complex = %.3f Watt",Sts_real)
printf(" \n Sts_imag = Imaginary Part of Power Complex = %.3f Var",Sts_imag)
|