jogano, Вам, конечно же, огромное-преогромное спасибо за помощь и проделанную работу!

Сообщение от jogano Откройте маткадовский файл

Наверное, это богохульство, но я не умею в маткад

Сообщение от jogano самое время перефутболить вашу тему в Матлаб

В целом, ведь не имеет значения ни код, ни способ, ни метод решения; на сам вопрос темы мне ответили ещё в 3 сообщении.
Решение же обсуждаемой нами моей проблемы состоит в правильном написании исчерпывающей системы уравнений аналитической геометрии, поэтому матлаб-то, логически, не причём. Код матлаба я привёл, чтобы Вы смогли увидеть какие уравнения и ограничения вида равенств я ставлю к системе, отправляющейся в решальщик для нахождения её корней, и, возможно, нашли бы ошибку, которая и ведёт к некоторой "несанкционированности" работы "чёрного ящика".

0

Вот вам картинки этого файла. Картинки нумерованы. Разбирайтесь.

Вложения

Photos.rar (1.13 Мб, 3 просмотров)

1

Сообщение от Ромуальд_7 радиус обеих окружностей принят равным 14 и является константой всегда;

Что изначальная задача? Выглядит как “забивать гвозди микроскопом”.

0

jogano, Большое спасибо!

Добавлено через 11 минут
Excalibur921,

Сообщение от Ромуальд_7 Решение же обсуждаемой нами моей проблемы состоит в правильном написании исчерпывающей системы уравнений аналитической геометрии

-
по изначальной задаче я просил помощи с нахождением правильного угла поворота, который решился вместе с указанием на ошибку в логике расчётов; таким образом, задача плавно перетекла в поиск этой ошибки, ведь уважаемые авторы двух сообщений с указанием на эту самую ошибку так и не изъявили желания любезно на неё указать. Конкретики от них не было, а мои эксперименты показали, что ошибки нет, но тогда почему система выдаёт 6 решений, а не 4?
Я не предлагал разбирать мои методы решения и создавать эксперименты, я интересуюсь, как же выглядит правильно написанная исчерпывающая СНлАУ для решения задачи касания (любой задаваемой) окружности с прямой, проходящей через точку (любую задаваемую), что является чистейшей задачей ангема.
Это более чем конкретная задача, и уж точно она не выглядит как

Сообщение от Excalibur921 “забивать гвозди микроскопом”

.

Добавлено через 8 минут
В целом, я вообще не хотел никого обременять созданием вариантов решений и разбирательством в моей программе; мне всего-то нужно было (да и сейчас нужно), чтобы 1 из 2 человек, что нашли в моей СНлАУ №2 ошибку чуть громче об этом сказали (в том смысле, что конкретики бы ещё добавить не помешало), и всё! тему можно будет считать закрытой ; и никаких доп вопросов и траты не принадлежащего мне времени.
Повторю - я не хочу тратить ничьё время, но если человек сказал, что есть ошибка, то пусть он не руководствуется логикой, что он самый умный, а я - сам догадайся он это как-то проаргументирует, поскольку я же мучаюсь, можно сказать, пытаясь найти ошибку, которой, может, и нет; полагаю, недельный срок достаточен для обдумывания, можно и вскрыть карты уже (прошу прощения заранее - накипело уже из-за этой молчанки).

0

Ромуальд_7, вы ещё мучаетесь, после того, что я вам сбросил на рисунках?! Ну... тогда, как говорят патологоанатомы, "я сделал всё, что мог".

Сообщение от Ромуальд_7 прошу прощения заранее - накипело уже из-за этой молчанки

А что вы хотите? Ваша задача (особенно вторая, с 4-мя окружностями) достаточно громоздкая, а вы находитесь в таком разделе форума, что другим вкладываться в это не очень-то нужно. У меня ушло на эти две версии задачи часов 8, не меньше. И я её решил. Но вы продолжаете чего-то ждать... Поэтому ждать вы будете длительно, скорее всего. Или вы меняете вашу программу, убирая решатель для нахождения центров окружностей, а делая это прямыми формулами (1-я страница маткадовского файла) и дальше используя не один решатель для поиска точек касания, а четыре, или будете "мучиться" и дальше.

0

jogano,

Сообщение от jogano А что вы хотите?

Вы это серьёзно? Снова этот вопрос?

Сообщение от jogano достаточно громоздкая

Ничего подобного - цикл на два шага в цикле на два шага и решальщик в нём (1 строчка кода, по сути), который работает прекрасно, но нужна правильная математика и всего-то.

Сообщение от jogano Или вы меняете вашу программу, ...

Правильно ли я Вас понял, что Вы считаете, что нет иных разумных методов, кроме предложенного Вами?
Ваш понятен и прост, но это не делает его единственным, верно? Возможно, я субъективен, но ведь действительно проще составить систему уравнений и всего-то её решить; "проще, чем что?" - спросите Вы, "проще, чем всё остальное" - отвечу я.
Разумеется, я могу создать ещё одну тему с конкретным вопросом, чтобы, видимо, соблюсти все формальности, и получить не иной метод решения, а исправленный существующий, но зачееееем?? Это просто лишнее.

0

Ромуальд_7, "Что вы хотите?" это был риторический вопрос...

Сообщение от Ромуальд_7 Правильно ли я Вас понял, что Вы считаете, что нет иных разумных методов, кроме предложенного Вами?

Я считаю, что второму решальщику нужно подавать на вход две окружности и должно быть известно направления их вращения, потому что применяются разные критерии отсева (конкретнее - знаки определителей, если вы открывали картинки), а не бросать в него скопом 4 центра окружностей и чтобы решальщик сам варился в этом. А для окружностей противоположного направления вращения применять решальщик только если расстояние между их центрами больше суммы их радиусов (28), иначе есть пересечение, а значит нет внутренних касательных. Т.е. внешнее условие if..., а внутри решальщик для пар окружностей противоположного направления вращения.
Комплексные решения можете не учитывать сразу.

1

jogano, Вы неправильно понимаете как работает мой алгоритм.
Распишу поэтапно:
1. по введённым координатам векторов ищутся центры четырёх окружностей - по 2 около каждого вектора (используется всё тот же решальщик, но, естественно, ему даётся другая геометрия);
2. с помощью цикла попарно берутся окружности у каждого из векторов (V1_1 с V2_1, V1_1 с V2_2, V1_2 с V2_1 и V1_2 с V2_2) и координаты их центров подаются в решальщик, иными словами,

Сообщение от jogano второму решальщику нужно подавать на вход две окружности

(ровно как Вы говорите оно и делается, всем скопом центры не передаются, это и не выходит сделать (я поначалу пробовал));
3. решальщик, используя геометрию, заданную ему в виде уравнений с участием этих двух центров, выдаёт точные координаты точек касания окружностей для всех прямых (если окружности являются пересекающимися то система даст комплексное решение - своеобразная автопроверка, выходит). Геометрия говорит, что их всегда 4, если окружности не пересекаются, 2 - если пересекаются и 0 - если совпадают;
4.0. из полученных координат выбираются и отбрасываются комплексные решения, как недостойные рассмотрения;
4. по полученным вектор-столбцам с координатами точек касания для каждой окружности (итого 8 пар координат) происходит отбор по признаку равенства знаков векторов векторного произведения радиусов окружностей (по очереди) и векторов, коллинеарных касательным, сначала по первой окружности ("правильная физика входа в касательную"), а затем - по второй ("правильная физика выхода из неё") для каждой пары окружностей в цикле их попарного перебора;
5. Удовлетворяющие всем условиям и прошедшие все этапы отбора решения сохраняются в массивы.

Всё работает великолепно, за исключением, единственно, того, что решений находится всегда 6, а не 4 в пункте 2, причём всегда чисто комплексных, и в некоторых случаях не находятся решения, которые на самом деле есть.
То есть, повторюсь, система принимает не совсем корректный набор уравнений геометрии процесса, вот и всё, и вот именно в этом я, по-прежнему, и прошу помощи.

0

Сообщение от Ромуальд_7 по изначальной задаче я просил помощи с нахождением правильного угла поворота,

Что изначальная задача? Для чего это вообще? Выглядит как игра преследования где точка поворачивает минимум по окружности радиус 14 преследуя цель. Это траектория 2д ракеты для игры?

0

Дабы закрыть тему окончательно отвечу сам себе (простите, что так долго отсутствовал - был в отпуске):
новообразовавшаяся проблема с СнЛАУ оказалось ложной - уравнения полностью верны (что, впрочем, я и ранее заявлял), но теперь это подтвердилось окончательно. Проблема исчезла, когда область расчётов сдвинулись от нулевых значений координат настолько, что даже значения десятков единиц x и y перестали фигурировать в расчётах. Таким образом выяснено, что алгоритм неустойчив, когда в расчётах принимают участие околонулевые координаты, во всех остальных же случаях он работает "идеально"; я думаю, что это связано с тем, что матлабовский "решальщик" СнЛАУ использует слишком большие (с точки зрения принимаемых на вход координат) интервалы поиска корней и регулярно через них "перешагивал", что и влекло за собой появление с последующим накоплением ошибки вычисления (но это всего лишь я так думаю).
В общем, изначальная и образовавшаяся проблема решена полностью, тему можно закрывать.
Всем помогавшим - огромное спасибо и удачи!!!

0