分享一个自己编写的SVD函数

雅凝馨子

Aires 发表于 2017-9-6 15:27
你可以选取该流域降水的指数（比如这个区域的空间平均），然后和全球每一个格点的海温求相关，得到一个相 ...

您说的区域的空间平均，是这个流域空间内所有降水的平均值吗？假设我要求月平均的，那这样得到的是12个元素的一维数组，全球每一个格点的月平均海温也是12个元素的一维数组，然后进行SVD计算吗？

Aires

雅凝馨子 发表于 2017-9-6 16:27
您说的区域的空间平均，是这个流域空间内所有降水的平均值吗？假设我要求月平均的，那这样得到的是12个元 ...

你不是求相关系数么？那还SVD干什么？请先自己把相关和SVD这些方法搞明白。

雅凝馨子

Aires 发表于 2017-9-7 07:38
你不是求相关系数么？那还SVD干什么？请先自己把相关和SVD这些方法搞明白。

抱歉，是我的语言表述不对。
求的是两个场之间的相关系数，一个大区域和一个小区域，计算方法是用SVD方法。
在实现的过程中，遇到的问题是，两个场大小不一样，不知道该怎么处理？之前您说用指数，然后求它和格点求相关，可这并没有用到SVD。
现在我想请教的是如何用SVD处理这样两个大小不一致的空间场，得到相关系数。

Aires

雅凝馨子 发表于 2017-9-7 14:57
抱歉，是我的语言表述不对。
求的是两个场之间的相关系数，一个大区域和一个小区域，计算方法是用SVD方 ...

我不知道。你可以用SVD得到两个场主模态对应的时间系数，这两个时间系数的相关在一定程度上也能代表两个场的耦合强度的。

雅凝馨子

Aires 发表于 2017-9-7 16:42
我不知道。你可以用SVD得到两个场主模态对应的时间系数，这两个时间系数的相关在一定程度上也能代表两个 ...

因为是非专业，看了很久的资料也不懂主模态，不知道怎么找到两个场的主模态，尤其是大小范围还不一样的两个场

夏夜

Aires 发表于 2016-1-3 15:33
You can calculate the autocorrelation and estimate the decorrelation timescale (i.e. the first tim ...

楼主，这段解释计算DOF的有没有参考文章呀，引用下

雅凝馨子

老师，你好。我根据您的程序，算了一些系数结果，有两个问题：
一：最后的各种系数，都是一维数组，这样是对的吗？
二：另外请教您，我画系数相关图，要怎么画呢？因为我是reshape原始数据（180*89，12），而且去除了无效值，数据个数减少了，最后一维数组是（1*19088），不能再reshape回去类似原始数据的经纬度分布了，这样子，要怎么做图呢？

雅凝馨子

@Aires 雅凝馨子 发表于 2017-9-11 23:49
老师，你好。我根据您的程序，算了一些系数结果，有两个问题：
一：最后的各种系数，都是一维数组，这样是 ...

@mofangbao @Aires

Aires

雅凝馨子 发表于 2017-9-11 23:49
@mofangbao @Aires

二、没必要挑出NaN的值，做回归或者相关的时候算出来的值就是NaN，出来的就是矩阵直接画图就好了。

最后一个匈奴

感谢楼主分享！！！