% Gerak Peluru Pada Bidang Datar
disp(['=================================================='])
disp(['PERSAMAAN GERAK PELURU'])
disp(['Y= Yo + Voy*t - 0.5*g*t^2'])
disp(['X=Xo+Vox*t'])
disp(['Keterangan:'])
disp([' Y merupakan ketinggian suatu peluru pada saat t sekon'])
disp([' Yo merupakan ketinggian suatu peluru pada saat t=0 sekon pada
posisi di titik nol'])
disp([' X merupakan jarak suatu peluru pada saat t sekon'])
disp([' Xo merupakan jarak suatu peluru pada saat t=0 sekon pada posisi di
titik nol'])
teta=input('Sudut awal peluru ditembakkan (dalam derajat) = ');
Vo=input('Kecepatan awal peluru ditembakkan (dalam m/s) = ');
sudut=(teta*pi())/180;
g=9.78 %m/s^2
Voy=Vo*sin(sudut);
Vox=Vo*cos(sudut);
t_akhir=(2*Voy)/g;
t=0:0.01:t_akhir
disp(['t akhir merupakan waktu pada saat ketinggian akhir nol meter adalah
',num2str(t_akhir),' sekon'])
%rumus ketinggian dan jarak peluru pada saat t sekon
Y=Voy*t-(0.5*g*t.^2);
X=+(Vox*t);
%untuk membuat plot untuk ketinggian terhadap jarak
subplot(3,1,1)
plot(X,Y,'b')
grid on
title('Grafik Y terhadap X')
h=get(gca,'title');
set(h,'fontsize',13)
xlabel('Jarak (m)')
ylabel('Ketinggian (m) ')
subplot(3,1,2)
plot(t,Y,'r')
grid on
title('Grafik Y terhadap t')
h=get(gca,'title');
set(h,'fontsize',13)
xlabel('Waktu (sekon)')
ylabel('Ketinggian (m)')
subplot(3,1,3)
plot(t,X,'g')
grid on
title('Grafik X terhadap t')
h=get(gca,'title');
set(h,'fontsize',13)
xlabel('Waktu (sekon)')
ylabel('Jarak (m)')
disp(['Terima kasih telah menggunakan program pemodelan ini'])
disp(['created by: yusia '])
Untuk mencobanya silakan copy script di atas, masukkan ke M-file. Berikut beberapa pengetesan dengan menggunakan Matlab untuk gerak parabola:

0 comments:
Post a Comment