задача на экстремум

Тема в разделе "Университет", создана пользователем Львица, 19 фев 2006.

Статус темы:
Закрыта.
  1. TopicStarter Overlay

    Львица Учаcтник

    • Участник
    Рег.:
    19.02.2006
    Сообщения:
    24
    Симпатии:
    0
    Репутация:
    0
    Оффлайн
    привет всем!как найти минимальное расстояние между двумя параболами?у меня все сводится либо к нахождению экстремума через решение кубического уравнения,что плохо,либо к нахождению минимума суммы квадратов линейной и квадратичной формы,что,кажется,более перспективно в общем виде.помоготе!срочно!
  2. Ivrin Модератор

    • Команда форума
    Рег.:
    07.02.2006
    Сообщения:
    387
    Симпатии:
    4
    Репутация:
    0
    Оффлайн
    не в бровь, а в глаз! :D
    это лихо. периодические числа нервно курят в сторонке...
  3. Grigoriy Старожил

    • Участник
    • Старожил
    Рег.:
    11.02.2006
    Сообщения:
    4.120
    Симпатии:
    87
    Репутация:
    5
    Оффлайн
    Общая схема:
    Пусть есть 2 кривые, заданные ур-ями
    Ф1(х,у) = 0
    Ф2(х,у) = 0
    Найти минимум расстояний между ними.
    Если они пересекаются, то очевидно ето точки пересечения и расстояние там 0.
    Если нет, то отрезок, соединяющий точки, на которых достигается минимум, перпендикулярен к соответствующим кривым в етих точках.
    Т е если обозначить искомые точки А1 =(х1,у1), А2 = (х2,у2), то вектор
    (х1-х2, у1-у2) параллелен градиентам Ф1 и Ф2 в точках А1 и А2 соответственно( и, соответственно, ети градиенты параллельны между собой).
    В нашем конкретном случае имеем 2 квадратичные функции
    Ф1 = а1*х**2 + ... + ф1
    Ф2 = ...
    Просматриваются 2 пути решения задачи:
    1.
    Для любого данного вектора направления находим точки на параболах, где касательная совпадает с этим векторам, и смотрим, перпендикулярна ли прямая, соединяющая ети точки данному вектору.
    2. Для любой данной точки на, скажем, первой параболе, проводим прямую параллельно градиенту до пересечения со 2-ой, и смотрим, параллелен ли ли градиент 2-ой функции в получившейся точке исходноу.
    Понятно. что способы д б еквивалентны, но какой удобнее - надо считать.
    Также заметим, что я не пользовался тем, что кривые - именно параболы. Но Вы и не сказали, что они и расположены канонически, т е что уравнение имеет вид
    у = ах**2 + 2вх + с
    а тогда и непонятно, как использовать, что это параболы
  4. atoku Модератор

    • Команда форума
    Рег.:
    05.02.2006
    Сообщения:
    2.949
    Симпатии:
    9
    Репутация:
    0
    Адрес:
    USA
    Оффлайн
    Все проще должно быть. Надо таки использовать что это параболы и не надо использовать мат анализ. Это задача из аналитической геометрии.

    Я правда ее не решил. Но уверен в этом.
  5. TopicStarter Overlay

    Львица Учаcтник

    • Участник
    Рег.:
    19.02.2006
    Сообщения:
    24
    Симпатии:
    0
    Репутация:
    0
    Оффлайн
    То,что пишет Григорий,я делала,но там вылезает уравнение третьей степени,которое не раздраконишь в общем виде,везде гнусно.а параболы канонические.мне тоже,как Атоку,кажется,что все проще.но как?!
  6. Guest

    Рег.:
    Сообщения:
    0
    Симпатии:
    0
    Репутация:
    0
    Оффлайн
    У меня проще не получилось. За исключением случая, когда они пересекаются :)
    А что страшного с решением кубического уравнения? Есть же формулы для вычисления корней. Тем более, что если параболу представить в виде: y=a*(x-b)^2+c то в кубическом уравнении вторая степень уже будет исключена. Едиственная тонкость, из трех корней два могут ложными, соотвествовать либо локальному максимуму, либо минимуму, но тоже локальному.
  7. Муркенштейн Гастролёр

    • Участник
    Рег.:
    20.02.2006
    Сообщения:
    1.794
    Симпатии:
    15
    Репутация:
    2
    Адрес:
    Nowhere
    Оффлайн
    Зачастую эти формулы представляют корни ну в совсем уж жутком виде, в то время, как на самом деле они могут быть даже натуральными числами. И самое неприятное здесь - невозможность конструктивного приведения их из "трёхэтажного" состояния в более естественное. Остаются только искуственные методы, но сами понимаете, они далеко не всегда работают... :(

    В качестве яркого примера вышенаписанного приведу такое уравнение: 8x^3 - 4x^2 - 4x + 1 = 0.
  8. Grigoriy Старожил

    • Участник
    • Старожил
    Рег.:
    11.02.2006
    Сообщения:
    4.120
    Симпатии:
    87
    Репутация:
    5
    Оффлайн
    К сожалению, 3-я степень видимо по существу, и от неё не убежать :( Мне-то казалось, что там квадратное уравнение, но это очевидно не так. Возьмём стандартную параболу у = х**2 и немного её приподымем. Сразу видно, что есть минимумы - видимо по одному на ветвях и локальный экстремум (но не минимум)в вершинах. Но уравнение то должно давать и его!
    В общем же случае для квадратичной формы, как показывает пример эллипса, должно получаться уравнение 4-ой степени.
  9. TopicStarter Overlay

    Львица Учаcтник

    • Участник
    Рег.:
    19.02.2006
    Сообщения:
    24
    Симпатии:
    0
    Репутация:
    0
    Оффлайн
    Третья степень по существу вылезает,если решать через нормаль или шарахать матаном по функции,формула Кордано-Тартальи только в теории хороша,но если обойтись без поиска экстремума,а найти способ отыскать сразу наименьшее значение,неравенством Коши например(у меня не получается) или как нибудь еще.?Ладно,выкладываю конкретную задачу: найти расстояние между параболами у=х^2-2х и у=-х^2+18х-41(МИФИ).Но даже если здесь удасться отыскать корни в кубическом уравнении,которое меня не вдохнавляет на подвиги,то все равно интересно,что делать в общем случае?
  10. atoku Модератор

    • Команда форума
    Рег.:
    05.02.2006
    Сообщения:
    2.949
    Симпатии:
    9
    Репутация:
    0
    Адрес:
    USA
    Оффлайн
    Я не хочу ввязываться в решение: но у меня возник по этому поводу вопрос. А является ли кратчайшее расстояние каким-либо образом инвариантным относительно проективных преобразований? Я бы сразу захотел из одной параболы круг сделать.
  11. TopicStarter Overlay

    Львица Учаcтник

    • Участник
    Рег.:
    19.02.2006
    Сообщения:
    24
    Симпатии:
    0
    Репутация:
    0
    Оффлайн
    Вообще говоря,расстояние в проективном преобразовании,конечно,не инвариант,но если точки,дающие минимум оставить на месте(кто бы знал,где они),то расстояние сохранится,если остальные точки ближе не приползут
  12. Guest

    Рег.:
    Сообщения:
    0
    Симпатии:
    0
    Репутация:
    0
    Оффлайн
    Так они вроде пересекаются, поэтому ответ очевиден - 0.
  13. PP Заблокирован

    • Заблокирован
    Рег.:
    21.02.2006
    Сообщения:
    1.023
    Симпатии:
    148
    Репутация:
    5
    Нарушения:
    31
    Оффлайн
    Не хочу вас огорчать, но они пересекаются в точках х=5 +/- 3/sqrt(2)
  14. TopicStarter Overlay

    Львица Учаcтник

    • Участник
    Рег.:
    19.02.2006
    Сообщения:
    24
    Симпатии:
    0
    Репутация:
    0
    Оффлайн
    прошу прощения,ошибка при выкладывании,вторая функция у=-2х^2+18х-41
  15. atoku Модератор

    • Команда форума
    Рег.:
    05.02.2006
    Сообщения:
    2.949
    Симпатии:
    9
    Репутация:
    0
    Адрес:
    USA
    Оффлайн
    Так они обе соосные... Фи. Надо бы было лучше для любых парабол.
  16. Guest

    Рег.:
    Сообщения:
    0
    Симпатии:
    0
    Репутация:
    0
    Оффлайн
    Я не перепроверял, но там кубическое уравнение очень простое, типа (x-1)^3 = 7/6.
  17. TopicStarter Overlay

    Львица Учаcтник

    • Участник
    Рег.:
    19.02.2006
    Сообщения:
    24
    Симпатии:
    0
    Репутация:
    0
    Оффлайн
    Нет,кубическое урвнение гораздо хуже,но если толком нарисовать картинку,то видна нормаль к графикам,она проходит через (2;0) на первой и (4;-1) на второй.Убиться веником.Но что делать в общем случае?Похоже,что ничего.Спасибо всем.:)
  18. Ivrin Модератор

    • Команда форума
    Рег.:
    07.02.2006
    Сообщения:
    387
    Симпатии:
    4
    Репутация:
    0
    Оффлайн
    я могу спросить Шихановича... =)
  19. TopicStarter Overlay

    Львица Учаcтник

    • Участник
    Рег.:
    19.02.2006
    Сообщения:
    24
    Симпатии:
    0
    Репутация:
    0
    Оффлайн
    Хорошая идея про Шихановича,пусть тоже на форуме тусуется:)
  20. Ivrin Модератор

    • Команда форума
    Рег.:
    07.02.2006
    Сообщения:
    387
    Симпатии:
    4
    Репутация:
    0
    Оффлайн
    не, ты что! я тут кину какую-нибудь его задачку, а он заметит :o
    к тому же, у него вроде нет инета...)
  21. Guest

    Рег.:
    Сообщения:
    0
    Симпатии:
    0
    Репутация:
    0
    Оффлайн
    OK, пока дочка играет в бадминтон, у меня есть время все перепроверить и оформить решение.
    Итак у нас две параболы:
    y1 = f(x1) = (x1-1)^2 - 1
    y2 = g(x2) = -2*(x2-4.5)^2 -0.5
    Для справки, производные функций f и g равны:
    f'(x1) = 2*(x1-1)
    g'(x2) = -4*(x2-4.5)

    Нам надо минимизировать
    D^2 = (f(x1)-g(x2))^2 + (x1-x2)^2

    Для минимума необходимо, чтобы обе частные производные по x1 и x2 были равны нулю:
    1. 2*(f(x1)-g(x2))*f'(x1) + 2*(x1-x2) = 0
    2. 2*(g(x2)-f(x1))*g'(x2) + 2*(x2-x1) = 0

    Складываем уравнения 1. и 2. и получаем:
    2*(f(x1)-g(x2))*(f'(x1)-g'(x2))=0
    то есть, либо f(x1)=g(x2), но тогда в совокупности с уравнением 1. и x1=x2 и параболы пересекаются, а наши параболы не пересекаются, значит - f'(x1)=g'(x2) (это как раз то, что Григорий говорил о паралельных касательных) или
    3. (x1-1) = -2*(x2-4.5)
    Введем переменную t=x1-1 и выразим уравнение 1. через t с учетом равенства 3. ((x2-4.5)=-0.5*t):
    4. 2*(t^2-1 + 2*0.25*t^2 + 0.5)*2*t +2*(t+1 + 0.5*t - 4.5) = 0
    или
    5. 6*t^3 + t - 7 =0
    Вот этот довесочек в t я и прозевал вчера, но и с ним решение угадывается очень легко, обратите внимание, что сумма коэффициентов равна нулю, и значит t=1 является корнем, а других вещественных корней нет. По уму надо конечно посчитать вторые производные и убедиться, что Якобиан положительно определен, но по сути понятно, что ничего кроме минимума у нас быть в данном случае не может.
    Ну а чему равно D?
    x1=t+1 = 2,
    x2 = 4.5 - 0.5*t = 4,
    f(x1) = 0,
    g(x2) = -1,
    D = sqrt((0+1)^2 + (2-4)^2) = sqrt(5)

Поделиться этой страницей

Статус темы:
Закрыта.