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

Aires

雅凝馨子 发表于 2017-8-9 22:27
想请教楼主，那个相关系数的等值线图是怎么绘制的呀，最近一直被这困扰。。。。

contour

雅凝馨子

没有贡献啊啊啊┭┮﹏┭┮

雅凝馨子

楼主能发我一份吗？没有贡献，下不了，一直被困扰。。

雅凝馨子

Aires 发表于 2017-8-9 23:29
contour

我用contour画出来的很奇怪，因为不知道怎么处理那些nan，又是非专业的，一直不成功。想请问下楼主，那些nan值怎么处理呢？

Aires

雅凝馨子 发表于 2017-8-11 10:45
我用contour画出来的很奇怪，因为不知道怎么处理那些nan，又是非专业的，一直不成功。想请问下楼主，那些 ...

去查一下Matlab的书吧。你做SVD之前可以用isnan函数把不是NaN（即有效的数据）挑出来再做SVD。

雅凝馨子

Aires 发表于 2017-8-11 11:28
去查一下Matlab的书吧。你做SVD之前可以用isnan函数把不是NaN（即有效的数据）挑出来再做SVD。

用了isnan函数，剔除后，不知道怎么重组数组，用reshape的话，会把原来的经纬度和对应的数据顺序都弄乱呀。
另外有个疑问，想问下老师，海温场，降水场都是三维的，我看论坛里很多svd程序都是二维的场，请问三维是怎么变成二维的呢？假设是180*89*12，直接用reshape变成16020*12吗？这样的话，顺序不是都乱了吗？

Aires

雅凝馨子 发表于 2017-8-11 11:33
用了isnan函数，剔除后，不知道怎么重组数组，用reshape的话，会把原来的经纬度和对应的数据顺序都弄乱呀 ...

SVD或者EOF本身就是找一个场的主模态。你用reshape将[lon lat time]的场变换为[lon*lat time]，才能对二维的矩阵做SVD或EOF，做出来的模态（特征向量）也会是lon*lat，你再用reshape变为[lon lat]就可以了，并没有对数据分析有任何影响。你可以想象原先每个时刻是一个二维场，用reshape只是变为一维场处理了而已，画图时再变为二维矩阵画就是了。

雅凝馨子

Aires 发表于 2017-8-11 11:44
SVD或者EOF本身就是找一个场的主模态。你用reshape将[lon lat time]的场变换为[lon*lat time]，才能对二 ...

我的理解是：假如数据中有一处 sst(12,3,1)=nan, 剔除它后，重组数据，那样sst (36,1)处会有一个非nan的数，这样子，是相当于在该处填补了一个新的数，如果这样的话，就不对吧

Aires

雅凝馨子 发表于 2017-8-11 12:47
我的理解是：假如数据中有一处 sst(12,3,1)=nan, 剔除它后，重组数据，那样sst (36,1)处会有一个非nan的 ...

你应该reshape成sst(36,1)后再挑出NaN的位置。

雅凝馨子

Aires 发表于 2017-8-11 14:05
你应该reshape成sst(36,1)后再挑出NaN的位置。

老师，前几天一直在忙别的事，忘记看回复了。如果重组后再剔除，根据SVD的算法，要求两个矩阵中一个行数和另一个列数相等，两个场的nan不同，重组剔除后，怎么计算呢？