矩阵分解——谱分解

  • 28
    点赞
  • 134
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
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 = [];% 构造系数矩阵 for i = 1:m     for j = 1:n         X(i,j) = xt(n i-j);     end end beta = inv(X'*X)*X'*Y'; 复制代码 beta即为用最小二乘法估计出的模型参数。 此外,还有估计AR模型参数的Yule-Walker方程法、基于线性预测理论的Burg算法和修正的协方差算法等[26]。相应的参数估计方法在MATLAB中都有现成的函数,比如aryule、arburg以及arcov等。 4.3.3 AR模型阶次的选择及实验设计 文献[26]中介绍了五种不同的AR模型定阶准则,分别为矩阵奇异值分解(Singular Value Decomposition, SVD)定阶法、最小预测定误差阶准则(Final Prediction Error Criterion, FPE)、AIC定阶准则(Akaika’s Information theoretic Criterion, AIC)、MDL定阶准则以及CAT定阶准则。文献[28]中还介绍了一种BIC定阶准则。SVD方法是对Yule-Walker方程中的自相关矩阵进行SVD分解来实现的,在MATLAB工具箱中arorder函数就是使用的该算法。其他五种算法的基本思想都是建立目标函数,阶次估计的标准是使目标函数最小化。 以上定阶准则在MATLAB中也可以方便的实现,下面是本文实现FPE、AIC、MDL、CAT定阶准则的程序(部分): for m = 1:N-1    ……       % 判断是否达到所选定阶准则的要求    if strcmp(criterion,'FPE')        objectfun(m 1) = (N (m 1))/(N-(m 1))*E(m 1);    elseif strcmp(criterion,'AIC')        objectfun(m 1) = N*log(E(m 1)) 2*(m 1);    elseif strcmp(criterion,'MDL')        objectfun(m 1) = N*log(E(m 1)) (m 1)*log(N);    elseif strcmp(criterion,'CAT')        for index = 1:m 1            temp = temp (N-index)/(N*E(index));        end        objectfun(m 1) = 1/N*temp-(N-(m 1))/(N*E(m 1));    end        if objectfun(m 1) >= objectfun(m)        orderpredict = m;        break;    end end 复制代码 orderpredict变量即为使用相应准则预测的AR模型阶次。 (注:以上代码为结合MATLAB工具箱函数pburg,arburg两个功率估计函数增加而得,修改后的pburg等函数会在附件中示意,名为pburgwithcriterion) 登录/注册后可看大图 程序1.JPG (35.14 KB, 下载次数: 20352) 下载附件  保存到相册 2009-8-28 20:54 上传 登录/注册后可看大图 程序2.JPG (51.78 KB, 下载次数: 15377) 下载附件  保存到相册 2009-8-28 20:54 上传
非负矩阵分解(Non-negative Matrix Factorization,NMF)是一种在机器学习和数据挖掘领域广泛应用的算法。其核心思想是将一个非负矩阵分解为两个非负的低秩矩阵的乘积,从而发现原始矩阵的潜在结构和特征。 在MATLAB中,可以通过调用内置函数或开源工具箱来实现非负矩阵分解。一种常用的方法是使用MATLAB自带的函数`nnmf`。这个函数可以通过多种不同的算法实现非负矩阵分解,如乘法更新法、反射法等。 使用`nnmf`函数可以指定要分解的原始矩阵以及分解后的矩阵的秩。此外,还可以选择不同的代价函数和优化算法以获得更好的分解结果。函数的输出为分解后的两个非负矩阵及相应的误差。 除了`nnmf`函数外,MATLAB中还有其他一些工具箱和函数可以实现非负矩阵分解,如NMF_ARD、NMF库等。这些工具箱和函数提供了更多的参数选项和算法选择,可以根据实际需求进行调整和优化。 在使用非负矩阵分解时,需要注意选择合适的参数和算法,以及对分解结果进行合理的解释和应用。同时,也要注意数据的预处理和归一化,以提高非负矩阵分解的效果和可解释性。 总之,MATLAB提供了多种实现非负矩阵分解的函数和工具箱,可以根据需求选择合适的方法进行分解和应用。非负矩阵分解在机器学习、图像处理、文本挖掘等领域有广泛的应用前景。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值