MATLAB做EOF分析，距平没有NaN，距平求得的协方差出现NaN是为什么啊

默尘

aeof=double(ap*ap'/nt);    nt是时间长度，不为零，ap是距平，已经剔除缺测了，但是做协方差时求得的aeof中又出现缺测是为什么呢，求指导

772615260

程序贴多一点，按道理，ap非Nan或inf，aeof也不会出现，你的eig求的是哪个矩阵

默尘

772615260 发表于 2018-5-28 22:49
程序贴多一点，按道理，ap非Nan或inf，aeof也不会出现，你的eig求的是哪个矩阵

是的，我仔细看了，我用
……
z=find(~isnan(ac(:,1)));
ap=ac(z,:);
……
剔除缺测值,但是没有剔除干净，还是有一些NaN，怎么办呢

默尘

772615260 发表于 2018-5-28 22:49
程序贴多一点，按道理，ap非Nan或inf，aeof也不会出现，你的eig求的是哪个矩阵

eig求的是aeof:     [V,D]=eig(aeof);

772615260

……
z=find(~isnan(ac(:,1)));
ap=ac(z,:);
……
你这句只是针对ac的第一列是否有nan进行了判断，看你ap的公式又是针对ac的所有列，因为不清楚你的数据的具体格式，如果只是单纯去掉所有nan的行的话，可以：
[z,z1]=find(isnan(ac)==1);
ap=ac;ap(m,:)=[];

默尘

772615260 发表于 2018-5-31 17:33
……
z=find(~isnan(ac(:,1)));
ap=ac(z,:);

ap=ac;ap(m,:)=[];这句看不懂，m是什么？

默尘

772615260 发表于 2018-5-31 17:33
……
z=find(~isnan(ac(:,1)));
ap=ac(z,:);

ap=ac;ap(m,:)=[];这句看不懂，m是什么？

默尘

772615260 发表于 2018-5-31 17:33
……
z=find(~isnan(ac(:,1)));
ap=ac(z,:);

ap=ac;ap(m,:)=[];这句看不懂，m是什么？

m0932

最好可以把测试的数据和程序贴上来，这都是在猜

oyhx_cas

感谢楼主！