分享一个自己编写的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，出来的就是矩阵直接画图就好了。

最后一个匈奴

感谢楼主分享！！！

雅凝馨子

Aires 发表于 2017-9-12 07:42
二、没必要挑出NaN的值，做回归或者相关的时候算出来的值就是NaN，出来的就是矩阵直接画图就好了。

如果不挑出NaN,  svd 就无法分解，这要怎么办呢

yingnars

感谢楼主分享，很好用，{:5_196:}{:5_196:}{:5_196:}

Aires

雅凝馨子 发表于 2017-9-12 09:46
如果不挑出NaN,  svd 就无法分解，这要怎么办呢

做SVD的时候是要把NaN挑出来啊。SVD之后得到的时间系数和原场做相关或回归的时候就不需要了。