matlab实现多元高斯分布概率密度计算
最近要应用多元高斯分布进行异常检测,使用matlab实现。
1.使用matlab自带函数mvnpdf
mvnpdf(data,miu,COV)
在这里,data是需要计算多元高斯分布概率密度的m*n维数组。m是数据的条数,n是每条数据的特征数。miu是属性的平均值组成的向量(是行向量),COV为协方差矩阵。
2.用原始公式
多元高斯分布的概率密度公式如下。使用矩阵计算实现
i=17;%特征个数
m=1;
while m<=30687x=data(m,:)';%带入的数据p(m)=(1/((2*3.1415926)^(i/2)*det(COV)^(1/2)))*(exp(-((x-miu)'*inv(COV)*(x-miu))/2));m=m+1;
end
其中i为特征个数,COV为协方差矩阵,miu为每个属性的平均值组成的向量(是列向量)。此处每一条数据也要变成列向量的形式(x=data(m,:)';)。
inv()用于矩阵求逆,det()用于求方阵的行列式,cov()函数用于求协方差(小写)。
最终两种方式的结果近似相同。
此处要注意的就是平均值向量是行向量还是列向量。如果报错就试试转置一下
matlab实现多元高斯分布概率密度计算相关推荐
- 绘制卡方分布的概率密度函数 matlab,MATLAB如何使用chi2pdf函数计算卡方分布的概率密度...
MATLAB如何使用chi2pdf函数计算卡方分布的概率密度 [语法说明] Y=chi2pdf(X,V):函数计算X中的元素在参数V指定的卡方分布下的概率密度函数值.如果输入参数中有一个为标量,则将其 ...
- matlab卡方分布函数,MATLAB如何使用ncx2pdf函数计算非中心卡方分布的概率密度
MATLAB如何使用ncx2pdf函数计算非中心卡方分布的概率密度 [语法说明] Y=ncx2pdf(X,V,DELTA):计算X中各元素在参数V.DELTA指定下的非中心卡方分布概率密度函数值.X与 ...
- matlab 求概率密度,MATLAB如何使用pdf函数计算指定分布的概率密度函数
MATLAB如何使用pdf函数计算指定分布的概率密度函数 [语法说明] Y=pdf('name',X,A) Y=pdf('name',X,A,B) Y=pdf('name',X,A,B,C) 字符串n ...
- matlab韦伯分布,MATLAB如何使用wblpdf函数计算韦伯分布的概率密度
MATLAB如何使用wblpdf函数计算韦伯分布的概率密度 [语法说明] Y=wblpdf(X,A,B):计算X中的元素在参数A.B指定的韦伯分布下的概率密度函数值.Y 是与 X.A.B 同型的数组, ...
- matlab伽玛分布如何表示,MATLAB如何使用gampdf函数计算伽马分布的概率密度
MATLAB如何使用gampdf函数计算伽马分布的概率密度 [语法说明] Y=gampdf(X,A,B):计算X中的元素在参数A.B指定的伽马分布下的概率密度函数值.X与A.B是同型的数组,如果输入参 ...
- 多元高斯分布及多元条件高斯分布
高斯那些公式 已知 DD 维向量 x\mathbf {x},其高斯概率分布为: N(x|μ,Σ)==1(2π)D/21|Σ|1/2exp(−12(x−μ)TΣ−1(x−μ))1|Σ|(2π)D−−−− ...
- 第二课.多元高斯分布与其几何特征
目录 多元高斯分布 一元高斯与多元高斯 多元高斯分布的参数 二元高斯分布 多元高斯分布的几何特征 几何特征 实例演示 多元高斯分布 一元高斯与多元高斯 在第一课的一元高斯分布中,处理的是一组样本X=( ...
- 多元高斯分布(Multivariate Gaussian Distribution)
from:https://www.jianshu.com/p/d6c8ca915f69 还是对计算机的监测,我们发现CPU负载和占用内存之间,存在正相关关系. CPU负负载增加的时候占用内存也会增加: ...
- 多元线性回归分析matlab实验报告,利用MATLAB进行多元线性回归.ppt
<利用MATLAB进行多元线性回归.ppt>由会员分享,可在线阅读,更多相关<利用MATLAB进行多元线性回归.ppt(15页珍藏版)>请在人人文库网上搜索. 1.2.线性回归 ...
- 多元高斯分布是非参_(二)多元高斯分布与概率图条件独立性假设
Author: Pan Date: 2020/7/15 首先,我们通过随机向量的介绍来引出多元高斯分布,通过多元高斯分布,我们将介绍概率图的条件独立性假设. 先快速过一遍随机向量: 1.随机向量 ...
最新文章
- 市政管网检测机器人收费标准_管道无损CCTV检测:复杂的地下管网工作,其实可以交给机器人来做...
- 独家 | 将时间信息编码用于机器学习模型的三种编码时间信息作为特征的三种方法...
- 云原生时代,CNStack 如何解决企业数字化转型难题?
- pycharm(windows)安装及其设置中文菜单
- remote addr 取到内网ip_内网穿透工具frp
- Java局域网对战游戏、天气预报项目
- CentOS7 基于http服务搭建本地yum仓库
- 【采访】腾讯社交广告高校算法大赛决赛第二周周冠军——nju_newbiew比赛经验及心得分享
- 服务器运行多个jdk版本_如何使用中央管理服务器运行多个查询
- 微信秘笈之--微信多开
- 线程的同步互斥之事件对象(Event)
- 「2020目标一」项目管理PMP
- Nordic-nRF52840环境搭建及学习
- [BFS]愿天下有情人都是失散多年的兄妹
- 【组队学习】【38期】组队学习内容详情!
- 前端的一些实用算法题
- Knockout.js-111
- Node总结及核心模块fs
- python基础课程学习笔记-第二节课
- spider-flow初步使用