爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4603|回复: 0

[源程序] MATLAB中前N个特征值对应的时间系数是哪个变量呢?

[复制链接]

新浪微博达人勋

发表于 2014-5-12 12:29:04 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 栤唥de嶶笑 于 2014-5-12 12:31 编辑


我要找的是前几个特征值对应的时间系数,但哪个才是呢?
方法一;
运用svd函数进行分解
[u,s,v]=svd(x)
其中 u为x的空间函数
pc=s*v'为x的时间函数
s为x分解的对角阵
在计算主分量方差奉献时需要X*X'的对角阵。
因此还需要再运行一次svd函数
[u1,ss,v1]=svd(x*x')
则ss为x*x'的对角阵,可用来计算方差奉献。

方法二;
运用eig函数进行分解
[eof,e]=eig(x*x')
其中eof为x的空间函数
eof'*x为x的时间函数
pc=eof'*x
e为x的特征值组成的对角阵
这里需要注意的是e的对角阵不是按大到小排列的。因此若要使对角阵e的对角元素从大到小排列的话。时间场和空间场都应该有相应的变化。具体变化如下
eof=fliplr(eof)      最后的空间场      掉换矩阵左右两侧的数据
pc=flipud(pc)掉换矩阵上下两端的数据
e=fliplr(flipud(e))  最后按对角元素大到小排列的对角阵与eof分解的完全一致
pc=flipud(pc)        最后的时间场
其中e为*x*x'的对角阵,eof为x的空间分量,pc为时间分量。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

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