Simulasi ini dapat dijalankan menggunakan ‘Octave’ maupun ‘Matlab’
Parameter yang dapat diubah adalah R, L, C, dan t-akhir simulasi.
hasil simulasi ditampilkan dalam bentuk grafik
***
clear all;
clc;
R=50;
L=0.005;
C=0.000001;
vs=5;
a=[0 1;-1/(L*C) -R/L]
b=[0;1]
x=[0;1];
t=0; dt=0.000001; takhir=0.004;
u=vs/(L*C);
tampil=[];
vekt=[];
vekx1=[];
vekx2=[];
kk=[];
while t<=takhir
k1=dt*(a*x+b*u);
k2=dt*(a*(x+0.5*k1)+b*u);
k3=dt*(a*(x+0.5*k2)+b*u);
k4=dt*(a*(x+k3)+b*u);
tampil=[tampil; t x(1) x(2)];
vekt=[vekt;t];
vekx1=[vekx1;x(1)];
vekx2=[vekx2;x(2)];
kk=[kk; k1 k2 k3 k4];
x=x+(k1+2*k2+2*k3+k4)/6;
t=t+dt;
end
%vekt
%vekx1
%vekx2
%tampil
%kk
plot (vekt,vekx1);
legend (‘vc, volt’,’dvc/dt’ );
grid on;
title (‘Simulasi Tegangan kapasitor pada rangkaian RLC seri Menggunakan Metode Runge-Kutta’);
xlabel (‘t, detik’,’Fontsize’,12);
ylabel (‘vc’,’Fontsize’,12);