0
  • 聊天消息
  • 系统消息
  • 评论与回复
登录后你可以
  • 下载海量资料
  • 学习在线课程
  • 观看技术视频
  • 写文章/发帖/加入社区
会员中心
创作中心

完善资料让更多小伙伴认识你,还能领取20积分哦,立即完善>

3天内不再提示

基于Matlab实验仿真对经典功率谱估计性能分析

西西 来源:博客园 作者:J博士 2020-10-08 17:46 次阅读

功率谱估计在分析平稳各态遍历随机信号频率成分领域被广泛使用,并且已被成功应用到雷达信号处理、故障诊断等实际工程中。本文给出了经典功率谱估计的几类方法,并通过Matlab的实验仿真对经典功率谱估计方法性能进行了分析,最后说明了经典功率谱估计法的局限性和造成这种局限性的原因。

1.引言

给定一个标准的正弦信号,我们可以通过傅里叶变换来分析它的频率成分。然而,实际工程应用中,由于存在着各种干扰、噪声,我们得到的信号往往不是理想的,如图1-1这种信号,具有不确定性,幅度不能预知,非周期,但往往服从一定的统计特性,这种信号叫作随机信号。需要注意的是,本文所说的随机信号是指平稳各态遍历的随机信号,关于非平稳随机信号的分析方法[1]本文不予讨论。

基于Matlab实验仿真对经典功率谱估计性能分析

图1-1 一种随机信号时域形式

对于图1-1的随机信号,我们可以通过功率谱来分析它的频率成分,如图1-2所示为图1-1随机信号的功率谱。实际过程中,我们只能获得随机信号的一些离散数据点(假设为N个),本文将讨论如何利用这N个数据点,来得到一个"非精确"的功率谱来对真实随机信号的功率谱进行估计,并讨论如何更好的估计,即在下一章要讲述的几个经典的功率谱估计法。

基于Matlab实验仿真对经典功率谱估计性能分析

图1-2 上图所示的随机信号功率谱

2.经典功率谱估计法

上一章我们已经知道功率谱估计法是通过利用已经获得的N个数据点,来得到一个"非精确"的功率谱对真实随机信号的功率谱进行估计,所以在给出具体的方法之前,如何来评价我们得出的这个"非精确"的功率谱的好坏呢?

评价功率谱性能好坏的标准有很多,本文只给出两个影响最大的标准:分辨率和方差。分辨率即功率谱上能够区分的最小相邻频率成分,分辨率越高,我们观察信号的频率成分越清晰;方差大小则反映到功率谱波动性的大小,如果方差太大,功率谱波动性大,则很容易造成有用的频率成分被噪声淹没。所以,我们希望得到的这个"非精确"的功率谱,分辨率越高越好,方差越小越好。

2.1.2周期图法性能(Matlab仿真)

上一小节我们已经给出了周期图法的原理。本节将通过Matlab仿真给出数据点数N对功率谱性能好坏的影响,正如上文所述,将通过对所得功率谱的分辨率和方差两方面进行分析。

图2-1 实验所用的随机信号

当数据点数N分别为128、256、512和1024时,得到的功率谱分别如图2-2、图2-3、图2-4和图2-5所示。分辨率能够直观的通过功率谱图形看出,方差的数值由表2-1给出。

表2-1 不同N值得到功率谱的方差值

N 128 256 512 1024
方差 92.7108 130.9109 160.9187 483.5894

通过上面实验结果的比较,我们很容易发现,周期图法得到的功率谱随着数据点数N的增大,分辨率变大、方差变也大。

2.1.3平均周期图法

周期图法得到的功率谱与我们所期望的"分辨率大、方差小"是矛盾的。为了进一步降低方差,将N个观测样本数据点分为L段,每段数据长度为M, 分别对每段数据求周期图功率谱估计,然后求平均值,这种方法称平均周期图法。

2.1.4平均周期图法性能(Matlab仿真)

当数据点数N为1024,分段数分别为8、4、2时,平均周期图法得到的功率谱分别如图2-6、图2-7、图2-8所示。分辨率能够直观的通过功率谱图形看出,方差的数值由表2-2给出。

表2-2 不同L值得到功率谱的方差值

L 8 4 2 1
方差 96.3756 190.9647 400.6464 483.5894

L=1时,平均周期图法退化为周期图法。通过上面实验结果的比较,我们很容易发现,平均周期图法得到的功率谱随着分段数L变大,方差变小,但分辨率变小。

当观测样本序列数据个数N固定时,要降低方差需要增加分段数L。当N不大时分段长度M取值较小,则功率谱分辨率降低到较低的水平。若分段数L固定时,增加分辨率需要分段长度M,则需要采集到更长的检测数据序列。实际中恰恰是检测样本序列长度不足。

2.1.5修正的平均周期图法

上一节已经提到实际中检测样本序列长度是有限的。对现有数据长度N,如果能获得更多的段数分割,将会得到更小的方差。允许数据段间有重叠部分,来得到更多的段数。对段间重叠长度的选取,最简单是取为段长度M的一半。由式(2-5)可知更多的段数可以进一步降低方差。

数据截断的过程中相当于数据加矩形窗,矩形窗幅度较大的旁瓣会造成"频谱泄漏"。我们分段时采取的窗函数更为多样(三角窗,海明窗等), 以减小截断数据(加矩形窗)窗函数带来的影响[2]

2.1.6修正的平均周期图法性能(Matlab仿真)

利用修正平均周期图法,分别使用矩形窗、Blackman窗和Hamming窗得到的功率谱如图2-9所示。

图2-9 不同窗函数的修正平均周期图法得到的功率谱

可以发现,矩形窗的分辨率最高,但是方差也最大,这是由于矩形窗频谱主瓣最窄,分辨率因此最高,旁瓣也高,导致频谱泄漏最严重,方差最大。

2.1.7总结

周期图法获得的功率谱随着样本点数越多,分辨率越大、方差越大;平均周期图法以牺牲分辨率来进一步改善方差;修正的平均周期图法允许段的重叠来进一步增大分段数、或者分段数相同,每段样本点数变多。无论是哪种方法都没有彻底结局方差与分辨率之间的矛盾。

2.2相关功率谱估计法-BT法

正如我们之前介绍的,要提高功率谱估计的分辨率,必须增加数据序列的长度N,但是较长的数据序列,由噪声引起的随机性得到更为充分的体现-较大的方差。事实上,当N无穷大时,方差为一非零常数。即周期图法无法实现功率谱的一致估计。而这节讲述的相关功率谱估计法(下文称作BT法),是一致估计。

2.2.1 BT法的原理

维纳辛钦定理指出,随机信号的相关函数与它的功率谱是一对傅里叶变换对。BT法就是基于这个原理。先由观测数据估计出自相关函数,然后求自相关函数的傅立叶变换,以此变换作为对功率谱的估计,也称为间接法。BT法要求信号长度N以外的信号为零,这也造成BT法的局限性。

2.2.2 BT法的性能(Matlab仿真)

数据点数N分别为128、256、512和1024的BT法,得到的功率谱如图2-10、图2-11、图2-12和图2-13所示。

图2-13 N=1024时,BT法得到的功率谱

由上面实验可以发现,M随着N的增大而增大时,分辨率提高,方差变大。BT法仍然没有解决分辨率与方差之间的矛盾,但是BT法得到的功率谱当N为无穷大时,方差会趋向于零,即为一致估计[2]。

2.2.3 周期图法与BT法的关系

结 论

本文通过Matlab仿真,以一个具体的随机信号为例,简单介绍了周期图法、平均周期图法、修正的平均周期图法以及BT法的基本原理,并对这些方法的性能进行分析。可以看出,无论是周期图法及其改进算法还是BT法都没有从根本上解决分辨率与方差的矛盾。经典功率谱估计是利用傅里叶变换估计功率谱,而我们之前分析随机信号不满足傅里叶变换的条件,所以经典功率谱估计方法不得不从无限长数据点截取有限长数据点,加入限制条件(周期图法实际上假定N点外数据周期重复、BT法假定N点外数据为零)来"强制"作傅里叶变换,这也是造成它局限性的原因。

参考资料

[1]朱哲,钟宏伟. 非平稳随机信号分析处理方法研究[J] 安徽电子信息技术学院学报2008.6:28-28

[2]皇甫堪.现代数字信号处理[M].电子工业出版社

声明:本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人,不代表电子发烧友网立场。文章及其配图仅供工程师学习之用,如有内容侵权或者其他违规问题,请联系本站处理。 举报投诉
  • 功率谱
    +关注

    关注

    0

    文章

    13

    浏览量

    8345
  • 雷达信号处理

    关注

    2

    文章

    13

    浏览量

    9299
收藏 人收藏

    评论

    相关推荐

    安泰功率放大器在电机特性实验研究中的应用

    形式进行优化设计,提出了三种不同的电机形式,分别对电机进行理论模型分析,运用有限元进行仿真,对电机进行整体结构设计并制作样机,进行性能实验研究。 测试设备:
    的头像 发表于 10-31 11:18 149次阅读
    安泰<b class='flag-5'>功率</b>放大器在电机特性<b class='flag-5'>实验</b>研究中的应用

    这些电源常用仿真软件,你都知道吗?

    ,便于工程师验证分析。 电源仿真为电源系统开发带来很多便利 : ▍电源的设计与分析: 在电源方案制定过程中,可利用电源仿真手段辅助电源电路设计与
    发表于 10-25 14:20

    Matlab/Simulink/Stateflow建模开发及仿真测试

    matlab 模拟仿真 熟悉Matlab/Simulink/Stateflow建模开发及仿真测试,熟悉V模型开发流程。 熟悉自动代码生成,能够编写或者配置自动代码生成脚本。
    发表于 10-24 17:23

    PPEC-HIL 三相整流逆变实时仿真测试

    研/教学解决方案,全面覆盖电力电子学科的教学和科研需求,通过虚拟仿真与实物实验相结合,有效解决了教学内容复杂、理论与实际脱节的问题。教学平台涵盖几乎所有经典电力电子实验,配备完备的课程
    发表于 10-12 15:40

    宽带功率放大器在可调谐激光器解调实验的应用

    实验名称:可调谐激光器解调实验研究方向:光纤光栅传感、可调谐激光实验目的:基于使用光谱仪解调的Buneman频率估计解调算法,编写了一种适用于可调谐激光器解调法的解调程序作为解调结果的
    的头像 发表于 08-26 15:25 868次阅读
    宽带<b class='flag-5'>功率</b>放大器在可调谐激光器解调<b class='flag-5'>实验</b>的应用

    常用的仿真软件有哪些

    在当今的工程和技术领域,仿真软件已经成为了一种非常重要的工具。它们可以帮助工程师和研究人员在实际制造和测试之前,对系统或产品进行模拟和分析。以下是一些常用的仿真软件及其特点: MATLAB
    的头像 发表于 08-22 09:18 2778次阅读

    教学验证篇丨PPEC+HIL 单相逆变仿真验证

    —————————————————EasyGo电力电子科研/教学解决方案,全面覆盖电力电子学科的教学和科研需求,通过虚拟仿真与实物实验相结合,有效解决了教学内容复杂、理论与实际脱节的问题。教学平台涵盖几乎所有经典电力电子
    发表于 08-09 10:25

    如何利用Matlab进行神经网络训练

    Matlab作为一款强大的数学计算软件,广泛应用于科学计算、数据分析、算法开发等领域。其中,Matlab的神经网络工具箱(Neural Network Toolbox)为用户提供了丰富的函数和工具
    的头像 发表于 07-08 18:26 1866次阅读

    功率分析仪在高校科研、实验室的应用

    理论模型、进行精确的能效分析实验数据的可靠性至关重要。02应用需求1.准确的功率测量对于验证理论模型、进行精确的能效分析实验数据的可靠性
    的头像 发表于 06-04 08:33 345次阅读
    <b class='flag-5'>功率</b><b class='flag-5'>分析</b>仪在高校科研、<b class='flag-5'>实验</b>室的应用

    基于MATLAB的信号处理系统与分析

    在信息技术日新月异的今天,信号处理与分析技术在通信、控制、生物医学、图像处理等领域的应用日益广泛。MATLAB作为一款功能强大的数学软件,为信号处理与分析提供了强大的工具和支持。本文旨在探讨
    的头像 发表于 05-17 14:24 1067次阅读

    MATLAB电路仿真时间要考虑那些因素

    MATLAB电路仿真是一种常用的工具,用于模拟和分析电子电路的性能。在进行电路仿真时,需要考虑各种因素,以确保
    的头像 发表于 04-21 09:27 739次阅读

    MATLAB电路仿真中能讲电流汇总的元件

    MATLAB电路仿真中,可以使用许多元件来汇总电流。以下是一些常见的元件和它们在电路中的作用: 电阻(Resistor):电阻是电路中最基本的元件之一,用于限制电流的流动。在MATLAB电路
    的头像 发表于 04-21 09:23 3234次阅读

    Matlab仿真】单相交交变频器研究与实现,Simulink建模

    数量多使之应用受到了一定的限制,因此它广泛应用于大功率交流电动机调速传动系统。 单相交交变频器框架 仿真分析 (1)当给定信号输出频率f=10Hz 时,触发角为30度 时,电压电流波形如下。 (2
    发表于 03-26 17:20

    DSP软件 MATLAB仿真实验报告

    电子发烧友网站提供《DSP软件 MATLAB仿真实验报告.pdf》资料免费下载
    发表于 03-24 09:49 13次下载

    如何使用MATLABMATLAB Parallel Server扩展整车仿真呢?

    使用 MATLAB 和 Simulink 设计和构建 Rivian 车辆仿真界面平台帮助我们实现了关键目标。我们为工程师和非工程师创建了统一平台,用于运行整车仿真、后处理结果和创建报告。
    的头像 发表于 01-10 18:22 1215次阅读
    如何使用<b class='flag-5'>MATLAB</b>和<b class='flag-5'>MATLAB</b> Parallel Server扩展整车<b class='flag-5'>仿真</b>呢?