201
y3 = 0; y4 = vo*sind(alpha(i));
y0 = [y1 y2 y3 y4]'; % vector chua dieu kien dau
[t,y] = ode45('dandao',[0: 0.05: t_end],y0);
ymax(i) = max(y(:,3));
% tim thoi gian chuyen dong
for k = 10:length(t)
if y(k,3)>=0 && y(k+1,3)<0
t_ground(i) = t(k); % thoi gian tiep dat
ki = k+1;
end
end
tamxa(i) = y(ki, 1); % Xac dinh tam xa va tam cao
tamcao(i) = max(y(1:ki, 3));
plot(y(1:ki,1),y(1:ki,3),'k-', 'linewidth',1.5),
xlabel('x [m]'), ylabel('y [m]'), grid on
hold on
end
axis equal
title('Qui dao chuyen dong khi bo qua luc can khong khi')
disp(['goc ban (do) ', 'Thoi gian bay (s)' , ' tam xa (m)', '
tam cao (m)'])
for i = 1:length(alpha)
disp([sprintf('%4d',alpha(i)), sprintf(' %6.3f',t_ground(i)),
sprintf(' %6.3f',tamxa(i)), sprintf(' %6.3f',tamcao(i))])
end
0
1000
2000
3000
4000
5000
6000
7000
8000
9000
0
500
1000
1500
2000
2500
3000
3500
4000
4500
x [m]
y
[
m
]
Qui dao chuyen dong khi bo qua luc can khong khi
Hình 9-14. Quĩ ñạo chuyển ñộng của viên ñạn khi không có cản
với các góc bắn khác nhau 10, 20, . … , 90
o
