- 1、本文档共90页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
高等数字信号处理第3章 高阶谱估计
第三章 高阶谱估计 3.1 累量及高阶谱 3.2 高阶谱估计 3.3 有色噪声背景下的频率估计 3.4 高阶谱的应用 3.1 累量与高阶谱 (Cumulants and Higher Order Spectral 简记:HOS) 3.1.1、累量的定义 1、随机变量的特征函数和矩函数 对于高斯分布的随机变量, 矩函数的定义 2、累量的定义 高阶矩与高阶累量的关係(M-C公式): 3、平稳随机过程的累量 对于零均值实平稳随机过程{x(n)},其k阶矩(k阶相关函数)和k阶累量分别为: 4、高斯过程的累量 3.1.3、高阶谱 1、定义:假定随机过程{x(n)}的k阶累量是绝对可和的,则其k阶谱是k阶累量的(k-1)维富里叶变換,即 3.2 高阶谱估计 从己知一段样本序列{x(1),x(2),…….,x(N)}出发,进行高阶谱估计的方法,与功率谱估计类似,也可分为非参数法和参数法两大类。 3.2.1、非参数法谱估计 1、基本思路: 假定n=0或n=N+1范围内,样本值x(n)=0, 由高阶谱的定义直接构造谱估计式。 3.3 有色噪声中的频率估计 3.4 高阶谱的应用 应用高阶谱的动机大致有以下几点: 3.4.1、 从非高斯信号中提取信息。 这是基于累量描述了信号与高斯分布偏离的程度。实际上,任何周期信号、准周期信号都是非高斯的。例如,复杂的机械系统自身“辐射”的信号都是非高斯的。 3、无相位耦合过程的累量 奇阶次累量恒为零; 为非零整数 不为零的四阶累量可定义为: 4、二次相位耦合过程的累量: , 谐波过程中各分量由于相位耦合而相关。 最后可得三阶累量的定义为 由于二次相位耦合, 3.3.2、高斯噪声背景下的谐波恢复 1、线性予测法: 利用高斯噪声的四阶累量为零来抑制噪声的影响; 利用谐波过程为退化的AR过程的原理,将谐波恢复过程转化为ARMA参数辨识过程。 若令 2、谐波过程的四阶累量与二阶累量的关系 则有下式成立: 3、基于特征值分解的谐波分析法 基本思想是利用四阶累量与二阶累量(相关函数)的关係,将功率谱估计中的方法进行推广。 MUSIC法 (Multiple Signal Classification ) 包含p个复谐波的谐波信号x(n)的相关矩阵 的秩为P,所以有 前m个特征向量张成信号子空间S, 后(M-p+1)个特征值张成与S正交的噪声 子空间G。 由S与G正交得: 基于高阶累量的MUSIC法的关键是如何由观 测信号y(n)来获取它的相关矩阵的特征矢量 该信号与x(n)的频率完全相同。因此,可通过计 算y(n)的四阶累量来获取它的自相关矩阵及其特 征矢量。 频率分辨率:在非参数法高阶谱估计中,其富里叶变换都是用DFT实现的。因此,最后得到的高阶谱谱线间的距离(频率分辨率)必然与所用的样本序列的长度成反比。即用于计算DFT的时间序列长度越长,则频率分辨率越高。 估计方差:可以证明,非参数法高阶谱估计是渐近 无偏的,但一般存在较大的估计方差。为减少估计 方差,可采用时域平滑或频域平滑的方法,但平滑 的结果必然使频率分辨率下降。 因此,估计方差与频率分辨率之间的矛盾是非参数法谱估计的固有矛盾。 3、确定性信号的高阶谱 4、主要方法: 平滑周期图法(直接法) MATLAB实现: [bspec,waxis]=bispecd(x,nfft,wind,samp_seg,overlap) x:时域信号; nfft:FFT的长度; wind:Rao最优窗函数的长度; samp_seg:每个分段的长度;overlap:每段重迭长度; bspec:等高线显示的直接法双谱;waxis:频率点矩阵; 间接法:先估计高阶累量,再进行DFT。 MATLAB实现: [bspec,waxis]=bispeci(x,nlag, samp_seg, overlap , flag,nfft,wind,) x:时域信号; nfft:FFT的长度;wind:窗函数类型; samp_seg:每个分段的长度;overlap:每段重迭长度; nlag:计算累积量的最大延迟;flag:是否有偏; bspec:等高线显
您可能关注的文档
最近下载
- 第八章、特殊时期施工进度及保证措施.docx VIP
- 物理光学作业参考答案4.pdf
- 国家开放大学(山东)《公务员制度讲座》形考任务1-3+终结性考试参考答案.pdf
- 博世 gsh16-30 分解图 转子 定子 开关 碳刷 零配件.pdf
- 突发公共卫生事件心理干预研究进展.docx VIP
- 皇帝的新装省名师优质课赛课获奖课件市赛课一等奖课件.pptx VIP
- 湖北省2024年中考:《生物》考试真题与参考答案.pdf
- OptiXOSN9800系统硬件介绍-20130702-B概述.pptx
- 国开山东电大《素质与思想品德教育》形成性考核1-3+终结性考试答案.docx
- 借用资质合同协议书范本.doc
文档评论(0)