SCRIP PERHITUNGAN MATLAB
clear all clc a = cosd(120)+1i*sind(120);A = [1 1 1;a^2 a 1;a a^2 1]; Vf = 1; dasar_tegangan = 13.8; %kV dasar_daya = 300000; %VA dasar_arus = dasar_daya / (sqrt(3)*13.8); Z1_saluran = 1i*0.1815; Z2_saluran = Z1_saluran; Z0_saluran = 1i*0.5445; Z1_gen = 1i*0.2; Z2_gen = Z1_gen; Z0_gen = 1i*0.05; Z1_trafo_1 = 1i*0.0857; Z2_trafo_1 = Z1_trafo_1; Z0_trafo_1 = Z1_trafo_1; Z1_trafo_2 = 1i*0.0915; Z2_trafo_2 = Z1_trafo_2; Z0_trafo_2 = Z1_trafo_2; Z1_mot_1 = 1i*0.2745; Z2_mot_1 = Z1_mot_1; Z0_mot_1 = 1i*0.0686; Z1_mot_2 = 1i*0.549; Z2_mot_2 = Z1_mot_2; Z0_mot_2 = 1i*0.1372; XN_gen = 1i*0.3; XN_mot_1 = 1i*0.63;
Y_pos(1,1) = 1/(Z1_gen) + (1/Z1_trafo_1); Y_pos(1,2) = -1/Z1_trafo_1; Y_pos(2,1) = Y_pos(1,2); Y_pos(1,2); Y_pos(2,2) = (1/Z1_trafo_1) + (1/Z1_saluran); Y_pos(1,3) = 0; Y_pos(3,1) = Y_pos(1,3); Y_pos(1,3);
Y_pos(3,3) = (1/Z1_trafo_2) + (1/Z1_saluran); Y_pos(1,4) = 0; Y_pos(4,1) = Y_pos(1,4); Y_pos(2,4) = 0; Y_pos(4,2) = Y_pos(2,4); Y_pos(3,4) = -1/Z1_trafo_2; Y_pos(4,3) = Y_pos(3,4); Y_pos(4,4) = (1/Z1_trafo_2) + ((Z1_mot_1*Z1_mot_2)/(Z1_mot_1+Z1_mot_2)) Y_neg = Y_pos
Y_nol(1,1) = 1/(Z0_gen + XN_gen); Y_nol(2,2) = (1/Z1_trafo_1) + (1/Z1_saluran); Y_nol(3,3) = (1/Z1_trafo_2) + (1/Z1_saluran); Y_nol(4,4) = 1/(Z0_mot_1+ XN_mot_1) Y_nol(1,2) = -1/Z0_saluran; Y_nol(2,1) = Y_nol(1,2); Y_nol(1,3) = 0; Y_nol(3,1) = Y_nol(1,3); Y_nol(1,4) = 0; Y_nol(4,1) = Y_nol(4,1); Y_nol(2,4) = 0; Y_nol(4,2) = Y_nol(2,4); Y_nol(3,4) = 0; Y_nol(4,3) = Y_nol(3,4);
Z_pos = inv(Y_pos) Z_neg = inv(Y_neg) Z_nol = inv(Y_nol)
Z_positif_44 = Z_pos(4,4) Z_negatif_44 = Z_neg(4,4) Z_nol_44 = Z_nol(4,4)
I_a1 = Vf/((Z_positif_44)+(Z_negatif_44)+(Z_nol_44)) I_a2 = I_a1
I_a0 = I_a1 If = I_a1 + I_a2 + I_a0 Arus_gangguan = If*dasar_arus
V_a1_4 = Vf - (Z_positif_44 * I_a1) V_a2_4 = -(Z_negatif_44*I_a2) V_a0_4 = -(Z_nol_44*I_a0)
VA = V_a1_4 + V_a2_4 + V_a0_4 VB = a^2*V_a1_4 + a*V_a2_4 + V_a0_4 VC = a*V_a1_4 + a^2*V_a2_4 + V_a0_4
Magnitude_VA = abs(VA) % PU Magnitude_VB = abs(VB) % PU Magnitude_VC = abs(VC) % PU Teg_fasa_sehat_VA = abs(VA) * dasar_tegangan % kV Teg_fasa_sehat_VB = abs(VB) * dasar_tegangan % kV Teg_fasa_sehat_VC = abs(VC) * dasar_tegangan % kV Sudut_VA = angle(VA)*180/3.14 Sudut_VB = angle(VB)*180/3.14 Sudut_VC = angle(VC)*180/3.14
FP = (Teg_fasa_sehat_VB/dasar_tegangan)*100 % dalam persentase
