Функция Гаусса - Gaussian function
Эта статья нужны дополнительные цитаты для проверка.Август 2009 г.) (Узнайте, как и когда удалить этот шаблон сообщения) ( |
В математика, а Функция Гаусса, часто называемый просто Гауссовский, это функция формы
для произвольных настоящий константы а, б и ненулевой c. Назван в честь математика. Карл Фридрих Гаусс. В график гауссиана является характеристикой симметричной "кривая колокола "форма. Параметр а - высота пика кривой, б это положение центра пика и c (в стандартное отклонение, иногда называемый гауссовой RMS width) контролирует ширину «колокола».
Гауссовские функции часто используются для представления функция плотности вероятности из нормально распределенный случайная переменная с ожидаемое значение μ = б и отклонение σ2 = c2. В этом случае гауссиан имеет вид:
Функции Гаусса широко используются в статистика описать нормальные распределения, в обработка сигналов определять Гауссовы фильтры, в обработка изображений где двумерные гауссианы используются для Размытие по Гауссу, а в математике решить уравнения теплопроводности и уравнения диффузии и определить Преобразование Вейерштрасса.
Характеристики
Гауссовы функции возникают при составлении экспоненциальная функция с вогнутый квадратичная функция:
куда:
Таким образом, гауссовские функции - это те функции, логарифм - вогнутая квадратичная функция.
Параметр c относится к полная ширина на половине максимальной (FWHM) пика согласно
Затем функция может быть выражена через FWHM, представленную ш:
В качестве альтернативы параметр c можно интерпретировать, сказав, что два точки перегиба функции происходят в Икс = б − c и Икс = б + c.
В полная ширина в десятых долях от максимума (FWTM) для гауссиана может представлять интерес и является
Гауссовы функции аналитический, и их предел в качестве Икс → ∞ равно 0 (для указанного выше случая б = 0).
Гауссовы функции относятся к числу тех функций, которые элементарный но не хватает элементарного первообразные; в интеграл функции Гаусса является функция ошибки. Тем не менее их несобственные интегралы по всей действительной прямой можно вычислить точно, используя Гауссовский интеграл
и получается
Этот интеграл равен 1 тогда и только тогда, когда (в нормализующая константа ), и в этом случае гауссовский функция плотности вероятности из нормально распределенный случайная переменная с ожидаемое значение μ = б и отклонение σ2 = c2:
Эти гауссианы показаны на прилагаемом рисунке.
Гауссовские функции с центром в нуле минимизируют Фурье принцип неопределенности.
Произведение двух гауссовых функций является гауссовой, а свертка двух гауссовских функций также является гауссовой, причем дисперсия представляет собой сумму исходных дисперсий: . Однако произведение двух гауссовых функций плотности вероятности (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]
- Шум в измеряемом профиле либо i.i.d. Гауссов, или шум Распределенный по Пуассону.
- Расстояние между каждой выборкой (т. Е. Расстояние между пикселями, измеряющими данные) одинаково.
- Пик является «хорошо отобранным», так что менее 10% площади или объема под пиком (площадь, если гауссиан 1D, объем, если гауссиан 2D), лежит за пределами области измерения.
- Ширина пика намного больше, чем расстояние между точками выборки (т.е. пиксели детектора должны быть по крайней мере в 5 раз меньше, чем гауссова FWHM).
Когда эти предположения выполняются, следующие ковариационная матрица K применяется для параметров профиля 1D , , и под i.i.d. Гауссов шум и под шумом Пуассона:[8]
куда ширина пикселей, используемых для выборки функции, - квантовая эффективность детектора, а обозначает стандартное отклонение шума измерения. Таким образом, индивидуальные отклонения параметров в случае гауссовского шума равны
а в случае пуассоновского шума
Для параметров 2D профиля, дающих амплитуду , позиция , и ширина профиля применяются следующие ковариационные матрицы:[9]
где дисперсии отдельных параметров даны диагональными элементами ковариационной матрицы.
Дискретный гауссовский
Можно попросить дискретный аналог гауссианы; это необходимо в дискретных приложениях, в частности цифровая обработка сигналов. Простой ответ - выбрать непрерывный гауссиан, получив выбранное ядро Гаусса. Однако эта дискретная функция не имеет дискретных аналогов свойств непрерывной функции и может приводить к нежелательным эффектам, как описано в статье. реализация масштабного пространства.
Альтернативный подход - использовать дискретное гауссово ядро:[10]
куда обозначает модифицированные функции Бесселя целочисленного порядка.
Это дискретный аналог непрерывного гауссиана в том смысле, что он является решением дискретной уравнение диффузии (дискретное пространство, непрерывное время), точно так же, как непрерывный гауссиан является решением уравнения непрерывной диффузии.[11]
Приложения
Гауссовские функции появляются во многих контекстах в естественные науки, то социальные науки, математика, и инженерное дело. Вот некоторые примеры:
- В статистика и теория вероятности, Функции Гаусса появляются как функция плотности нормальное распределение, что является ограничивающим распределение вероятностей сложных сумм, согласно Центральная предельная теорема.
- Гауссовы функции - это Функция Грина для (однородного и изотропного) уравнение диффузии (и к уравнение теплопроводности, что одно и то же), a уравнение в частных производных который описывает эволюцию плотности массы во времени при распространение. В частности, если массовая плотность во время т= 0 задается Дельта Дирака, что по сути означает, что сначала масса сосредоточена в одной точке, а затем распределение массы во времени т будет задана функцией Гаусса с параметром а линейно связаны с 1 /√т и c линейно связаны с √т; этот изменяющийся во времени гауссиан описывается тепловое ядро. В более общем смысле, если начальная массовая плотность равна φ (Икс), то плотность массы на более поздних временах получается путем свертка функции φ с функцией Гаусса. Свертка функции с гауссианом также известна как Преобразование Вейерштрасса.
- Функция Гаусса - это волновая функция из основное состояние из квантовый гармонический осциллятор.
- В молекулярные орбитали используется в вычислительная химия возможно линейные комбинации функций Гаусса, называемых Гауссовские орбитали (смотрите также базисный набор (химия) ).
- Математически производные функции Гаусса можно представить с помощью Функции Эрмита. В п-я производная от функции Гаусса - это сама функция Гаусса, умноженная на п-й Многочлен Эрмита, до масштаба.
- Следовательно, гауссовские функции также связаны с состояние вакуума в квантовая теория поля.
- Гауссовы пучки используются в оптических системах, микроволновых системах и лазерах.
- В масштабное пространство представления, функции Гаусса используются как сглаживающие ядра для генерации многомасштабных представлений в компьютерное зрение и обработка изображений. В частности, производные от гауссианов (Функции Эрмита ) используются как основа для определения большого количества типов визуальных операций.
- Функции Гаусса используются для определения некоторых типов искусственные нейронные сети.
- В флуоресцентная микроскопия 2D-функция Гаусса используется для аппроксимации Диск Эйри, описывающий распределение интенсивности, создаваемое точечный источник.
- В обработка сигналов они служат для определения Гауссовы фильтры, например, в обработка изображений где 2D гауссианы используются для Размытие по Гауссу. В цифровая обработка сигналов, используется дискретное гауссово ядро, который может быть определен путем выборки по Гауссу или другим способом.
- В геостатистика они использовались для понимания различий между паттернами сложного тренировочный образ. Они используются с методами ядра для кластеризации шаблонов в пространстве функций.[12]
Смотрите также
Рекомендации
- ^ Сквайрс, Г. Л. (30 августа 2001 г.). Практическая физика (4-е изд.). Издательство Кембриджского университета. Дои:10.1017 / cbo9781139164498. ISBN 978-0-521-77940-1.
- ^ Вайсштейн, Эрик В. «Преобразование Фурье - Гауссово». MathWorld. Получено 19 декабря 2013.
- ^ Наури, Николай. "Berechnung von Kovarianzellipsen" (PDF). Получено 14 августа 2019.
- ^ Родитель А., М. Морин и П. Лавин. «Распространение супергауссовских распределений поля». Оптическая и квантовая электроника 24.9 (1992): S1071-S1079.
- ^ "Руководство по командам оптического программного обеспечения GLAD, Запись по команде GAUSSIAN" (PDF). Прикладные оптические исследования. 2016-12-15.
- ^ Каруана, Ричард А .; Searle, Roger B .; Хеллер, Томас .; Шупак, Саул И. (1986). «Быстрый алгоритм разрешения спектров». Аналитическая химия. Американское химическое общество (ACS). 58 (6): 1162–1167. Дои:10.1021 / ac00297a041. ISSN 0003-2700.
- ^ а б Хунвэй Го, "Простой алгоритм подбора функции Гаусса", IEEE Sign. Proc. Mag. 28 (9): 134-137 (2011).
- ^ а б Н. Хаген, М. Купинский и Э. Л. Дереняк, "Оценка гауссовского профиля в одном измерении", Прил. Опт. 46: 5374–5383 (2007).
- ^ а б Н. Хаген и Е. Л. Дерениак, "Оценка гауссовского профиля в двух измерениях", Прил. Опт. 47: 6842–6851 (2008).
- ^ Линдеберг, Т., "Масштабное пространство для дискретных сигналов", ПАМИ (12), № 3, март 1990 г., стр. 234–254.
- ^ Кэмпбелл, Джей, 2007, Модель SMM как краевая задача с использованием дискретного уравнения диффузии, Theor Popul Biol. 2007 декабрь; 72 (4): 539–46.
- ^ Хонархах, М. и Каерс, Дж., 2010 г., Стохастическое моделирование паттернов с использованием дистанционного моделирования паттернов, Математические науки о Земле, 42: 487–517.
внешняя ссылка
- Mathworld, включает доказательство связи между c и FWHM
- «Интеграция колоколообразной кривой». MathPages.com.
- Реализация распределения Гаусса на Haskell, Erlang и Perl
- Бенсимхаун Майкл, N-Кумулятивная функция измерений и другие полезные факты о гауссианах и нормальной плотности (2009 г.)
- Код для подгонки гауссиан в ImageJ и Фиджи.