Исследование полета снаряда

m = 6.3; % масса снаряда

cal = 76; % калибр (в мм)

S = pi * ( cal / 1000 ) ^ 2 % площадь поперечного сечения снаряда (мидель)

Cx = 0.2; % аэродинамический(баллистический) коэффициент

Pv = 1.225; % плотность воздуха (на уровне моря)

NU = [0; 0]; % координаты в начальный момент времени

t = 100; % время, в течение которого ведется расчет

colour = ['r' 'g' 'b' 'c' 'm'] %

цвета графиков

kV0 = [340 580 680 780 1360] % начальная скорость

ka = [25 35 45 55 65] % угол бросания (возвышения орудия)

kCx = [0.1 0.15 0.2 0.25 0.3] % аэродинамический(баллистический) коэффициент

km = [3.15 5.3 6.3 7.3 12.6] % масса снаряда

% Зависимость от Vo

figure

hold

for i = 1:5

V0 = kV0(i)

[T2,Y2]=ode45(@Complex, [0 t], NU);

plot(Y2(:, 1), Y2(:, 2), colour(i))

end

title 'Зависимость от Vo'

xlabel 'Дальность'

ylabel 'Высота'

legend 'Vo = 340' 'Vo = 580' 'Vo = 680' 'Vo = 780' 'Vo = 1360'

axis([0 10000 0 10000])

grid

V0 = kV0(3)

% Зависимость от a

figure

hold

for i = 1:5

a = ka(i)

[T2,Y2]=ode45(@Complex, [0 t], NU);

plot(Y2(:, 1), Y2(:, 2), colour(i))

end

title 'Зависимость от a'

xlabel 'Дальность'

ylabel 'Высота'

legend 'a = 25' 'a = 35' 'a = 45' 'a = 55' 'a = 65'

axis([0 10000 0 10000])

grid

a = ka(3)

% Зависимость от Cx

figure

hold

for i = 1:5

Cx = kCx(i)

[T2,Y2]=ode45(@Complex, [0 t], NU);

plot(Y2(:, 1), Y2(:, 2), colour(i))

end

title 'Зависимость от Cx'

xlabel 'Дальность'

ylabel 'Высота'

legend 'Cx = 0.1' 'Cx = 0.15' 'Cx = 0.2' 'Cx = 0.25' 'Cx = 0.3'

axis([0 10000 0 10000])

grid

Cx = kCx(3)

% Зависимость от m

figure

hold

for i = 1:5

m = km(i)

[T2,Y2]=ode45(@Complex, [0 t], NU);

plot(Y2(:, 1), Y2(:, 2), colour(i))

end

title 'Зависимость от m'

xlabel 'Дальность'

ylabel 'Высота'

legend 'm = 3.15' 'm = 5.3' ' m = 6.3' 'm = 7.3' 'm = 12.6'

axis([0 10000 0 10000])

grid

function Complex = Complex(t, x)

global Cx S Pv m a V0 g

Complex = [ ( 2*m*V0*cosd(a) ) / ( Cx*S*Pv*V0*cosd(a).*t + 2*m ) ; ( 2*m*V0*sind(a) ) / ( Cx*S*Pv*V0*sind(a).*t + 2*m ) - g.*t ];

Размещено на Allbest.ru

 Скачать реферат