MathCAD   КМ   ТПОИ   3GL   к экономической информатике   4GL - визуальным средам

Решение задачи Кеплера в пакете Mathcad

Поршнев С.В.

Проверка второго закона Кеплера

Для проверки второго закона Кеплера необходимо провести сравнение площади, заметаемой радиус-вектором за равные промежутки времени, используя значения кинематических характеристик движения тела в гравитационном поле, полученные численным решением системы ДУ (24), (25), возвращенные в матрицу Z. Напомним, что размерность матрицы Z, содержащей решения системы ДУ, -- N* 5; нумерация столбцов и строк матрицы начинается с нулевого значения; значения t, содержатся соответственно в 0-ом, 1-ом, 2-ом, 3-ем, 4-ом столбцах матрицы.

   Так как решение уравнений движения проводится на равномерной временной сетке, то для задания начального и конечного значений временного интервала, в течение которого вычисляется площадь сектора, заметаемого радиус-вектором, достаточно указать номер соответствующей строки матрицы Z и количество используемых точек (т.е. длину временного интервала).

При достаточно малом шаге интегрирования D t уравнений движения (определяемым количеством точек N, в которых ищется решение системы ДУ, ) площадь, заметаемая радиус вектором за время D t, примерно равна площади треугольника с вершинами в точках (0,0), , (рис. 8).

рис. 8

Для вычисления площади треугольника, у которого заданы координаты вершин, можно воспользоваться формулой Герона

, (28)

где p - полупериметр треугольника, изображенного на рис. 8, a, b, c - длины его сторон:

, (29)

, (30)

. (31)

Однако программа вычисления площади сектора, заметаемого радиус-вектором, получается более компактной, если использовать встроенную в пакет Mathcad функцию, вычисляющую векторное произведение трехмерных векторов. Как известно из аналитической геометрии, площадь треугольника, изображенного на рис. 8, равна

, (32)

поэтому алгоритм вычисления площади сектора, заметаемого радиус-вектором, реализуется следующей последовательностью действий:

                  1. Задание начального момента времени и длины временного интервала.

                  2. Присвоение координатам вектора соответствующих значений координат орбиты () в момент времени.

                 3. Присвоение координатам вектора соответствующих значений координат орбиты () в момент времени .

                 4. Вычисление площади треугольника , заметаемого радиус-вектором на временном интервале .

                 5. Нахождение площади, заметаемой радиус-вектором на временном интервале , суммированием известной площади, заметаемой радиус-вектором на интервале , и D S.

   Для реализации этого алгоритма документ, описанный в предыдущем разделе, необходимо дополнить следующими строками.

                  1) Задание номера начальной точки, начиная с которой производится вычисление площади, заметаемой радиус-вектором

;

                   2) Задание количества точек, по которым вычисляется площадь, заметаемая радиус-вектором

;

                   3) Программа для вычисления площади, реализующая описанный выше алгоритм

   Далее, задавая различные начальные значения номера начальной точки, можно убедиться в том, что значение переменной S с определенной точностью остается.

                                                                                   Рис.9

 

   Представляет определенный интерес не только вычислить площадь сектора, заметаемого радиус-вектором, но и изобразить данный сектор на рисунке вместе с траекторией тела. Однако анализ этой задачи обнаруживает известную ограниченность графических возможностей Mathcad, который позволяет изображать на рисунках только функциональные зависимости. Средств, позволяющих нарисовать замкнутую фигуру произвольной формы и залить ее выбранным цветом, разработчики Mathcad, к сожалению, не предусмотрели. Однако выход из данной ситуации все-таки существует. Он состоит в том, чтобы последовательно нарисовать каждый из секторов, заметаемых радиус-вектором за время Dt. Для построения последовательности секторов необходимо сформировать двумерный массив S1, содержащий координаты вершин последовательных треугольников, аппроксимирующих сектор, заметаемый радиус вектором (рис. 9), и затем построить график зависимости, заданной в виде таблицы. Ниже мы приводим текст соответствующей программы, который следует добавить в ранее созданный документ.

   Затем построить на одном чертеже орбиту, отметить фокус эллипса, и нарисовать сектор, заметаемый радиус-вектором (рис. 10).

                                                                                           Рис.10

на верх   MathCAD   КМ   ТПОИ   3GL   к экономической информатике   4GL - визуальным средам

Знаете ли Вы, что такое мысленный эксперимент, gedanken experiment?
Это несуществующая практика, потусторонний опыт, воображение того, чего нет на самом деле. Мысленные эксперименты подобны снам наяву. Они рождают чудовищ. В отличие от физического эксперимента, который является опытной проверкой гипотез, "мысленный эксперимент" фокуснически подменяет экспериментальную проверку желаемыми, не проверенными на практике выводами, манипулируя логикообразными построениями, реально нарушающими саму логику путем использования недоказанных посылок в качестве доказанных, то есть путем подмены. Таким образом, основной задачей заявителей "мысленных экспериментов" является обман слушателя или читателя путем замены настоящего физического эксперимента его "куклой" - фиктивными рассуждениями под честное слово без самой физической проверки.
Заполнение физики воображаемыми, "мысленными экспериментами" привело к возникновению абсурдной сюрреалистической, спутанно-запутанной картины мира. Настоящий исследователь должен отличать такие "фантики" от настоящих ценностей.

Релятивисты и позитивисты утверждают, что "мысленный эксперимент" весьма полезный интрумент для проверки теорий (также возникающих в нашем уме) на непротиворечивость. В этом они обманывают людей, так как любая проверка может осуществляться только независимым от объекта проверки источником. Сам заявитель гипотезы не может быть проверкой своего же заявления, так как причина самого этого заявления есть отсутствие видимых для заявителя противоречий в заявлении.

Это мы видим на примере СТО и ОТО, превратившихся в своеобразный вид религии, управляющей наукой и общественным мнением. Никакое количество фактов, противоречащих им, не может преодолеть формулу Эйнштейна: "Если факт не соответствует теории - измените факт" (В другом варианте " - Факт не соответствует теории? - Тем хуже для факта").

Максимально, на что может претендовать "мысленный эксперимент" - это только на внутреннюю непротиворечивость гипотезы в рамках собственной, часто отнюдь не истинной логики заявителя. Соответсвие практике это не проверяет. Настоящая проверка может состояться только в действительном физическом эксперименте.

Эксперимент на то и эксперимент, что он есть не изощрение мысли, а проверка мысли. Непротиворечивая внутри себя мысль не может сама себя проверить. Это доказано Куртом Гёделем.

Понятие "мысленный эксперимент" придумано специально спекулянтами - релятивистами для шулерской подмены реальной проверки мысли на практике (эксперимента) своим "честным словом". Подробнее читайте в FAQ по эфирной физике.

НОВОСТИ ФОРУМА

Форум Рыцари теории эфира


Рыцари теории эфира
 10.11.2021 - 12:37: ПЕРСОНАЛИИ - Personalias -> WHO IS WHO - КТО ЕСТЬ КТО - Карим_Хайдаров.
10.11.2021 - 12:36: СОВЕСТЬ - Conscience -> РАСЧЕЛОВЕЧИВАНИЕ ЧЕЛОВЕКА. КОМУ ЭТО НАДО? - Карим_Хайдаров.
10.11.2021 - 12:36: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от д.м.н. Александра Алексеевича Редько - Карим_Хайдаров.
10.11.2021 - 12:35: ЭКОЛОГИЯ - Ecology -> Биологическая безопасность населения - Карим_Хайдаров.
10.11.2021 - 12:34: ВОЙНА, ПОЛИТИКА И НАУКА - War, Politics and Science -> Проблема государственного терроризма - Карим_Хайдаров.
10.11.2021 - 12:34: ВОЙНА, ПОЛИТИКА И НАУКА - War, Politics and Science -> ПРАВОСУДИЯ.НЕТ - Карим_Хайдаров.
10.11.2021 - 12:34: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Вадима Глогера, США - Карим_Хайдаров.
10.11.2021 - 09:18: НОВЫЕ ТЕХНОЛОГИИ - New Technologies -> Волновая генетика Петра Гаряева, 5G-контроль и управление - Карим_Хайдаров.
10.11.2021 - 09:18: ЭКОЛОГИЯ - Ecology -> ЭКОЛОГИЯ ДЛЯ ВСЕХ - Карим_Хайдаров.
10.11.2021 - 09:16: ЭКОЛОГИЯ - Ecology -> ПРОБЛЕМЫ МЕДИЦИНЫ - Карим_Хайдаров.
10.11.2021 - 09:15: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Екатерины Коваленко - Карим_Хайдаров.
10.11.2021 - 09:13: ВОСПИТАНИЕ, ПРОСВЕЩЕНИЕ, ОБРАЗОВАНИЕ - Upbringing, Inlightening, Education -> Просвещение от Вильгельма Варкентина - Карим_Хайдаров.
Bourabai Research - Технологии XXI века Bourabai Research Institution