Помогите с кодом. Как видоизменить или что добавить к формуле, что бы величины менялись от высоты или времени в каждом цикле ? while y(i)>=0 % цикл продолжается, пока высота полета материальной точки положительна или равно 0 (м.т. приземлилась); Vx(i+1) = Vx(i)+ax(i)*dt; % обновление скорости вдоль оси x по дискретному аналогу первого порядка точности; Vy(i+1) = Vy(i)+ay(i)*dt; % обновление скорости вдоль оси y по дискретному аналогу первого порядка точности; x(i+1) = x(i)+ Vx(i)*dt; % обновление координаты x по дискретному аналогу первого порядка точности; y(i+1) = y(i)+Vy(i)*dt; % обновление координаты y по дискретному аналогу первого порядка точности; t(i+1) = t(i) + dt; % временная ось; X(i+1) = (ro*(Vx(i+1)^2+ Vy(i+1)^2)*S*Cx)/2; % обновление силы лобового сопротивления; A(i+1) = atan(Vy(i+1)/Vx(i+1))*180/pi; % угол наклона вектора начальной скорости, градусы; ax(i+1) = -X(i+1)* cosd(A(i+1))/m; % обновление ускорения вдоль оси x; ay(i+1) = -X(i+1)* sind(A(i+1))/m-g; % обновление ускорения вдоль оси y; i = i+1; % обновление счетчика итераций; end if Lmax<x(i) % условие, если текущая дальность полета больше максимальной дальности; Lmax=x(i); % обновление максимальной дальности; Aopt=A(1)*180/pi; % присваивание оптимального угла наклона вектора начальной скорости; end end
Любые из этих величин должны меняться либо по времени, либо по высоте
Centurio
Просветленный
(33562)
Марат Суспенский, нужна формула зависимости, не от фонаря же зависимость выдумывать. Где вы взяли эту задачу, как она полностью выглядит?
while y(i)>=0 % цикл продолжается, пока высота полета материальной точки положительна или равно 0 (м.т. приземлилась);
Vx(i+1) = Vx(i)+ax(i)*dt; % обновление скорости вдоль оси x по дискретному аналогу первого порядка точности;
Vy(i+1) = Vy(i)+ay(i)*dt; % обновление скорости вдоль оси y по дискретному аналогу первого порядка точности;
x(i+1) = x(i)+ Vx(i)*dt; % обновление координаты x по дискретному аналогу первого порядка точности;
y(i+1) = y(i)+Vy(i)*dt; % обновление координаты y по дискретному аналогу первого порядка точности;
t(i+1) = t(i) + dt; % временная ось;
X(i+1) = (ro*(Vx(i+1)^2+ Vy(i+1)^2)*S*Cx)/2; % обновление силы лобового сопротивления;
A(i+1) = atan(Vy(i+1)/Vx(i+1))*180/pi; % угол наклона вектора начальной скорости, градусы;
ax(i+1) = -X(i+1)* cosd(A(i+1))/m; % обновление ускорения вдоль оси x;
ay(i+1) = -X(i+1)* sind(A(i+1))/m-g; % обновление ускорения вдоль оси y;
i = i+1; % обновление счетчика итераций;
end
if Lmax<x(i) % условие, если текущая дальность полета больше максимальной дальности;
Lmax=x(i); % обновление максимальной дальности;
Aopt=A(1)*180/pi; % присваивание оптимального угла наклона вектора начальной скорости;
end
end