Русская Википедия:Преобразование Вейерштрасса
В математике преобразование Вейерштрасса[1] функции Шаблон:Math, названное в честь Карла Вейерштрасса, представляет собой «сглаженную» версию Шаблон:Math, полученную путём усреднения значений Шаблон:Mvar, взвешенных с помощью гауссиана с центром в точке x.
В частности, это функция Шаблон:Mvar, определённая формулой
- <math>F(x)=\frac{1}{\sqrt{4\pi}}\int_{-\infty}^\infty f(y) \; e^{-\frac{(x-y)^2}{4}} \; dy = \frac{1}{\sqrt{4\pi}}\int_{-\infty}^\infty f(x-y) \; e^{-\frac{y^2}{4}} \; dy~,</math>
то есть свёртка Шаблон:Mvar с функцией Гаусса
- <math>\frac{1}{\sqrt{4\pi}} e^{-x^2/4}~.</math>
Коэффициент 1/√(4π) выбран так, чтобы общий интеграл гауссиана был равен 1, вследствие чего постоянные функции не изменяются преобразованием Вейерштрасса.
Вместо Шаблон:Math также пишется Шаблон:Math. Заметим, что Шаблон:Math не обязательно должно существовать для каждого действительного числа Шаблон:Mvar, когда определяющий интеграл не сходится.
Преобразование Вейерштрасса тесно связано с уравнением теплопроводности (или, что эквивалентно, уравнению диффузии с постоянным коэффициентом диффузии). Если функция Шаблон:Mvar описывает начальную температуру в каждой точке бесконечно длинного стержня, который имеет постоянную теплопроводность, равную 1, то распределение температуры стержня через единицу времени t = 1 будет задано функцией F. Используя значения t, отличные от 1, мы можем определить обобщённое преобразование Вейерштрасса для Шаблон:Mvar.
Обобщённое преобразование Вейерштрасса предоставляет средство для сколь угодно хорошей аппроксимации заданной интегрируемой функции Шаблон:Mvar аналитическими функциями.
Названия
Вейерштрасс использовал это преобразование в своём первоначальном доказательстве теоремы Вейерштрасса — Стоуна. Оно также известно как преобразование Гаусса или преобразование Гаусса-Вейерштрасса в честь Карла Гаусса и как преобразование Хилле в честь Шаблон:Iw, который его тщательно изучал. Обобщение преобразования Вейерштрасса, упомянутое ниже, известно в обработке сигналов как фильтр Гаусса, а в обработке изображений (при реализации на R2) как размытие по Гауссу.
Преобразования некоторых важных функций
Как упоминалось выше, каждая постоянная функция является собственным преобразованием Вейерштрасса. Преобразование Вейерштрасса любого многочлена является многочленом той же степени и фактически тем же ведущим коэффициентом (асимптотическое разложение не меняется). Действительно, если Шаблон:Math обозначает многочлены Эрмита степени n, тогда преобразование Вейерштрасса Шаблон:Math(Шаблон:Mvar/2) — это просто Шаблон:Math. Это можно показать, используя тот факт, что производящая функция последовательности для многочленов Эрмита тесно связана с ядром Гаусса, используемым в определении преобразования Вейерштрасса.
Преобразованием Вейерштрасса функции eax (где a является произвольной константой) является ea2 eax. Функция eax таким образом, является Шаблон:Iw преобразования Вейерштрасса. (Это, на самом деле, в более общем случае верно для всех свёрточных преобразований.)
Установив a=bi, где i является мнимой единицей, и применив тождество Эйлера, можно увидеть, что преобразование Вейерштрасса функции cos(bx) равно e−b2 cos(bx), а преобразование Вейерштрасса функции sin(bx) равна e−b2 sin(bx).
Преобразование Вейерштрасса функции eax2 имеет вид
- <math>\frac{1}{\sqrt{1-4a}}e^{\frac{ax^2}{1-4a}}</math>, если a < 1/4, и не определено, если a ≥ 1/4.
В частности, при выборе отрицательного значения очевидно, что преобразование Вейерштрасса гауссовой функции снова является гауссовой функцией, но «более широкой».
Общие свойства
Преобразование Вейерштрасса ставит в соответствие каждой функции f новую функцию F; это соответствие является линейным. Оно также является трансляционно-инвариантным, что означает, что преобразованием функции f(x + a) является F(x + a). Оба эти факта в более общем случае справедливы для любого интегрального преобразования, определённого с помощью свёртки.
Если преобразование F(x) существует для действительных чисел x = a и x = b, тогда оно также существует для всех промежуточных действительных значений и образует там аналитическую функцию; кроме того, F(x) будет существовать для всех комплексных чисел x с a ≤ Re(x) ≤ b и образует голоморфную функцию на этой полоске комплексной плоскости. Это формальное утверждение о «гладкости» F, упомянутой выше.
Если функция f интегрируема по всей вещественной оси (то есть f ∈ L1(R)), то это справедливо и для её преобразования Вейерштрасса F, и если, кроме того, f(x) ≥ 0 для всех x, тогда также F(x) ≥ 0 для всех x и интегралы от f и F равны. Это выражает физический факт, что полная тепловая энергия или тепло сохраняется уравнением теплопроводности или что общее количество диффундирующего вещества сохраняется уравнением диффузии.
Используя вышесказанное, можно показать, что для 0 < p ≤ ∞ и f ∈ Lp(R), мы имеем F ∈ Lp(R) и ||F||p ≤ ||f||p. Следовательно, преобразование Вейерштрасса даёт ограниченный оператор W : Lp(R) → Lp(R).
Если f является достаточно гладкой, тогда преобразование Вейерштрасса k-й производной функции f равно k-ой производной преобразования Вейерштрасса от f.
Существует формула, связывающая преобразование Вейерштрасса W и двустороннее преобразование Лапласа L. Если мы определяем
- <math>g(x)=e^{-\frac{x^2}{4}} f(x),</math>
тогда
- <math>W[f](x)=\frac{1}{\sqrt{4\pi}} e^{-x^2/4} L[g]\left(-\frac{x}{2}\right).</math>
Фильтр нижних частот
Шаблон:DetailsМы видели выше, что преобразованием Вейерштрасса cos(bx) является e−b2 cos(bx), и аналогично для sin(bx). В терминах обработки сигналов это говорит о том, что если сигнал f содержит частоту b (то есть содержит слагаемое, представляющее собой комбинацию sin(bx) и cos(bx)), тогда преобразованный сигнал F будет содержать ту же частоту, но с амплитудой, умноженной на коэффициент e−b2. Это приводит к тому, что более высокие частоты уменьшаются больше, чем более низкие, и, таким образом, преобразование Вейерштрасса действует как фильтр нижних частот. Это также можно показать с помощью преобразования Фурье. Преобразование Фурье анализирует сигнал с точки зрения его частот, преобразует свёртки в произведения и преобразует функцию Гаусса в функцию Гаусса. Преобразование Вейерштрасса является свёрткой с гауссианом и, следовательно, произведением преобразования Фурье сигнала и гауссиана, с последующим применением обратного преобразования Фурье. Это умножение на гауссиан в частотном пространстве смешивает высокие частоты, что является другим способом описания свойства «сглаживания» преобразования Вейерштрасса.
Обратное преобразование
Относительно легко установить следующую формулу, тесно связанную с преобразованием Лапласа функции Гаусса и действительным аналогом Шаблон:Iw:
- <math>e^{u^2}=\frac{1}{\sqrt{4\pi}} \int_{-\infty}^\infty e^{-uy} e^{-y^2/4}\;dy.</math>
Теперь заменим в ней u оператором формального дифференцирования D = d/dx и используем Шаблон:Iw
- <math>e^{-yD}f(x)=f(x-y)</math>,
(следствие формулы ряда Тейлора и определения показательной функции), чтобы получить
- <math>
\begin{align} e^{D^2}f(x) & = \frac{1}{\sqrt{4\pi}} \int_{-\infty}^\infty e^{-yD}f(x) e^{-y^2/4}\;dy \\ & =\frac{1}{\sqrt{4\pi}} \int_{-\infty}^\infty f(x-y) e^{-y^2/4}\;dy=W[f](x) \end{align} </math>
чтобы таким образом получить следующее формальное выражение для преобразования Вейерштрасса W,Шаблон:Equation box 1где оператор справа следует понимать как действующий на функцию f(x) как
- <math>e^{D^2} f(x) = \sum_{k=0}^\infty \frac{D^{2k}f(x)}{k!}~.</math>
Приведённый выше формальный вывод скрывает детали сходимости, а формула W = eD2, таким образом, не является универсальной; существует несколько функций f, которые имеют чётко определённое преобразование Вейерштрасса, но для которых eD2f(x) не может быть определено осмысленно.
Тем не менее, правило всё ещё довольно полезно и может, например, использоваться для получения преобразований Вейерштрасса полиномов, экспоненциальных и тригонометрических функций, упомянутых выше.
Таким образом, формальное обратное преобразование Вейерштрасса задаётся формулой
- <math>W^{-1}=e^{-D^2} ~. </math>
Эта формула не является универсальной, но может служить руководством. Можно показать, что она корректна для определённых классов функций, если правильно определён оператор правой части[2].
В качестве альтернативы можно попытаться обратить преобразование Вейерштрасса немного другим способом: учитывая аналитическую функцию
- <math>F(x)=\sum_{n=0}^\infty a_n x^n ~,</math>
применить W−1 для получения
- <math>f(x)=W^{-1}[F(x)]=\sum_{n=0}^\infty a_n W^{-1}[x^n]=\sum_{n=0}^\infty a_n H_n(x/2),</math>
ещё раз используя фундаментальное свойство многочленов Эрмита Шаблон:Math.
Опять же, эта формула для f(x) в лучшем случае формальна, поскольку не проверялось, сходится ли конечный ряд. Но если, например, f ∈ L2(R), то знание всех производных от F в x = 0 достаточно для получения коэффициентов an; таким образом, можно представить Шаблон:Mvar в виде ряда многочленов Эрмита.
Третий метод обращения преобразования Вейерштрасса использует его связь с преобразованием Лапласа, упомянутым выше, и хорошо известной формулой обращения преобразования Лапласа. Результат приведён ниже для распределений.
Обобщения
Мы можем использовать свёртку с ядром Гаусса <math display="inline">\frac{1}{\sqrt{4\pi t}} e^{-\frac{x^2}{4t}}</math> (для Шаблон:Math) вместо <math display="inline">\frac{1}{\sqrt{4\pi}} e^{-\frac{x^2}{4}}</math>, таким образом определяя оператор Шаблон:Math обобщённого преобразования Вейерштрасса.
При малых значениях Шаблон:Math очень близко к Шаблон:Math, но гладко. Чем больше Шаблон:Math, тем больше этот оператор усредняет и изменяет Шаблон:Math. Физически, Шаблон:Math соответствует уравнению теплопроводности (или диффузии) для единицы времени Шаблон:Math и аддитивно:<math display="block">W_s \circ W_t = W_{s+t},</math>то есть "диффузия для единицы времени Шаблон:Math, а затем единицы времени Шаблон:Math, эквивалентна диффузии для единицы времени Шаблон:Math ". Можно распространить это на Шаблон:Math, установив Шаблон:Math тождественным оператором (то есть свёрткой дельта-функции) и таким образом получив полугруппу операторов.
Ядро <math display="inline">\frac{1}{\sqrt{4\pi t}} e^{-\frac{x^2}{4t}},</math> используемое для обобщённого преобразования Вейерштрасса, иногда называют ядром Гаусса-Вейерштрасса, это функция Грина для уравнения диффузии <math>(\partial_t -D^2) (e^{tD^2} f(x))=0 </math> на Шаблон:Math.
Шаблон:Math может быть вычислен из Шаблон:Math: задана функция Шаблон:Math, определим новую функцию Шаблон:Math; тогда Шаблон:Math, что является следствием правила замены переменной.
Преобразование Вейерштрасса также может быть определено для обобщённой функции[3]. Например, преобразование Вейерштрасса дельта-функции является гауссианом<math display="inline">\frac{1}{\sqrt{4\pi}} e^{-x^2/4}</math>.
В этом контексте могут быть доказаны строгие формулы обращения, например,<math display="block">f(x)=\lim_{r\to\infty}\frac{1}{i\sqrt{4\pi}} \int_{x_0-ir}^{x_0+ir} F(z)e^{\frac{(x-z)^2}{4}}\;dz</math>где Шаблон:Math — любое фиксированное действительное число, для которого Шаблон:Math существует, интеграл простирается по вертикальной линии в комплексной плоскости с вещественной частью Шаблон:Math, и предел следует принимать в смысле распределений.
Кроме того, преобразование Вейерштрасса может быть определено для вещественных (или комплексных) функций (или распределений), определённых на Шаблон:Math. Мы используем ту же формулу свёртки, что и выше, но интерпретируем интеграл как распространяющийся на всё Шаблон:Math и выражение Шаблон:Math как квадрат евклидовой длины вектора Шаблон:Math; множитель перед интегралом должен быть скорректирован так, чтобы общий интеграл гауссиана был равен 1.
В более общем смысле, преобразование Вейерштрасса может быть определено на любом Римановом пространстве: уравнение теплопроводности может быть сформулировано там (используя оператор Лапласа — Бельтрами), и преобразование Вейерштрасса Шаблон:Math получается путём решения уравнения теплопроводности за единицу времени, начиная с начального «распределения температуры» Шаблон:Math.
Связанные преобразования
Если рассматривать свёртку с ядром Шаблон:Math вместо гауссиана, получается ядро Пуассона, которое сглаживает и усредняет заданную функцию способом, аналогичным преобразованию Вейерштрасса.
См. также
Примечания
- ↑ Ahmed I. Zayed, Handbook of Function and Generalized Function Transformations, Chapter 18. CRC Press, 1996.
- ↑ G. G. Bilodeau, «The Weierstrass Transform and Hermite Polynomials». Duke Mathematical Journal 29 (1962), p. 293—308
- ↑ Yu A. Brychkov, A. P. Prudnikov. Integral Transforms of Generalized Functions, Chapter 5. CRC Press, 1989
