生物医学信号处理--随机信号的数字特征

前言

概率密度函数完整地表现了随机变量和随机过程的统计性质。但是信号经处理后再求其概率密度函数往往较难,而且往往也并不需要完整地了解随机变量或过程的全部统计性质只要了解其某些特定方面即可。这时就可以引用几个数值来表示该变量或过程在这几方面的特征。
例如,只要知道某随机过程在某一时刻的平均值和实际值相对于这个平均值的分散程度,这时只要了解它的均值和方差即可。

均值

E ( x ) = ∫ − ∞ + ∞ x p ( x ; t ) d x E(x)=\int_{-\infty}^{+\infty}xp(x;t)dx E(x)=+xp(x;t)dx
相当于 t 时刻所有样本 x i ( t ) x_i(t) xi(t)总体均值 lim ⁡ N → + ∞ 1 N ∑ i = 1 N x i ( t ) { \lim_{N \to +\infty} \frac{1}{N}\sum_{i=1}^N x_i(t) } N+limN1i=1Nxi(t)当过程平稳时,均值与时间无关是常数。
如果平稳过程又是各态遍历的,总体均值将等于时间均值。此时均值可根据单样本来求
E ( x ) = m x = lim ⁡ T → ∞ 1 2 T ∫ − T + T x ( t ) d t E(x)=m_x={ \lim_{T \to \infty} \frac{1}{2T}\int_{-T}^{+T} x(t) }dt E(x)=mx=Tlim2T1T+Tx(t)dt
时间均值和总体均值 两者区别:
时间均值:单一样本在全部时间历程内的统计结果
总体均值:全部样本集合在固定瞬间的统计结果。

均方

E ( x 2 ) = ∫ − ∞ + ∞ x 2 p ( x ; t ) d x E(x^2)=\int_{-\infty}^{+\infty}x^2p(x;t)dx E(x2)=+x2p(x;t)dx
均方即全部样本集合在固定时刻的平均平方值 lim ⁡ N → + ∞ 1 N ∑ i = 1 N x i 2 ( t ) { \lim_{N \to +\infty} \frac{1}{N}\sum_{i=1}^N x^2_i(t) } N+limN1i=1Nxi2(t)
如果平稳过程又是各态遍历的,它等于时间均方。
E ( x 2 ) = D x = lim ⁡ T → ∞ 1 2 T ∫ − T + T x 2 ( t ) d t E(x^2)=D_x={ \lim_{T \to \infty} \frac{1}{2T}\int_{-T}^{+T} x^2(t) }dt E(x2)=Dx=Tlim2T1T+Tx2(t)dt
它反映的是随机信号的平均功率。

方差

反映随机变量在其均值附近的集中程度。方差愈小,随机变量取值愈集中。
E [ x − E ( x ) ] 2 = ∫ − ∞ + ∞ [ x − E ( x ) ] 2 p ( x ; t ) d x E[x-E(x)]^2=\int_{-\infty}^{+\infty}[x-E(x)]^2p(x;t)dx E[xE(x)]2=+[xE(x)]2p(x;t)dx

如果平稳过程又是各态遍历的,它也等于
E [ x − E ( x ) ] 2 = σ x 2 = lim ⁡ T → ∞ 1 2 T ∫ − T + T [ x ( t ) − m x ] 2 d t E[x-E(x)]^2=\sigma_x^2={ \lim_{T \to \infty} \frac{1}{2T}\int_{-T}^{+T} [x(t)-m_x]^2}dt E[xE(x)]2=σx2=Tlim2T1T+T[x(t)mx]2dt

自相关

数字特征表示单一时刻随机变量的特征。
自相关函数则表征信号在不同时刻取值间的关联程度,在确定性情况下,
对于非周期信号: R x ( τ ) = ∫ − ∞ + ∞ x ( t ) x ( t + τ ) p ( x t , x t + τ ; t , t + τ ) d t R_x(\tau)=\int_{-\infty}^{+\infty}x(t)x(t+\tau)p(x_t,x_{t+\tau};t,t+\tau)dt Rx(τ)=+x(t)x(t+τ)p(xt,xt+τ;t,t+τ)dt
对于周期信号: R x ( τ ) = lim ⁡ T → ∞ 1 T ∫ T 2 T 2 x ( t ) x ( t + τ ) d t R_x(\tau)={ \lim_{T \to \infty} \frac{1}{T}\int_{\frac{T}{2}}^{\frac{T}{2}} x(t)x(t+\tau) }dt Rx(τ)=TlimT12T2Tx(t)x(t+τ)dt
随机信号情况下:设 t 1 t_1 t1时刻的随机变量是 x 1 x_1 x1 t 2 t_2 t2时刻的随机变量是 x 2 x_2 x2
在这里插入图片描述
对于平稳过程,因为联合概率密度函数只是 τ = t 1 − t 2 \tau=t1-t2 τ=t1t2差值的函数,所以所得结果只与 τ \tau τ有关,与时间起点的取法无关。
在这里插入图片描述

过程各态遍历时,它又等于单一样本的自相关函数

在这里插入图片描述
对于信号是非周期的情况,一般 τ → + ∞ \tau→+∞ τ+时, R x ( τ ) → 0 R_x(\tau)→0 Rx(τ)0。波形的相邻值愈是独立无关, R x ( τ ) R_x(\tau) Rx(τ) τ \tau τ衰减得愈快。例如,白噪声的特点是不同时刻取值完全独立,因此,它的自相关函数只在 τ = 0 \tau=0 τ=0时有值, τ ≠ 0 \tau≠0 τ=0时, R x ( τ ) = 0 R_x(\tau)=0 Rx(τ)=0,即 R ( z ) R(z) R(z) δ \delta δ函数型。

平稳信号的自相关函数的性质

在这里插入图片描述

自协方差函数

进行信号处理时,往往先把均值(直流分量)除去再做剩余部分的相关函数。所得结果为自协方差, m x m_x mx为直流分量:
在这里插入图片描述

随机信号不同时刻的关系

在这里插入图片描述

平稳随机信号的功率谱

对于功率信号 s ( t ) s(t) s(t)其功率P定义为 P = lim ⁡ T → ∞ 1 2 T ∫ − T T s 2 ( t ) d t P={ \lim_{T \to \infty} \frac{1}{2T}\int_{{-T}}^{T} s^2(t)}dt P=Tlim2T1TTs2(t)dt
在这里插入图片描述
s ( t ) s(t) s(t)功率的频域表示:
P = 1 2 π ∫ − ∞ + ∞ lim ⁡ T → ∞ 1 2 T ∣ F s T ( w ) ∣ 2 d w P={\frac{1}{2\pi}\int_{-\infty}^{+\infty} \lim_{T \to \infty} \frac{1}{2T}|F_{s_T}(w)|^2}dw P=2π1+Tlim2T1FsT(w)2dw

周期性信号可以在一个周期范围内做傅里叶级数展开,因此有离散的功率谱,而非周期信号能在全时间范围内做傅里叶变换,因此有连续的能量谱。随机信号虽然是非周期的,但持续时间却往往是无限的,因此性质上属于功率信号。
单个样本 x i ( t ) x_i(t) xi(t) 的功率谱密度为: S x i ( w ) = lim ⁡ T → ∞ 1 2 T ∣ X i T ( w ) ∣ 2 S_{x_i}(w)={ \lim_{T \to \infty} \frac{1}{2T}|X_{i_T}(w)|^2} Sxi(w)=Tlim2T1XiT(w)2

其中:(这里的 X i T ( w ) X_{i_T}(w) XiT(w) F s T ( w ) F_{s_T}(w) FsT(w)同个意思)在这里插入图片描述
平稳随机信号𝑋(𝑡) 的功率谱密度(简称功率谱)定义为:
S X ( w ) = E [ S x i ( w ) ] S_X(w)=E[S_{x_i}(w)] SX(w)=E[Sxi(w)]
经过推导可得:
S X ( w ) = ∫ − ∞ ∞ R x ( τ ) e − j w τ d τ S_X(w)={\int_{{-\infty}}^{\infty} R_x(\tau)e^{-jw\tau}}d\tau SX(w)=Rx(τ)ejwτdτ
其中 R x ( τ ) R_x(\tau) Rx(τ)为自相关函数

维纳-辛钦定理

在这里插入图片描述

功率谱密度的性质

在这里插入图片描述

例题

随机信号有下述三个样本函数组成: x 1 t = 1 , x 2 t = s i n t , x 3 t = c o s t x_1t = 1,x_2 t = sin t, x_3 t = cos t x1t=1,x2t=sint,x3t=cost,且等概率发生,计算信号的均值和自相关函数
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

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

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

LabVIEW数据库管理系统

LabVIEW数据库管理系统(DBMS)是一种集成了数据库技术与数据采集、控制系统的解决方案。通过LabVIEW的强大图形化编程环境,结合数据库的高效数据存储与管理能力,开发人员可以实现高效的数据交互、存储、查询、更新和报告生成。LabV…

合并模型带来的更好性能

研究背景与问题提出 在人工智能领域,当需要处理多个不同任务时,有多种方式来运用模型资源。其中,合并多个微调模型是一种成本效益相对较高的做法,相较于托管多个专门针对不同任务设计的模型,能节省一定成本。然而&…

Virgo:增强慢思考推理能力的多模态大语言模型

每周跟踪AI热点新闻动向和震撼发展 想要探索生成式人工智能的前沿进展吗?订阅我们的简报,深入解析最新的技术突破、实际应用案例和未来的趋势。与全球数同行一同,从行业内部的深度分析和实用指南中受益。不要错过这个机会,成为AI领…

本地缓存:Guava Cache

这里写目录标题 一、范例二、应用场景三、加载1、CacheLoader2、Callable3、显式插入 四、过期策略1、基于容量的过期策略2、基于时间的过期策略3、基于引用的过期策略 五、显示清除六、移除监听器六、清理什么时候发生七、刷新八、支持更新锁定能力 一、范例 LoadingCache<…

Android adb shell GPU信息

Android adb shell GPU信息 先 adb shell 进入控制台。 然后&#xff1a; dumpsys | grep GLES Android adb shell命令捕获systemtrace_android 抓trace-CSDN博客文章浏览阅读2.5k次&#xff0c;点赞2次&#xff0c;收藏8次。本文介绍了如何使用adbshell命令配合perfetto工…

ElasticSearch | Elasticsearch与Kibana页面查询语句实践

关注&#xff1a;CodingTechWork 引言 在当今大数据应用中&#xff0c;Elasticsearch&#xff08;简称 ES&#xff09;以其高效的全文检索、分布式处理能力和灵活的查询语法&#xff0c;广泛应用于各类日志分析、用户行为分析以及实时数据查询等场景。通过 ES&#xff0c;用户…

RK3588平台开发系列讲解(系统篇)Linux Kconfig的语法

文章目录 一、什么是Kconfig二、config模块三、menuconfig四、menu 和 endmenu五、choice 和 endchoice六、source七、depends on八、default九、help十、逻辑表达式沉淀、分享、成长,让自己和他人都能有所收获!😄 一、什么是Kconfig Kconfig的语法及代码结构非常简单。本博…

STM32 USB组合设备 MSC CDC

STM32 USB组合设备 MSC CDC实现 教程 教程请看大佬niu_88 手把手教你使用USB的CDCMSC复合设备&#xff08;基于stm32f407&#xff09; 大佬的教程很好&#xff0c;很详细&#xff0c;我调出来了&#xff0c;代码请见我绑定的资源 注意事项 值得注意的是&#xff1a; 1、 cu…

深入学习RabbitMQ的Direct Exchange(直连交换机)

RabbitMQ作为一种高性能的消息中间件&#xff0c;在分布式系统中扮演着重要角色。它提供了多种消息传递模式&#xff0c;其中Direct Exchange&#xff08;直连交换机&#xff09;是最基础且常用的一种。本文将深入介绍Direct Exchange的原理、应用场景、配置方法以及实践案例&a…

Node.js——path(路径操作)模块

个人简介 &#x1f440;个人主页&#xff1a; 前端杂货铺 &#x1f64b;‍♂️学习方向&#xff1a; 主攻前端方向&#xff0c;正逐渐往全干发展 &#x1f4c3;个人状态&#xff1a; 研发工程师&#xff0c;现效力于中国工业软件事业 &#x1f680;人生格言&#xff1a; 积跬步…

【Verdi实用技巧-Part2】

Verdi实用技巧-Part2 2 Verdi实用技巧-Part22.1 Dump波形常用的task2.1.1 Frequently Used Dump Tasks2.1.2 Demo 2.2 提取波形信息小工具--FSDB Utilities2.3 Debug in Source code view2.3.1 Find Scopes By Find Scope form 2.3.2 Go to line in Souce code View2.3.3 Use B…

web-前端小实验4

实现以上图片中的内容 代码 <!DOCTYPE html> <html lang"en"> <head><meta charset"UTF-8"><meta name"viewport" content"widthdevice-width, initial-scale1.0"><title>用户注册</title&…

NLP项目实战——基于Bert模型的多情感评论分类(附数据集和源码)

在当今数字化的时代&#xff0c;分析用户评论中的情感倾向对于了解产品、服务的口碑等方面有着重要意义。而基于强大的预训练语言模型如 Bert 来进行评论情感分析&#xff0c;能够取得较好的效果。 在本次项目中&#xff0c;我们将展示如何利用 Python 语言结合transformers库&…

TAS测评倍智题库 | 益丰大药房2025年中高层测评BA商业推理测评真题考什么?

您好&#xff01;您已被邀请参加360评估。您的评估与反馈将有助于被评估人更深入地了解个人情况&#xff0c;发现个人优势和潜在风险。请您秉持公正、开放的心态进行评估。请尽快完成评估&#xff0c;在此衷心感谢您的配合与支持&#xff01; ​ 相关事宜&#xff1a; 请您在…

优秀的大模型会不会做坏事?

主要围绕大型语言模型&#xff08;LLMs&#xff09;在特定情境下可能出现的欺骗行为及相关研究展开&#xff0c;具体如下&#xff1a; 研究背景与核心发现&#xff1a;研究发现即使在用户无意激励的情况下&#xff0c;LLMs 也可能说谎&#xff0c;而能使用工具的模型更易被诱导…

fiscoBcos落盘加密介绍

落盘加密 落盘加密是在机构内部进行的&#xff0c;每个机构独立管理自己硬盘数据的安全。内网中&#xff0c;每个节点的硬盘数据是被加密的。所有加密数据的访问权限&#xff0c;通过Key Manager来管理。Key Manager是部署在机构内网内&#xff0c;专门管理节点硬盘数据访问秘…

完全二叉树的删除

&#xff08;1&#xff09;删除叶子节点 找到要删除的节点 targetNode找到要删除节点的父节点parent&#xff08;父节点是否存在&#xff09;要删除的节点是父节点的左子树还是右子树如果是左子树&#xff0c;则parent.leftnull;如果是右子树则parent.rightnull。 &#xff08;…

ModuleNotFoundError: No module named ‘setuptools_rust‘ 解决方案

大家好,我是爱编程的喵喵。双985硕士毕业,现担任全栈工程师一职,热衷于将数据思维应用到工作与生活中。从事机器学习以及相关的前后端开发工作。曾在阿里云、科大讯飞、CCF等比赛获得多次Top名次。现为CSDN博客专家、人工智能领域优质创作者。喜欢通过博客创作的方式对所学的…

【算法】时间复杂度以及O(N^2)的排序

目录 1.常数时间的操作 2.时间复杂度 2.1.以选择排序为例 2.2.O(n^2)从何而来 2.3.冒泡排序 2.3.1.抑或运算 2.4.插入排序 3.二分法 3.1.局部最小 4.递归 4.1.递归行为时间复杂度的估计 1.常数时间的操作 一个操作如果和样本的数据量无关&#xff0c;每次都是固定时…

uni app 写的 小游戏,文字拼图?文字拼写?不知道叫啥

从下方的偏旁部首中选在1--3个组成上面文章中的文字&#xff0c;完成的文字标红 不喜勿喷 《满江红》 其中用到了两个文件 strdata.json parameters.json 这两个文件太大 放到资源中了 资源文件 <template><view class"wenzi_page_main"><view c…