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

海小宝

学长 我想问你第一对左右奇向量图，同性相关系数图和异性相关系数图分别代表什么物理意义

Aires

海小宝 发表于 2015-12-27 15:10
学长 我想问你第一对左右奇向量图，同性相关系数图和异性相关系数图分别代表什么物理意义

你想想相关的含义就知道了

dreamdong

找标准化方法的时候看到了这个帖子，很厉害呢。学习了，等我有了贡献一定来下。

mikeflm

楼主，你好！
我想问一下，我使用你给的程序时，rXU_jian/rYV_jian和rXV_jian/rYU_jian的检验为什么总是零，我用的资料是1982~2011年1月的SST和风场做相关性分析。
希望能帮小弟解答一下。

Aires

mikeflm 发表于 2016-1-3 10:34
楼主，你好！
我想问一下，我使用你给的程序时，rXU_jian/rYV_jian和rXV_jian/rYU_jian的检验为什么总是零 ...

你看看你相关系数大概有多少吧，小的话一般就不痛过检验，还有你的自由度，也有可能你SST和风场关系太弱

mikeflm

Aires 发表于 2016-1-3 11:12
你看看你相关系数大概有多少吧，小的话一般就不痛过检验，还有你的自由度，也有可能你SST和风场关系太弱

能再帮我具体看看吗？
[L,R,U,V,SCF,r_alpha,rUV,rUV_jian,rXU,rXU_jian,rYV,rYV_jian,...
    rXV,rXV_jian,rYU,rYU_jian] = SVD_method(x,y,-9999,0.8,28,0.05);
这是我设置的，
r_alpha=0.36；rUV=[0.885   0.684   0.841];
30年的数据自由度不是30-28吗？

Aires

mikeflm 发表于 2016-1-3 12:11
能再帮我具体看看吗？
[L,R,U,V,SCF,r_alpha,rUV,rUV_jian,rXU,rXU_jian,rYV,rYV_jian,...
    rXV,rXV ...

30年的数据自由度不是28啊，你要看你数据的decorrelation timescale，你看看我程序有没有问题？你的结果应该通过检验的

mikeflm

学长，我也是海大的；看在同门师弟的份上，能说一下自由度怎么算吗？
比如30年1月的月平均数据。
谢谢。

Aires

mikeflm 发表于 2016-1-3 14:54
学长，我也是海大的；看在同门师弟的份上，能说一下自由度怎么算吗？
比如30年1月的月平均数据。
谢谢。

You can calculate the autocorrelation and estimate the decorrelation timescale (i.e. the first time that the correlation coefficient drops to zero). Decorrelation timescale is roughly represented as the 1/4 period. Half period represents one degree of freedom. Thus, two times of decorrelation timescale represents one DOF. Therefore, you can based on this timescale to estimate the DOF of your time series. Notice that the DOF of decadal variability in the ocean-atmosphere system will significantly decrease due to their long periods.

mikeflm

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

非常感谢。