使用
bsxfun创建正确的索引应该肯定会有所帮助：
ind = bsxfun(@plus, 1:W, (0:X:numel(T)-W).');
out = T(ind);
创建正确的索引是第一步,由第一行代码描述.这段代码的作用是创建一个2D矩阵,其中每一行都是每个感兴趣的窗口访问的元素.如果你想直观了解代码如何生成索引,请特别注意第一种情况,即X = 1;和W = 3;
我们可以看到第一行包括访问元素1,2,3.第二行包括访问元素2,3,4 ……直到最后一行,即5,6,7.我们可以看到我们必须访问窗口中的相邻元素,因此基本索引需要从1,2,3或一般从1到W.我们现在需要偏移这些索引,以便它们以正确的元素为中心每个窗口T.第一个窗口的偏移量仅为0,第二个窗口的下一个偏移量仅为1,直到最后一行为3.我们看到,对于每一行,随着行的增加,我们再向基本索引添加1个.因此,我们为第二行的每个基本索引添加1,然后为第三行中的每个基本索引添加2,依此类推.如果使用偏移索引添加基本索引,则最终获得正确的索引以访问T中的正确元素.
同样,如果X = 2;并且W = 3;,我们看到我们仍然有1,2,3的基本索引.但是,现在访问的正确元素是第一行的1,2,3,然后是第二行的3,4,5然后是第5排,第6排,第7排.对于每一行,我们现在将基本索引偏移2而不是1.因此,第二行我们为每个基本索引添加2,然后我们为第三行的每个基本索引添加4,依此类推.
通常,使用向量1：W创建基本索引,并且使用向量0：X：numel(T)-W创建偏移索引.需要减去W,以便在按照要求对信号进行采样时不会超出范围.为了创建我们刚才谈到的这些索引,bsxfun为我们处理了这个问题.
我们创建一个1：W的行向量,它对应于基本索引和列向量(0：X：numel(T)-W).它对应于每个窗口的偏移量.请注意,第一个偏移从0开始,然后我们增加X个量以确保计算正确的中心以将基本索引置于.我们停止直到我们击中了数字(T)-W元素,这是你说过的条件.通过使用bsxfun,创建了两个临时2D矩阵,其中行向量被复制的行数与列向量中的行一样多,并且列向量与行向量中的列一样多.将这两个矩阵一起添加后,就可以得到结果索引矩阵.
运行W = 3的代码;和X = 1;得到：
>> T = [1, 5, 6, 8, 10, 14, 22];
>> X = 1;
>> W = 3;
>> ind = bsxfun(@plus, 1:W, (0:X:numel(T)-W).')
ind =
1 2 3
2 3 4
3 4 5
4 5 6
5 6 7
同样,如果W = 3;和X = 2;我们还得到：
>> T = [1, 5, 6, 8, 10, 14, 22];
>> X = 2;
>> W = 3;
>> ind = bsxfun(@plus, 1:W, (0:X:numel(T)-W).')
ind =
1 2 3
3 4 5
5 6 7
您可以自己验证这些索引是否与T中的正确元素相对应,以便在这种情况下创建所需的矩阵.
我们最终使用它来索引我们的矩阵以获取正确的元素：
out = T(ind);
为X = 1执行此操作;和W = 3;得到：
>> out = T(ind)
out =
1 5 6
5 6 8
6 8 10
8 10 14
10 14 22
同样,X = 2;和W = 3;得到：
>> out = T(ind)
out =
1 5 6
6 8 10
10 14 22