深入理解Allan方差:用体重数据分析误差的时间尺度与稳定性
文章目录
- 1. 什么是Allan方差?
- Allan方差的特点
- 2. Allan方差与传统方差的区别
- 3. 用体重数据举例分析波动性
- 场景A:体重变化较平稳
- 场景B:体重变化波动较大
- 4. Allan方差的计算公式与详细步骤
- 5. 不同时间块长度下的Allan方差计算
- 场景A的Allan方差计算(块长度为1天)
- 场景B的Allan方差计算(块长度为1天)
- 6. Allan偏差与块长度的双对数曲线
- 双对数曲线的解读
- 7. 总结
在日常生活和科学研究中,体重、温度、金融市场波动等数据往往包含复杂的噪声和趋势成分。传统的统计方法(如均值和标准差)在揭示数据动态变化方面存在局限。Allan方差(Allan Variance)作为一种分析时间序列数据稳定性的方法,能够帮助我们深入探究数据在不同时间尺度下的波动特性。本文将详细介绍Allan方差的概念及其计算过程,并通过实际的体重数据示例分析其短期和长期波动性。我们将通过双对数曲线展示Allan方差在不同块长度下的变化,从而全面揭示体重数据的动态特性。
1. 什么是Allan方差?
Allan方差最初被应用于分析频率标准中的噪声特性,特别适用于分析时间序列的长期和短期稳定性。它能够揭示数据在不同时间尺度下的波动特性,从而帮助我们理解数据的动态变化。
Allan方差的特点
与传统的方差不同,Allan方差能够根据不同时间尺度(块长度)分析数据的波动。这意味着我们可以通过调整块长度,分别观察短期和长期的波动特性。Allan方差的这种特性使其在分析数据的稳定性时具有优势。
2. Allan方差与传统方差的区别
传统方差衡量的是数据整体的波动情况,无法分辨短期波动与长期趋势。而Allan方差则通过计算不同块长度下的数据差异,揭示了不同时间尺度上的稳定性。这一特性使得Allan方差在体重管理、温度分析、金融市场波动等领域表现得尤为出色。
3. 用体重数据举例分析波动性
为了更好地理解Allan方差的计算,我们使用一周的体重数据进行分析,假设每天早晨记录体重,并分为两个场景:
场景A:体重变化较平稳
| 天数 | 体重(kg) |
|---|---|
| 第1天 | 70.0 |
| 第2天 | 70.2 |
| 第3天 | 70.1 |
| 第4天 | 70.3 |
| 第5天 | 70.4 |
| 第6天 | 70.3 |
| 第7天 | 70.5 |
场景B:体重变化波动较大
| 天数 | 体重(kg) |
|---|---|
| 第1天 | 70.0 |
| 第2天 | 70.6 |
| 第3天 | 70.1 |
| 第4天 | 70.5 |
| 第5天 | 70.0 |
| 第6天 | 70.6 |
| 第7天 | 70.1 |
场景A中体重变化较为平稳,而场景B中则显示出更显著的波动。
4. Allan方差的计算公式与详细步骤
Allan方差的计算公式如下:
σ y 2 ( τ ) = 1 2 ( N − 1 ) ∑ i = 1 N − 1 ( y ˉ i + 1 − y ˉ i ) 2 \sigma_y^2(\tau) = \frac{1}{2(N - 1)} \sum_{i=1}^{N-1} (\bar{y}_{i+1} - \bar{y}_i)^2 σy2(τ)=2(N−1)1i=1∑N−1(yˉi+1−yˉi)2
其中:
- σ y 2 ( τ ) \sigma_y^2(\tau) σy2(τ) 表示块长度为 τ \tau τ 的Allan方差;
- τ \tau τ 表示时间块长度(如1天、3天等);
- N N N 为数据块的总数;
- y ˉ i \bar{y}_i yˉi 表示第 i i i 个块的平均值。
计算步骤:
- 选择时间块长度 τ \tau τ。
- 按块长度将数据分组,计算每组的平均值。
- 计算相邻时间块的均值差的平方。
- 求和并除以 2 ( N − 1 ) 2(N-1) 2(N−1)。
5. 不同时间块长度下的Allan方差计算
通过选择不同的时间块长度(如1天、3天、7天等),可以分析不同时间尺度下的波动性。例如,当块长度为1天时,主要分析短期波动;块长度为7天时,观察更长时间尺度的趋势变化。
以下为基于场景A和场景B数据的详细计算:
场景A的Allan方差计算(块长度为1天)
-
计算每个时间块的平均值:
- y ˉ 1 = 70.0 \bar{y}_1 = 70.0 yˉ1=70.0, y ˉ 2 = 70.2 \bar{y}_2 = 70.2 yˉ2=70.2, y ˉ 3 = 70.1 \bar{y}_3 = 70.1 yˉ3=70.1, y ˉ 4 = 70.3 \bar{y}_4 = 70.3 yˉ4=70.3, y ˉ 5 = 70.4 \bar{y}_5 = 70.4 yˉ5=70.4, y ˉ 6 = 70.3 \bar{y}_6 = 70.3 yˉ6=70.3, y ˉ 7 = 70.5 \bar{y}_7 = 70.5 yˉ7=70.5
-
计算相邻时间块的均值差的平方:
- ( y ˉ 2 − y ˉ 1 ) 2 = 0.04 (\bar{y}_2 - \bar{y}_1)^2 = 0.04 (yˉ2−yˉ1)2=0.04, ( y ˉ 3 − y ˉ 2 ) 2 = 0.01 (\bar{y}_3 - \bar{y}_2)^2 = 0.01 (yˉ3−yˉ2)2=0.01, ( y ˉ 4 − y ˉ 3 ) 2 = 0.04 (\bar{y}_4 - \bar{y}_3)^2 = 0.04 (yˉ4−yˉ3)2=0.04, ( y ˉ 5 − y ˉ 4 ) 2 = 0.01 (\bar{y}_5 - \bar{y}_4)^2 = 0.01 (yˉ5−yˉ4)2=0.01, ( y ˉ 6 − y ˉ 5 ) 2 = 0.01 (\bar{y}_6 - \bar{y}_5)^2 = 0.01 (yˉ6−yˉ5)2=0.01, ( y ˉ 7 − y ˉ 6 ) 2 = 0.04 (\bar{y}_7 - \bar{y}_6)^2 = 0.04 (yˉ7−yˉ6)2=0.04
-
计算Allan方差:
σ y 2 ( τ ) = 1 2 × 6 ( 0.04 + 0.01 + 0.04 + 0.01 + 0.01 + 0.04 ) = 0.15 12 = 0.0125 \sigma_y^2(\tau) = \frac{1}{2 \times 6} (0.04 + 0.01 + 0.04 + 0.01 + 0.01 + 0.04) = \frac{0.15}{12} = 0.0125 σy2(τ)=2×61(0.04+0.01+0.04+0.01+0.01+0.04)=120.15=0.0125 -
计算Allan偏差:
σ y ( τ ) = 0.0125 ≈ 0.1118 kg \sigma_y(\tau) = \sqrt{0.0125} \approx 0.1118 \text{ kg} σy(τ)=0.0125≈0.1118 kg
场景B的Allan方差计算(块长度为1天)
同样的方法计算得出:
- Allan方差为 0.1358 0.1358 0.1358;
- Allan偏差为 σ y ( τ ) ≈ 0.3687 kg \sigma_y(\tau) \approx 0.3687 \text{ kg} σy(τ)≈0.3687 kg。
6. Allan偏差与块长度的双对数曲线
为了直观展示波动性,我们通常绘制Allan偏差与块长度的双对数曲线。假设我们分别计算了1天、2天、3天等不同块长度的Allan方差,并开方得到Allan偏差。
双对数曲线的解读
- 负斜率:块长度增加时,Allan偏差下降,表示短期波动显著,随着时间延长趋于平稳。
- 正斜率:块长度增加时,Allan偏差上升,说明数据可能存在长期趋势或周期性变化。
- 水平段:波动性变化较小,数据的波动相对随机且无显著趋势。
7. 总结
Allan方差能够帮助我们揭示数据在不同时间尺度下的波动特性。通过计算和分析体重数据的Allan方差,我们展示了其在实际应用中的效果。进一步,通过双对数曲线的绘制,我们能更深入理解和应用Allan方差,在更广泛的数据分析中发挥作用。
应用实例:
- 体重管理:Allan方差可以帮助判断体重变化的稳定性,识别出短期波动与长期趋势。
- 金融数据分析:通过不同时间尺度的Allan方差,可以分析市场的短期波动性和长期趋势。
- 温度监测:Allan方差可用于分析温度数据的波动特性,识别季节性变化或长期趋势。
Allan方差不仅在科研中具有重要的应用价值,在日常生活中同样能提供有力的数据分析支持。通过深入理解和应用Allan方差,我们可以更精准地分析和解释数据的动态变化,提升数据分析的深度与广度。
相关文章:
深入理解Allan方差:用体重数据分析误差的时间尺度与稳定性
文章目录 1. 什么是Allan方差?Allan方差的特点 2. Allan方差与传统方差的区别3. 用体重数据举例分析波动性场景A:体重变化较平稳场景B:体重变化波动较大 4. Allan方差的计算公式与详细步骤5. 不同时间块长度下的Allan方差计算场景A的Allan方差…...
Elasticsearch 解析:倒排索引机制/字段类型/语法/常见问题
Elasticsearch 是一个分布式的开源搜索引擎,广泛用于全文搜索、分析和数据存储。它基于 Apache Lucene 构建,支持 RESTful 风格的 API,使得开发者能够高效地存储和检索数据。本文将详细讲解 Elasticsearch 的基本原理,特别是其倒排…...
数字后端零基础入门系列 | Innovus零基础LAB学习Day5
###Module 12 RC参数提取和时序分析 数字后端零基础入门系列 | Innovus零基础LAB学习Day4 数字后端零基础入门系列 | Innovus零基础LAB学习Day3 数字后端零基础入门系列 | Innovus零基础LAB学习Day2 数字后端零基础入门系列 | Innovus零基础LAB学习Day1 ###LAB12-1 这个章节…...
Redis 内存回收策略小结
Redis 内存回收策略 及时回收内存中不需要的数据,能有效地保持性能和防止内存溢出。Redis内存回收主要有两种场景 删除过期的键值对内存使用达到maxmemory时触发回收策略 删除过期的键值对 惰性删除: 在查询时如果发现 该键值对已经过期则执行删除操作…...
React常用前端框架合集
React 是 Facebook 开发的一款用于构建用户界面的 JavaScript 库。由于其高效、灵活的特性,React 成为了目前最流行的前端框架之一。为了帮助开发者更好地利用 React 构建应用,市场上涌现了许多优秀的辅助工具和框架。本文将详细介绍几个常用的 React 前…...
python对文件的读写操作
任务:读取文件夹下的批量txt数据,并将其写入到对应的word文档中。 txt文件中包含:编号、报告内容和表格数据。写入到word当中:编号、报告内容、表格数据、人格雷达图以及对应的详细说明(详细说明是根据表格中的标识那一列中的加号…...
Redis工具类(解决缓存穿透、缓存击穿)
文章目录 前言IBloomFilterObjectMapUtilsCacheClient使用示例具体业务的布隆过滤器控制层服务层 前言 该工具类包含以下功能: 1.将任意对象存储在 hash 类型的 key 中,并可以设置 TTL 2.将任意对象存储在 hash 类型的 key 中,并且可以设置…...
Air780E量产binpkg文件的获取方法
Air780E量产binpkg文件如何获取呢?操作方法如下。 一、背景 最近luatos开发客户增多,客户在量产烧录的时候需要binpkg文件,但是有些客户不知道binpkg文件是什么,在哪里获取,是否可以用soc文件提取出来,使…...
C++STL之stack
1.stack的使用 函数说明 接口说明 stack() 构造空的栈 empty() 检测 stack 是否为空 size() 返回 stack 中元素的个数 top() 返回栈顶元素的引用 push() 将元素 val 压入 stack 中 pop() 将 stack 中尾部的元素弹出 2.stack的模拟实现 #include<vector> namespace abc { …...
git的学习之远程进行操作
1.代码托管GitHub:充当中央服务器仓库的角色 2.git远程进行操作 3.配置本地服务器的公钥 4.推送 5.git远程操作 pull .gitignore 6.给命令配置别名 git config --global alias.st status 7.标签管理 git tag -a [name] -m "XXX" [commit_id] 操作标签…...
蓝桥杯普及题
[蓝桥杯 2024 省 B] 好数 题目描述 一个整数如果按从低位到高位的顺序,奇数位(个位、百位、万位……)上的数字是奇数,偶数位(十位、千位、十万位……)上的数字是偶数,我们就称之为“好数”。 给定一个正整数 N N N,请计算从 1 1...
Spreadsheet导出excel
记录下常用的方法 数字转字符:Coordinate::stringFromColumnIndex(27); 输出 AA字符转数字:Coordinate::columnIndexFromString(AA); 输出27设置单元格式 eg:(设置为保留两位小数点) $sheet->getStyle($columnLetter)->getNumberFormat()->set…...
Leetcode|454.四数相加II ● 383. 赎金信 ● 15. 三数之和 ● 18. 四数之和
15.三数之和 哈希解法: 用俩个for循环求出,所需的a和b,再用哈希表,判断剩余的那个c是否在数组 class Solution { public:vector<vector<int>> threeSum(vector<int>& nums) {vector<vector<int>…...
使用ceph-csi把ceph-fs做为k8s的storageclass使用
背景 ceph三节点集群除了做为对象存储使用,计划使用cephfs替代掉k8s里面现有的nfs-storageclass。 思路 整体实现参考ceph官方的ceph csi实现,这套环境是arm架构的,即ceph和k8s都是在arm上实现。实测下来也兼容。 ceph-fs有两种两种挂载方…...
太速科技-212-RCP-601 CPCI刀片计算机
RCP-601 CPCI刀片计算机 一、产品简介 RCP-601是一款基于Intel i7双核四线程的高性能CPCI刀片式计算机,同时,将CPCI产品的欧卡结构及其可靠性、可维护性、可管理性与计算机的抗振动、抗冲击、抗宽温环境急剧变化等恶劣环境特性进行融合。产品特别…...
【解决 Windows 下 SSH “Bad owner or permissions“ 错误及端口转发问题详解】
使用 Windows 连接远程服务器出现 Bad owner or permissions 错误及解决方案 在 Windows 系统上连接远程服务器时,使用 SSH 可能会遇到以下错误: Bad owner or permissions on C:\Users\username/.ssh/config这个问题通常是由于 SSH 配置文件 .ssh/con…...
使用预训练的BERT进行金融领域问答
获取更多完整项目代码数据集,点此加入免费社区群 : 首页-置顶必看 1. 项目简介 本项目旨在开发并优化一个基于预训练BERT模型的问答系统,专注于金融领域的应用。随着金融市场信息复杂性和规模的增加,传统的信息检索方法难以高效…...
ReactOS系统中MM_REGION结构体的声明
ReactOS系统中MM_REGION结构体的声明 ReactOS系统中MM_REGION结构体的声明 文章目录 ReactOS系统中MM_REGION结构体的声明MM_REGION MM_REGION typedef struct _MM_REGION {ULONG Type;//MEM_COMMIT,MEM_RESERVEULONG Protect;//PAGE_READONLYY,PAGE_READ_WRITEULONG Length;…...
web相关知识学习笔记
一, web安全属于网络信息安全的一个分支,www即全球广域网,也叫万维网,是一个分布式图形信息系统 二, 1.①安全领域,通常将用户端(浏览器端)称为前端,服务器端称为后端 ②…...
App测试环境部署
一.JDK安装 参考以下AndroidDevTools - Android开发工具 Android SDK下载 Android Studio下载 Gradle下载 SDK Tools下载 二.SDK安装 安装地址:https://www.androiddevtools.cn/ 解压 环境变量配置 变量名:ANDROID_SDK_HOME 参考步骤: A…...
Oracle查询表空间大小
1 查询数据库中所有的表空间以及表空间所占空间的大小 SELECTtablespace_name,sum( bytes ) / 1024 / 1024 FROMdba_data_files GROUP BYtablespace_name; 2 Oracle查询表空间大小及每个表所占空间的大小 SELECTtablespace_name,file_id,file_name,round( bytes / ( 1024 …...
QMC5883L的驱动
简介 本篇文章的代码已经上传到了github上面,开源代码 作为一个电子罗盘模块,我们可以通过I2C从中获取偏航角yaw,相对于六轴陀螺仪的yaw,qmc5883l几乎不会零飘并且成本较低。 参考资料 QMC5883L磁场传感器驱动 QMC5883L磁力计…...
Cesium1.95中高性能加载1500个点
一、基本方式: 图标使用.png比.svg性能要好 <template><div id"cesiumContainer"></div><div class"toolbar"><button id"resetButton">重新生成点</button><span id"countDisplay&qu…...
汽车生产虚拟实训中的技能提升与生产优化
在制造业蓬勃发展的大背景下,虚拟教学实训宛如一颗璀璨的新星,正发挥着不可或缺且日益凸显的关键作用,源源不断地为企业的稳健前行与创新发展注入磅礴强大的动力。就以汽车制造企业这一极具代表性的行业主体为例,汽车生产线上各类…...
Web 架构之 CDN 加速原理与落地实践
文章目录 一、思维导图二、正文内容(一)CDN 基础概念1. 定义2. 组成部分 (二)CDN 加速原理1. 请求路由2. 内容缓存3. 内容更新 (三)CDN 落地实践1. 选择 CDN 服务商2. 配置 CDN3. 集成到 Web 架构 …...
适应性Java用于现代 API:REST、GraphQL 和事件驱动
在快速发展的软件开发领域,REST、GraphQL 和事件驱动架构等新的 API 标准对于构建可扩展、高效的系统至关重要。Java 在现代 API 方面以其在企业应用中的稳定性而闻名,不断适应这些现代范式的需求。随着不断发展的生态系统,Java 在现代 API 方…...
在鸿蒙HarmonyOS 5中使用DevEco Studio实现指南针功能
指南针功能是许多位置服务应用的基础功能之一。下面我将详细介绍如何在HarmonyOS 5中使用DevEco Studio实现指南针功能。 1. 开发环境准备 确保已安装DevEco Studio 3.1或更高版本确保项目使用的是HarmonyOS 5.0 SDK在项目的module.json5中配置必要的权限 2. 权限配置 在mo…...
深入理解 React 样式方案
React 的样式方案较多,在应用开发初期,开发者需要根据项目业务具体情况选择对应样式方案。React 样式方案主要有: 1. 内联样式 2. module css 3. css in js 4. tailwind css 这些方案中,均有各自的优势和缺点。 1. 方案优劣势 1. 内联样式: 简单直观,适合动态样式和…...
FTXUI::Dom 模块
DOM 模块定义了分层的 FTXUI::Element 树,可用于构建复杂的终端界面,支持响应终端尺寸变化。 namespace ftxui {...// 定义文档 定义布局盒子 Element document vbox({// 设置文本 设置加粗 设置文本颜色text("The window") | bold | color(…...
欢乐熊大话蓝牙知识17:多连接 BLE 怎么设计服务不会乱?分层思维来救场!
多连接 BLE 怎么设计服务不会乱?分层思维来救场! 作者按: 你是不是也遇到过 BLE 多连接时,调试现场像网吧“掉线风暴”? 温度传感器连上了,心率带丢了;一边 OTA 更新,一边通知卡壳。…...
