声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 4433|回复: 3

[编程技巧] 求Matlab中实现双谱估计!!

[复制链接]
发表于 2007-3-26 09:53 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
求Matlab中实现双谱估计!!


哪位高人知道  谢谢!
matlab7版本的   双谱估计函数
回复
分享到:

使用道具 举报

发表于 2009-9-23 16:27 | 显示全部楼层

回复 楼主 qingyi 的帖子

双谱估计和传统的功率谱估计一样,有三种方法:
1、直接法,使用函数bispecd(sig);
2、间接法,使用函数bispeci(sig);
3、参数模型估计,这里面主要有三种模型,AR 模型、MA模型和ARMA模型。

而参数模型的估计方法需要确定模型的阶次和模型的参数。
  arorder函数,用来确定AR模型的阶次,用p=arorder(y);%其中P返回估计的AR模型的阶次。y为观测序列信号。
  maorder函数,用来确定MA模型的阶次,用p=maorder(y)%其中p返回估计的MA模型的阶次。y为观测序列信号。
  arrcest函数,用来确定AR模型的参数,用avec=arrcest(y,p);%其中avec返回AR模型的参数,y为观测序列信号,P为AR模型的阶次。
  maest用来估计MA模型的参数,用bvec=maest(y,q);%其中bvec返回MA模型的参数,y为观测序列信号,q为MA模型的阶次。
  armaqs函数用q-slice算法来估计ARMA模型参数,用[avec,bvec]=armaqs(y,p,q),含义同上。
  armarts函数用残差时间序列估计ARMA模型参数,用[avec,bvec]=armarts(y,p,q),含义同上。
  bispect(ma,ar,128)函数用来实现双谱估计,其中若ma为1时,ar是由函数arrcest产生的是基于AR模型的双谱估计,ar为1,ma是由maest函数产生的时是基于ma模型的双谱估计,如果ar 和ma都是由armaqs或者是armarts产生的,则是基于ARMA模型的双谱估计。

[ 本帖最后由 ChaChing 于 2009-9-25 09:47 编辑 ]

评分

1

查看全部评分

发表于 2009-9-26 11:15 | 显示全部楼层

回复 沙发 我爱计算机 的帖子

AR模型定阶的函数arorder     matlab中没有呢?你自己编的嘛还是工具箱里的?能否提供
发表于 2009-9-29 17:02 | 显示全部楼层

回复 板凳 limingaini 的帖子

本身就有,你只要安装一下高阶谱分析工具箱就可以了!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-9-22 06:48 , Processed in 0.073527 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表