求解EOF所需的距平场和标准化场

shijian0819

言深深 发表于 2014-4-16 11:09
貌似乱码了···
你想做什么呢？对数据进行eof分解，还是分析结果？

是这个样子的，比如我导出一个SST的原始场数据，同一个程序，我发现原始场和距平能运行程序，而如果ks=1，也就是标准化的时候，会提示有run time error，sqrt之类的错误。。也就是会报错，应该是jacobi方法求特征值和特征向量那一步出错了~~~~  不知道为何，同一个数据。。之前也发现了。。选不同的ks值有的能出来结果，有的出不来~

shijian0819

言深深 发表于 2014-4-16 11:09
貌似乱码了···
你想做什么呢？对数据进行eof分解，还是分析结果？

是这个样子的，我刚才做了一个实验。我用导出的SST原始场的数据，选择不同的ks值，ks=0，和-1的时候都能瞬间出结果。可是ks=1，也就是要标准化原始场的时候，程序报错了。具体如图。麻烦亲指导一下~~~~

shijian0819

言深深 发表于 2014-4-16 11:09
貌似乱码了···
你想做什么呢？对数据进行eof分解，还是分析结果？

不知道为什么点高级模式总是看不到我回复您的帖子。那就简单回复一下。我利用我导出的原始场的SST数据，选择不同的ks值，选择ks=-1和ks=0都可以秒出结果。但是选择标准化的时候，ks=1，就提示在jacobi那一步有一个run time error，好像是sqrt函数报错了~~ 应该是求出的A矩阵出现了错误。。表示不理解，求指导~~~

littlen9ne

言深深 发表于 2012-11-28 20:31
做eof场距平是但看空间的，也就是你行······标准化亦然

如果是对单一的时间序列进行标准化/距平， ...

你好，我想问一下，EOF的距平场是空间的距平场还是时间的距平场？
我做了原场的标准化EOF，就是ks=1的情况，但是我将k改成0就陷入死循环里了。
所以我后来想直接把原场数据换成距平场数据，来做标准化EOF，可是出来的规律跟原场差不多相似。

看了这个帖子之后，我想说我的距平是时间的距平。

所以EOF的距平场是空间的距平场吗？


还有一个问题，就是我的EOF之后标准化序列是可以出结果的，另外两个ks=0和ks=-1是初步来结果的。

言深深

littlen9ne 发表于 2014-4-16 14:07
你好，我想问一下，EOF的距平场是空间的距平场还是时间的距平场？
我做了原场的标准化EOF，就是ks=1的情 ...

你问的很好，这个问题是很多人模棱两可都带过去的。
距平场，这里说的是时间的距平，你看一下对ks判断并距平处理的程序段就能看出来了。
计算时间距平的目的是去除气候态，计算空间距平在这里是没有意义的。

言深深

littlen9ne 发表于 2014-4-16 14:07
你好，我想问一下，EOF的距平场是空间的距平场还是时间的距平场？
我做了原场的标准化EOF，就是ks=1的情 ...

应该是时间距平。
ks选取决定程序运算情况不同是不应该的，你的情况比较罕见，一般都是能出来的。出不来的原因我估计是你的矩阵出现奇异。不知道你是否存在缺测数据，并对你的缺测数据经行剔除？这个程序是不支持缺测数据运算的···

言深深

shijian0819 发表于 2014-4-16 12:17
不知道为什么点高级模式总是看不到我回复您的帖子。那就简单回复一下。我利用我导出的原始场的SST数据， ...

ks选择不同，应该都是能够出来结果的。出不来的原因可能是你的矩阵出现奇异，检查一下是否有缺测，并对缺测进行剔除。

言深深

shijian0819 发表于 2014-4-16 12:14
是这个样子的，我刚才做了一个实验。我用导出的SST原始场的数据，选择不同的ks值，ks=0，和-1的时候都能 ...

单单从这个是看不出来的，这里的程序意思如你所言，是说求解过程中出现根号为负一类的情况。我在另一个帖子回复说可能出现的缺测情况，请思考一下看看···

littlen9ne

言深深 发表于 2014-4-16 15:40
单单从这个是看不出来的，这里的程序意思如你所言，是说求解过程中出现根号为负一类的情况。我在另一个帖 ...

数据上没有缺省值，但是怎么说呢，换成矩阵的话，第一列都是0。
数据确实很小，精度上确实有可能被忽略掉。那我先试着放大数据做着试试。

言深深

littlen9ne 发表于 2014-4-16 16:25
数据上没有缺省值，但是怎么说呢，换成矩阵的话，第一列都是0。
数据确实很小，精度上确实有可能被忽略 ...

是的，应该考虑到这种清理。已经下载并查看你的数据，确实太小了···