示例代码-Matlab绘制高斯分布曲面图(1)

 

高斯分布

\quad在数据建模时,经常会用到高斯分布模型[1,2],下面我们就使用Matlab来绘制高斯分布曲面。
 

1. 首先,需要得到 xoyxoy 平面每个交叉点的坐标。这里,我们需要使用 matlab 中的 meshgrid 函数。
2. 然后,根据高斯分布的均值以及协方差矩阵来计算每个点的概率密度。这里,我们需要使用 matlab 中的 mvnpdf 函数。
3. 最后,将计算出的交叉点的概率绘制出来,形成一个曲面。这里,我们需要使用 matlab 中的 surf 函数。
 

代码如下:

clear;
clc;

mu = zeros(1,2);            			%该高斯分布的均值向量
sigma = [0.9 0.5; 0.5 0.9]; 			%该高斯分布的协方差矩阵

% 第1步.
% 使用 meshgrid 函数得到xoy平面中每个交叉点的坐标,横坐标矩阵 X 、和纵坐标矩阵 Y
[X,Y] = meshgrid(-1.5:0.05:1.5); 
[rows,cols] = size(X); 
x = reshape(X,rows*cols,1);
y = reshape(Y,rows*cols,1);
p = [x,y];                  			%将每个点的横坐标矩阵 X 和纵坐标矩阵 Y 合并到一个矩阵中

% 第2步.
% 使用 mvnpdf 函数得到每个点的概率密度
z = mvnpdf(p,mu,sigma);     			%求得网格点上的概率密度
Z = reshape(z,rows,cols);

% 第3步.
% 使用 surf 函数画出高斯分布的曲面
figure
surf(X,Y,Z);                			%画曲面 

 
此段代码的效果图如下:

在这里插入图片描述
 
此外绘制高斯曲面还可以只用 mesh, surfc 和 meshc函数,其中最后两个函数可以显示等高线。可自行运行如下代码观察效果:

surf(X,Y,Z);                %画曲面 
mesh(X,Y,Z);                %画网格

surfc(X,Y,Z);               %画曲面 ,可以显示等高线
meshc(X,Y,Z);               %画网格 ,可以显示等高线

 
示例代码地址: https://github.com/Kai-Xuan/MyNote/blob/master/Matlab/Gaussian/gaussian_surf.m/ [link]

 
参考:


1. Chen K X, Ren J Y, Wu X J, et al. Covariance Descriptors on a Gaussian Manifold and their Application to Image Set Classification[J]. Pattern Recognition, 2020: 107463. [link]
2. https://ww2.mathworks.cn/help/stats/multivariate-normal-distribution-1.html [link]

版权声明:本文为u013515929原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/u013515929/article/details/108246728

智能推荐

CentOS学习之路1-wget下载安装配置

参考1: https://blog.csdn.net/zhaoyanjun6/article/details/79108129 参考2: http://www.souvc.com/?p=1569 CentOS学习之路1-wget下载安装配置 1.wget的安装与基本使用 安装wget yum 安装软件 默认安装保存在/var/cache/yum ,用于所有用户使用。 帮助命令 基本用法 例子:下载...

深入浅出Spring的IOC容器,对Spring的IOC容器源码进行深入理解

文章目录 DispatcherServlet整体继承图 入口:DispatcherServlet.init() HttpServletBean.init() FrameworkServlet.initServletBean() 首先大家,去看Spring的源码入口,第一个就是DispatcherServlet DispatcherServlet整体继承图 入口:DispatcherServlet....

laravel框架的课堂知识点概总

1. MVC 1.1 概念理解 MVC全名是Model View Controller,是模型(model)-视图(view)-控制器(controller)的缩写,一种软件设计典范,用一种业务逻辑、数据、界面显示分离的方法组织代码,将业务逻辑聚集到一个部件里面,在改进和个性化定制界面及用户交互的同时,不需要重新编写业务逻辑 MVC 是一种使用 MVC(Model View Controller ...

Unity人物角色动画系统学习总结

使用动画系统控制人物行走、转向、翻墙、滑行、拾取木头 混合树用来混合多个动画 MatchTarget用来匹配翻墙贴合墙上的某一点,人物以此为支点翻墙跳跃 IK动画类似于MatchTarget,控制两只手上的两个点来指定手的旋转和位置,使得拾取木头时更逼真 创建AnimatorController: 首先创建一个混合树,然后双击 可以看到该混合树有五种状态机,分别是Idle、WalkForward、...

Composer 安装 ThinkPHP6 问题

Composer 安装 ThinkPHP6 问题 先说说问题 一.运行环境要求 二.配置 参考: ThinkPHP6.0完全开发手册 先说说问题 执行ThinkPHP6的安装命令 遇到问题汇总如下: 看提示是要更新版本,执行命令更新。 更新之后,再次安装ThinkPHP,之后遇到如下问题。 尝试了很多方法,依然不能解决。其中包括使用https://packagist.phpcomposer.com...

猜你喜欢

Spring Boot 整合JDBC

今天主要讲解一下SpringBoot如何整合JDBC,没啥理论好说的,直接上代码,看项目整体结构 看一下对应的pom.xml 定义User.java 定义数据源配置,这里使用druid,所以需要写一个配置类 上面指定druid的属性配置,和用户登录的账号信息以及对应的过滤规则: 下面定义数据访问接口和对应的实现: 数据访问层很简单,直接注入JdbcTemplate模板即可,下面再看对应的servi...

html鼠标悬停显示样式

1.显示小手:     在style中添加cursor:pointer 实现鼠标悬停变成小手样式     实例:         其他参数: cursor语法: cursor : auto | crosshair | default | hand | move | help | wait | tex...

Yupoo(又拍网)的系统架构

Yupoo!(又拍网) 是目前国内最大的图片服务提供商,整个网站构建于大量的开源软件之上。以下为其使用到的开源软件信息: 操作系统:CentOS、MacOSX、Ubuntu 服务器:Apache、Nginx、Squid 数据库:MySQLmochiweb、MySQLdb 服务器监控:Cacti、Nagios、 开发语言:PHP、Python、Erlang、Java、Lua 分布式计算:Hadoop...

创建一个Servlet项目流程(入门)

版本 IDEA 2020.2 JDK1.8 apache-tomcat-9.0.36 项目流程 一、IDEA中新建JaveEE项目 项目起名,选择项目存放地址,点击finish创建成功 进入项目后,右键选择项目,选择add Framework Support 选择Web Application,点击OK 此时项目文件夹 在WEB-INF下创建两个目录classes和lib 按ctrl+alt+sh...

Docker部署SpringCloud ELK+RabbitMQ日志

Docker部署SpringCloud ELK+RabbitMQ日志  Im_Coder 原文:https://www.jianshu.com/p/f773f23096a9 一、效果图 image.png 二、ELK是什么? ELK由ElasticSearch、Logstash和Kiabana三个开源工具组成。 其中Elasticsearch是个开源分布式搜索引擎,它的特点有:分布式,索...