MATLAB

%=== menentukan kurva ayunan selama 2 detik dan waktu pemutus

gangguan ===
%===pada 0. detik melalui penyelesaian persamaan numerik dengan
metode==
%== rungu-kutta orde ke 4 pada mesin majemuk
===
%== Di PT Inracawood MFG
===
format short g
delta_t = 0.01;
fo = 50;
% data pada generator:
Xd1 = 1i*0.21;
H1 = 8.2;
Xd2 = 1i*0.22;
H2 = 8.0;
Xd3 = 1i*0.21;
H3 = 8.2;
%data saluran data transformator:
%=========================================================
%== Impendansi ! Y-shunt ! bus ke bus ==
%== (per unit) ! (per unit) ! ==
%==========================================================
Z_15 =0.060+1i*0.229; Ysh_15=1i*0.022; %Saluran 1-5 %
Z_25 =0.060+1i*0.229; Ysh_25=1i*0.022; %Saluran 2-5 %
Z_26 = 0.060+1i*0.229; Ysh_26=1i*0.022; %Saluan 2-6 %
Z_34 = 0.060+1i*0.229; Ysh_34=1i*0.022; %Saluan 3-4 %
Z_36 = 0.060+1i*0.229; Ysh_36=1i*0.022; %Saluran 3-6 %
Z_46 = 0.060+1i*0.229; Ysh_46=1i*0.022; %Saluran 4-6 %
Z_56 = 0.060+1i*0.229; Ysh_56=1i*0.022; %Saluran 5-6 %
%==========================================================
% data harga harga nilai aliran daya sebelum gangguan :
%===============================================================
=========
%== tegangan ! pembangkit ! beban
! !
%== (per unit) ! (per unit) ! (per unit)
! bus !
 %==-------------------------------------------------------------
-------==
%== magnitude ! sudut ! P ! Q ! P ! Q
! !
%===============================================================
=========
V1 =1; sd_V1 = 14.07; P1=
1.319;Q1=1i*0.023;Pb1=0.76;Qb1=1i*0.013; % 1 %
V2 =1; sd_V2 = 14.07; P2= 3.88;
Q2=1i*0.067;Pb2=2.218;Qb2=1i*0.038;% 2 %
V3 =1; sd_V3 = 14.07; P3=
1.319;Q5=1i*0.023;Pb3=0.88;Qb3=1i*0.011; % 3 %
v4 =1; sd_V4 =0;
%4%
%===============================================================
==========
% arus yang mengalir ke jala jala dari bus 1 dan 2
I1 = (Pb1 -
Qb1)/conj(V1*(cos(sd_V1*pi/180)+1i*sin(sd_V1*pi/180)))

I1 =
0.74036 + 0.17215i

I2 = (Pb2 - Qb2)/conj
(V2*(cos(sd_V2*pi/180)+1i*sin(sd_V2*pi/180)))

I2 =
2.1607 + 0.50235i

I3 = (Pb3 - Qb3)/conj
(V3*(cos(sd_V3*pi/180)+1i*sin(sd_V3*pi/180)))

I3 =
0.85627 + 0.20326i

% tegangan pada masing masing generator dan infine bus :

E1 = V1*(cos(sd_V1*pi/180)+1i*sin (sd_V1*pi/180))+Xd1*I1

E1 =
0.93385 + 0.39858i

E2 = V2*(cos(sd_V2*pi/180)+1i*sin (sd_V2*pi/180))+Xd2*I2
E2 =
0.85948 + 0.71846i

E3 = V3*(cos(sd_V3*pi/180)+1i*sin (sd_V3*pi/180))+Xd3*I3

E3 =
0.92731 + 0.42292i

% kondisi awal :
delta1 = 180/pi*atan2(imag(E1) ,real (E1));
delta2 = 180/pi*atan2(imag(E2) ,real (E2));
delta3 = 180/pi*atan2(imag(E3) ,real (E3));
omega1 = 0;
omega2 = 0;
omega3 = 0;
disp('==========================================================
========')

==================================================================

disp('== harga delta dan omega setiap mesin melalui

penyelesaian ==')

== harga delta dan omega setiap mesin melalui penyelesaian ==

disp('== numerik dengan metode runge-kutta orde 4

==')

== numerik dengan metode runge-kutta orde 4 ==

disp('==========================================================
========')

==================================================================

disp('== waktu ! delta1 ! delta 2 ! delta 3 !

omega1 ! omega2 ! omega 3 ==')

== waktu ! delta1 ! delta 2 ! delta 3 ! omega1 ! omega2

! omega 3 ==

disp('== (detik) ! (derajat)! (derajat) ! (derajat)

!(radian/detik)!(radian/detik)!(radian/detik)==')

== (detik) ! (derajat)! (derajat) ! (derajat)

!(radian/detik)!(radian/detik)!(radian/detik)==
 disp('==========================================================
========')

==================================================================

tp = input (' waktu pemutus (detik)=');

disp('==========================================================
========')

==================================================================

% harga delta dan omega selama gangguan :

for t = 0 : 0.01 : (tp -0.01)
K1m1 = omega1*delta_t
K1m2 = omega2*delta_t
K1m3 = omega3*delta_t
L1m1 = (pi*fo/H1)*(P1)*delta_t
L1m2 =
(pi*fo/H2)*(P2-(2.051-3.561*cos((-74.99-delta2)*pi/180)))*delta_
t
L1m3 =
(pi*fo/H3)*(P3-(0.607-3.861*cos((-74.99-delta3)*pi/180)))*delta_
t
K2m1 = (omega1 + 0.5 *L1m1)*delta_t
K2m2 = (omega2 + 0.5 *L1m2)*delta_t
K2m3 = (omega3 + 0.5 *L1m3)*delta_t
L2m1 = (pi*fo/H1)*(P1)*delta_t
L2m2 = (pi*fo/H2)*(P2-(2.051-3.561*cos((-74.99-
(delta2+(0.5*K1m2*180/pi)))*pi/180)))*delta_t
L2m3 = (pi*fo/H3)*(P3-(0.607-3.861*cos((-74.99-
(delta3+(0.5*K1m3*180/pi)))*pi/180)))*delta_t
K3m1 = (omega1 + 0.5 *L2m1)*delta_t
K3m2 = (omega2 + 0.5 *L2m2)*delta_t
K3m3 = (omega3 + 0.5 *L2m3)*delta_t
L3m1 = (pi*fo/H1)*(P1)*delta_t
L3m2 = (pi*fo/H2)*(P2-(2.051-3.561*cos((-74.99-
(delta2+(0.5*K2m2*180/pi)))*pi/180)))*delta_t
L3m3 = (pi*fo/H3)*(P3-(0.607-3.861*cos((-74.99-
(delta3+(0.5*K2m3*180/pi)))*pi/180)))*delta_t
K4m1 = (omega1 + 0.5 *L3m1)*delta_t
K4m2 = (omega2 + 0.5 *L3m2)*delta_t
K4m3 = (omega3 + 0.5 *L3m3)*delta_t
L4m1 = (pi*fo/H1)*(P1)*delta_t
 L4m2 = (pi*fo/H2)*(P2-(2.051-3.561*cos((-74.99-(delta2+
(K3m2*180/pi)))*pi/180)))*delta_t
L4m3 = (pi*fo/H3)*(P3-(0.607-3.861*cos((-74.99-(delta3+
(K3m3*180/pi)))*pi/180)))*delta_t
fprintf('!%7.2f !%10.3f !%10.3f !%10.3f !%9.3f !%10.3f
!%10.3f! \n',t,delta1,delta2,delta3,omega1,omega2,omega3)
delta1 = delta1 + ((1/6)*((K1m1 + 2*K2m1 + 2*K3m1 +
K4m1)*180/pi));
delta2 = delta2 + ((1/6)*((K1m2 + 2*K2m2 + 2*K3m2 +
K4m2)*180/pi));
delta3 = delta3 + ((1/6)*((K1m3 + 2*K2m3 + 2*K3m3 +
K4m3)*180/pi));
plot(t,delta1,'x')
plot(t,delta2,'x')
plot(t,delta3,'x')
grid on
title ('kurva ayunan')
xlabel('waktu, detik')
ylabel ('delta, derajat')
hold on
omega1=omega1 + ((1/6)*(L1m1 + 2*L2m1 + 2*L3m1 + L4m1));
omega2=omega2 + ((1/6)*(L1m2 + 2*L2m2 + 2*L3m2 + L4m2));
omega3=omega3 + ((1/6)*(L1m3 + 2*L2m3 + 2*L3m3 + L4m3));
end

K1m1 =
0
K1m2 =
0
K1m3 =
0
L1m1 =
0.25267
L1m2 =
0.064923
L1m3 =
0.014236
K2m1 =
0.0012633
K2m2 =
0.00032461
K2m3 =
7.1181e-05
L2m1 =
0.25267
L2m2 =
0.064923
L2m3 =
0.014236
K3m1 =
0.0012633
K3m2 =
0.00032461
K3m3 =
7.1181e-05
L3m1 =
0.25267
L3m2 =
0.06482
L3m3 =
0.01421
K4m1 =
0.0012633
K4m2 =
0.0003241
K4m3 =
7.1051e-05
L4m1 =
0.25267
L4m2 =
0.064717
L4m3 =
0.014184
! 0.00 ! 23.114 ! 39.893 ! 24.517 ! 0.000 !
0.000 ! 0.000!
K1m1 =
0.0025267
K1m2 =
0.00064854
K1m3 =
0.00014219
L1m1 =
0.25267
L1m2 =
0.064751
L1m3 =
0.014193
K2m1 =
0.00379
K2m2 =
0.0009723
K2m3 =
0.00021315
L2m1 =
0.25267
L2m2 =
0.064546
L2m3 =
0.014141
K3m1 =
0.00379
K3m2 =
0.00097127
K3m3 =
 0.00021289
L3m1 =
0.25267
L3m2 =
0.064443
L3m3 =
0.014115
K4m1 =
0.00379
K4m2 =
0.00097076
K4m3 =
0.00021276
L4m1 =
0.25267
L4m2 =
0.064135
L4m3 =
0.014038
! 0.01 ! 23.174 ! 39.909 ! 24.520 ! 0.253 !
0.065 ! 0.014!
K1m1 =
0.0050534
K1m2 =
0.0012933
K1m3 =
0.00028343
L1m1 =
0.25267
L1m2 =
0.064169
L1m3 =
0.014046
K2m1 =
0.0063167
K2m2 =
0.0016142
K2m3 =
0.00035366
L2m1 =
0.25267
L2m2 =
0.06376
L2m3 =
0.013943
K3m1 =
0.0063167
K3m2 =
0.0016121
K3m3 =
0.00035314
L3m1 =
0.25267
L3m2 =
0.063658
L3m3 =
0.013917
K4m1 =
0.0063167
K4m2 =
0.0016116
K4m3 =
0.00035301
L4m1 =
0.25267
L4m2 =
0.063148
L4m3 =
0.013789
! 0.02 ! 23.379 ! 39.961 ! 24.531 ! 0.505 !
0.129 ! 0.028!
K1m1 =
0.0075801
K1m2 =
0.0019302
K1m3 =
0.00042268
L1m1 =
0.25267
L1m2 =
0.063181
L1m3 =
0.013797
K2m1 =
0.0088434
K2m2 =
0.0022461
K2m3 =
0.00049167
L2m1 =
0.25267
L2m2 =
0.06257
L2m3 =
0.013643
K3m1 =
0.0088434
K3m2 =
0.0022431
K3m3 =
0.0004909
L3m1 =
0.25267
L3m2 =
0.06247
L3m3 =
0.013618
K4m1 =
0.0088434
K4m2 =
 0.0022426
K4m3 =
0.00049077
L4m1 =
0.25267
L4m2 =
0.061761
L4m3 =
0.013439
! 0.03 ! 23.729 ! 40.050 ! 24.551 ! 0.758 !
0.193 ! 0.042!
K1m1 =
0.010107
K1m2 =
0.0025553
K1m3 =
0.00055894
L1m1 =
0.25267
L1m2 =
0.061793
L1m3 =
0.013447
K2m1 =
0.01137
K2m2 =
0.0028642
K2m3 =
0.00062618
L2m1 =
0.25267
L2m2 =
0.060985
L2m3 =
0.013243
K3m1 =
0.01137
K3m2 =
0.0028602
K3m3 =
0.00062516
L3m1 =
0.25267
L3m2 =
0.060887
L3m3 =
0.013219
K4m1 =
0.01137
K4m2 =
0.0028597
K4m3 =
0.00062504
L4m1 =
0.25267
L4m2 =
0.059984
L4m3 =
0.012991
! 0.04 ! 24.224 ! 40.176 ! 24.579 ! 1.011 !
0.256 ! 0.056!
K1m1 =
0.012633
K1m2 =
0.0031645
K1m3 =
0.00069121
L1m1 =
0.25267
L1m2 =
0.060016
L1m3 =
0.012999
K2m1 =
0.013897
K2m2 =
0.0034646
K2m3 =
0.00075621
L2m1 =
0.25267
L2m2 =
0.059016
L2m3 =
0.012747
K3m1 =
0.013897
K3m2 =
0.0034596
K3m3 =
0.00075495
L3m1 =
0.25267
L3m2 =
0.058921
L3m3 =
0.012723
K4m1 =
0.013897
K4m2 =
0.0034591
K4m3 =
0.00075483
L4m1 =
0.25267
L4m2 =
0.057831
L4m3 =
0.012448
! 0.05 ! 24.863 ! 40.337 ! 24.614 ! 1.263 !
0.316 ! 0.069!
K1m1 =
0.01516
K1m2 =
0.003754
K1m3 =
0.00081853
L1m1 =
0.25267
L1m2 =
0.057861
L1m3 =
0.012456
K2m1 =
0.016423
K2m2 =
0.0040433
K2m3 =
0.00088081
L2m1 =
0.25267
L2m2 =
0.056677
L2m3 =
0.012157
K3m1 =
0.016423
K3m2 =
0.0040374
K3m3 =
0.00087931
L3m1 =
0.25267
L3m2 =
0.056586
L3m3 =
0.012135
K4m1 =
0.016423
K4m2 =
0.0040369
K4m3 =
0.0008792
L4m1 =
0.25267
L4m2 =
0.055316
L4m3 =
0.011815
! 0.06 ! 25.647 ! 40.533 ! 24.656 ! 1.516 !
0.375 ! 0.082!
K1m1 =
0.017687
K1m2 =
 0.0043202
K1m3 =
0.00093995
L1m1 =
0.25267
L1m2 =
0.055345
L1m3 =
0.011822
K2m1 =
0.01895
K2m2 =
0.0045969
K2m3 =
0.00099906
L2m1 =
0.25267
L2m2 =
0.053985
L2m3 =
0.011479
K3m1 =
0.01895
K3m2 =
0.0045901
K3m3 =
0.00099735
L3m1 =
0.25267
L3m2 =
0.053898
L3m3 =
0.011458
K4m1 =
0.01895
K4m2 =
0.0045897
K4m3 =
0.00099724
L4m1 =
0.25267
L4m2 =
0.052457
L4m3 =
0.011095
! 0.07 ! 26.576 ! 40.761 ! 24.706 ! 1.769 !
0.432 ! 0.094!
K1m1 =
0.020213
K1m2 =
0.0048595
K1m3 =
0.0010546
L1m1 =
0.25267
L1m2 =
0.052484
L1m3 =
0.011101
K2m1 =
0.021477
K2m2 =
0.0051219
K2m3 =
0.0011101
L2m1 =
0.25267
L2m2 =
0.050958
L2m3 =
0.010717
K3m1 =
0.021477
K3m2 =
0.0051143
K3m3 =
0.0011082
L3m1 =
0.25267
L3m2 =
0.050876
L3m3 =
0.010697
K4m1 =
0.021477
K4m2 =
0.0051138
K4m3 =
0.0011081
L4m1 =
0.25267
L4m2 =
0.049274
L4m3 =
0.010294
! 0.08 ! 27.650 ! 41.022 ! 24.763 ! 2.021 !
0.486 ! 0.105!
K1m1 =
0.02274
K1m2 =
0.0053685
K1m3 =
0.0011616
L1m1 =
0.25267
L1m2 =
0.0493
L1m3 =
0.0103
K2m1 =
 0.024003
K2m2 =
0.005615
K2m3 =
0.0012131
L2m1 =
0.25267
L2m2 =
0.047618
L2m3 =
0.0098763
K3m1 =
0.024003
K3m2 =
0.0056066
K3m3 =
0.001211
L3m1 =
0.25267
L3m2 =
0.047541
L3m3 =
0.0098576
K4m1 =
0.024003
K4m2 =
0.0056062
K4m3 =
0.0012109
L4m1 =
0.25267
L4m2 =
0.04579
L4m3 =
0.0094171
! 0.09 ! 28.868 ! 41.313 ! 24.826 ! 2.274 !
0.537 ! 0.116!
K1m1 =
0.025267
K1m2 =
0.0058442
K1m3 =
0.0012603
L1m1 =
0.25267
L1m2 =
0.045813
L1m3 =
0.0094226
K2m1 =
0.02653
K2m2 =
0.0060733
K2m3 =
0.0013074
L2m1 =
0.25267
L2m2 =
0.043988
L2m3 =
0.0089635
K3m1 =
0.02653
K3m2 =
0.0060641
K3m3 =
0.0013051
L3m1 =
0.25267
L3m2 =
0.043917
L3m3 =
0.0089463
K4m1 =
0.02653
K4m2 =
0.0060638
K4m3 =
0.001305
L4m1 =
0.25267
L4m2 =
0.042029
L4m3 =
0.0084718
! 0.10 ! 30.231 ! 41.632 ! 24.895 ! 2.527 !
0.584 ! 0.126!
K1m1 =
0.027794
K1m2 =
0.0062836
K1m3 =
0.0013498
L1m1 =
0.25267
L1m2 =
0.042049
L1m3 =
0.0084767
K2m1 =
0.029057
K2m2 =
0.0064939
K2m3 =
0.0013922
L2m1 =
0.25267
L2m2 =
0.040093
L2m3 =
 0.0079851
K3m1 =
0.029057
K3m2 =
0.0064841
K3m3 =
0.0013897
L3m1 =
0.25267
L3m2 =
0.040028
L3m3 =
0.0079696
K4m1 =
0.029057
K4m2 =
0.0064837
K4m3 =
0.0013897
L4m1 =
0.25267
L4m2 =
0.038015
L4m3 =
0.0074644
! 0.11 ! 31.739 ! 41.977 ! 24.969 ! 2.779 !
0.628 ! 0.135!
K1m1 =
0.03032
K1m2 =
0.0066841
K1m3 =
0.0014296
L1m1 =
0.25267
L1m2 =
0.038034
L1m3 =
0.0074687
K2m1 =
0.031584
K2m2 =
0.0068743
K2m3 =
0.0014669
L2m1 =
0.25267
L2m2 =
0.03596
L2m3 =
0.0069482
K3m1 =
0.031584
K3m2 =
0.0068639
K3m3 =
0.0014643
L3m1 =
0.25267
L3m2 =
0.035901
L3m3 =
0.0069346
K4m1 =
0.031584
K4m2 =
0.0068636
K4m3 =
0.0014642
L4m1 =
0.25267
L4m2 =
0.033779
L4m3 =
0.0064024
! 0.12 ! 33.392 ! 42.347 ! 25.049 ! 3.032 !
0.668 ! 0.143!
K1m1 =
0.032847
K1m2 =
0.0070433
K1m3 =
0.0014989
L1m1 =
0.25267
L1m2 =
0.033795
L1m3 =
0.006406
K2m1 =
0.03411
K2m2 =
0.0072123
K2m3 =
0.001531
L2m1 =
0.25267
L2m2 =
0.031617
L2m3 =
0.0058603
K3m1 =
0.03411
K3m2 =
0.0072014
K3m3 =
0.0015282
L3m1 =
0.25267
L3m2 =
 0.031565
L3m3 =
0.0058487
K4m1 =
0.03411
K4m2 =
0.0072012
K4m3 =
0.0015282
L4m1 =
0.25267
L4m2 =
0.029346
L4m3 =
0.0052934
! 0.13 ! 35.190 ! 42.739 ! 25.132 ! 3.285 !
0.704 ! 0.150!
K1m1 =
0.035374
K1m2 =
0.0073592
K1m3 =
0.0015575
L1m1 =
0.25267
L1m2 =
0.029361
L1m3 =
0.0052963
K2m1 =
0.036637
K2m2 =
0.007506
K2m3 =
0.001584
L2m1 =
0.25267
L2m2 =
0.027094
L2m3 =
0.0047295
K3m1 =
0.036637
K3m2 =
0.0074947
K3m3 =
0.0015811
L3m1 =
0.25267
L3m2 =
0.027049
L3m3 =
0.0047199
K4m1 =
0.036637
K4m2 =
0.0074944
K4m3 =
0.0015811
L4m1 =
0.25267
L4m2 =
0.024749
L4m3 =
0.0041456
! 0.14 ! 37.132 ! 43.150 ! 25.219 ! 3.537 !
0.736 ! 0.156!
K1m1 =
0.0379
K1m2 =
0.0076299
K1m3 =
0.0016047
L1m1 =
0.25267
L1m2 =
0.024761
L1m3 =
0.0041478
K2m1 =
0.039164
K2m2 =
0.0077537
K2m3 =
0.0016254
L2m1 =
0.25267
L2m2 =
0.02242
L2m3 =
0.0035639
K3m1 =
0.039164
K3m2 =
0.007742
K3m3 =
0.0016225
L3m1 =
0.25267
L3m2 =
0.022382
L3m3 =
0.0035564
K4m1 =
0.039164
K4m2 =
0.0077418
K4m3 =
0.0016225
L4m1 =
 0.25267
L4m2 =
0.020016
L4m3 =
0.0029672
! 0.15 ! 39.219 ! 43.578 ! 25.310 ! 3.790 !
0.763 ! 0.160!
K1m1 =
0.040427
K1m2 =
0.0078538
K1m3 =
0.0016403
L1m1 =
0.25267
L1m2 =
0.020026
L1m3 =
0.0029687
K2m1 =
0.04169
K2m2 =
0.0079539
K2m3 =
0.0016551
L2m1 =
0.25267
L2m2 =
0.017627
L2m3 =
0.0023721
K3m1 =
0.04169
K3m2 =
0.007942
K3m3 =
0.0016522
L3m1 =
0.25267
L3m2 =
0.017596
L3m3 =
0.0023667
K4m1 =
0.04169
K4m2 =
0.0079418
K4m3 =
0.0016521
L4m1 =
0.25267
L4m2 =
0.01518
L4m3 =
0.001767
! 0.16 ! 41.451 ! 44.021 ! 25.403 ! 4.043 !
0.785 ! 0.164!
K1m1 =
0.042954
K1m2 =
0.0080299
K1m3 =
0.001664
L1m1 =
0.25267
L1m2 =
0.015187
L1m3 =
0.0017677
K2m1 =
0.044217
K2m2 =
0.0081058
K2m3 =
0.0016728
L2m1 =
0.25267
L2m2 =
0.012745
L2m3 =
0.0011626
K3m1 =
0.044217
K3m2 =
0.0080936
K3m3 =
0.0016698
L3m1 =
0.25267
L3m2 =
0.012722
L3m3 =
0.0011594
K4m1 =
0.044217
K4m2 =
0.0080935
K4m3 =
0.0016698
L4m1 =
0.25267
L4m2 =
0.010271
L4m3 =
0.00055348
! 0.17 ! 43.828 ! 44.476 ! 25.497 ! 4.295 !
0.803 ! 0.166!
K1m1 =
0.04548
K1m2 =
 0.0081572
K1m3 =
0.0016756
L1m1 =
0.25267
L1m2 =
0.010275
L1m3 =
0.00055345
K2m1 =
0.046744
K2m2 =
0.0082086
K2m3 =
0.0016784
L2m1 =
0.25267
L2m2 =
0.0078062
L2m3 =
-5.5609e-05
K3m1 =
0.046744
K3m2 =
0.0081963
K3m3 =
0.0016753
L3m1 =
0.25267
L3m2 =
0.0077907
L3m3 =
-5.6615e-05
K4m1 =
0.046744
K4m2 =
0.0081962
K4m3 =
0.0016753
L4m1 =
0.25267
L4m2 =
0.0053199
L4m3 =
-0.00066437
! 0.18 ! 46.349 ! 44.939 ! 25.593 ! 4.548 !
0.816 ! 0.168!
K1m1 =
0.048007
K1m2 =
0.0082352
K1m3 =
0.001675
L1m1 =
0.25267
L1m2 =
0.0053214
L1m3 =
-0.00066514
K2m1 =
0.04927
K2m2 =
0.0082618
K2m3 =
0.0016717
L2m1 =
0.25267
L2m2 =
0.0028411
L2m3 =
-0.0012738
K3m1 =
0.04927
K3m2 =
0.0082494
K3m3 =
0.0016687
L3m1 =
0.25267
L3m2 =
0.0028331
L3m3 =
-0.0012726
K4m1 =
0.04927
K4m2 =
0.0082494
K4m3 =
0.0016687
L4m1 =
0.25267
L4m2 =
0.00035842
L4m3 =
-0.0018778
! 0.19 ! 49.015 ! 45.408 ! 25.689 ! 4.801 !
0.824 ! 0.168!

disp('-------------------waktu
pemutus--------------------------')

-------------------waktu pemutus--------------------------

for t = (tp + 0.01) :0.002: 0.25

 K1m1 = omega1*delta_t
K1m2 = omega2*delta_t
K1m3 = omega3*delta_t
L1m1 =
(pi*fo/H1)*(P1-(1.377-9.871*cos((26.58-delta1)*pi/180)))*delta_t
L1m2 =
(pi*fo/H2)*(P2-(0.051-4.974*cos((-80.71-delta2)*pi/180)))*delta_
t
L1m3 =
(pi*fo/H3)*(P3-(1.115-3.954*cos((-73.32-delta3)*pi/180)))*delta_
t
K2m1 = (omega1 + 0.5 *L1m1)*delta_t
K2m2 = (omega2 + 0.5 *L1m2)*delta_t
K2m3 = (omega3 + 0.5 *L1m3)*delta_t
L2m1 = (pi*fo/H1)*(P1-(1.377-9.871*cos((26.58-
(delta1+(0.5*K1m1*180/pi)))*pi/180)))*delta_t
L2m2 = (pi*fo/H2)*(P2-(0.051-4.974*cos((-80.71-
(delta2+(0.5*K1m2*180/pi)))*pi/180)))*delta_t
L2m3 = (pi*fo/H3)*(P3-(1.115-3.954*cos((-73.32-
(delta3+(0.5*K1m3*180/pi)))*pi/180)))*delta_t
K3m1 = (omega1 + 0.5 *L2m1)*delta_t;
K3m2 = (omega2 + 0.5 *L2m2)*delta_t;
K3m3 = (omega3 + 0.5 *L2m3)*delta_t;
L3m1 = (pi*fo/H1)*(P1-(1.377-9.871*cos((26.58-
(delta1+(0.5*K2m1*180/pi)))*pi/180)))*delta_t;
L3m2 = (pi*fo/H2)*(P2-(0.051-4.974*cos((-80.71-
(delta2+(0.5*K2m2*180/pi)))*pi/180)))*delta_t;
L3m3 = (pi*fo/H3)*(P3-(1.115-3.954*cos((-73.32-
(delta3+(0.5*K2m3*180/pi)))*pi/180)))*delta_t;
K4m1 = (omega1 + L3m1)*delta_t
K4m2 = (omega2 + L3m2)*delta_t
K4m3 = (omega3 + L3m3)*delta_t
L4m1 = (pi*fo/H1)*(P1-(1.377-9.871*cos((26.58-(delta1+
(K3m1*180/pi)))*pi/180)))*delta_t
L4m2 = (pi*fo/H2)*(P2-(0.051-4.974*cos((-80.71-(delta2+
(K3m2*180/pi)))*pi/180)))*delta_t
L4m3 = (pi*fo/H3)*(P3-(1.115-3.954*cos((-73.32-(delta3+
(K3m3*180/pi)))*pi/180)))*delta_t
fprintf('!%7.2f !%10.3f !%10.3f !%9.3f !%10.3f !%10.3f
!%10.3f! \n',t,delta1,delta2,delta3,omega1,omega2,omega3)
delta1 = delta1+ ((1/6)*((K1m1 + 2*K2m1 + 2*K3m1 +
K4m1)*180/pi));
 delta2 = delta2+ ((1/6)*((K1m2 + 2*K2m3 + 2*K3m2 +
K4m2)*180/pi));
delta3 = delta3+ ((1/6)*((K1m3 + 2*K2m2 + 2*K3m3 +
K4m3)*180/pi));
plot(t,delta1,'x')
plot(t,delta2,'x')
plot(t,delta3,'x')
omega1=omega1 + ((1/6)*(L1m1 + 2*L2m1 + 2*L3m1 + L4m1));
omega2=omega2 + ((1/6)*(L1m2 + 2*L2m2 + 2*L3m2 + L4m2));
omega3=omega3 + ((1/6)*(L1m3 + 2*L2m3 + 2*L3m3 + L4m3));
end

K1m1 =
0.050534
K1m2 =
0.0082636
K1m3 =
0.0016623
L1m1 =
1.6992
L1m2 =
0.16965
L1m3 =
-0.080778
K2m1 =
0.05903
K2m2 =
0.0091118
K2m3 =
0.0012584
L2m1 =
1.6783
L2m2 =
0.16641
L2m3 =
-0.081399
K4m1 =
0.06728
K4m2 =
0.0099244
K4m3 =
0.00084983
L4m1 =
1.6487
L4m2 =
0.16254
L4m3 =
-0.081717
! 0.21 ! 51.826 ! 45.881 ! 25.785 ! 5.053 ! 0.826 !
0.166!
K1m1 =
0.06729
K1m2 =
0.0099255
K1m3 =
0.00084933
L1m1 =
1.6487
L1m2 =
0.16457
L1m3 =
-0.083675
K2m1 =
0.075533
K2m2 =
0.010748
K2m3 =
0.00043096
L2m1 =
1.6173
L2m2 =
0.16071
L2m3 =
-0.083992
K4m1 =
0.083423
K4m2 =
0.011529
K4m3 =
1.0976e-05
L4m1 =
1.5758
L4m2 =
0.15624
L4m3 =
-0.083995
! 0.21 ! 55.204 ! 46.253 ! 26.007 ! 6.729 ! 0.993 !
0.085!
K1m1 =
0.083432
K1m2 =
0.011531
K1m3 =
1.0457e-05
L1m1 =
1.5757
L1m2 =
0.1589
L1m3 =
-0.086566
K2m1 =
0.091311
K2m2 =
0.012325
K2m3 =
-0.00042237
L2m1 =
 1.5314
L2m2 =
0.15443
L2m3 =
-0.08657
K4m1 =
0.098704
K4m2 =
0.013072
K4m3 =
-0.00085362
L4m1 =
1.4756
L4m2 =
0.14939
L4m3 =
-0.08625
! 0.21 ! 59.525 ! 46.671 ! 26.228 ! 8.343 ! 1.153 !
0.001!
K1m1 =
0.098713
K1m2 =
0.013073
K1m3 =
-0.00085416
L1m1 =
1.4755
L1m2 =
0.15266
L1m3 =
-0.089423
K2m1 =
0.10609
K2m2 =
0.013836
K2m3 =
-0.0013013
L2m1 =
1.4161
L2m2 =
0.14763
L2m3 =
-0.089104
K4m1 =
0.11283
K4m2 =
0.014546
K4m3 =
-0.0017435
L4m1 =
1.3438
L4m2 =
0.14207
L4m3 =
-0.088453
! 0.22 ! 64.748 ! 47.132 ! 26.448 ! 9.871 ! 1.307 !
-0.085!
K1m1 =
0.11284
K1m2 =
0.014547
K1m3 =
-0.0017441
L1m1 =
1.3437
L1m2 =
0.14592
L1m3 =
-0.092218
K2m1 =
0.11956
K2m2 =
0.015277
K2m3 =
-0.0022052
L2m1 =
1.2672
L2m2 =
0.14036
L2m3 =
-0.091568
K4m1 =
0.12546
K4m2 =
0.015948
K4m3 =
-0.0026581
L4m1 =
1.1773
L4m2 =
0.13431
L4m3 =
-0.090576
! 0.22 ! 70.815 ! 47.635 ! 26.662 ! 11.284 ! 1.455 !
-0.174!
K1m1 =
0.12547
K1m2 =
0.015949
K1m3 =
-0.0026586
L1m1 =
1.1771
L1m2 =
0.13872
L1m3 =
-0.094921
K2m1 =
0.13136
K2m2 =
 0.016643
K2m3 =
-0.0031332
L2m1 =
1.0825
L2m2 =
0.13268
L2m3 =
-0.09393
K4m1 =
0.13625
K4m2 =
0.017273
K4m3 =
-0.0035962
L4m1 =
0.97494
L4m2 =
0.12618
L4m3 =
-0.092589
! 0.22 ! 77.650 ! 48.175 ! 26.870 ! 12.547 ! 1.595 !
-0.266!
K1m1 =
0.13626
K1m2 =
0.017274
K1m3 =
-0.0035968
L1m1 =
0.9747
L1m2 =
0.13113
L1m3 =
-0.097502
K2m1 =
0.14113
K2m2 =
0.01793
K2m3 =
-0.0040843
L2m1 =
0.86256
L2m2 =
0.12464
L2m3 =
-0.096162
K4m1 =
0.14484
K4m2 =
0.018518
K4m3 =
-0.0045566
L4m1 =
0.7389
L4m2 =
0.11773
L4m3 =
-0.094463
! 0.22 ! 85.157 ! 48.750 ! 27.068 ! 13.626 ! 1.727 !
-0.360!
K1m1 =
0.14485
K1m2 =
0.018519
K1m3 =
-0.0045572
L1m1 =
0.73858
L1m2 =
0.12319
L1m3 =
-0.099929
K2m1 =
0.14854
K2m2 =
0.019135
K2m3 =
-0.0050568
L2m1 =
0.611
L2m2 =
0.1163
L2m3 =
-0.098232
K4m1 =
0.15093
K4m2 =
0.01968
K4m3 =
-0.0055376
L4m1 =
0.47458
L4m2 =
0.10903
L4m3 =
-0.096168
! 0.22 ! 93.222 ! 49.356 ! 27.255 ! 14.485 ! 1.852 !
-0.456!
K1m1 =
0.15094
K1m2 =
0.019681
K1m3 =
-0.0055383
L1m1 =
0.4742
L1m2 =
0.11497
L1m3 =
 -0.10217
K2m1 =
0.15331
K2m2 =
0.020256
K2m3 =
-0.0060491
L2m1 =
0.33503
L2m2 =
0.10772
L2m3 =
-0.10011
K4m1 =
0.15426
K4m2 =
0.020756
K4m3 =
-0.0065375
L4m1 =
0.19074
L4m2 =
0.10013
L4m3 =
-0.097675
! 0.23 ! 101.708 ! 49.989 ! 27.428 ! 15.094 ! 1.968 !
-0.554!
K1m1 =
0.15427
K1m2 =
0.020757
K1m3 =
-0.0065381
L1m1 =
0.19031
L1m2 =
0.10653
L1m3 =
-0.1042
K2m1 =
0.15522
K2m2 =
0.021289
K2m3 =
-0.0070591
L2m1 =
0.044826
L2m2 =
0.098953
L2m3 =
-0.10177
K4m1 =
0.15471
K4m2 =
0.021744
K4m3 =
-0.0075538
L4m1 =
-0.10141
L4m2 =
0.091094
L4m3 =
-0.098955
! 0.23 ! 110.465 ! 50.645 ! 27.584 ! 15.427 ! 2.076 !
-0.654!
K1m1 =
0.15471
K1m2 =
0.021745
K1m3 =
-0.0075545
L1m1 =
-0.10186
L1m2 =
0.097911
L1m3 =
-0.10598
K2m1 =
0.15421
K2m2 =
0.022235
K2m3 =
-0.0080844
L2m1 =
-0.24755
L2m2 =
0.090063
L2m3 =
-0.10317
K4m1 =
0.15224
K4m2 =
0.022644
K4m3 =
-0.0085843
L4m1 =
-0.38954
L4m2 =
0.081973
L4m3 =
-0.099979
! 0.23 ! 119.331 ! 51.322 ! 27.721 ! 15.471 ! 2.175 !
-0.755!
K1m1 =
0.15225
K1m2 =
0.022645
K1m3 =
-0.0085849
L1m1 =
 -0.38999
L1m2 =
0.08918
L1m3 =
-0.10749
K2m1 =
0.1503
K2m2 =
0.023091
K2m3 =
-0.0091224
L2m1 =
-0.52978
L2m2 =
0.0811
L2m3 =
-0.1043
K4m1 =
0.14697
K4m2 =
0.023454
K4m3 =
-0.009626
L4m1 =
-0.66186
L4m2 =
0.07282
L4m3 =
-0.10072
! 0.23 ! 128.139 ! 52.016 ! 27.838 ! 15.225 ! 2.264 !
-0.858!
K1m1 =
0.14697
K1m2 =
0.023455
K1m3 =
-0.0096267
L1m1 =
-0.66228
L1m2 =
0.080388
L1m3 =
-0.1087
K2m1 =
0.14366
K2m2 =
0.023857
K2m3 =
-0.01017
L2m1 =
-0.79085
L2m2 =
0.072117
L2m3 =
-0.10512
K4m1 =
0.13909
K4m2 =
0.024175
K4m3 =
-0.010676
L4m1 =
-0.90865
L4m2 =
0.063689
L4m3 =
-0.10115
! 0.23 ! 136.723 ! 52.722 ! 27.931 ! 14.697 ! 2.345 !
-0.963!
K1m1 =
0.13909
K1m2 =
0.024175
K1m3 =
-0.010677
L1m1 =
-0.90901
L1m2 =
0.071584
L1m3 =
-0.10958
K2m1 =
0.13454
K2m2 =
0.024533
K2m3 =
-0.011224
L2m1 =
-1.0225
L2m2 =
0.063163
L2m3 =
-0.10561
K4m1 =
0.1289
K4m2 =
0.024806
K4m3 =
-0.011731
L4m1 =
-1.1232
L4m2 =
0.054627
L4m3 =
-0.10125
! 0.24 ! 144.930 ! 53.437 ! 27.999 ! 13.909 ! 2.418 !
-1.068!
K1m1 =
0.1289
K1m2 =
 0.024806
K1m3 =
-0.011731
L1m1 =
-1.1235
L1m2 =
0.062819
L1m3 =
-0.11011
K2m1 =
0.12328
K2m2 =
0.02512
K2m3 =
-0.012282
L2m1 =
-1.2197
L2m2 =
0.054287
L2m3 =
-0.10575
K4m1 =
0.11674
K4m2 =
0.025348
K4m3 =
-0.012787
L4m1 =
-1.3025
L4m2 =
0.045679
L4m3 =
-0.10099
! 0.24 ! 152.617 ! 54.159 ! 28.039 ! 12.890 ! 2.481 !
-1.173!
K1m1 =
0.11673
K1m2 =
0.025349
K1m3 =
-0.012788
L1m1 =
-1.3027
L1m2 =
0.054135
L1m3 =
-0.11026
K2m1 =
0.11022
K2m2 =
0.025619
K2m3 =
-0.013339
L2m1 =
-1.3811
L2m2 =
0.04553
L2m3 =
-0.10551
K4m1 =
0.10297
K4m2 =
0.025803
K4m3 =
-0.013841
L4m1 =
-1.4463
L4m2 =
0.036886
L4m3 =
-0.10036
! 0.24 ! 159.663 ! 54.882 ! 28.051 ! 11.673 ! 2.535 !
-1.279!
K1m1 =
0.10296
K1m2 =
0.025804
K1m3 =
-0.013841
L1m1 =
-1.4465
L1m2 =
0.045577
L1m3 =
-0.11001
K2m1 =
0.095727
K2m2 =
0.026032
K2m3 =
-0.014391
L2m1 =
-1.5079
L2m2 =
0.036933
L2m3 =
-0.10486
K4m1 =
0.087923
K4m2 =
0.026172
K4m3 =
-0.014888
L4m1 =
-1.5572
L4m2 =
0.028286
L4m3 =
-0.099323
! 0.24 ! 165.963 ! 55.604 ! 28.032 ! 10.296 ! 2.580 !
-1.384!
K1m1 =
0.087915
K1m2 =
0.026173
K1m3 =
-0.014889
L1m1 =
-1.5573
L1m2 =
0.037181
L1m3 =
-0.10934
K2m1 =
0.080128
K2m2 =
0.026359
K2m3 =
-0.015435
L2m1 =
-1.6037
L2m2 =
0.028531
L2m3 =
-0.10381
K4m1 =
0.071918
K4m2 =
0.026457
K4m3 =
-0.015925
L4m1 =
-1.6393
L4m2 =
0.019913
L4m3 =
-0.097877
! 0.24 ! 171.437 ! 56.322 ! 27.980 ! 8.791 ! 2.617 !
-1.489!
K1m1 =
0.071909
K1m2 =
0.026458
K1m3 =
-0.015925
L1m1 =
-1.6393
L1m2 =
0.028983
L1m3 =
-0.10824
K2m1 =
0.063713
K2m2 =
 0.026603
K2m3 =
-0.016466
L2m1 =
-1.6728
L2m2 =
0.020358
L2m3 =
-0.10231
K4m1 =
0.055218
K4m2 =
0.026661
K4m3 =
-0.016946
L4m1 =
-1.6971
L4m2 =
0.011797
L4m3 =
-0.096004
! 0.25 ! 176.019 ! 57.032 ! 27.895 ! 7.191 ! 2.646 !
-1.593!
K1m1 =
0.055209
K1m2 =
0.026661
K1m3 =
-0.016947
L1m1 =
-1.6972
L1m2 =
0.021013
L1m3 =
-0.10668
K2m1 =
0.046723
K2m2 =
0.026766
K2m3 =
-0.01748
L2m1 =
-1.7201
L2m2 =
0.012442
L2m3 =
-0.10038
K4m1 =
0.038042
K4m2 =
0.026786
K4m3 =
-0.017949
L4m1 =
-1.7352
L4m2 =
0.0039625
L4m3 =
-0.093691
! 0.25 ! 179.664 ! 57.732 ! 27.775 ! 5.521 ! 2.666 !
-1.695!
K1m1 =
0.038032
K1m2 =
0.026786
K1m3 =
-0.01795
L1m1 =
-1.7352
L1m2 =
0.013299
L1m3 =
-0.10467
K2m1 =
0.029356
K2m2 =
0.026852
K2m3 =
-0.018473
L2m1 =
-1.7497
L2m2 =
0.004806
L2m3 =
-0.097987
K4m1 =
0.020567
K4m2 =
0.026834
K4m3 =
-0.018927
L4m1 =
-1.7572
L4m2 =
-0.0035674
L4m3 =
-0.090931
! 0.25 ! 182.337 ! 58.419 ! 27.620 ! 3.803 ! 2.679 !
-1.795!

disp('==========================================================
======')

================================================================

text (0.23,157,'mesin1')
text (0.23 ,69, 'mesin2')
text (0.23,28,'mesin3')