Функция Гаусса - Gaussian function

В математика, а Функция Гаусса, часто называемый просто Гауссовский, это функция формы

для произвольных настоящий константы а, б и ненулевой c. Назван в честь математика. Карл Фридрих Гаусс. В график гауссиана является характеристикой симметричной "кривая колокола "форма. Параметр а - высота пика кривой, б это положение центра пика и cстандартное отклонение, иногда называемый гауссовой RMS width) контролирует ширину «колокола».

Гауссовские функции часто используются для представления функция плотности вероятности из нормально распределенный случайная переменная с ожидаемое значение μ = б и отклонение σ2 = c2. В этом случае гауссиан имеет вид:

[1]

Функции Гаусса широко используются в статистика описать нормальные распределения, в обработка сигналов определять Гауссовы фильтры, в обработка изображений где двумерные гауссианы используются для Размытие по Гауссу, а в математике решить уравнения теплопроводности и уравнения диффузии и определить Преобразование Вейерштрасса.

Характеристики

Гауссовы функции возникают при составлении экспоненциальная функция с вогнутый квадратичная функция:

куда:

Таким образом, гауссовские функции - это те функции, логарифм - вогнутая квадратичная функция.

Параметр c относится к полная ширина на половине максимальной (FWHM) пика согласно

Затем функция может быть выражена через FWHM, представленную ш:

В качестве альтернативы параметр c можно интерпретировать, сказав, что два точки перегиба функции происходят в Икс = б − c и Икс = б + c.

В полная ширина в десятых долях от максимума (FWTM) для гауссиана может представлять интерес и является

Гауссовы функции аналитический, и их предел в качестве Икс → ∞ равно 0 (для указанного выше случая б = 0).

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

и получается

Этот интеграл равен 1 тогда и только тогда, когда нормализующая константа ), и в этом случае гауссовский функция плотности вероятности из нормально распределенный случайная переменная с ожидаемое значение μ = б и отклонение σ2 = c2:

Эти гауссианы показаны на прилагаемом рисунке.

Нормализованный Гауссовы кривые с ожидаемое значение μ и отклонение σ2. Соответствующие параметры: , б = μ и c = σ.

Гауссовские функции с центром в нуле минимизируют Фурье принцип неопределенности.

Произведение двух гауссовых функций является гауссовой, а свертка двух гауссовских функций также является гауссовой, причем дисперсия представляет собой сумму исходных дисперсий: . Однако произведение двух гауссовых функций плотности вероятности (PDF), как правило, не является гауссовской PDF.

Принимая Преобразование Фурье (унитарная, угловая частота) функции Гаусса с параметрами а = 1, б = 0 и c дает другую функцию Гаусса с параметрами , б = 0 и .[2] Так, в частности, гауссовские функции с б = 0 и фиксируются преобразованием Фурье (они собственные функции преобразования Фурье с собственным значением 1) .Физическая реализация дифракционная картина: например, a фотографический слайд чей коэффициент пропускания имеет гауссовскую вариацию, также является гауссовой функцией.

Тот факт, что функция Гаусса является собственной функцией непрерывного преобразования Фурье, позволяет нам получить следующие интересные[требуется разъяснение ] личность из Формула суммирования Пуассона:

Интеграл от функции Гаусса

Интеграл от произвольной гауссовой функции равен

Альтернативная форма -

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

Связь со стандартным гауссовским интегралом

Интегральный

для некоторых настоящий константы a, b, c> 0 можно вычислить, представив его в виде Гауссовский интеграл. Во-первых, постоянная а можно просто вычесть из интеграла. Далее переменная интегрирования изменяется с Икс к у = Икс - б.

а затем в

Затем, используя Гауссовское интегральное тождество

у нас есть

Двумерная функция Гаусса

Гауссова кривая с двумерной областью

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

Частным примером двумерной функции Гаусса является

Здесь коэффициент А это амплитуда, Иксо, yо центр, а σИкс, σу являются Икс и у распространение капли. Рисунок справа был создан с использованием А = 1, Иксо = 0, уо = 0, σИкс = σу = 1.

Объем под функцией Гаусса определяется выражением

В общем, двумерная эллиптическая функция Гаусса выражается как

где матрица

является положительно определенный.

Используя эту формулировку, рисунок справа можно создать с помощью А = 1, (Иксо, уо) = (0, 0), а = c = 1/2, б = 0.

Значение параметров для общего уравнения

Для общего вида уравнения коэффициент А - высота пика и (Иксоуо) является центром капли.

Если мы установим

затем мы поворачиваем каплю на угол по часовой стрелке (для вращения против часовой стрелки переверните знаки в б коэффициент).[3] Это можно увидеть на следующих примерах:

Используя следующие Октава кода, легко увидеть эффект изменения параметров

А = 1;x0 = 0; y0 = 0;sigma_X = 1;sigma_Y = 2;[Икс, Y] = сетка(-5:.1:5, -5:.1:5);за тета = 0:число Пи/100:число Пи    а = потому что(тета)^2/(2*sigma_X^2) + грех(тета)^2/(2*sigma_Y^2);    б = -грех(2*тета)/(4*sigma_X^2) + грех(2*тета)/(4*sigma_Y^2);    c = грех(тета)^2/(2*sigma_X^2) + потому что(тета)^2/(2*sigma_Y^2);    Z = А*exp( - (а*(Икс-x0).^2 + 2*б*(Икс-x0).*(Y-y0) + c*(Y-y0).^2));серфить(Икс,Y,Z);затенение интерп;Посмотреть(-36,36)подождите, нажмите кнопкуконец

Такие функции часто используются в обработка изображений и в вычислительных моделях зрительная система функция - см. статьи на масштабное пространство и аффинный shn.

Также см многомерное нормальное распределение.

Гауссова или супергауссова функция высшего порядка

Более общая формулировка функции Гаусса с плоской вершиной и спадом по Гауссу может быть взята путем возведения содержания показателя в степень, :

Эта функция известна как супергауссова функция и часто используется для формулировки гауссова пучка.[4] В двумерной постановке функция Гаусса вдоль и можно комбинировать с потенциально разными и чтобы сформировать эллиптическое гауссово распределение, или прямоугольное распределение Гаусса, .[5]

Многомерная функция Гаусса

В -мерном пространстве функцию Гаусса можно определить как

куда столбец координаты, это положительно определенный матрица и обозначает транспонирование матрицы.

Интеграл от этой гауссовой функции по всей -мерное пространство задается как

Его легко вычислить, диагонализуя матрицу и заменяя переменные интегрирования на собственные векторы .

В более общем смысле смещенная функция Гаусса определяется как

куда - вектор сдвига и матрица можно считать симметричным, , и положительно-определенный. Следующие интегралы с этой функцией могут быть вычислены с помощью того же метода:

Оценка параметров

Ряд полей, таких как звездная фотометрия, Гауссов пучок характеристика, и спектроскопия эмиссионных / абсорбционных линий работать с выборочными функциями Гаусса и вам необходимо точно оценить параметры высоты, положения и ширины функции. Есть три неизвестных параметра для одномерной функции Гаусса (а, б, c) и пять для двумерной гауссовой функции .

Наиболее распространенный метод оценки гауссовских параметров - это логарифм данных и соответствовать параболе в результирующий набор данных.[6][7] Хотя это обеспечивает простой подгонка кривой процедуры, результирующий алгоритм может быть искажен из-за чрезмерного взвешивания малых значений данных, что может привести к большим ошибкам в оценке профиля. Частично эту проблему можно компенсировать за счет взвешенный метод наименьших квадратов оценка, уменьшая вес небольших значений данных, но это тоже может быть смещено, позволяя хвосту гауссианы преобладать при подгонке. Чтобы устранить предвзятость, вместо этого можно использовать методом наименьших квадратов с повторным взвешиванием процедура, в которой веса обновляются на каждой итерации.[7]Также возможно выполнение нелинейная регрессия непосредственно на данных, без привлечения логарифмическое преобразование данных; Дополнительные параметры см. аппроксимация распределения вероятностей.

Точность параметра

Если у кого-то есть алгоритм для оценки параметров функции Гаусса, также важно знать, как точный эти оценки. Любой наименьших квадратов Алгоритм оценки может предоставить числовые оценки дисперсии каждого параметра (т. е. дисперсии расчетной высоты, положения и ширины функции). Также можно использовать Граница Крамера – Рао теории, чтобы получить аналитическое выражение для нижней границы дисперсии параметров, учитывая определенные предположения о данных.[8][9]

  1. Шум в измеряемом профиле либо i.i.d. Гауссов, или шум Распределенный по Пуассону.
  2. Расстояние между каждой выборкой (т. Е. Расстояние между пикселями, измеряющими данные) одинаково.
  3. Пик является «хорошо отобранным», так что менее 10% площади или объема под пиком (площадь, если гауссиан 1D, объем, если гауссиан 2D), лежит за пределами области измерения.
  4. Ширина пика намного больше, чем расстояние между точками выборки (т.е. пиксели детектора должны быть по крайней мере в 5 раз меньше, чем гауссова FWHM).

Когда эти предположения выполняются, следующие ковариационная матрица K применяется для параметров профиля 1D , , и под i.i.d. Гауссов шум и под шумом Пуассона:[8]

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

а в случае пуассоновского шума

Для параметров 2D профиля, дающих амплитуду , позиция , и ширина профиля применяются следующие ковариационные матрицы:[9]

где дисперсии отдельных параметров даны диагональными элементами ковариационной матрицы.

Дискретный гауссовский

В дискретное гауссово ядро (сплошной), по сравнению с выбранное ядро ​​Гаусса (пунктир) для шкалы

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

Альтернативный подход - использовать дискретное гауссово ядро:[10]

куда обозначает модифицированные функции Бесселя целочисленного порядка.

Это дискретный аналог непрерывного гауссиана в том смысле, что он является решением дискретной уравнение диффузии (дискретное пространство, непрерывное время), точно так же, как непрерывный гауссиан является решением уравнения непрерывной диффузии.[11]

Приложения

Гауссовские функции появляются во многих контекстах в естественные науки, то социальные науки, математика, и инженерное дело. Вот некоторые примеры:

Смотрите также

Рекомендации

  1. ^ Сквайрс, Г. Л. (30 августа 2001 г.). Практическая физика (4-е изд.). Издательство Кембриджского университета. Дои:10.1017 / cbo9781139164498. ISBN  978-0-521-77940-1.
  2. ^ Вайсштейн, Эрик В. «Преобразование Фурье - Гауссово». MathWorld. Получено 19 декабря 2013.
  3. ^ Наури, Николай. "Berechnung von Kovarianzellipsen" (PDF). Получено 14 августа 2019.
  4. ^ Родитель А., М. Морин и П. Лавин. «Распространение супергауссовских распределений поля». Оптическая и квантовая электроника 24.9 (1992): S1071-S1079.
  5. ^ "Руководство по командам оптического программного обеспечения GLAD, Запись по команде GAUSSIAN" (PDF). Прикладные оптические исследования. 2016-12-15.
  6. ^ Каруана, Ричард А .; Searle, Roger B .; Хеллер, Томас .; Шупак, Саул И. (1986). «Быстрый алгоритм разрешения спектров». Аналитическая химия. Американское химическое общество (ACS). 58 (6): 1162–1167. Дои:10.1021 / ac00297a041. ISSN  0003-2700.
  7. ^ а б Хунвэй Го, "Простой алгоритм подбора функции Гаусса", IEEE Sign. Proc. Mag. 28 (9): 134-137 (2011).
  8. ^ а б Н. Хаген, М. Купинский и Э. Л. Дереняк, "Оценка гауссовского профиля в одном измерении", Прил. Опт. 46: 5374–5383 (2007).
  9. ^ а б Н. Хаген и Е. Л. Дерениак, "Оценка гауссовского профиля в двух измерениях", Прил. Опт. 47: 6842–6851 (2008).
  10. ^ Линдеберг, Т., "Масштабное пространство для дискретных сигналов", ПАМИ (12), № 3, март 1990 г., стр. 234–254.
  11. ^ Кэмпбелл, Джей, 2007, Модель SMM как краевая задача с использованием дискретного уравнения диффузии, Theor Popul Biol. 2007 декабрь; 72 (4): 539–46.
  12. ^ Хонархах, М. и Каерс, Дж., 2010 г., Стохастическое моделирование паттернов с использованием дистанционного моделирования паттернов, Математические науки о Земле, 42: 487–517.

внешняя ссылка