二维高斯函数的两种形式

news2024/11/16 2:47:07

第一种形式很常见

多元正态分布

多元正态分布(Multivariate Normal Distribution),也称为多变量正态分布或多维正态分布,是统计学中一种重要的概率分布,用于描述多个随机变量的联合分布。

假设有 n n n 个随机变量 X 1 , X 2 , … , X n X_1, X_2, \ldots, X_n X1,X2,,Xn,它们服从多元正态分布,则这 n n n 个随机变量的联合概率密度函数(PDF)可以表示为:

f ( x 1 , x 2 , … , x n ) = 1 ( 2 π ) n ∣ Σ ∣ exp ⁡ ( − 1 2 ( x − μ ) ⊤ Σ − 1 ( x − μ ) ) f(x_1, x_2, \ldots, x_n) = \frac{1}{\sqrt{(2\pi)^n |\Sigma|}} \exp\left(-\frac{1}{2} (x - \mu)^\top \Sigma^{-1} (x - \mu)\right) f(x1,x2,,xn)=(2π)n∣Σ∣ 1exp(21(xμ)Σ1(xμ))

其中:

  • x = ( x 1 , x 2 , … , x n ) ⊤ x = (x_1, x_2, \ldots, x_n)^\top x=(x1,x2,,xn) 是一个 n × 1 n \times 1 n×1 的列向量,表示 n n n 个随机变量的观测值。
  • μ = ( μ 1 , μ 2 , … , μ n ) ⊤ \mu = (\mu_1, \mu_2, \ldots, \mu_n)^\top μ=(μ1,μ2,,μn) 是一个 n × 1 n \times 1 n×1 的列向量,表示 n n n 个随机变量的均值。
  • Σ \Sigma Σ 是一个 n × n n \times n n×n 的对称正定矩阵,称为协方差矩阵,其元素 σ i j \sigma_{ij} σij 表示 X i X_i Xi X j X_j Xj 之间的协方差。
  • ∣ Σ ∣ |\Sigma| ∣Σ∣ 表示协方差矩阵 Σ \Sigma Σ 的行列式。
  • Σ − 1 \Sigma^{-1} Σ1 是协方差矩阵 Σ \Sigma Σ 的逆矩阵。
  • ( x − μ ) ⊤ (x - \mu)^\top (xμ) ( x − μ ) (x - \mu) (xμ) 的转置。

二维高斯函数(也称为二维正态分布函数)是统计学和物理学中常用的一个函数,用于描述二维空间中的随机变量的分布。这里我们讨论两种形式的二维高斯函数,其中第二种形式通过几何变换(如旋转和平移)从第一种形式推导出来。

二元正态分布

在二元正态分布中,如果已知两个变量 X X X Y Y Y的相关系数 ρ \rho ρ以及它们的标准差 σ X \sigma_X σX σ Y \sigma_Y σY,那么协方差矩阵 Σ \Sigma Σ可以表示为:

Σ = ( σ X 2 ρ σ X σ Y ρ σ X σ Y σ Y 2 ) \Sigma = \begin{pmatrix} \sigma_X^2 & \rho \sigma_X \sigma_Y \\ \rho \sigma_X \sigma_Y & \sigma_Y^2 \end{pmatrix} Σ=(σX2ρσXσYρσXσYσY2)

这里, σ X 2 \sigma_X^2 σX2 σ Y 2 \sigma_Y^2 σY2分别是 X X X Y Y Y的方差,而 ρ σ X σ Y \rho \sigma_X \sigma_Y ρσXσY是它们的协方差。

接下来,需要计算这个协方差矩阵的逆矩阵 Σ − 1 \Sigma^{-1} Σ1。首先,计算协方差矩阵的行列式 ∣ Σ ∣ |\Sigma| ∣Σ∣

∣ Σ ∣ = σ X 2 σ Y 2 − ( ρ σ X σ Y ) 2 = σ X 2 σ Y 2 ( 1 − ρ 2 ) |\Sigma| = \sigma_X^2 \sigma_Y^2 - (\rho \sigma_X \sigma_Y)^2 = \sigma_X^2 \sigma_Y^2 (1 - \rho^2) ∣Σ∣=σX2σY2(ρσXσY)2=σX2σY2(1ρ2)

然后,利用行列式的值和协方差矩阵的元素来构造逆矩阵。逆矩阵 Σ − 1 \Sigma^{-1} Σ1的元素可以表示为:

Σ − 1 = 1 ∣ Σ ∣ ( σ Y 2 − ρ σ X σ Y − ρ σ X σ Y σ X 2 ) \Sigma^{-1} = \frac{1}{|\Sigma|} \begin{pmatrix} \sigma_Y^2 & -\rho \sigma_X \sigma_Y \\ -\rho \sigma_X \sigma_Y & \sigma_X^2 \end{pmatrix} Σ1=∣Σ∣1(σY2ρσXσYρσXσYσX2)

将行列式的值代入上式,得到:

Σ − 1 = 1 σ X 2 σ Y 2 ( 1 − ρ 2 ) ( σ Y 2 − ρ σ X σ Y − ρ σ X σ Y σ X 2 ) \Sigma^{-1} = \frac{1}{\sigma_X^2 \sigma_Y^2 (1 - \rho^2)} \begin{pmatrix} \sigma_Y^2 & -\rho \sigma_X \sigma_Y \\ -\rho \sigma_X \sigma_Y & \sigma_X^2 \end{pmatrix} Σ1=σX2σY2(1ρ2)1(σY2ρσXσYρσXσYσX2)

进一步化简,得到:

Σ − 1 = 1 1 − ρ 2 ( 1 σ X 2 − ρ σ X σ Y − ρ σ X σ Y 1 σ Y 2 ) \Sigma^{-1} = \frac{1}{1 - \rho^2} \begin{pmatrix} \frac{1}{\sigma_X^2} & -\frac{\rho}{\sigma_X \sigma_Y} \\ -\frac{\rho}{\sigma_X \sigma_Y} & \frac{1}{\sigma_Y^2} \end{pmatrix} Σ1=1ρ21(σX21σXσYρσXσYρσY21)

这个逆矩阵代入多元正态分布就是第一种形式。

第一种形式:标准二维高斯函数

标准二维高斯函数(也称为二维正态分布的概率密度函数)的表达式为:

f ( x , y ) = 1 2 π σ x σ y 1 − ρ 2 exp ⁡ ( − 1 2 ( 1 − ρ 2 ) [ ( x − μ x σ x ) 2 − 2 ρ ( x − μ x σ x ) ( y − μ y σ y ) + ( y − μ y σ y ) 2 ] ) f(x, y) = \frac{1}{2\pi\sigma_x\sigma_y\sqrt{1-\rho^2}} \exp\left(-\frac{1}{2(1-\rho^2)}\left[\left(\frac{x-\mu_x}{\sigma_x}\right)^2 - 2\rho\left(\frac{x-\mu_x}{\sigma_x}\right)\left(\frac{y-\mu_y}{\sigma_y}\right) + \left(\frac{y-\mu_y}{\sigma_y}\right)^2\right]\right) f(x,y)=2πσxσy1ρ2 1exp(2(1ρ2)1[(σxxμx)22ρ(σxxμx)(σyyμy)+(σyyμy)2])

其中,

  • μ x \mu_x μx μ y \mu_y μy 分别是 x x x y y y 的均值(位置参数)。
  • σ x \sigma_x σx σ y \sigma_y σy 分别是 x x x y y y 的标准差(尺度参数)。
  • ρ \rho ρ x x x y y y 的相关系数(形状参数)。

第二种出现在Gabor滤波器中,第二种又有多种变形,认清它是高斯函数。

Gabor函数

Gabor滤波器是一种基于Gabor函数的特定频率和方向选择性滤波器,是一种用于图像纹理分析和特征提取的线性滤波器,由Dennis Gabor于1946年提出,广泛应用于图像处理和计算机视觉领域。在空间域中,一个二维Gabor滤波器可以看作是一个正弦平面波和高斯核函数的乘积。

g ( x , y ; λ , θ , ψ , σ , γ ) = exp ⁡ ( − x ′ 2 + γ 2 y ′ 2 2 σ 2 ) cos ⁡ ( 2 π x ′ λ + ψ ) g(x, y; \lambda, \theta, \psi, \sigma, \gamma) = \exp\left(-\frac{x'^2 + \gamma^2 y'^2}{2\sigma^2}\right) \cos\left(2\pi\frac{x'}{\lambda} + \psi\right) g(x,y;λ,θ,ψ,σ,γ)=exp(2σ2x′2+γ2y′2)cos(2πλx+ψ)

其中,

  • ( x , y ) (x, y) (x,y) 是图像中像素的坐标。
  • λ \lambda λ 是正弦函数的波长。
  • θ \theta θ 是Gabor核函数的方向,表示Gabor函数平行条纹的方向。
  • ψ \psi ψ 是相位偏移,通常设置为0或 π 2 \frac{\pi}{2} 2π
  • σ \sigma σ 是高斯函数的标准差,决定了Gabor滤波器的带宽。
  • γ \gamma γ 是空间纵横比,决定了Gabor滤波器的椭圆率。当 γ = 1 \gamma = 1 γ=1时,滤波器形状是圆的;当 γ < 1 \gamma < 1 γ<1时,滤波器形状沿平行于 θ \theta θ的方向拉长。

x ′ x' x y ′ y' y 是通过旋转坐标轴得到的坐标,它们与原始坐标 ( x , y ) (x, y) (x,y)的关系为:

( x ′ y ′ ) = ( cos ⁡ θ sin ⁡ θ − sin ⁡ θ cos ⁡ θ ) ( x y ) \begin{pmatrix} x' \\ y' \end{pmatrix} = \begin{pmatrix} \cos\theta & \sin\theta \\ -\sin\theta & \cos\theta \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix} (xy)=(cosθsinθsinθcosθ)(xy)

这个旋转矩阵用于将Gabor滤波器的方向调整到 θ \theta θ指定的角度。

Gabor滤波器的一个重要特性是它能够捕捉图像中特定频率和方向的信息。通过调整 λ \lambda λ θ \theta θ σ \sigma σ γ \gamma γ等参数,可以生成一系列具有不同频率和方向选择性的Gabor滤波器,从而用于图像的纹理分析和特征提取。

在实际应用中,通常会使用多个Gabor滤波器(具有不同的参数组合)对图像进行滤波,然后分析滤波后的图像以提取有用的特征信息。

第二种形式高斯函数

1

f ( x , y ) = A exp ⁡ ( − x ′ 2 + γ 2 y ′ 2 2 σ 2 ) f(x, y) =A\exp\left(-\frac{x'^2 + \gamma^2 y'^2}{2\sigma^2}\right) f(x,y)=Aexp(2σ2x′2+γ2y′2)是高斯函数,可由标准正态分布通过几何变换推出。

2

也可以写为
f ( x , y ) = A e − 1 2 ( x ˊ 2 σ x 2 + x ˊ 2 σ y 2 ) f(x, y) =A{\rm e} ^{-{ \tfrac {1}{2} \left( \tfrac { \acute{x} ^{2} } {\sigma_x ^{2}} + \tfrac { \acute{x} ^{2} } { \sigma_y ^{2}} \right) }} f(x,y)=Ae21(σx2xˊ2+σy2xˊ2)

σ x \sigma_x σx σ y \sigma_y σy x x x y y y方向上高斯函数的标准差, γ = σ x / σ y \gamma = \sigma_x /\sigma_y γ=σx/σy表示高斯函数的椭圆率。

滤波器在 x x x方向上的“宽度”可以由标准差 σ x \sigma_x σx控制,而在 y y y方向上的“宽度”可以由另一个标准差 σ y \sigma_y σy控制。但在实际应用中,我们并不直接设置这两个标准差,而是通过设置一个主标准差 σ \sigma σ和一个空间纵横比 γ \gamma γ来间接控制。

空间纵横比 γ \gamma γ定义了滤波器在 x x x y y y方向上的相对“宽度”。当 γ = 1 \gamma = 1 γ=1时,滤波器是圆形的;当 γ < 1 \gamma < 1 γ<1时,滤波器在 y y y方向上的“宽度”比 x x x方向小,形成椭圆形,且椭圆的长轴与 x x x轴平行;当 γ > 1 \gamma > 1 γ>1时,情况相反。

在Gabor滤波器的标准定义中,通常只关注一个标准差 σ \sigma σ和一个空间纵横比 γ \gamma γ。这两个参数共同决定了滤波器在不同方向上的形状和带宽。

3

f ( x , y ) f(x, y) f(x,y)写成关于 x x x y y y的函数,这是一个指数项有交叉项的形式。

z = A e − [ a ( x − μ x ) 2 + 2 b ( x − μ x ) ( y − μ y ) + c ( y − μ y ) 2 ] z = A e^{-\left[a(x - \mu_x)^2 + 2b(x - \mu_x)(y - \mu_y) + c(y - \mu_y)^2\right]} z=Ae[a(xμx)2+2b(xμx)(yμy)+c(yμy)2]

其中,

  • a = cos ⁡ 2 θ 2 δ x 2 + sin ⁡ 2 θ 2 δ y 2 a = \frac{\cos^2\theta}{2\delta_x^2} + \frac{\sin^2\theta}{2\delta_y^2} a=2δx2cos2θ+2δy2sin2θ
  • b = − sin ⁡ ( 2 θ ) 4 δ x 2 + sin ⁡ ( 2 θ ) 4 δ y 2 b = -\frac{\sin(2\theta)}{4\delta_x^2} + \frac{\sin(2\theta)}{4\delta_y^2} b=4δx2sin(2θ)+4δy2sin(2θ)
  • c = sin ⁡ 2 θ 2 δ x 2 + cos ⁡ 2 θ 2 δ y 2 c = \frac{\sin^2\theta}{2\delta_x^2} + \frac{\cos^2\theta}{2\delta_y^2} c=2δx2sin2θ+2δy2cos2θ
    在这里插入图片描述
    在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/2108761.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

自己设计的QT系统,留个档

注册登录 主界面展示 天气预报 音乐播放

卷积神经网络与小型全连接网络在MNIST数据集上的对比

卷积神经网络&#xff08;CNN&#xff09; 深度卷积神经网络中&#xff0c;有如下特性 很多层&#xff08;Compositionality&#xff0c;组合性&#xff09;: 深度卷积神经网络通常由多层卷积和非线性激活函数组成。这种多层结构使得网络能够逐步提取和组合低层次的特征&…

shell 学习笔记:数组

目录 1. 定义数组 2. 读取数组元素值 3. 关联数组 4. 在数组前加一个感叹号 ! 可以获取数组的所有键 5. 在数组前加一个井号 # 获取数组的长度 6. 数组初始化的时候&#xff0c;也可以用变量 7. 循环输出数组的方法 7.1 for循环输出 7.2 while循环输出 7.2.1 …

大数据-120 - Flink Window 窗口机制-滑动时间窗口、会话窗口-基于时间驱动基于事件驱动

点一下关注吧&#xff01;&#xff01;&#xff01;非常感谢&#xff01;&#xff01;持续更新&#xff01;&#xff01;&#xff01; 目前已经更新到了&#xff1a; Hadoop&#xff08;已更完&#xff09;HDFS&#xff08;已更完&#xff09;MapReduce&#xff08;已更完&am…

Redis实战宝典:开发规范与最佳实践

目录标题 Key命名设计&#xff1a;可读性、可管理性、简介性Value设计&#xff1a;拒绝大key控制Key的生命周期&#xff1a;设定过期时间时间复杂度为O(n)的命令需要注意N的数量禁用命令&#xff1a;KEYS、FLUSHDB、FLUSHALL等不推荐使用事务删除大key设置合理的内存淘汰策略使…

Java | Leetcode Java题解之第387题字符串中的第一个唯一字符

题目&#xff1a; 题解&#xff1a; class Solution {public int firstUniqChar(String s) {Map<Character, Integer> position new HashMap<Character, Integer>();Queue<Pair> queue new LinkedList<Pair>();int n s.length();for (int i 0; i …

【python因果推断库8】工具变量回归与使用 pymc 验证工具变量1

目录 工具变量回归与使用 pymc 验证工具变量 回归机制与局部平均处理效应 旁白&#xff1a;从多元正态分布中采样 import arviz as az import daft import matplotlib.pyplot as plt import numpy as np import pandas as pd import pymc as pm import scipy from matplotli…

如何阅读PyTorch文档及常见PyTorch错误

如何阅读PyTorch文档及常见PyTorch错误 文章目录 如何阅读PyTorch文档及常见PyTorch错误阅读PyTorch文档示例常见Pytorch错误Tensor在不同设备上维度不匹配cuda内存不足张量类型不匹配 参考 PyTorch文档查看https://pytorch.org/docs/stable/ torch.nn -> 定义神经网络 torc…

红队攻防 | 利用GitLab nday实现帐户接管

在一次红队任务中&#xff0c;目标是一家提供VoIP服务的公司。该目标拥有一些重要的客户&#xff0c;如政府组织&#xff0c;银行和电信提供商。该公司要求外部参与&#xff0c;资产测试范围几乎是公司拥有的每一项互联网资产。 第一天是对目标进行信息收集。这一次&#xff0…

结构开发笔记(七):solidworks软件(六):装配摄像头、摄像头座以及螺丝,完成摄像头结构示意图

若该文为原创文章&#xff0c;转载请注明原文出处 本文章博客地址&#xff1a;https://hpzwl.blog.csdn.net/article/details/141931518 长沙红胖子Qt&#xff08;长沙创微智科&#xff09;博文大全&#xff1a;开发技术集合&#xff08;包含Qt实用技术、树莓派、三维、OpenCV…

成功之路:如何获得机器学习和数据科学实习机会

一年内获得两份实习机会的数据科学家的建议和技巧 欢迎来到雲闪世界。在当今竞争激烈的就业市场中&#xff0c;获得数据科学实习机会可以成为您在科技领域取得成功的门票。 但申请者如此之多&#xff0c;你该如何脱颖而出呢&#xff1f; 无论您是学生、应届毕业生还是想要转行…

IDEA2024.2最新工具下载

​软件使用 1、解压缩包 2、打开如图第三个 3、运行过十来秒等待提示以下信息即可

Ubuntu 无法全局安装 node 包

Anchor: $: cat /etc/lsb* DISTRIB_IDUbuntu DISTRIB_RELEASE22.04 DISTRIB_CODENAMEjammy DISTRIB_DESCRIPTION"Ubuntu 22.04.4 LTS" $: node -v v20.17.0 $: npm -v 10.8.2Question: $: npm install -g docsify-cli结果&#xff1a;超时或者如下图 Answer: 有…

【Python 千题 —— 算法篇】字符串替换

Python 千题持续更新中 …… 脑图地址 &#x1f449;&#xff1a;⭐https://twilight-fanyi.gitee.io/mind-map/Python千题.html⭐ 题目背景 在日常编程中&#xff0c;我们经常会遇到需要对字符串中的特定字符或子串进行替换的需求。比如&#xff0c;替换文本中的敏感词汇、…

html初体验标准标签

<!DOCTYPE html> <html lang"en"> <head><meta charset"UTF-8"><title>Document</title> </head> <body><input type"text"> </body> </html> 内容展示

C和Java实现杨辉三角

C: #include <stdio.h> #define N 15 int main() {int arr[N][N] { 0 }; //初始化int i 0;//行数int j 0;//列数for (i 0; i < N; i){arr[i][0] 1; //每行首元素是 1for (j 0; j < i; j)//为啥j<i呢&#xff1f;因为每一行的个数&#xff08;每一列多少个…

爵士编曲:如何编写爵士钢琴

排列&#xff08;Voicing&#xff09; 由于爵士和声组成音较为复杂&#xff0c;故此衍声排列内容 密集排列&#xff1a;直接堆一起或者左手低音右手和弦音 。 开放排列&#xff1a;各个声部大于等于纯四度&#xff0c;小于八度&#xff0c;符合两只手能弹的情况 混合排列&a…

Python | Leetcode Python题解之第388题文件的最长绝对路径

题目&#xff1a; 题解&#xff1a; class Solution:def lengthLongestPath(self, input: str) -> int:ans, i, n 0, 0, len(input)level [0] * (n 1)while i < n:# 检测当前文件的深度depth 1while i < n and input[i] \t:depth 1i 1# 统计当前文件名的长度l…

深入CSS 布局——WEB开发系列29

CSS 页面布局技术允许我们拾取网页中的元素&#xff0c;并且控制它们相对正常布局流、周边元素、父容器或者主视口/窗口的位置。 一、正常布局流&#xff08;Normal Flow&#xff09; CSS的布局基础是“正常流”&#xff0c;也就是页面元素在没有特别指定布局方式时的默认排列…

文件操作详解:fgetc,fputc,fgets,fputs,fscanf,,fprintf,fread,fwrite的使用和例子 C语言

前言 在日常应用中&#xff0c;我们为了持续的使用一些数据&#xff0c;为了让数据可以在程序退出后可以保存并正常使用&#xff0c;引入了文件的概念和操作。本文分享了一些常用的文件操作函数的使用方法和各自的区别。 一、常用文件顺序读写函数 下面例程所使用的VS工程代码…