删除满足特定条件的单元格

Question

我有一个单元阵列，每个单元是一个n×n矩阵。 我想删除inv(cell{i}'*cell{i})发出矩阵接近单数的警告的inv(cell{i}'*cell{i}) 。 谢谢

Answer 1

通常，删除元素是很容易的部分。 如果C是您的数组，则可以通过以下方式删除向量idx索引指定的单元格：

C(idx) = {};

关于你的具体问题，检查矩阵是否“几乎”单数或不是rcond （如果结果接近零，它可能是单数）。 要将其应用于所有单元格，您可以通过以下方式使用cellfun ：

idx = cellfun(@(x)(rcond(x' * x) < 1e-12), C);

根据您的喜好调整阈值。 得到的idx是一个逻辑数组，在奇异矩阵的位置有1 s。 使用idx从C删除这些元素，如上所示。

Answer 2

创建一个检查您的条件的函数：

function state = CheckElement(element)

if ([condition])
   state = 1;
else
   state = 0;
end

end

然后对所有单元格数组元素执行cellfun，如下所示：

indices = cellfun(@CheckElement,myCellArray);
cellArray(indices ) = [];

Answer 3

假设您已经定义了一个定义的issingular函数，您可以使用cellfun来获取包含要删除的矩阵的单元格的索引。

c; % cell arry

singularIdx = cellfun((@x) issingular( inv(x' * x)), c);  %' added single quote for SO syntax

cFiltered = c(~singluarIdx);

您可能需要编写自己的函数来检查奇点，有关详细信息，请参阅此问题： 检查矩阵是否为单数的快速方法？ （不可逆，det = 0）