现代信号处理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

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

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

相关文章

【OpenHarmony】TypeScript 语法 ③ ( 条件语句 | if else 语句 | switch case 语句 )

文章目录 一、条件语句1、if else 语句2、switch case 语句 参考文档 : <HarmonyOS第一课>ArkTS开发语言介绍 一、条件语句 1、if else 语句 TypeScript 中的 if 语句 / if else 语句 用法 , 与 JavaScript 语言中的 if 语句 / if else 语句 语法 基本相同 ; if else 语…

项目质量管理

目录 1.概述 2.三个关键过程 2.1.规划质量管理&#xff08;Plan Quality Management&#xff09; 2.2.管理质量&#xff08;Manage Quality&#xff09; 2.3.控制质量&#xff08;Control Quality&#xff09; 3.应用场景 3.1.十个应用场景 3.2.产品设计与开发 4.小结…

使用PyCharm 开发工具创建工程

一. 简介 前面学习了 安装 python解释器。如何安装python的一种开发工具 PyCharm。 本文来简单学习一下&#xff0c;如何使用 PyCharm 开发工具创建一个简单的 python工程。 二. PyCharm 开发工具创建一个工程 1. 首先&#xff0c;首先打开PyCharm 开发工具。选择 创建一…

Docker部署SiYuan笔记-Unraid

使用unraid的docker部署SiYuan笔记&#xff0c;简单记录 笔记说明 Siyuan笔记是一款基于markdown语法的笔记工具&#xff0c;具有活跃的社区和多设备支持。大部分功能都是免费&#xff0c;源代码开源&#xff0c;支持插件安装&#xff0c;具有很不错的使用体验。 Docker地址&a…

linux---生产者和消费者模型

生产者消费者模式就是通过一个容器来解决生产者和消费者的强耦合问题。生产者和消费者彼此之间不直接通讯&#xff0c;而通过阻塞队列来进行通讯&#xff0c;所以生产者生产完数据之后不用等待消费者处理&#xff0c;直接扔给阻塞队列&#xff0c;消费者不找生产者要数据&#…

2024年海南省三支一扶报名指南,照片要求

2024年海南省三支一扶报名指南&#xff0c;照片要求 一、考试时间安排&#xff1a; 报名时间&#xff1a;6月1日8:00至6月7日18:00 准考证打印时间&#xff1a;6月17日8:00 考试时间&#xff1a;6月22日 二、招聘人数 海南省计划招募390名高校毕业生

STM32_IIC

1、IIC简介 I2C&#xff0c;即Inter IC Bus。是由Philips公司开发的一种串行通用数据总线&#xff0c;主要用于近距离、低速的芯片之间的通信&#xff1b;有两根通信线&#xff1a;SCL&#xff08;Serial Clock&#xff09;用于通信双方时钟的同步、SDA&#xff08;Serial Data…

友善RK3399v2平台利用rkmpp实现硬件编解码加速

测试VPU 编译mpp sudo apt update sudo apt install gcc g cmake make cd ~ git clone https://github.com/rockchip-linux/mpp.git cd mpp/build/linux/aarch64/ sed -i s/aarch64-linux-gnu-gcc/gcc/g ./arm.linux.cross.cmake sed -i s/aarch64-linux-gnu-g/g/g ./arm.lin…

如何学习ai agent?

如何学习Agent&#xff0c;推荐阅读《动手做AI Agent》这本书。 推荐理由&#xff1a; 1&#xff1a;一本书能够全方位了解并探索Agent的奥秘&#xff01; &#xff08;1&#xff09;Agent的发展进程。 &#xff08;2&#xff09;可以帮我们做哪些事&#xff1a;自动办公&am…

码随想录算法训练营第二十四天| 77. 组合

77. 组合 - 力扣&#xff08;LeetCode&#xff09; class Solution {ArrayList<Integer> path new ArrayList<>();ArrayList<List<Integer>> result new ArrayList<>();public List<List<Integer>> combine(int n, int k) {if(n &…

15-通过JS代码处理窗口滚动条

selenium并不是万能的&#xff0c;页面上有些操作无法实现时&#xff0c;就需要借助JS代码来完成了。selenium提供了一个方法&#xff1a;execute_script()&#xff0c;可以执行JS脚本代码。 比如&#xff1a;当页面上的元素超过一屏后&#xff0c;想操作屏幕下方的元素&#x…

[leetcode hot 150]第一百零八题,将有序数组转换为二叉搜索树

题目&#xff1a;给你一个整数数组 nums &#xff0c;其中元素已经按 升序 排列&#xff0c;请你将其转换为一棵 平衡二叉搜索树。 给定一个有序的整数数组,我们需要构建一棵平衡的二叉搜索树。平衡二叉树是指任意一个节点的左右子树的高度差不超过1。 由于给定的数组是有序的…

Prime1 - 提权的另一种解法,彻底搞懂OpenSSL解密渗透提权,超强思路版。

提权枚举 现在我们直接从低权限用户开始&#xff1b;我们先按照提权步骤&#xff0c;简单的系统枚举 虽然我们知道可以利用系统版本低进行内核提权&#xff0c;内核提权虽然比较快比较方便&#xff0c;但也比较暴力&#xff0c;缺点非常明显&#xff1b;很容易导致系统服务中…

【云原生】Kubernetes----POD控制器

目录 引言 一、Pod控制器概述 二、Pod控制器的种类 &#xff08;一&#xff09;ReplicaSet &#xff08;二&#xff09;Deployment &#xff08;三&#xff09;StatefulSet &#xff08;四&#xff09;DaemonSet &#xff08;五&#xff09;Job 三、使用POD控制器 &a…

【Seafile】Seafile容器版文件删除后存储空间不释放问题解决

Seafile是一款非常优秀的网盘系统&#xff0c;我们可以根据官方文档&#xff0c;在本地虚拟机研究Seafile免费版的安装和使用&#xff0c;安装建议采用使用docker容器的方式。 不过在使用过程中&#xff0c;刚接触的小伙伴可能会遇到这样的问题&#xff1a; 删除网盘里面的文…

C++设计模式-状态模式

文章目录 28. 状态模式 运行在VS2022&#xff0c;x86&#xff0c;Debug下。 28. 状态模式 状态模式让一个对象的行为随着内部状态的改变而改变&#xff0c;而该对象也像换了类一样。应用&#xff1a;如在游戏开发中&#xff0c;游戏有不同场景&#xff0c;如主菜单、开始、战斗…

Houdini的PythonScript基本使用

前言 Houdini内置了Python脚本和相应的编辑器, 很多时候想灵活的制作各种Houdini工具, 基本是必须用到 Python。Houdini官方的python提供了非常完善的接口, 比如可以创建各种节点&#xff0c;连接各种节点&#xff0c;遍历节点各种数据&#xff0c;遍历节点参数等等。 Houdin…

word避免画质画质模糊方法

问题描述&#xff1a;   近期写文档时会高频率贴图&#xff0c;粘图过程中发现Word会自动压缩图片画质&#xff0c;而且压缩得很严重&#xff0c;下面是一幅图被压缩前后的画质对比 &#xff08;图片压缩前&#xff09; &#xff08;图片压缩后&#xff09; 解决方案&#x…

基于JSP的九宫格日志网站

你好呀&#xff0c;我是学长猫哥&#xff01;如果有需求可以文末加我。 开发语言&#xff1a;Java 数据库&#xff1a;MySQL 技术&#xff1a;JSP技术 工具&#xff1a;浏览器/服务器&#xff08;B/S&#xff09;结构 系统展示 首页 管理员功能模块 用户功能模块 摘要 本…

GPT-4o VS GPT-3.5 完胜

前言&#xff1a; 最近&#xff0c;GPT-4o已经限时免费开放了&#xff0c;试了一下&#xff0c;然后&#xff0c;说我的时间到了&#xff0c;然后&#xff0c;有给我转到3.5&#xff0c;正好遇到一个问题做一下对吧&#xff0c;感觉4O完胜啊。3.5还是很好胡诌&#xff0c;也就…