Matlab:循环数组段

Question

我有2 m x n个矩阵，并且想要计算给定窗口长度内的分段的典型相关性。例如，如果我的窗口长度是100，我希望得到

canoncorr(X(1:100,:),Y(1:100,:))
canoncorr(X(101:200,:),Y(101:200,:))
canoncorr(X(201:300,:),Y(201:300,:))
...

以此类推，所有这些都累积到一个矩阵中。我只对相关系数r感兴趣。

我正在尝试以下操作：

win=100;
r=[];
for i=1:win:size(X,1)-win-2
    [A,B,r(i,:)] = canoncorr(X(i:i+win,:),Y(i:i+win,:));
end

但是，我得到的矩阵不仅保存了行1、101、201等的值，而且还用零填充了1到101之间的行。

如果我尝试

 [A,B,r(i:i+win,:)] = canoncorr(X(i:i+win,:),Y(i:i+win,:));

则下标分配维度不匹配。

我做错了什么？

Answer 1

i从1,101,201，...

所以，请检查第101行，看看它们是否为零。

您还可以使用spy(r)查看矩阵元素。