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

雅凝馨子 · 发表于 2017-9-12 17:17:21

本帖最后由雅凝馨子于 2017-9-12 17:27 编辑

Aires 发表于 2017-9-12 13:09
做SVD的时候是要把NaN挑出来啊。SVD之后得到的时间系数和原场做相关或回归的时候就不需要了。

现在的困惑是 A=X*Y',[U，S,V]=svd(A),Matlab 报错说svd分解不能含有NaN和Inf。然后我挑出了X和Y里的NaN和Inf，之后得到了一些相关系数，想作图，可是因为挑出了nan, 作图的话，不知道怎么做。
原本经纬度数据是180*89，现在系数个数是（1*19088），怎么把它和经纬度相对应起来作图呢？

Aires · 发表于 2017-9-13 09:16:40

雅凝馨子发表于 2017-9-12 17:17
现在的困惑是 A=X*Y',=svd(A),Matlab 报错说svd分解不能含有NaN和Inf。然后我挑出了X和Y里的NaN和Inf， ...

你做SVD的时候当然是要挑出NaN了，但是你画图的时候要把NaN位置的数据原封不动的放进去啊。

雅凝馨子 · 发表于 2017-9-13 11:47:52

Aires 发表于 2017-9-13 09:16
你做SVD的时候当然是要挑出NaN了，但是你画图的时候要把NaN位置的数据原封不动的放进去啊。

哦哦，懂啦，谢谢老师耐心解答，麻烦了，非常感谢您

云雾里 · 发表于 2017-9-13 14:59:37

请教楼主一个检验问题，自己看了一些资料后还是不太明白，做模态检验的时候比如说SVD做了十个模态，那么蒙特卡洛检验就有1000个值，要怎么确定哪个值是显著性临界值呢？

Aires · 发表于 2017-9-13 19:12:03

云雾里发表于 2017-9-13 14:59
请教楼主一个检验问题，自己看了一些资料后还是不太明白，做模态检验的时候比如说SVD做了十个模态，那么蒙 ...

抱歉，我没有做过蒙特卡洛检验。

想个名字 · 发表于 2017-10-7 08:01:02

{:eb502:}{:eb502:}{:eb502:}{:eb502:}

liuyg123 · 发表于 2017-11-3 12:53:09

感谢楼主分享

liuyg123 · 发表于 2017-11-7 21:29:54

楼主您好想问一个比较幼稚的问题。。您在读取海温的时候为什么要从第1402个开始读取？

553825467 · 发表于 2017-11-28 00:09:32

本帖最后由 553825467 于 2017-11-28 00:29 编辑

老师您好，借用您的代码，我分析了9年108个月GPP对CO2 的影响，通过去趋势后进行SVD分析，得到的第一模态的相关系数如图所示。我参考了一些文献对时间系数的分析，认为如果时间系数的正负发生转变则说明在这这个时间点上发生了突变，但是这种分析方式对于去趋势后的数据是不是也成立呢？时间系数图像中的正负代表的是什么含义呢？最大小值是否代表的就是极值点呢？针对这个时间系数图我是否可以这么分析呢？
（1）CO2变化滞后于GPP的变化两个月左右
（2）CO2的极值点在4、8月份，而GPP在3、8月份取得极值
（3）CO2与GPP年、季循环的特征
谢谢老师了

yf9527 · 发表于 2017-12-11 16:20:03

下载没反应啊

		自动登录	找回密码
密码			立即注册

[源程序] 分享一个自己编写的SVD函数