现代信号处理12_谱估计的4种方法(CSDN_20240602)

Slepian Spectral Estimator(1950)

做谱估计的目标是尽可能看清楚信号功率谱在某一个频率上的情况,假设我们想了解零频时的分布,最理想的情况是滤波器的传递函数H(ω) 是一个冲激函数,这样就没有旁瓣,也就没有泄漏;其次,主瓣宽度为零,分辨率极好。然而在现实中,理想的冲激函数是无法实现的,所以,只能允许H(ω) 有一定通带(假设为-βπ,βπ )。另外,由于滤波器是有限长度的,所以H(ω) 不可避免地会有泄漏,但我们可以要求泄漏尽可能地小,即做下面的优化

时,有最优解,其中U是矩阵B的特征向量,而h正是B最大特征值对应的那个特征向量。

小结

与传统的周期图谱估计方法相比,Slepian Spectral Estimator将着眼点转移到要观察的频率点附近,其目标是,设计一个对信号功率谱进行估计的滤波器,希望信号通过滤波器后,想要的频谱分量能够有效保存下来,而无关的谱分量尽可能被抑制。       这种估计存在的问题:

  1. 仅考虑滤波器自身的响应,而没有考虑信号,不同的信号通过滤波器会有不同的表现,特别是随机信号会特别复杂。这种谱分析是信号无关的,不管信号是什么样的,都使用同一个滤波器进行估计。
  2. 没有用到统计的观念,所有的信息都是确定的,只依赖于频带的宽度。然而对随机信号进行谱估计是,不能不用到统计。

Capon Spectral Estimator (1969)

与一般谱估计的对比

小结

        传统的谱分析:为了获得信号在某个频率点上能量情况,就极大化信号在这个频率点上的响应,Slepian方法就是典型代表:在要分析的频率的附近划一个区域,然后极大化信号在这个区域上的响应,尽可能抑制其它频率的响应。而Capon方法的思路是,在要观察的频率点上给定一个约束条件,使信号在这个频率上的响应得到保证,然后在此基础上提出新的要求:极小化信号在其它频率的响应。

MUSIC(Multiple Signal Classification)

MUSIC

至此,我们可以看到MUSIC方法至少在3个方面进行了创造性的工作:

·1. 对Y求相关阵,充分考虑了噪声

·2. 在解方程很困难的情况下,放弃了直接求解方程的方法,而是分析等式两端的秩,这样虽然不能完全得到方程的解,但是我们感兴趣的信息(信号频率)可以得到。

·3. 得到了信号子空间和方向矢量构成的子空间相同,进而得到方向矢量张成的子空间与噪声子空间正交的结论,这也是MUSIC方法的核心所在。

小结

MUSIC方法提出了“超分辨率”的理念,首次呈现了分辨率极高的谱图。在之前的谱分析中,分辨率取决于信号的长度,信号越长,分辨率越高。而在MUSIC方法中,分辨率和信号长度没有太大关系,以很短的数据就可以获得很高的分辨率,能得到这种效果的原因主要有2个:

·1. MUSIC是一种非线性方法,它脱离了过去用线性滤波提取频率分量的传统路线,因此它达到的水平是Slepian和Capon等方法无法相比的。

·2. MUSIC得到的是一种“伪谱”,谱峰的位置代表信号的频率,但高度并不代表信号在这个频率上的能量大小,而且这种方法也没有给出谱峰高度与信号能量之间的关系。

MUSIC谱是一种“伪谱”,只反映了方向矢量与噪声子空间之间的正交关系的良好程度。在理想情况下,它们之间应该是严格正交的,但是因为噪声的存在,谱图只反映了它们之间关系的一种估计。

       在上面的5步中,最困难的是第3步,因为有些时候无法准确地对特征值分组,如果分组不准,会对MUSIC性能造成致命的影响,因为特征分组错误,就意味着子空间估计错误,信号频率分量的个数也会估计错误。因此MUSIC是一种很脆弱的方法,这种方法对信噪比的要求特别高,只有在高信噪比的条件下,才能准确估计子空间的维数。在信噪比的条件下,周期图估计效果比较好,通常我将MUSIC和周期图两种方法结合使用:先用周期图法确定信号个数,然后在用MUSIC方法获取信号频率分量准确位置。

ROOT MUSIC

Min-Norn MUSIC

MUSIC的核心是方向矢量子空间与噪声子空间正交。从正交性出发,容易导致“伪峰”的出现。因为在寻找正交点的过程中,频率点在整个频率轴上滑动,方向矢量与噪声子空间的关系在不断变换,在某一频点两者正交时,我们判定这个频点就是我们希望得到的频点之一,即这个频点包含在信号中。然而我们仅仅知道正交时,该频点在信号中,那么不正交时有什么意义呢?或者说两者之间的夹角的大小有什么意义呢?两者之间的夹角为89度时是不是比60度更像信号呢?MUSIC方法并没有给出上面问题的答案。

那么能不能构造一种新的方法,既能利用MUSIC带来的超分辨率,又能引入误差的概念,作为频点与信号符合程度的判定依据呢?答案就是使用Min-Norm MUSIC方法。

小结

MUSIC方法最核心的思想是:方向矢量张成的子空间和信号矢量张成的子空间是同一个子空间,判断2个子空间是否相同时很困难的,尤其是两个子空间对应的基不相同的时候。在这里,信号子空间的基使彼此正交的,但是方向子空间的基并不正交,为了得到信号频率的信息,必须让两个子空间产生关系,所以前面提出了很多必要条件,Min-Norm MUSIC 对谱峰的高度做出了解释。

MUSIC与Capon方法比较

ESPRIT(Estimation of Signal Parameters with Rotation Invariant Technique)用旋转不变技术估计参数

ESPRIT方法可以一次求出我们想要的所有信号频点。

------------------------------------------------------------------------------------------------
因为文档中公式较多,不方便编辑,所以本文使用截图的方式展现。如需电子版文档,可以通过下面的链接进行下载。

链接icon-default.png?t=N7T8http://generatelink.xam.ink/change/makeurl/changeurl/11779

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

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

相关文章

回溯算法 -- 216. 组合总和 III

目录 一.题目描述 二.解题思路 三.回溯三部曲 3.1确定递归函数的参数 3.2确认递归的终止条件 3.3确定单层循环逻辑 四.具体的代码 一.题目描述 找出所有相加之和为 n 的 k 个数的组合,且满足下列条件: 只使用数字1到9每个数字 最多使用一次 返…

经典的滑动窗口的题目 力扣 2799. 统计完全子数组的数目(面试题)

给你一个由 正 整数组成的数组 nums 。 如果数组中的某个子数组满足下述条件,则称之为 完全子数组 : 子数组中 不同 元素的数目等于整个数组不同元素的数目。 返回数组中 完全子数组 的数目。 子数组 是数组中的一个连续非空序列。 示例 1&#xff1…

【1】AI介绍

迎接 AGI 时代 AGI(Artificial General Intelligence),人工通用智能,AGI是一种可以执行复杂任务的人工智能,能够完全模仿人类智能的行为。应用领域涉及医疗、交通、智能家居等多个与人类活动密切相关的领域。 AGI 多久会到来? 乐观预测:明年(未来已来)主流预测:3-5…

【云原生】Docker Compose 使用详解

目录 一、前言 二、Docker Compose 介绍 2.1 Docker Compose概述 2.2 Docker Compose特点 2.3 Docker Compose使用场景 三、Docker Compose 搭建 3.1 安装docker环境 3.2 Docker Compose安装方式一 3.2.1 下载最新版/如果不是最新可替换最新版本 3.2.2 设置权限 3.2.…

Zigbee +PC上位机 无线控制二维云台开发笔记

今日尝试开发一款简单好学的PC上位机无线控制二维云台的小试验品: 主要开发环境与工具介绍: 单片机 STM32F103C8T6 使用标准库函数编程 Visual Studio 2022软件C# Winform 开发 上位机控制软件 DL_20 无线串口模块 + USB-TTL 模块 实现无线通…

使用 Scapy 库编写 Ping of Death 攻击脚本

一、介绍 1.1 概述 Ping of Death(PoD)攻击是一种历史悠久的拒绝服务(DoS)攻击,攻击者通过发送特制的畸形ICMP Echo请求数据包,导致目标系统无法正确处理,从而导致系统崩溃、重启或无法响应正…

用于相似图片搜索引擎的Python OpenCV图像直方图

图像直方图 那么,图像直方图到底是什么? 图片 图像的构成是由像素点构成的,每个像素点的值代表着该点的颜色(灰度图或者彩色图)。所谓直方图就是对图像的中的这些像素点的值进行统计,得到一个统一的整体的…

手眼标定学习笔记

目录 标定代码: 手眼标定原理学习 什么是手眼标定 手眼标定的目的 eye in hand eye to hand AXXB问题的求解 标定代码: GitHub - pumpkin-ws/HandEyeCalib 推荐博文: https://zhuanlan.zhihu.com/p/486592374 手眼标定原理学习 参…

YOLOv8: 标注石头、识别边缘及计算面积的方案

YOLOv8: 标注石头、识别边缘及计算面积的方案 引言 YOLO(You Only Look Once)是一种非常有效的实时目标检测算法,自其首次发布以来就受到了广泛的关注和应用。YOLOv8 是这一系列算法的最新版本,继承了之前版本的高效性和准确性&a…

如何在一台电脑上安装多个版本的JDK并且切换使用?

如何在一台电脑上安装多个版本的JDK并且切换使用? 文章目录 如何在一台电脑上安装多个版本的JDK并且切换使用?1.目录管理2.下载JDK3.配置环境变量4. 验证 1.目录管理 我们需要首先对不同版本的JDK进行版本管理,如下所示,我在C盘路…

vue3状态管理,pinia的使用

状态管理 我们知道组件与组件之间可以传递信息,那么我们就可以将一个信息作为组件的独立状态(例如,单个组件的颜色)或者共有状态(例如,多个组件是否显示)在组件之传递,这样的话我们希…

STM32作业实现(七)OLED显示数据

目录 STM32作业设计 STM32作业实现(一)串口通信 STM32作业实现(二)串口控制led STM32作业实现(三)串口控制有源蜂鸣器 STM32作业实现(四)光敏传感器 STM32作业实现(五)温湿度传感器dht11 STM32作业实现(六)闪存保存数据 STM32作业实现(七)OLED显示数据 STM32作业实现(八)触摸按…

最新h5st(4.7.2)参数分析与纯算法还原(含算法源码)

文章目录 1. 写在前面2. 加密分析3. 算法还原 【🏠作者主页】:吴秋霖 【💼作者介绍】:擅长爬虫与JS加密逆向分析!Python领域优质创作者、CSDN博客专家、阿里云博客专家、华为云享专家。一路走来长期坚守并致力于Python…

HALCON-从入门到入门-最常用的算子-二值化

1.废话 图像处理中的二值化是一种将灰度图像转换为只有两种可能值(通常是0和255,分别代表黑色和白色)的过程。这个过程在数字图像处理中非常常见,因为它可以简化图像数据,突出图像的主要特征,并降低后续处…

5.25.12 数字组织病理学的自我监督对比学习

无监督学习可以弥补标记数据集的稀缺性。 无监督学习的一个有前途的子类是自监督学习,其目的是使用原始输入作为学习信号来学习显著特征。在这项工作中,我们解决了在没有任何监督的情况下学习领域特定特征的问题,以提高数字组织病理学界感兴…

【vue实战项目】通用管理系统:作业列表

目录 目录 1.前言 2.后端API 3.前端API 4.组件 5.分页 6.封装组件 1.前言 本文是博主前端Vue实战系列中的一篇文章,本系列将会带大家一起从0开始一步步完整的做完一个小项目,让你找到Vue实战的技巧和感觉。 专栏地址: https://blog…

【Linux】Linux工具——gcc/g++

1.使用vim更改信用名单——sudo 我们这里来补充sudo的相关知识——添加信任白名单用户 使用sudo就必须将使用sudo的那个账号添加到信用名单里,而且啊,只有超级管理员才可以添加 信用名单在/etc/sudoers里 我们发现它的权限只是可读啊,所以…

MD中 面料的物理属性参数

该图片是Marvelous Designer软件中"Fabric Physical Properties"(面料物理属性)面板的截图,用于调整面料在弯曲、折叠时的硬度(Buckling Stiffness)。 目标部分解释了调整Buckling Stiffness的作用:通过调整该百分比值来决定面料角落处的硬度。进入80%的Buckling St…

Linux网络编程:应用层协议|HTTPS

目录 1.预备知识 1.1.加密和解密 1.2.常见加密方式 1.2.1.对称加密 1.2.2.非对称加密 ​编辑 1.3.数据摘要(数据指纹)和数据签名 1.4.证书 1.4.1.CA认证 1.4.2.证书和数字签名 2.HTTPS协议 2.1.自行设计HTTPS加密方案 2.1.1.只使用对称加密 …

java调用科大讯飞离线语音合成SDK --内附完整项目

科大讯飞语音开放平台基础环境搭建 1.用户注册 注册科大讯飞开放平台账号 2.注册好后先创建一个自己的应用 创建完成后进入应用选择离线语音合成(普通版)可以看到我们开发需要的SDK,选择windows MSC点击下载。 3.选择你刚刚创建的应用,选择…