有许多方法可以在没有循环的情况下完成此操作.大多数解决方案涉及将向量x和y扩展为更大的索引矩阵,并且可能使用函数
REPMAT,
BSXFUN或
SUB2IND中的一个或多个.可以在
here找到用于矩阵索引的良好教程.
但是,既然你要求一个优雅的解决方案,这里有一点不寻常.它使用anonymous functions以及函数ARRAYFUN和CAT:
indexFcn = @(r,c) cam1(r:(r+99),c:(c+99));
result = arrayfun(indexFcn,x,y,'UniformOutput',false);
result = cat(3,result{:});
说明:
第一行创建一个匿名函数.这是一个简单的单行函数,可以在运行时创建,而无需将其放在m文件中.该函数定义了两个输入r和c,用于从cam1中提取100×100的子矩阵.变量indexFcn存储用于调用该函数的function handle.请注意,匿名函数使用的cam1值是静态的.即使变量cam1中的值发生更改,匿名函数仍会使用创建函数时cam1中的值.
第二行调用ARRAYFUN,它将函数应用于数组的每个元素. ARRAYFUN循环遍历x和y中的每个条目,将值传递给indexFcn.输出存储在结果中,这是一个15元素的单元阵列,其中每个单元格包含一个100乘100的矩阵.
第三行使用CAT函数将100×100矩阵连接成100×100×15的矩阵.