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

553825467

斑竹您好，借助您的代码，我想分析GPP对CO2的影响，于是我对每月的GPP与CO2进行了分析，其时间系数如下第一个图所示，接着我对co2进行了延迟处理，结果出来的时间系数比较好的表明了CO2逐年增加的特点，如下第二个图所示，可SVD不是需要进行去趋势的处理吗？为什么会出现延迟后co2趋势相反的现象呢？谢谢版主

Aires

553825467 发表于 2017-7-26 18:56
斑竹您好，借助您的代码，我想分析GPP对CO2的影响，于是我对每月的GPP与CO2进行了分析，其时间系数如下第一 ...

不知道为什么你的GPP是周期振荡的，你有没有去掉季节循环？你的CO2没有去趋势？

553825467

没有呢，就是直接用您的代码，计算2003年到2011年每个月的GPP对co2的影响，如何去掉季节循环以及去趋势呢？去季节循环是不是就是选择同一个季节进行计算呢？去趋势是不是就是先做拟合然后减去拟合后的值呢？谢谢您了。

Aires

553825467 发表于 2017-7-27 16:47
没有呢，就是直接用您的代码，计算2003年到2011年每个月的GPP对co2的影响，如何去掉季节循环以及去趋势呢？ ...

如果研究两个场自然变率（扣除人类强迫的影响）的相互作用，应该要去掉季节循环和趋势的。季节循环可是是每年的1月减去多年1月的平均，每年2月减2月的平均，以此类推。你说的去趋势的做法基本是对的。

553825467

奥，好的呢，我再做个去趋势和去季节。斑竹，对于风（U，V）如何做SVD呢？谢谢斑竹的耐心讲解。

Aires

553825467 发表于 2017-7-27 22:20
奥，好的呢，我再做个去趋势和去季节。斑竹，对于风（U，V）如何做SVD呢？谢谢斑竹的耐心讲解。

用U, V做和用你之前的两个变量又有什么区别呢？

553825467

Aires 发表于 2017-7-28 14:17
用U, V做和用你之前的两个变量又有什么区别呢？

因为风除了风速以外还有风向，这些都会影响CO2的空间分布，如何将风速和风向合并到一个场中去呢？

Aires

553825467 发表于 2017-7-28 16:20
因为风除了风速以外还有风向，这些都会影响CO2的空间分布，如何将风速和风向合并到一个场中去呢？

你可以试着做u, v, CO2的联合EOF，就是把这三个场都放到同一个矩阵里去做EOF，然后看看主模态是什么，如何解释。或者你可以先分析下是U还是V对CO2的影响大，然后用U或者V去和CO2做SVD，然后把得到的时间序列回归到U，V场和CO2场。你可以先初步分析下，再选取合适的统计方法进行处理。

553825467

Aires 发表于 2017-7-28 19:09
你可以试着做u, v, CO2的联合EOF，就是把这三个场都放到同一个矩阵里去做EOF，然后看看主模态是什么，如 ...

谢谢斑竹的耐心答复，我按照您的方法用eof先分析一下，再看看后面该怎么做

雅凝馨子

很想可是没贡献，呜呜呜呜