第 4 卷 第 11 期 中 国 水 运 Vol.4 No.11 2006 年 11 月 China Water Transport Novembdr 2006 收稿日期:2006-9-16 作者简介:张丹丹 武汉理工大学土木工程与建筑学院 (430070) 基于 MATLAB 的 AR 模型参数估计 张丹丹 徐 振 孙希宁 万 平 夏亚伟 摘 要:分析在地震作用下建筑物扭转振动效应的模型参数估计,根据系统振动加速度作为时间序列,来建立 AR ( n )模型,并利用时序模型参数与模态参数之间的关系,求解振动系统的模态参数。 关键词:时间序列 AR模型 参数估计 MATLAB 中图分类号:TP311.51 文献标识码:A 文章编号:1006-7973(2006)11-0109-02 一、时间序列的建模 1.ARMA 模型的定义 对于一个振动系统,某测点的输出(响应)为{} tx ,可 有 n 阶自回归—m 阶滑动平均混合时序模型来描述这个振动 系统,记为 AR MA( n , m )。即: 1 1 1 1t t i t i n t n t t j t j m t mx x x x a a a a ϕ ϕ ϕ θ θ θ − − − − − − − − − − − = − − − − 式中, i ϕ 为自回归系数,( i =1,2,⋯, n) ; j θ 为滑动平 均系数,( j =1,2,⋯, m ); t a 为白噪声序列, 2 (0, ) ta a NID σ ∼ , 2 a σ 为白噪声方差。 用线性后移算子 B,上式可表示为: ( ) ( ) ttB x B a ϕθ = 式中 1 ()1 in inB B B B ϕϕϕϕ = − − − − −, 1()1 im im B B B B θθθθ = − − − − −, 特殊地,若 n =0,模型称为纯滑动平远近模型,记为 MA ( m );若 m =0,模型称为纯自回归模型,记为 AR ( n );若 n = m =0,模型(1)退化为 tt x a= 即{} tx 为白噪声序列。 2.数据的采集与标准化处理 时序模型的建立需要离散的时间序列{} tx 。在本文所涉及的算例分析和试验分析中,所获取的加速度信号均为离散 信号,且信号的时间间隔Δ=0.02 s ,信号的持续时间为 16 s 。 对于时间序列{} tx ,当其取值过大或过小时,为保证计算精度、减少误差、避免溢出,可对{} tx 进行标准化处理。 记所得的时序为{} tx ,当{} tx 满足均值为 ˆx µ 、方差为 2ˆ x σ 的正态分布时,对{} tx 中各数据进行如下标准化处理: ˆ ˆ tx t x x y µ σ − = 注意:为了方便起见,本文将经过预处理的时间 序列仍然记为{} tx 。 3.模型阶次的确定 对时间序列的{} tx ( t =1,2,⋯ N ),首先要进行相关 性分析。相关性分析的任务是计算序列{} tx 的样本自相关函数 和样本偏自相关函数,并由他们的截尾性来进行模型类别的判断。可根据表1进行模型结构的初选。 表1 ARMA(n,m)模型的序列特征表 AR( n ) MA( m ) ARMA( n , m ) 自相关函数 拖尾 截尾 k = m 处 拖尾 偏自相关函数 截尾 k = n 处 拖尾 拖尾 ARMA( n , m )中的 n , m 参数的确定: n , m 并不能直接确定,而是需要先假定一组值,一般是 从(1,1)开始,建立模型,然后逐步升高 n , m 的值,求 出一系列模型

matlab 模型参数估计值,基于MATLAB的AR模型参数估计相关推荐

  1. matlab传递函数参数辨识,基于matlab/Simulink的参数辨识

    基于Simulink的辨识 Simulink自带Parameter Estimation功能.可以对Simulink模型中的参数进行估计.MATLAB的Parameter Estimation官方说明 ...

  2. matlab 线性回归 参数显著性,基于Matlab的栅格数据一元线性回归及显著性检验(slope趋势分析)...

    %by yinlichang3064@163.com 在进行长时间序列的栅格数据分析时,如NDVI,fvc,LAI,NPP,GPP,需要知道每个格点的长期趋势. 如果再arcgis中进行一元回归计算, ...

  3. matlab 赋权法,基于MATLAB的非线性曲线拟合赋权法

    成都航空职业技术学院学报 Journal of Chengdu Aeronautic Polytechnic 2018 年 12 月第 4 期(总第 117 期) Vol.34 No.4(Serial ...

  4. Matlab之DNN:基于Matlab利用神经网络模型(epochs=10000000)预测勒布朗詹姆斯的2018年总决赛(骑士VS勇士)第一场得分、篮板、助攻

    Matlab之DNN:基于Matlab利用神经网络模型(epochs=10000000)预测勒布朗詹姆斯的2018年总决赛(骑士VS勇士)第一场得分.篮板.助攻 目录 输出结果 预测勒布朗詹姆斯的20 ...

  5. matlab功率谱的峰值,基于Matlab功率谱密度估计方法.doc

    基于Matlab功率谱密度估计方法 基于Matlab功率谱密度估计方法要在实际情况下, 许多平稳信号无法导出数学表达式, 要准确获取这些信号的功率谱密度存在一定的困难.根据维纳辛钦 (Wiener K ...

  6. 为什么MATLAB用ar模型预测误差很大,关于AR模型预测的参数问题

    本帖最后由 zhuyuanzhao 于 2017-4-19 09:27 编辑 问题1:最近利用AR模型仿真电力线背景噪声,利用示波器录制了一段波形,进行计算,已经求得AR模型的参数以及预测方差值,如何 ...

  7. 距离矢量算法matlab实现,一种基于最小费用距离模型的城市生态网络构建方法与流程...

    本发明涉及生态网络构建技术领域,特别是涉及一种城市网络的构建方法. 背景技术: 最小费用距离是网络分析的一种计算方法,这种方法被用于物种保护.自然保护区功能规划.动物栖息地的确定.区域生态安全格局设计 ...

  8. matlab 随机骨料程序,基于matlab的混凝土三维圆形骨料模型随机投放方法.pdf

    一试验研究 斑楚遽1村 2Ol2年 基于matlab的混凝土三维圆形骨料模型 随机投放 方法 张海波 ,.何军拥. (1.广州航海高等专科学校,广东 广州 510330:2.广东工贸职业技术学院,广东 ...

  9. matlab 随机骨料程序,基于matlab的混凝土三维圆形骨料模型随机投放方法

    2012 年■ 试验研究 基于 matlab 的混凝土三维圆形骨料模型随机投放方法 张海波 1,何军拥 2 (1.广州航海高等专科学校,广东广州 510330: 2.广东工贸职业技术学院,广东广州 5 ...

最新文章

  1. 【蓝桥java】递归基础之振兴中华
  2. 使用JDBC进行批处理
  3. Java 取整数的个、十、百位
  4. C++模板中的函数式参数
  5. ArcIMS 投影问题 如何向地图配置文件中添加投影信息(转载)
  6. python 读取鼠标选中文本_python怎么读取文本文件
  7. 学生选课系统代码-2view视图层代码【MVC--v】代码
  8. VS code 尝试在目标目录创建文件时发生一个错误
  9. 谷歌chrome浏览器突然不记录登录状态,怎么办
  10. [源码和文档分享]基于SSM的网上购物系统的设计与开发
  11. 人生辣么多的谎言,没必要一个个试一下
  12. 仿京东轮播的广告展示栏
  13. [读书][笔记]WINDOWS PE权威指南《零》PE基础
  14. Automated_bounty_Hunter全自动漏洞赏金猎人使用场景一
  15. 转:别把赌注压在明天:如何了解自己,提高自控力
  16. centos中nginx使用
  17. 电梯黑板签到的测试用例
  18. proteus实现电子钟计时
  19. Numpy和Pandas的简单使用
  20. java5.0下载_java虚拟机下载 v5.0 官方免费版

热门文章

  1. 大数据传输为何这么慢?如何提高大数据传输效率
  2. 匿名飞控TI版_PID部分,串级PID,微分先行,前馈控制
  3. Python数据分析--平民最强宝可梦
  4. jzxx1087歌德巴赫猜想
  5. stm32、gd32等芯片固件升级思路
  6. GAD计算机辅助诊断,GAD-7-PHQ-9可打印.doc
  7. Linux-Mysql 日志
  8. Packet Tracer has stopped working unexpectedly.
  9. 模型精度不降反升!飞桨是这样改进PACT量化算法的!
  10. 阿里云SSL证书到期怎么办