ny_quant: (Default)
ny_quant ([personal profile] ny_quant) wrote2011-02-25 07:53 pm

Требуется ортогональный базис на вещественной оси

Полиномы и триг. функции я знаю. А что еще бывает?

[identity profile] kdv2005.livejournal.com 2011-02-26 12:55 am (UTC)(link)
Э-э, тебе гладкие или разрывные тоже подойдут?

[identity profile] c0nfigure.livejournal.com 2011-02-26 01:08 am (UTC)(link)
если разрывные, то можно с wavelet basis начать

[identity profile] ny-quant.livejournal.com 2011-02-26 03:01 am (UTC)(link)
Нам бы гладких ...

[identity profile] ny-quant.livejournal.com 2011-02-26 03:06 am (UTC)(link)
Я бы хотел что-то неосциллирующее. Но мысль интересная.

[identity profile] rsokolov.livejournal.com 2011-02-26 03:09 am (UTC)(link)
А вам надо на отрезке или на всей оси? И функции должны быть "чисто ортогональными" или можно с весовой функцией?

[identity profile] ny-quant.livejournal.com 2011-02-26 03:16 am (UTC)(link)
В принципе на отрезке тоже могло бы подойти, т.к. я с отрезка некоторым "естественным" способом начинаю. Но мне могут понадобиться бесконечности на обеих концах, поэтому я больше думаю об отображении на ось. С весами или нет мне все равно. Задача чисто практическая - о приближении заданной функции комбинацией нескольких элементов базиса. Старо как мир, в сущности.

[identity profile] rsokolov.livejournal.com 2011-02-26 03:20 am (UTC)(link)
Просто не очень понятно, что значит "неосциллирующее", "ортогональное" и "интегрируемое на бесконечности".

А чем вас функции Эрмита не устраивают?

[identity profile] ny-quant.livejournal.com 2011-02-26 03:37 am (UTC)(link)
Ну, мне главное, что б приближение получилось хорошее. Исходя из природы аппроксимируемых функций, хорошо бы, что бы у элементов базиса было поменьше корней и экстремумов.

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

[identity profile] kdv2005.livejournal.com 2011-02-26 04:23 am (UTC)(link)
По-моему, так не бывает, чтобы приближающие функции не осциллировали или мало осциллировали, если они ортогональный базис образуют, хоть в пространстве с весом, хоть без оного. Ну сам посуди.

[identity profile] ny-quant.livejournal.com 2011-02-26 02:49 pm (UTC)(link)
В голове вспыхивают цитаты и кинофильмов, типа

- Это ловушка!

или

- Вы даёте неральные планы!

или даже

- Будем искать!

Соколов уже напомнил мне ниже теорему о корнях ортогонального базиса. Такое впечатление, что все же придется изобретать велосипед самому.

[identity profile] kdv2005.livejournal.com 2011-02-27 01:24 am (UTC)(link)
Я думаю, что лучше обсудить вживую.

[identity profile] ny-quant.livejournal.com 2011-02-27 02:13 am (UTC)(link)
Доживем до понедельника :)

[identity profile] ny-quant.livejournal.com 2011-02-26 03:37 am (UTC)(link)
Эрмит не устраивает тем, что многочлены.

[identity profile] rsokolov.livejournal.com 2011-02-26 04:06 am (UTC)(link)
Ну, что значит, "многочлены"? Функции-то неполиномиальные.

А так, любая функция раскладывается в ряд Тейлора или, на крайний случай, Лорана. Рациональными функциями еще можно аппроксимировать. Да чем угодно, на самом деле.

А если базис ортогональный, то у его n-ного элемента будет как минимум n-1 корней.

[identity profile] ny-quant.livejournal.com 2011-02-26 02:44 pm (UTC)(link)
Виноват, это был дурацкий ответ. Сейчас поясню.

Приближение мночленами Эрмита (или кого угодно) можно ведь потом переписать как самый обычный многочлен. Поскольку задача все равно решается численной оптимизацией невязки, какой бы полиномиальный базис я не использовал, результат получается один и тот же.

Свойство, от которого хочется избавиться это то, что все многочлены проходят через точку (1,1).

Рациональные функции, в классической формулировке, к сожалению, не подойдут, т.к. аппроксимируемая функция известна только в нескольких точках, да и степеней свободы у меня всего несколько. Я как рaз думал о том, чтобы определить "базис" как семейство функций типа (x^n+a_n)/(x^2+1). Но решил пока не изобретать велосипед, а спросить, что умные люди подскажут.

Теорему о корнях ортогонального базиса я забыл, конечно. Видимо, от ортогональности придется отказаться.

[identity profile] rsokolov.livejournal.com 2011-02-26 03:13 pm (UTC)(link)
(Функции Эрмита и полиномы Эрмита - вещи слегка разные, но не суть).

А у вас точка (x=0) чем-то отличается от прочих точек? Или же с точки зрения "физики" все точки более-менее эквивалентны?

[identity profile] ny-quant.livejournal.com 2011-02-26 04:05 pm (UTC)(link)
Да, отличается. Значение аппроксимируемой функции в нуле известно лучше, чем в других точках. Модель строится так, как будто значение в нуле известно точно. В принципе, это не обязательно, но так удобнее.

[identity profile] rsokolov.livejournal.com 2011-02-26 04:17 pm (UTC)(link)
Ну, тогда, наверное, надо подробнее разбираться в "физическом смысле" вашей функции, чем она могла бы быть в самом простом идеальном случае, какие поправки к ней имеют физический смысл, etc.

[identity profile] prof-yura.livejournal.com 2011-02-27 03:44 am (UTC)(link)
Полиномы, помноженные на убывающую экспоненту $\exp{-x^2/2}$.

О! Тяжелая артиллерия подтянулась!

[identity profile] ny-quant.livejournal.com 2011-02-27 04:32 pm (UTC)(link)
Видимо это будет не ортогональный базис, но и бог с ним. Это вполне может получиться, но мне надо будет подсунуть какую-то удачно подобранную последовательность в эспоненту, например записать энный член как

x^n * exp{-x^2/(2*n)}

На что бы я мог сослаться, говоря о свойствах такого разложения?

Спасибо!

Это будет ортогональный базис

[identity profile] prof-yura.livejournal.com 2011-02-27 04:43 pm (UTC)(link)
С весом $exp(-x^2)$.

Попробуйте посмотреть ссылки в http://en.wikipedia.org/wiki/Hermite_polynomials .

Re: Это будет ортогональный базис

[identity profile] ny-quant.livejournal.com 2011-02-27 05:47 pm (UTC)(link)
А, про полиномы Эрмита я и сам знал. Все же чему-то учили. Мне что-то другое нужно.