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
大图找小图程序在网上找不到现成的,不是收费就是不完整,或者不符合要求写了一个找图的子程序(开源),通过参数调用,可以自定义找图区域,定义透明色,定义图片相似程度精度为1效率最高,1024*768找图不会超过300毫秒。
精度越低效率越低,如果你有更好的优化方法请qq告知,感谢QQ83833303
2025/6/7 0:49:20 2.11MB 找图 找小图
1
XIU主题是一款Wordpress主题,支持百度熊掌号,适用于图片展示、多元化图片新闻展示、个人博客、资源分享站,扁平化设计、简洁风、全面SEO优化、多重列表展示方式,响应式布局支持电脑、平板和手机的完美展示。
功能展示:清新·扁平:前端展示小清新,扁平化设计风格SEO优化:你能想到的各种SEO优化,这里都做到了速度优化:排除你的主机慢和安装超多插件,它的速度绝对给力图片和头像提速:文章列表和头像都使用了lazyload做按需异步加载,你会发现访问速度上提升一大截2种可选框架:支持3栏和2栏13种可选颜色风格:基本上小清新的颜色都有了3种列表模式:多图展示、单图展示、无图展示7种连接符:设置你自己的很重要文章来源:后台编辑文章有来源,前台有来源展示,尊重原创很重要新窗口打开文章:有这个开关,说多了都是2分享:让百度分享和网站更和谐置顶焦点图:绝美展示,可自行设置置顶文章,且第一章大图可自定义近期热门排行榜:最近多少天的评论热门文章排行,可设置天数和显示数量独立页面:点赞最多的文章排行导航多级菜单:不管你听不听,我是建议你最多二级就行
2025/3/30 19:14:54 75B WordPress主题 模板 源码 网站模版
1
点击缩略图查看大图,支持滑动切换图片,缩放图片,封装成了一个工具类,基本两行代码实现上面需求。
代码很清晰,可以自已微调。
2025/2/19 0:09:20 376KB ViewPager 缩放图片 查看大图
1
仿微信朋友圈和QQ控件点击图片查看大图功能,传入图片的集合和显示的索引,一行代码就能实现,使用开源库Photo和Glide;
可在项目中使用。
请大家多多指点;
2025/1/1 17:39:53 18.76MB 图片浏览
1
hbuilder打包APP后根据版本号更新,以及简单的点击图片查看大图
2024/12/2 6:05:48 10.54MB mui
1
plist格式大图拆分小图工具,合成分解容量优化等功能~
2024/11/30 22:39:24 1.26MB plist 合图 拆图 大图拆分小图
1
拆图神器-通过自动计算透明边界,把打包好的大图自动拆分开
2024/5/31 11:28:51 106KB 拆图神器 png 拆分图集
1
可以免plist文件把大图集直接逆向拆分成小图,透明通道拆分方法。
使用方法,鼠标拖图片进去,在图片区域点右键导出
2024/2/27 0:55:29 153KB 拆图 TexturePacke 逆向 plist
1
淘价宝是时尚商城门下开发人员又一大作,这款软件秉承方便,快捷,价格透明的设计理念;
我们努力为各大买家争取到最大的利益。
不要质疑,真的可以帮你省钱的哦。
软件没有花哨的界面,只有简简单单的一个窗口,加上一个按钮,但是功能绝对是强大的。
下面是安装淘价宝后我们浏览商品的画面,大家看到在右下角会有个小窗口弹出。
11(点击查看大图)点击一键比价就会转到同款宝贝的列表。
22(点击查看大图)这样我们就可以挑选出价格最低的卖家了。
在淘宝,大家都知道,同一款商品在不同的卖家那的价格往往是会有所出入的,其实都是同一个商品,这款软件就是协助买家发掘出最优惠的宝贝了。
这样就给买家省到不少钱呢,就按照每个商品少五块钱,一个月我们购买五次东西,那样一个月最少可以省下二十多块,绝对是居家旅行,淘宝购物必备。
软件只有不到200kb,不占电脑空间,运行内存超小,完全不影响用户正常使用电脑。
软件可以自动更新,所以用户永远不用担心会落后啦。
2024/2/13 23:25:04 67KB 比价软件
1
共 38 条记录 首页 上一页 下一页 尾页
在日常工作中,钉钉打卡成了我生活中不可或缺的一部分。然而,有时候这个看似简单的任务却给我带来了不少烦恼。 每天早晚,我总是得牢记打开钉钉应用,点击"工作台",再找到"考勤打卡"进行签到。有时候因为工作忙碌,会忘记打卡,导致考勤异常,影响当月的工作评价。而且,由于我使用的是苹果手机,有时候系统更新后,钉钉的某些功能会出现异常,使得打卡变得更加麻烦。 另外,我的家人使用的是安卓手机,他们也经常抱怨钉钉打卡的繁琐。尤其是对于那些不太熟悉手机操作的长辈来说,每次打卡都是一次挑战。他们总是担心自己会操作失误,导致打卡失败。 为了解决这些烦恼,我开始思考是否可以通过编写一个全自动化脚本来实现钉钉打卡。经过一段时间的摸索和学习,我终于成功编写出了一个适用于苹果和安卓系统的钉钉打卡脚本。
2024-04-09 15:03 15KB 钉钉 钉钉打卡