方差,协方差及协方差矩阵的计算

1.方差

  • 方差是用来衡量一组数据的离散程度,数序表达式如下:
    σ 2 = 1 N ∑ i = 1 N ( x i − μ ) 2 \sigma^2=\frac1N\sum_{i=1}^N(x_i-\mu)^2 σ2=N1i=1N(xiμ)2

    • σ 2 σ^2 σ2表示样本的总体方差,
    • N N N 表示样本总数,
    • x i x _i xi是第 i 个样本,
    • μ μ μ 是数据集的平均值 ( μ = 1 N ∑ i = 1 N x i ) (\mu=\frac1N\sum_{i=1}^Nx_i) (μ=N1i=1Nxi)
  • 假设某班级有 5 名学生的数学成绩为: 85 , 74 , 63 , 95 , 99 85, 74, 63, 95, 99 85,74,63,95,99, 计算平均值为 μ = 85 + 74 + 63 + 95 + 99 5 = 83.2 μ = \frac{85 +74 + 63 + 95 + 99} 5 =83.2 μ=585+74+63+95+99=83.2 ,方差的计算如下:
    σ 成绩 2 = ( 85 − 83.5 ) 2 + ( 74 − 83.5 ) 2 + ( 63 − 83.5 ) 2 + ( 95 − 83.5 ) 2 + ( 99 − 83.5 ) 2 5 = 177.05 {σ_{成绩}}^2 = \frac{(85-83.5)^2 +(74-83.5)^2 + (63-83.5)^2 + (95-83.5)^2 + (99 - 83.5)^2}{5} = 177.05 σ成绩2=5(8583.5)2+(7483.5)2+(6383.5)2+(9583.5)2+(9983.5)2=177.05

2.协方差

  • 协方差是用来衡量两组样本之间的关系,数序表达式如下:
    σ x σ y = 1 N ∑ i = 1 N ( x i − μ x ) ( y i − μ y ) \sigma_x\sigma_y=\frac{1}{N}\sum_{i=1}^N(x_i-\mu_x)(y_i-\mu_y) σxσy=N1i=1N(xiμx)(yiμy)
    • σ x σ y σ_xσ_y σxσy表示两组样本的总体协方差,
    • N N N 是样本对的总数,
    • x i , y i x _i, y_i xi,yi 是第 i 对样本,
    • μ x , μ y μ_x,μ_y μx,μy分别表示两组样本的均值。
  • 例如,要计算如下 x , y x, y x,y 两个样本的协方差:
xy
500200
600250
451180
    • 计算 x x x的平均值 u x = 500 + 600 + 451 3 = 517 u_x = \frac{500 +600 +451}3 = 517 ux=3500+600+451=517
    • 计算 y y y的平均值 u y = 200 + 250 + 180 3 = 210 u_y = \frac{200 +250 +180}3 = 210 uy=3200+250+180=210
      σ x σ y = ( 500 − 517 ) ( 200 − 210 ) + ( 600 − 517 ) ( 250 − 210 ) + ( 451 − 517 ) ( 180 − 210 ) 3 = 1823.33... σ_xσ_y = \frac{(500-517)(200 - 210) + (600-517)(250-210) + (451-517)(180-210) }{3} = 1823.33... σxσy=3(500517)(200210)+(600517)(250210)+(451517)(180210)=1823.33...

3.协方差矩阵

协方差矩阵就是将多组数据的方差和协方差用矩阵的形式表达出来,例如如下三组样本,其协方差的矩阵排布如下,从下图可以看出,对角线上是各个样本的方差,两边则是样本之间的协方差

xyz
x 1 x_1 x1 y 1 y_1 y1 z 1 z_1 z1
x 2 x_2 x2 y 2 y_2 y2 z 2 z_2 z2
x 3 x_3 x3 y 3 y_3 y3 z 3 z_3 z3

c o v = [ σ x 2 σ x σ y σ x σ z σ y σ x σ y 2 σ y σ z σ z σ x σ z σ y σ x 2 ] cov = \left[\begin{array}{ccc}σ_x^2&σ_xσ_y&σ_xσ_z\\ σ_yσ_x&σ_y^2&σ_yσ_z\\ σ_zσ_x&σ_zσ_y&σ_x^2\end{array}\right] cov= σx2σyσxσzσxσxσyσy2σzσyσxσzσyσzσx2
在这里插入图片描述

4.通过矩阵运算求解一个矩阵的协方差矩阵

  • 将第三节使用的样本写成矩阵的形式,本节主要求解这个矩阵的协方差矩阵
    [ x 1 y 1 z 1 x 2 y 2 z 2 x 3 y 3 z 3 ] \left[\begin{array}{ccc}x_1&y_1&z_1\\ x_2&y_2&z_2\\ x_3&y_3&z_3\end{array}\right] x1x2x3y1y2y3z1z2z3
  • 求出过度矩阵 a a a,根据下面计算可知,过度矩阵 a a a 其实是求解了每个样本的 标准差
    a = [ x 1 y 1 z 1 x 2 y 2 z 2 x 3 y 3 z 3 ] − 1 3 [ 1 1 1 1 1 1 1 1 1 ] [ x 1 y 1 z 1 x 2 y 2 z 2 x 3 y 3 z 3 ] − − − − a = [ x 1 y 1 z 1 x 2 y 2 z 2 x 3 y 3 z 3 ] − 1 3 [ x 1 + x 2 + x 3 y 1 + y 2 + y 3 z 1 + z 2 + z 3 x 1 + x 2 + x 3 y 1 + y 2 + y 3 z 1 + z 2 + z 3 x 1 + x 2 + x 3 y 1 + y 2 + y 3 z 1 + z 2 + z 3 ] − − − − a = [ x 1 − 1 3 ( x 1 + x 2 + x 3 ) y 1 − 1 3 ( y 1 + y 2 + y 3 ) z 1 − 1 3 ( z 1 + z 2 + z 3 ) x 2 − 1 3 ( x 1 + x 2 + x 3 ) y 2 − 1 3 ( y 1 + y 2 + y 3 ) z 2 − 1 3 ( z 1 + z 2 + z 3 ) x 3 − 1 3 ( x 1 + x 2 + x 3 ) y 3 − 1 3 ( y 1 + y 2 + y 3 ) z 3 − 1 3 ( z 1 + z 2 + z 3 ) ] a = \left[\begin{array}{ccc}x_1&y_1&z_1\\x_2&y_2&z_2\\x_3&y_3&z_3\end{array}\right] - \frac{1}{3} \left[\begin{array}{ccc}1&1&1\\1&1&1\\1&1&1\end{array}\right] \left[\begin{array}{ccc}x_1&y_1&z_1\\x_2&y_2&z_2\\x_3&y_3&z_3\end{array}\right] \\ ----\\ a = \left[\begin{array}{ccc}x_1&y_1&z_1\\x_2&y_2&z_2\\x_3&y_3&z_3\end{array}\right] - \frac{1}{3} \left[\begin{array}{ccc}x_1+x_2+x_3&y_1+y_2+y_3&z_1+z_2+z_3\\x_1+x_2+x_3&y_1+y_2+y_3&z_1+z_2+z_3\\x_1+x_2+x_3&y_1+y_2+y_3&z_1+z_2+z_3\end{array}\right] \\ ---- \\ a = \left[\begin{array}{ccc}x_1 - \frac{1}{3} (x_1+x_2+x_3)&y_1 - \frac{1}{3} (y_1+y_2+y_3)&z_1 -\frac{1}{3} (z_1+z_2+z_3) \\x_2 - \frac{1}{3} (x_1+x_2+x_3)&y_2 - \frac{1}{3} (y_1+y_2+y_3)&z_2 -\frac{1}{3} (z_1+z_2+z_3)\\x_3 - \frac{1}{3} (x_1+x_2+x_3)&y_3 - \frac{1}{3} (y_1+y_2+y_3)&z_3 -\frac{1}{3} (z_1+z_2+z_3)\end{array}\right] a= x1x2x3y1y2y3z1z2z3 31 111111111 x1x2x3y1y2y3z1z2z3 a= x1x2x3y1y2y3z1z2z3 31 x1+x2+x3x1+x2+x3x1+x2+x3y1+y2+y3y1+y2+y3y1+y2+y3z1+z2+z3z1+z2+z3z1+z2+z3 a= x131(x1+x2+x3)x231(x1+x2+x3)x331(x1+x2+x3)y131(y1+y2+y3)y231(y1+y2+y3)y331(y1+y2+y3)z131(z1+z2+z3)z231(z1+z2+z3)z331(z1+z2+z3)
  • 最后求解协方差矩阵P
    P = 1 3 ∗ a T ∗ a P = \frac{1}{3} * a^T *a P=31aTa

a T = [ x 1 − 1 3 ( x 1 + x 2 + x 3 ) x 2 − 1 3 ( x 1 + x 2 + x 3 ) x 3 − 1 3 ( x 1 + x 2 + x 3 ) y 1 − 1 3 ( y 1 + y 2 + y 3 ) y 2 − 1 3 ( y 1 + y 2 + y 3 ) y 3 − 1 3 ( y 1 + y 2 + y 3 ) z 1 − 1 3 ( z 1 + z 2 + z 3 ) z 2 − 1 3 ( z 1 + z 2 + z 3 ) z 3 − 1 3 ( z 1 + z 2 + z 3 ) ] a^T = \left[\begin{array}{ccc}x_1 - \frac{1}{3} (x_1+x_2+x_3)&x_2 - \frac{1}{3} (x_1+x_2+x_3)&x_3 - \frac{1}{3} (x_1+x_2+x_3) \\ y_1 - \frac{1}{3} (y_1+y_2+y_3)&y_2 - \frac{1}{3} (y_1+y_2+y_3)&y_3 - \frac{1}{3} (y_1+y_2+y_3)\\ z_1 - \frac{1}{3} (z_1+z_2+z_3)&z_2 -\frac{1}{3} (z_1+z_2+z_3)&z_3 -\frac{1}{3}(z_1+z_2+z_3)\end{array}\right] aT= x131(x1+x2+x3)y131(y1+y2+y3)z131(z1+z2+z3)x231(x1+x2+x3)y231(y1+y2+y3)z231(z1+z2+z3)x331(x1+x2+x3)y331(y1+y2+y3)z331(z1+z2+z3)

P = 1 3 ∗ [ x 1 − 1 3 ( x 1 + x 2 + x 3 ) x 2 − 1 3 ( x 1 + x 2 + x 3 ) x 3 − 1 3 ( x 1 + x 2 + x 3 ) y 1 − 1 3 ( y 1 + y 2 + y 3 ) y 2 − 1 3 ( y 1 + y 2 + y 3 ) y 3 − 1 3 ( y 1 + y 2 + y 3 ) z 1 − 1 3 ( z 1 + z 2 + z 3 ) z 2 − 1 3 ( z 1 + z 2 + z 3 ) z 3 − 1 3 ∗ ( z 1 + z 2 + z 3 ) ] ∗ [ x 1 − 1 3 ( x 1 + x 2 + x 3 ) y 1 − 1 3 ( y 1 + y 2 + y 3 ) z 1 − 1 3 ( z 1 + z 2 + z 3 ) x 2 − 1 3 ( x 1 + x 2 + x 3 ) y 2 − 1 3 ( y 1 + y 2 + y 3 ) z 2 − 1 3 ( z 1 + z 2 + z 3 ) x 3 − 1 3 ( x 1 + x 2 + x 3 ) y 3 − 1 3 ( y 1 + y 2 + y 3 ) z 3 − 1 3 ( z 1 + z 2 + z 3 ) ] P = \frac{1}{3} * \left[\begin{array}{ccc}x_1 - \frac{1}{3} (x_1+x_2+x_3)&x_2 - \frac{1}{3} (x_1+x_2+x_3)&x_3 - \frac{1}{3} (x_1+x_2+x_3) \\ y_1 - \frac{1}{3} (y_1+y_2+y_3)&y_2 - \frac{1}{3} (y_1+y_2+y_3)&y_3 - \frac{1}{3} (y_1+y_2+y_3)\\ z_1 - \frac{1}{3} (z_1+z_2+z_3)&z_2 -\frac{1}{3} (z_1+z_2+z_3)&z_3 -\frac{1}{3} *(z_1+z_2+z_3)\end{array}\right] * \left[\begin{array}{ccc}x_1 - \frac{1}{3} (x_1+x_2+x_3)&y_1 - \frac{1}{3} (y_1+y_2+y_3)&z_1 -\frac{1}{3} (z_1+z_2+z_3) \\x_2 - \frac{1}{3} (x_1+x_2+x_3)&y_2 - \frac{1}{3} (y_1+y_2+y_3)&z_2 -\frac{1}{3} (z_1+z_2+z_3)\\x_3 - \frac{1}{3} (x_1+x_2+x_3)&y_3 - \frac{1}{3} (y_1+y_2+y_3)&z_3 -\frac{1}{3} (z_1+z_2+z_3)\end{array}\right] P=31 x131(x1+x2+x3)y131(y1+y2+y3)z131(z1+z2+z3)x231(x1+x2+x3)y231(y1+y2+y3)z231(z1+z2+z3)x331(x1+x2+x3)y331(y1+y2+y3)z331(z1+z2+z3) x131(x1+x2+x3)x231(x1+x2+x3)x331(x1+x2+x3)y131(y1+y2+y3)y231(y1+y2+y3)y331(y1+y2+y3)z131(z1+z2+z3)z231(z1+z2+z3)z331(z1+z2+z3)

  • 如此运算便可通过矩阵运算求出协方差矩阵

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

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【2025】基于python+django的慢性病健康管理系统(源码、万字文档、图文修改、调试答疑)

系统功能结构图如下 慢性病健康管理系统 课题背景 随着全球人口老龄化的加剧以及生活方式的改变,慢性病的发病率呈上升趋势,给个人健康和社会医疗资源带来了巨大压力。传统的慢性病管理模式存在信息不畅、患者参与度低、医疗资源分配不均等问题&#xf…

2.2 B/S架构和Tomcat服务器

本文介绍了B/S架构、Tomcat服务器及其与IDEA的整合。B/S架构是一种基于浏览器的网络计算模式,具有跨平台、易用性强的特点,适用于互联网应用。Tomcat是Apache开源的Web服务器,支持Java Web应用的部署和运行。文章通过实例演示了如何下载、安装…

QT非UI设计器生成界面的国际化

目的 UI设计器生成界面的国际化,比较容易实现些,因为有现成的函数可以调用,基本过程如下: void MainWindow::on_actLang_CN_triggered() {//中文界面qApp->removeTranslator(trans);delete trans;transnew QTranslator;trans…

Hackme靶机通关攻略

1,打开靶机和kali,在kali终端中扫描靶机ip,得到靶机ip为192.168.50.137 arp-scan -l 2,使用工具扫描出后台目录后访问login.php 3,注册后登陆发现有输入框,可以尝试使用sql注入来得到用户名和密码,密码需要…

国产编辑器EverEdit - 工具栏自定义及认识工具栏上的按钮

1 设置-高级-工具条 1.1 设置说明 1.1.1 工具条自定义 选择主菜单工具 -> 设置 -> 常规,在弹出的选项窗口中选择工具条分类,如下图所示: 左侧窗口是当前支持所有功能按钮列表(上图中居中栏),右侧的窗口是当前显示在工具栏…

docker安装rabbitmq

第一步直接dokce拉取rabbitmq镜像docker 利用docker直接拉取镜像最新版:docker search rabbitmq 运行mq: 需要注意的是-p 5673:5672 解释:-p 外网端口:docker的内部端口 ,你们可以改成自己的外网端口号,我这…

【实战ES】实战 Elasticsearch:快速上手与深度实践-8.2.2成本优化与冷热数据分离

👉 点击关注不迷路 👉 点击关注不迷路 👉 点击关注不迷路 文章大纲 8.2.2AWS OpenSearch Serverless 成本优化与冷热数据分离深度实践1. 成本构成分析与优化机会识别1.1 Serverless模式成本分布1.2 冷热数据特征分析数据特征矩阵 2. 冷热数据…

安卓edge://inspect 和 chrome://inspect调试移动设备上的网页

edge://inspect 和 chrome://inspect 是用于调试浏览器中运行的网页和移动设备上的网页的工具。这两个工具分别属于 Microsoft Edge 和 Google Chrome 浏览器。以下是它们的详细介绍: chrome://inspect 如果是直接使用数据线调试,则只需要执行下面的第一…

checkpoint机制

1、什么是checkpoint 将缓冲池中的脏页刷新到磁盘,并更新redo log的checkpoint位点,确保数据库在发生故障时可以快速恢复到一致的状态。 2、checkpoint执行过程 确保需要刷新的脏页:从缓冲池中选取一部分需要刷新的页数据页刷新&#xff1…

【微服务日志收集①】使用FileBeat+Logstash+ES搭建ELK日志系统

使用FileBeatLogstashES搭建ELK日志系统,架构图如下: 1、 使用docker快速创建ES服务和Kibana服务 前置条件:需要在linux上提前安装好docker和docker-compose 1.1、在linux创建好一个用于存放docker-compose配置文件的文件夹 我的目录是/app/…

Centos 7 安装达梦数据库

一、环境准备 1. 确认操作系统的版本和数据库的版本是否一致 cat /etc/redhat-release 2. 关闭防火墙 查看防火墙状态 firewall-cmd --state 停止firewall systemctl stop firewalld.service 禁止firewall开机启动 systemctl disable firewalld.service 3. 修改文件l…

仿“东方甄选”直播商城小程序运营平台

在公域直播流量红利趋于饱和、流量成本大幅攀升的当下,私域直播为企业开辟了新的流量聚集和转化渠道,特别是对于那些希望在私域流量领域取得突破的品牌商家来说,直播场景以其独特的高频互动氛围,相比其他运营方式,展现…

ZED X系列双目3D相机的耐用性与创新设计解析

在工业自动化和学术研究领域,高精度的视觉设备正成为提升效率和质量的关键。ZED X系列AI立体相机,凭借其先进的技术和耐用的设计,为这一领域带来了新的可能。 核心技术:深度感知与精准追踪 ZED X系列的核心技术之一是Neural Dept…

Cursor的使用感受,帮你使用好自动化编程工具,整理笔记

使用感受 说实话,我觉得cursor还是好用的,可能我刚开始使用,没有使用的非常的熟练,运用也没有非常的透彻,总体体验还是不错的,在使用它时,我优先考虑,前端页面功能复用的时候&#…

《C#上位机开发从门外到门内》3-5:基于FastAPI的Web上位机系统

文章目录 一、项目概述二、系统架构设计三、前后端开发四、数据可视化五、远程控制六、系统安全性与稳定性七、性能优化与测试八、实际应用案例九、结论 随着互联网技术的快速发展,Web上位机系统在工业自动化、智能家居、环境监测等领域的应用日益广泛。基于FastAPI…

vue3单独引用element-plus的Infinite Scroll无限滚动;vue3自定义指令

文章目录 1.正常单独使用element-plus其他功能组件2.引入类似与指令的插件3.自定义指令钩子 1.正常单独使用element-plus其他功能组件 引入即可使用 import { ElSelect, ElOption } from "element-plus"2.引入类似与指令的插件 需要先引入,再注册&…

CMake学习笔记(二):变量设值,源文件/文件查找

一_变量设值: 在上一节中我们知道了如何去链接起来多个源文件并且生成可执行文件&#xff0c;但是当我们的源文件过多的时候会导致我们在add_executable里面写很长的一串&#xff0c;所以我们可以使用变量来进行设值: set(<variable> <value>... [PARENT_SCOPE])…

【Function】Azure Function通过托管身份或访问令牌连接Azure SQL数据库

【Function】Azure Function通过托管身份或访问令牌连接Azure SQL数据库 推荐超级课程: 本地离线DeepSeek AI方案部署实战教程【完全版】Docker快速入门到精通Kubernetes入门到大师通关课AWS云服务快速入门实战目录 【Function】Azure Function通过托管身份或访问令牌连接Azu…

案例5_1:单位数码管显示0

文章目录 文章介绍效果图仿真图5_1放置单位数码管 代码5_1.c 文章介绍 效果图 仿真图5_1 复制案例1_2的仿真图&#xff0c;在此基础上修改 注意&#xff1a;栅格大小需要缩小 放置单位数码管 代码5_1.c #include <reg52.h>#define uchar unsigned char #define uint un…

helm部署metricbeat

背景 在Elastic Stack 7.5版本之前&#xff0c;系统默认采用内置服务进行监控数据采集&#xff08;称为内部收集机制&#xff09;&#xff0c;这种设计存在显著局限性&#xff1a; 当ES集群崩溃时自带的节点监控也会随之崩溃&#xff0c;直到集群恢复前&#xff0c;崩溃期间的…