@Fasterovich

Всем доброго времени суток. Запустил первый раз матлаб, чтобы выполнить задание. Задание следующее:

Без учета сил сопротивления уравнения движения тела, брошенного под углом L к горизонту с начальной скоростью V , задаются формулами:

x=(V*cos(L))*t;
y=((V*sin(L))*t)-(10.*(t.^2))/2;


Написать программу 1 расчета траектории движения тела x=F(y) при различных значениях угла L и начальной скорости V . Ввод угла (в градусах) и начальной скорости осуществлять с помощью команды input. Для построения графика использовать функцию plot(x,y).

Вроде все элементарно (для проверки не ввожу числа с клавиатуры):

Проблема в том, что график из себя представляет нечто странное, он должен быть положительным, от нуля и в положительную сторону по обеим осям, но этого нет. Я не понимаю как правильно задать время (t -я его задал в виде массива значений, может совершенно не верно), ведь оно растет от нуля к бесконечности, ну или хотя бы до тех пор, пока y не станет нулем. Запутался я после первого запуска.(

График в приложении.

Миниатюры

 

0

@Heidegger

Угол – в радианах (L0 = pi/3), временной промежуток – поменьше (t = 0:0.1:1.5), абсциссу – вперёд (plot(x, y)).

2

@Fasterovich

Очень благодарен. Теперь мой график и вправду похож на физический процесс. Просто меня смутило описание задания функции x=F(y), поэтому наоборот и строил график)

0

@Heidegger

Да, странная формулировка – для параболы в общем случае нельзя построить обратную функцию, поэтому для одного и того же значения y возможны 2 различных значения x. Возможно, просто опечатка.

1

@Fasterovich

У меня еще один вопрос:
Нашел максимум высоты полета:

ymax=max(y);

но никак не могу с помощью команды disp(ymax) вывести на экран максимальную высоту полета (но число выдает при отладке верное), ругается компилятор, а если вычислить длину полета xmax=max(x) - то выдает какое-то число в экспонициальной форме(.

0

@Heidegger

А что пишет компилятор? Лучше привести полный код.

Тег [MATLAB][/MATLAB] вверху справа под значком динамика.

1

@Fasterovich

Вот что выдает:

Subscript indices must either be real positive integers or logicals.

Error in First (line 17)
disp(ymax);

Видимо, что-то не нравится по типу

График теперь выглядит так:

Миниатюры

 

0

@Heidegger

Это точно весь текст программы? Здесь нет строки 17.

Время (t = 0:0.1:100) задано слишком большое – за это время снаряд уйдёт глубоко под землю и пророет длинный тоннель (отсюда и огромные значения xmax).

Параметры параболы лучше искать по явным формулам. Или применять более изощрённые методы.

В начало программы можно добавить clear all.

1

@Fasterovich

Строка 17 - это строка 13, просто я не стал вставлять верхние закомментированные строки

Добавлено через 8 минут
А есть ли возможность вычислить при y=0? Мне как раз и нужно это значение.

0

@Heidegger

Где-то (в тексте программы или в командном окне), видимо, проскочило присваивание вида (disp = ...) – и теперь Matlab принимает disp за переменную. Инструкция clear all (или просто clear) очистит память от переменных, и disp снова станет функцией вывода на экран.

Можно так:

Ищем последнее положительное значение y и выводим соответствующее значение x.

0

@Fasterovich

Если я правильно понимаю, то это будет предыдущая точка перед пересечением с абсциссой, т.к. в данном случае должен быть равен нулю. но он равен нулю и в точке 0:0

0

@Heidegger

Да.
Он может быть ни в какой больше точке не равен нулю – вычисления ведь идут с определённым шагом.

Можно уменьшить шаг по времени до необходимой точности (t = 0:0.01:2) или вычислять среднее:

1

@Fasterovich

Еще есть вопрос. В задании сказано, что нужно построить график Xmax(L0), где Xmax -длина полета, а L- это начальный угол броска. Это же две постоянные, каким образом можно построить график? Наверное, я не совсем верно понимаю задание( Ни у кого нет мыслей?

0

@Heidegger

Максимальная длина полёта Xmax меняется в зависимости от начального угла броска L, так что график вполне можно построить. Функцию можно задать либо явно, выразив из уравнений движения, либо на основе анализа параболы, как это было сделано выше.

1

@Fasterovich

То есть, мне нужно взять диапазон значений начального угла и находить максимальное значение длины полета в зависимости от этих значений угла?
Добавлено через 27 минут
Он так не считает, я не понимаю как сделать так, чтобы он мог считать это максимальное значение Xm, если у нас есть две не постоянные L1 и t1.

0

@Heidegger

Да, идея правильная, но реализацию нужно продумать – например, заполнять массив в цикле. И максимальное расстояние лучше искать в отдельной функции.

Но самый разумный вариант – рассчитывать по готовой формуле (если только заданием не предусмотрено непосредственное исследование параболы).

0

@Fasterovich

А как я могу получить готовую формулу, если у меня xmax вычисляется в самом конце и отдельно

0

@Heidegger

Я выше приводил ссылку на явные формулы:

Можно посчитать обоими методами и сравнить.

Если считать программно, нужно для каждого значения L0 рассчитать xmax и результаты поместить в матрицу.

0

@Том Ардер

Вся задача - яркая иллюстрация тезиса: компьютер - зло для неграмотных, поскольку внушает ложную уверенность в том, что теперь можно не думать головой.
Решение и полное исследование (задачка известна уже много веков) есть в любом школьном учебнике, стоило бы сначала туда посмотреть, прежде чем вычислять (с компьютером) бессмыслицу.
Практический совет: прежде, чем включить компьютер, нужно сделать расчет вручную (одну, две точки траектории, калькулятора достаточно). Многое стало бы понятным. Напр., необходимый интервал времени.

1

@Fasterovich

Я правильно понимаю, график Xmax(L0) будет прямая, параллельная оси абсцисс?

Добавлено через 7 минут
Ох нет, это бред

Добавлено через 43 минуты
Вопрос по поводу инициализации массива.

В цикле хочу заполнить массивы L2 и массив xmax1. Не получается, выдает ошибку.

Undefined function or variable 'L2'.

Error in Untitled4 (line 39)
plot(L2, xmax1);

0