ny_quant: (Default)
ny_quant ([personal profile] ny_quant) wrote2016-07-22 10:27 am
Entry tags:

Математическое


Рассмотрим задачу минимизации в Rn:

F(x) -> min

при ограничениях

G(x)=0 Edit: обсуждение показало, что это условие лишнее.
H(x)>=0

Хочется сформулировать такого типа теорему, что при разумных ограничениях на функции F,G,H (скажем F, видимо, должна быть выпуклой) решение задачи (Edit: под этим понимается точка, где достигается инфимум по допустимой области если таковая существует) либо совпадает с глобальным минимумом F(x) либо лежит на границе допустимой области.

Поскольку я это придумал сегодня по дороге на работу, я вижу два варианта. Либо это совсем неверно по каким-то очевидным причинам, которые мне с утра не пришли в голову. Либо это давно все знают и умные люди легко подскажут где найти соответствующую теорему.

[identity profile] ex-juan-gan.livejournal.com 2016-07-22 03:02 pm (UTC)(link)
Как вообще на многомерном пространстве определяется выпуклость? По-моему, если взять сферу, то минимум достигается на границе. Соответственно, какая бы ни была фигура, если минимум достигается внутри, то есть такая сфера, нарушающая условие.

[identity profile] ny-quant.livejournal.com 2016-07-22 03:23 pm (UTC)(link)

Выпуклость функции определяется точно также как и в одномерном случае: «надграфик», т.е. множество {F(X)>=a} является выпуклой фигурой.


Последнего не понял. Сфера вся состоит из границы. Вот если допустимой областью является шар, то минимум может достигаться и внутри. Например, в одномерном случае шар это отрезок, скажем [-1, 1]. Функция x^2 достигает минимума внутри при x=0, в точке глобального минимума. Но функция (x-2)^2 уже имеет глобальный минимум за пределами [-1, 1] и она уже достигает минимума на границе при x=1.

[identity profile] ex-juan-gan.livejournal.com 2016-07-22 05:00 pm (UTC)(link)
В смысле шар.

Начнем с выпуклости функции. Тогда (локальный) максимум один, это как бы очевидно (а то прямая между двумя локальными не входит в надграфик). Ну и взять этот глобальный максимум; он или внутри фигуры, или нет (логика у нас пока что булева). Если нет, то ни один открытый шар внутри фигуры не содержит максимума - доказывать? А раз нет, то максимум по фигуре будет на границе.

[identity profile] ny-quant.livejournal.com 2016-07-22 05:06 pm (UTC)(link)

Именно.

(deleted comment)

[identity profile] ex-juan-gan.livejournal.com 2016-07-23 11:53 pm (UTC)(link)
Я тоже так думаю.

[identity profile] fortran-only.livejournal.com 2016-07-22 03:23 pm (UTC)(link)
Давно известно, и тривиально для выпуклых функций и ограничений, там всегда локальный максимум есть глобальный максимум чисто по определению.

Для невыпуклых функций нет ничего.

[identity profile] ny-quant.livejournal.com 2016-07-22 03:40 pm (UTC)(link)

См выше: функция (x-2)^2 имеет глобальный минимум за пределами допустимой области. То что этот глобальный минимум также и единственный локальный действительно тривиально, но никак мне не помогает доказать, что в моей задаче минимум лежит на границе. Вернее, помогает, но это другая теорема.

Edited 2016-07-22 15:45 (UTC)

[identity profile] fortran-only.livejournal.com 2016-07-22 03:52 pm (UTC)(link)
Если допустимая область образует выпуклый сет, то локальный минимум также является глобальным минимумом в допустимой области.

Все. Насчет где именно лежит минимум на то есть теорема Каруша-Кун-Такера от 1939 года, которую знает каждый экономист буквально с пеленок.

[identity profile] ny-quant.livejournal.com 2016-07-22 04:35 pm (UTC)(link)

Первое тривиально, но это другая теорема.


Я не вижу чтоб теорема Куна-Таккера говорила о том, где находится точка минимума. Она говорит о существовании множителей Лагранжа с определёнными свойствами.


Либо моё утверждение верно, либо нет, либо накрайняк это новая теорема и никто не знает.

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-22 17:11 (UTC) - Expand

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-22 17:34 (UTC) - Expand
(deleted comment)

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-24 15:10 (UTC) - Expand
(deleted comment)

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-24 16:24 (UTC) - Expand

[identity profile] sovkista.livejournal.com 2016-07-23 08:25 pm (UTC)(link)
которую знает каждый экономист буквально с пеленок.
Блестящий образец ужас-совкизма, который мне мешает выделить в отдельный пост только лишь моя чрезвычайная загруженность.
(deleted comment)

(no subject)

[identity profile] sovkista.livejournal.com - 2016-07-23 22:28 (UTC) - Expand
(deleted comment)

(no subject)

[identity profile] sovkista.livejournal.com - 2016-07-23 23:22 (UTC) - Expand
(deleted comment)

(no subject)

[identity profile] sovkista.livejournal.com - 2016-07-24 00:28 (UTC) - Expand
(deleted comment)

(no subject)

[identity profile] sovkista.livejournal.com - 2016-07-24 02:02 (UTC) - Expand

(no subject)

[identity profile] cass1an.livejournal.com - 2016-07-23 23:39 (UTC) - Expand

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-24 15:20 (UTC) - Expand

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-24 15:06 (UTC) - Expand
(deleted comment)

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-24 16:29 (UTC) - Expand

[identity profile] misha-b.livejournal.com 2016-07-22 04:45 pm (UTC)(link)
A convex function has exactly one minimum on a convex domain. If that minimum is different from the global minimum (in R^n), то он очевидно лежит на границе области, т.к. у выпуклой функции есть только один локальный (он же глобальный) минимум.

[identity profile] ny-quant.livejournal.com 2016-07-22 05:06 pm (UTC)(link)

Я и говорю. А нужна ли выпуклость области?

[identity profile] misha-b.livejournal.com 2016-07-22 05:13 pm (UTC)(link)

Да. Простейший пример -- минимизировать y=x^2 по области
[-2,-1]U[1,2].

[identity profile] ny-quant.livejournal.com 2016-07-22 05:28 pm (UTC)(link)

Может всё же односвязности достаточно?

[identity profile] ny-quant.livejournal.com 2016-07-22 05:30 pm (UTC)(link)

Подождите. Ну, два минимума +1 и -1, оба на границе. Чего я не понял?

(no subject)

[identity profile] misha-b.livejournal.com - 2016-07-22 18:15 (UTC) - Expand

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-22 18:26 (UTC) - Expand

(no subject)

[identity profile] misha-b.livejournal.com - 2016-07-22 18:28 (UTC) - Expand

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-22 18:44 (UTC) - Expand

(no subject)

[identity profile] misha-b.livejournal.com - 2016-07-22 21:02 (UTC) - Expand

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-22 19:09 (UTC) - Expand

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-22 20:52 (UTC) - Expand

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-22 21:53 (UTC) - Expand

(no subject)

[identity profile] misha-b.livejournal.com - 2016-07-22 20:57 (UTC) - Expand

(no subject)

[identity profile] misha-b.livejournal.com - 2016-07-22 22:05 (UTC) - Expand

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-22 21:58 (UTC) - Expand

(no subject)

[identity profile] misha-b.livejournal.com - 2016-07-23 03:42 (UTC) - Expand

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-24 15:23 (UTC) - Expand

(no subject)

[identity profile] misha-b.livejournal.com - 2016-07-24 17:42 (UTC) - Expand

[identity profile] spamsink.livejournal.com 2016-07-22 11:46 pm (UTC)(link)
Будем в R^2. Возьмем параболоид вращения, это у нас выпуклая F(x), проведем на его поверхности какую-нибудь замкнутую извилистую загогулину, это у нас G(x)=0, и сотрем половину загогулины; стёртое - это у нас H(x) < 0, а оставшееся - H(x) >= 0. Понятно, что самая нижняя точка оставшейся части загогулины может с легкостью быть вовсе не на конце.

[identity profile] ny-quant.livejournal.com 2016-07-23 02:52 pm (UTC)(link)
Не понял что отсюда следует.

Но вот [livejournal.com profile] juan_gandhi предложил простое доказательство, в котором не используется выпуклость области. Что не так?

http://ny-quant.livejournal.com/595516.html?thread=4899900#t4899900

[identity profile] spamsink.livejournal.com 2016-07-23 03:11 pm (UTC)(link)
Раз не на конце, то, значит, не на границе допустимой области.

[identity profile] ny-quant.livejournal.com 2016-07-23 05:45 pm (UTC)(link)
Этого тоже не понял.

(no subject)

[identity profile] spamsink.livejournal.com - 2016-07-23 18:38 (UTC) - Expand

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-07-24 15:10 (UTC) - Expand

[identity profile] spamsink.livejournal.com 2016-07-23 04:31 pm (UTC)(link)
Ганди доказывал без учета G(x) = 0.

[identity profile] ny-quant.livejournal.com 2016-07-23 05:44 pm (UTC)(link)
G(x) = 0 снижает размерность допустимой области. Любая допустимая точка лежит на границе допустимой области, так что утверждение становится тривиальным. Не надо было мне это ограничение добавлять. Не подумал.
Edited 2016-07-23 17:47 (UTC)

[identity profile] f-b-dima.livejournal.com 2016-09-01 04:31 pm (UTC)(link)
выпуклая либо монотонна, либо единственный экстремум. Таким образом любой непрерывный интервал попадает либо на монотонный участок, либо на участок с экстремумом. В первом случае решение на краю, во втором случае, если экстремум максимум - тоже, а если минимум, то этот минимум и будет решением, т.е. то, что и утверждалось.

[identity profile] ny-quant.livejournal.com 2016-09-01 07:47 pm (UTC)(link)
Причем выпуклость области в этом рассуждении нигде не использовалась. Попробуйте это объяснить моему альтернативно одаренному гостю, который любит фортран.

[identity profile] f-b-dima.livejournal.com 2016-09-02 05:54 am (UTC)(link)
Использовались свойства выпуклости (монотонность, либо единственный экстремум) такая лемма была точно.

Если у любителя фортрана есть проблемы с установлением программным способом факта выпуклости, то могу помочь )

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-09-02 13:19 (UTC) - Expand

(no subject)

[identity profile] f-b-dima.livejournal.com - 2016-09-02 17:07 (UTC) - Expand

(no subject)

[identity profile] ny-quant.livejournal.com - 2016-09-02 18:11 (UTC) - Expand