Saturday, July 6, 2013

Plot gerak jatuh bebas di MATLAB

Buat teman-teman sekalian yang ingin belajar plot gerak jatuh bebas di MATLAB berikut saya berikan sedikit tutorial. Semoga bermanfaat:




% gerak jatuh bebas; author moh. fajar
fig = figure('menubar', 'none');
clc;
clear all;
v0 = 48; % kecepatan awal (m/s)
angle = 45; % sudut (derajat)
g = 9.8; % gravitasi (m/s)
t = linspace(0,10,60); % time step (0 sampai 10 detik)
x = v0.*t.*cosd(angle); % persamaan lintasan untuk x (horizontal)
y = v0.*t.*sind(angle) - 0.5*g.*t.^2; % persamaan lintasan untuk y (vertikal)
plot(x,y); % buat lintasan
hold on;
yy = 0;
plot(x,yy, 'b.', 'markersize', 10); % buat garis permukaan tanah
hold on;
p = plot(0,0);
for i=1:length(t),
if y(i) < 0, % jika nilai y lebih kecil dari nilai permukaan tanah, hentikan looping
disp('jarak maksimum:');
disp(x(i));
disp('tinggi maksimum:');
disp(max(y));
break;
end
if ishandle(p)==1, delete(p);end % jika obkek plot masih ada, maka hapus biar tidak ter-double
p = plot(x(i),y(i),'r.','markersize', 23); % buat objek plot
drawnow; % langsung gambar objek plot
pause(0.2); % pause pergantian gambar tiap 0.2 detik
end
disp('END');
view raw proyektil.m hosted with ❤ by GitHub