MATLAB中AR模型功率谱估计中AR阶次估计的实现-psd_my.rar(最近看了几个关于功率谱的问题,有关AR模型的谱估计,在此分享一下,希望大家不吝指正)(声明:本文内容摘自我的毕业论文——心率变异信号的预处理及功率谱估计)(按:AR模型功率谱估计是对非平稳随机信号功率谱估计的常用方法,但是其模型阶次的估计,除了HOSA工具箱里的arorder函数外,没有现成的函数可用,arorder函数是基于矩阵SVD分解的阶次估计方法,为了比较各种阶次估计方法的区别,下面的函数使用了'FPE','AIC','MDL','CAT'集中准则一并估计,并采用试验方法确定那一个阶次更好。
)………………………………以上省略……………………………………………………………………假设原始数据序列为x,那么n阶参数使用最小二乘估计在MATLAB中实现如下:Y=x;Y(1:n)=[];m=N-n;X=[];%构造系数矩阵fori=1:m  forj=1:n      X(i,j)=xt(ni-j);  endendbeta=inv(X'*X)*X'*Y';复制代码beta即为用最小二乘法估计出的模型参数。
此外,还有估计AR模型参数的Yule-Walker方程法、基于线性预测理论的Burg算法和修正的协方差算法等[26]。
相应的参数估计方法在MATLAB中都有现成的函数,比如aryule、arburg以及arcov等。
4.3.3AR模型阶次的选择及实验设计文献[26]中介绍了五种不同的AR模型定阶准则,分别为矩阵奇异值分解(SingularValueDecomposition,SVD)定阶法、最小预测定误差阶准则(FinalPredictionErrorCriterion,FPE)、AIC定阶准则(Akaika’sInformationtheoreticCriterion,AIC)、MDL定阶准则以及CAT定阶准则。
文献[28]中还介绍了一种BIC定阶准则。
SVD方法是对Yule-Walker方程中的自相关矩阵进行SVD分解来实现的,在MATLAB工具箱中arorder函数就是使用的该算法。
其他五种算法的基本思想都是建立目标函数,阶次估计的标准是使目标函数最小化。
以上定阶准则在MATLAB中也可以方便的实现,下面是本文实现FPE、AIC、MDL、CAT定阶准则的程序(部分):form=1:N-1  ……    %判断是否达到所选定阶准则的要求  ifstrcmp(criterion,'FPE')    objectfun(m1)=(N(m1))/(N-(m1))*E(m1);  elseifstrcmp(criterion,'AIC')    objectfun(m1)=N*log(E(m1))2*(m1);  elseifstrcmp(criterion,'MDL')    objectfun(m1)=N*log(E(m1))(m1)*log(N);  elseifstrcmp(criterion,'CAT')    forindex=1:m1        temp=temp(N-index)/(N*E(index));    end    objectfun(m1)=1/N*temp-(N-(m1))/(N*E(m1));  end    ifobjectfun(m1)>=objectfun(m)    orderpredict=m;    break;  endend复制代码orderpredict变量即为使用相应准则预测的AR模型阶次。
(注:以上代码为结合MATLAB工具箱函数pburg,arburg两个功率谱估计函数增加而得,修改后的pburg等函数会在附件中示意,名为pburgwithcriterion)登录/注册后可看大图程序1.JPG(35.14KB,下载次数:20352)下载附件 保存到相册2009-8-2820:54上传登录/注册后可看大图程序2.JPG(51.78KB,下载次数:15377)下载附件 保存到相册2009-8-2820:54上传下面本文使用3.2.1实验设计的输出结果即20例经预处理的HRV信号序列作为实验对象,分别使用FPE、AIC、MAL和CAT定阶准则预测AR模型阶次,图4.1(见下页)为其中一例典型信号使用不同预测准则其目标函数随阶次的变化情况。
从图中可以看出,使用FPE、AIC以及MDL定阶准则所预测的AR模型阶次大概位于10附近,即阶次10左右会使相应的目标函数最小化,符合定阶准则的要求,使用CAT定阶准则预测的阶次较小,在5~10之间。
图4.2(见下页)为另一例信号的阶次估计情况,从中也可以得到同样的结论。
(注,实验信号为实验室所得,没有上传)登录/注册后可看大图图片1.JPG(28.68KB,下载次数:5674)下载附件 保存到相册2009-8-2820:54上传
2025/6/27 16:08:25 6KB matlab
1
与visio开发有关的书籍非常的少,为了做一个visio二次开发的项目,买了一本,光盘中带有书籍的电子版。
书中更多的讲解了VBA和C++的开发。
不管用什么语言开发,书中有很多值得学习和借鉴的地方。
特此与大家共享,供大家收藏!介绍Visio环境和有关《开发MicrosoftVisio解决方案》的概念性信息。
使用公式设计模拟现实世界物体和行为的智能图形的详细信息。
使用MicrosoftVisualBasicforApplication(VBA)扩展Visio,或者将Visio用作自己应用程序的一个组件的信息、提示和技术。
介绍如何使用MicrosoftVisualBasic和C++编程语言来开发将Visio应用程序当作组件使用的程序。
2025/6/27 13:04:28 21.82MB visio 二次开发
1
《MapInfo9.5中文版标准教程》以最新版本的MapInfo9.5中文版为操作平台,全面介绍使用该软件绘制、编辑桌面式地图的方法和技巧。
全书共分13章,内容涉及MapInfo9.5的数据采集、矢量地图绘制和编辑、属性数据编辑(表操作)、地图查询分析、空间模型、地图装饰和MapBasic基础等内容。
《MapInfo9.5中文版标准教程》在讲解软件功能的同时,安排了丰富的课堂练习和上机练习,帮助读者巩固知识。
配套光盘附有多媒体语音视频教程和大量的图形文件,供读者学习和参考。
2025/6/27 3:26:37 49MB MapInfo 9.5 中文教程
1
该项目实现的功能如下:获取手机相册,点击每个相册之后进入该相册的图片列表界面,在图片列表界面可以实现图片多选,然后进入所选择的图片界面,在该界面内可以实现所选图片的上传等功能。
该项目最大特色:1、获取相册列表,目前网络上面介绍获取相册的项目很少,本文专门讲述相册的获取。
2、使用Android-Universal-Image-Loader集成框架-第三方jar包加载本地图片,熟悉这个jar的开发者肯定不陌生,该jar包十分强大,除了可以获取网络图片,本地图片也是可以的。
同时,通过引用第三方jar可以有效解决OOM异常问题。
博文地址:http://blog.csdn.net/u010156024/article/details/44136543欢迎访问!【握手】
2025/6/27 3:03:27 1.32MB android图片
1
详细介绍了FSR压力传感器的应用以及芯片管脚、时序图等
2025/6/26 20:44:44 964KB FSR压力传感器
1
《水文小波分析》简要介绍小波分析的基本理论和常用的主要方法,重点论述小波分析方法在水文学中的各种应用。
主要内容包括小波分析的基本理论、小波函数及其构造、小波快速算法、水文序列滤波与去噪、水文过程复杂性描述、水文系统多时间尺度分析、水文序列奇异性及趋势性分析、水文预测预报和水文随机模拟等方面。
该书为国内水文小波分析领域的第一本专著。
其特点是内容新颖,理论联系实际,深入浅出,便于理解和实际分析计算。
作者王文圣,丁晶,李跃清等。
2025/6/26 14:45:04 12.31MB 水文小波分析
1
在Linux操作系统中,`man`命令是不可或缺的一部分,它提供了系统的在线帮助文档,允许用户查阅各种命令、系统调用、库函数、配置文件等的详细信息。
本资源为"Linuxman中文手册",包含了丰富的中文解释,使得非英语环境的用户也能方便地学习和理解Linux系统操作。
`man`命令的使用方法非常简单。
在终端中输入`man`后跟需要查询的命令或函数名,例如`manls`将显示关于`ls`命令的使用手册。
手册通常分为多个章节,每个章节涵盖不同的主题。
章节号在手册页的顶部显示,例如1表示用户可执行的命令,2表示系统调用,3表示库函数等。
在安装这个中文手册前,你需要检查当前系统的语言环境,以确保手册显示为中文。
通过运行`locale`命令,你可以看到诸如`LC_ALL`、`LANG`等环境变量的设置,它们决定了系统显示语言。
如果希望显示中文,确保这些变量设置为支持中文的语言代码,如`zh_CN.UTF-8`。
在压缩包`man-pages-zh_CN-1.5`中,包含的是中文版的Linux手册页。
这些页面详细解释了各种Linux内核接口、系统调用、C库函数、shell命令以及系统管理工具的使用方法。
每个页面通常包含以下几个部分:1.**NAME**:简短介绍该功能的名称和用途。
2.**SYNOPSIS**:展示命令的基本语法和参数,或者函数的声明。
3.**DESCRIPTION**:详述命令或函数的工作原理、参数含义及使用场景。
4.**RETURNVALUE**(对于函数):说明函数执行后的返回值及其含义。
5.**CONFORMINGTO**:指出该功能遵循的标准化规范,如POSIX或UNIX标准。
6.**NOTES**:提供额外的注意事项或警告。
7.**BUGS**:列举已知的问题或局限性。
8.**EXAMPLES**:给出使用示例,帮助理解如何实际应用。
9.**SEEALSO**:推荐相关的命令、函数或文档供进一步阅读。
通过这个中文手册,无论是初学者还是经验丰富的系统管理员,都能更轻松地查找和理解Linux中的各种工具和功能。
它不仅涵盖了基本的命令行操作,还包括了系统管理和程序开发的相关知识,是Linux用户不可或缺的学习资源。
记得适时更新手册,以获取最新的信息和功能介绍。
2025/6/26 14:30:49 3.59MB
1
图像平滑或去噪就是为了抑制噪声,以达到改善图像质量的目的,既可以在空间域又可以频率域中实现,在数字图像处理中起着重要的作用。
本文将主要介绍空间域的几种平滑法的算法:邻点平均法、K个邻点平均法、最大均匀性平滑,其中操作平台是matlab7.1
1
介绍和效图:http://blog.csdn.net/geniusice18/article/details/8393674
2025/6/26 4:08:19 1.79MB QQ2012 界面 directUI
1
在性能基础之浅谈常见接口性能压测一文中我们有简单介绍常见的RPC接口,本文将单篇详细介绍RPC框架。
RPC(RemoteProcedureCall)—远程过程调用,它是一种通过网络从远程计算机程序上请求服务,而不需要了解底层网络技术的协议。
RPC协议假定某些传输协议的存在,如TCP或UDP,为通信程序之间携带信息数据。
在OSI网络通信模型中,RPC跨越了传输层和应用层。
其实简单点的说,就是像调用本地的类的方法样来调用服务器端的方法实现。
比如两个不同的服务A,B部署在两台不同的机器上,那么服务A如果想要调用服务B中的某个方法该怎么办呢?使用HTTP请求当然可以,但是可能会比较慢而且一些优化做的并
2025/6/25 21:09:32 359KB 性能基础之常见RPC框架浅析
1
共 1000 条记录 首页 上一页 下一页 尾页
在日常工作中,钉钉打卡成了我生活中不可或缺的一部分。然而,有时候这个看似简单的任务却给我带来了不少烦恼。 每天早晚,我总是得牢记打开钉钉应用,点击"工作台",再找到"考勤打卡"进行签到。有时候因为工作忙碌,会忘记打卡,导致考勤异常,影响当月的工作评价。而且,由于我使用的是苹果手机,有时候系统更新后,钉钉的某些功能会出现异常,使得打卡变得更加麻烦。 另外,我的家人使用的是安卓手机,他们也经常抱怨钉钉打卡的繁琐。尤其是对于那些不太熟悉手机操作的长辈来说,每次打卡都是一次挑战。他们总是担心自己会操作失误,导致打卡失败。 为了解决这些烦恼,我开始思考是否可以通过编写一个全自动化脚本来实现钉钉打卡。经过一段时间的摸索和学习,我终于成功编写出了一个适用于苹果和安卓系统的钉钉打卡脚本。
2024-04-09 15:03 15KB 钉钉 钉钉打卡