現代信號處理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.pswp.cn/bicheng/20770.shtml
繁體地址,請注明出處:http://hk.pswp.cn/bicheng/20770.shtml
英文地址,請注明出處:http://en.pswp.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 語…

使用Java構建RESTful API:實現靈活、可擴展的Web服務

RESTful API已經成為構建現代Web應用的標準之一&#xff0c;它通過簡單的HTTP協議進行通信&#xff0c;提供了一種輕量級、靈活、可擴展的方式來構建和管理Web服務。Java作為一種強大的編程語言&#xff0c;提供了許多框架和庫來幫助開發者構建高效的RESTful API。本文將探討如…

項目質量管理

目錄 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…

JVM之【執行引擎】

執行引擎 執行引擎是JVM的核心組件之一&#xff0c;它負責將Java字節碼文件轉換為機器指令并執行。這一過程涉及多個組成部分&#xff0c;各部分協同工作來完成字節碼到機器指令的轉換和執行。以下是執行引擎的主要組成部分及其作用&#xff1a; 1. 解釋器&#xff08;Interp…

vue.js框架快速入門

Vue.js是一個漸進式JavaScript框架&#xff0c;用于構建用戶界面和單頁應用程序。以下是Vue.js快速入門的基本步驟和概念&#xff1a; 1. 環境準備 確保你的計算機上安裝了Node.js&#xff0c;它包括npm&#xff08;Node Package Manager&#xff09;&#xff0c;用于管理項目…

友善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…

TypeScript 中的迭代器和生成器

1. 迭代器 迭代器是一種對象&#xff0c;它提供了一種統一的方式來訪問集合中的元素&#xff0c;而不暴露集合的內部結構。在 TypeScript 中&#xff0c;迭代器通過實現 Iterator 接口來定義。 interface Iterator<T> {next(): IteratorResult<T>; }interface It…

Liunx登錄時相關bash配置文件(登錄腳本)

profile類的文件&#xff1a;設定環境變量&#xff0c;運行命令或腳本 bashrc類的文件&#xff1a;定義命令別名 全局配置&#xff1a; /etc/profile /etc/profile.d/*.sh /etc/bashrc 個人配置文件: ~/.bash_profile ~/.bashrc 用戶登錄時加載bash配置文件的過程&#xff0c;…

碼隨想錄算法訓練營第二十四天| 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 &…

升級Jenkins從2.263.3到2.440.2

升級Jenkins從2.263.3到2.440.2 ###任何一次升級前&#xff0c;先做整體備份&#xff0c;同時最好對plugins目錄和config.xml單獨備份。 ###對于任何一次插件安裝&#xff0c;務必安裝前先備份當前的plugins目錄&#xff0c;這是血的教訓。升級過程 1、 升級2.263.3到2.263.4…

15-通過JS代碼處理窗口滾動條

selenium并不是萬能的&#xff0c;頁面上有些操作無法實現時&#xff0c;就需要借助JS代碼來完成了。selenium提供了一個方法&#xff1a;execute_script()&#xff0c;可以執行JS腳本代碼。 比如&#xff1a;當頁面上的元素超過一屏后&#xff0c;想操作屏幕下方的元素&#x…

MyBatis查詢功能

MyBatis的各種查詢功能 1、若查詢出的數據只有一條 &#xff08;1)可以通過實體類對象或者集合接收 (2)可以通過List集合接收 (3&#xff09;可以通過map集合接收 結果&#xff1a;{password123456, sex女, id8, age22, email1234qq.com, usernameadmin4} 2、若查詢出來的數據有…

[leetcode hot 150]第一百零八題,將有序數組轉換為二叉搜索樹

題目&#xff1a;給你一個整數數組 nums &#xff0c;其中元素已經按 升序 排列&#xff0c;請你將其轉換為一棵 平衡二叉搜索樹。 給定一個有序的整數數組,我們需要構建一棵平衡的二叉搜索樹。平衡二叉樹是指任意一個節點的左右子樹的高度差不超過1。 由于給定的數組是有序的…

阿里云計算學習筆記(一)

運維管理 運維管理&#xff08;Operation and Maintenance Management, 簡稱O&M管理&#xff09;是指通過科學的管理方法和技術手段&#xff0c;對IT系統和基礎設施進行監控、維護、優化和保障&#xff0c;以確保系統的高可用性、穩定性、安全性和性能。運維管理涵蓋了硬件…