2005iF 12月 笫 21卷第 6期 武警工程学 院学报 jOURNAL OF ENGG COLLEGE OF ARMED POLICE FORCE DeC.2005 Vo1.21 No.6 【计算机技术及应用】 基于 AR模型法功率谱估计 的 Matlab实现 周 渊 王 炳和2 (1武警上海总队,上海 200435;2武警工程学 院,陕西 西安 710086) 【摘 要】 本文根据现代谱 估计 中的线性预测 自回归模型法(AR模型法)估计功率谱的原理 ,讨论 了 Levinsion算法和 Burg算法计算 AR模型参数、估计功率谱,并利用 Matlab工具进行 了实例计算和分析。结 果表明:不加窗的 burg算法得到的随机信号功率谱较为真实,没有 明显的频率偏移和假峰 ,并且具有较高的 频 率分 辨率 。 【关键词】 功率谱估计 ;AR模型法;建模 ;Matlab 引 言 通常 ,由于有用信号与噪声的频谱特性不同,因此谱估计方法成 为一种在噪声背景下提取有用信号(如 正弦信号)的有效方法。谱估计的方法主要有非参数化方法和参数化方法 ,或称为经典谱估计方法和现代谱 估计方法。经典谱估计方法的优点是方法简便、计算效率高 ,其不足是频率分辨率低 。现代谱估计方法具有 频率分辨率高的优点 ,因此又被称为高分辨率谱估计法。近年来 ,现代谱估计理论和技术的研究一直十分活 跃。现代谱估计的方法主要有模型法、熵谱法 、最大似然法和特征分解法等四大类。 模型法包括 AR模型法 、MA模型法、ARMA模型法 ,由于三者之间可以相互表示 ,而对 AR模型参数的 估计得到的是线性方程 ,故 AR模型比 MA、ARMA模型具有在计算上简便 的优点 ,且实际物理系统往往是 全极点系统 ,所以研究有理分式传递函数的模型主要是讨论 AR模型⋯。 AR模型法做功率谱估计的原理是:假定所分析的信号 z(n)是 由一个均方误差为 a 的白噪声 (n)激 励一线性移不变系统 H(z)(即 AR模型)所得到 ,则分析信号的功率谱估计为 P ( ) P . (∞): d2 I H(2)I : H(2)= ———1 ———一 1.+∑akz 其中 a 为输入序列的方差 ,a。、a2、⋯⋯、a 为代估参数。 由上可见 AR模型法估计功率谱实质是求解模型系数 a1、a2、⋯⋯、a 和 d 的问题。本文采用 Levin— sion递推算法与 Burg递推算法求解模型系数 ,得出功率谱估计值 ,并进行 比较。 1 算 法简介 1.1 Levinsion递 推算 法 Levinsion递推算法是 Yule—Walker方程的快速解法。其计算过程为 :先从一阶的 Yule—Walker方程 11 l‰ (1)‰ (0)I l口 } L J 中解出口11和 a 分别为 口ll=一 (1)/ (0),di=(1一I口ll I ) (0)。再~,-j--阶的 Yule—Walker方程 收稿 日期 :2005—03—01 作者筒介 :心渊(1980一 ),男 ,2005年武警工程学院通信工程系硕士毕业 ,现为武警 上海 总队参谋 ;王炳 和,武 警工程学 院通信工程系 教授。 73 维普资讯 http://www.cqvip.com http://www.cqvip.com/ 武警 工程学院学报 2005年第 6期 (0) (1) (2) (1) (0) 工 (1) (2) (1) (0) i 0 0 利用一阶 Yu1e—Walker方程的解 n11和 di计算出二阶 Yule—Walker方程的解。

基于AR模型的功率谱matlab,基于AR模型法功率谱估计的Matlab实现相关推荐

  1. matlab 计算变异系数,变异系数法求权重matlab 代码

    利用matlab编程,很容易根据变异系数法,求得多指标综合评价模型的权重. 代码如果有不懂的地方,可以联系我. 变异系数法求权重matlab 代码 clear;clc; [data1,header1] ...

  2. 用matlab编程节点电压法求电路,MATLAB在电路中的应用

    <MATLAB在电路中的应用>由会员分享,可在线阅读,更多相关<MATLAB在电路中的应用(59页珍藏版)>请在人人文库网上搜索. 1.MATLAB应用(三) Matlab在电 ...

  3. matlab 计算变异系数,变异系数法求权重matlab代码

    <变异系数法求权重matlab代码>由会员分享,可在线阅读,更多相关<变异系数法求权重matlab代码(1页珍藏版)>请在读根文库上搜索. 1.变异系数 法求权重 matlab ...

  4. MATLAB幂法原点平移法,幂法及其MATLAB程序.doc

    幂法及其MATLAB程序 5.2 幂法及其MATLAB程序 5.2.2 幂法的MATLAB程序 用幂法计算矩阵的主特征值和对应的特征向量的MATLAB主程序 function [k,lambda,Vk ...

  5. matlab符号函数绘图法_《MATLAB符号运算及其应用》(黄忠霖著).pdf

    <MATLAB符号运算及其应用>(黄忠霖著).pdf [General Information] 书名=MATLAB符号运算及其应用 作者=黄忠霖,黄京编著 页数=425 出版社=北京市: ...

  6. matlab转换为部分分式,【MATLAB用部分分式展开法资讯】MATLAB用部分分式展开法足球知识与常识 - 足球百科 - 599比分...

    MATLAB中对多项式进行部分分式展开. 用MATLAB进行部分分式展开 MATLAB有1个命令用于求B(s)/A(s)的部分分式展开式. 设s的有理分式为 式中 (i=)和(j=)的某些值可能为零. ...

  7. MATLAB之牛顿下山法

    MATLAB之牛顿下山法 算法原理 matlab程序 算法原理 上一篇博客,我介绍了牛顿法迭代法,接下来我就们接着讲解一下什么是牛顿下山法. 一.迭代公式 在牛顿迭代过程中,若满足单调性|f(x(k+ ...

  8. ar模型matlab fpe,基于Matlab的AR模型参数估计.pdf

    基于Matlab的AR模型参数估计.pdf 维普资讯 2OO5年第39卷No4 39 基于Matlab的AR模型参数估计* 陈国强 赵俊伟 黄俊杰 刘万里 河南理工大学 摘 要:基于Matlab用时间 ...

  9. matlab 模型参数估计值,基于MATLAB的AR模型参数估计

    第 4 卷 第 11 期 中 国 水 运 Vol.4 No.11 2006 年 11 月 China Water Transport Novembdr 2006 收稿日期:2006-9-16 作者简介 ...

最新文章

  1. 世界各大天文台联合预警:今晚公布“引力波重要发现”
  2. python画二维散点图-python3怎样画二维点图
  3. linux下xampp(apache)中配置域名访问,以及遇到的问题
  4. Golang 使用Protocol Buffer 案例
  5. (转载)SpringBoot 发送邮件和附件
  6. navigator工具_Javascript常用工具类
  7. java中线程调度遵循的原则_深入理解Java多线程核心知识:跳槽面试必备
  8. Khronos EGL and Apple EAGL
  9. Office Timeline一款非常专业好用的PPT插件
  10. 常用企业管理工具介绍
  11. Java Initialization Order
  12. python代码画樱花-如何用Python代码实现樱花树效果
  13. git 回退远端master/develop分支版本
  14. 几行 Java 代码搞定图片提取文字功能
  15. H5 video标签播放MP4
  16. 神经系统体格检查概要
  17. POI-word模板替换
  18. 计算机应用基础作业一(答案),计算机应用基础作业一(答案).doc
  19. 口语 计算机 评分标准,高考英语口语考试说明及评分标准20163661676.doc
  20. 7个月吸粉13.7万人!泰禾广场Wi-Fi上的“智慧商业”

热门文章

  1. java ntp服务器ip地址_NTP服务器时间同步
  2. 记一次因tcp_tw_recycle误设置为1导致服务器无法正常建立TCP连接的问题排查
  3. Github md文件换行和标签转义
  4. FLV视频加密系统解决方案
  5. 基于ROS的机械臂手眼标定-基础使用
  6. 入职新公司,如何快速熟悉一个系统
  7. 人工智能数学基础6:极限、极限运算、ε-δ语言、ε-N语言、级数和函数连续性
  8. 企业大数据部署 Hadoop的新选择
  9. 雾霾治理力度持续加大,空净企业如何破局?
  10. 2021年Java软件开发有前途吗?