如何在Matlab中从有序一维双精度数组移除元素子集？

给定有序数组 A = [1,3,5,6,7] 和有序子集 a = [3,6]，要得到 A_a = [1,5,7]，以下是几种高效实现方式，避免低效的循环遍历：

一、通用简洁方法（支持有序/无序数组） #

直接使用ismember函数判断元素是否属于子集，取反筛选即可：

A = [1,3,5,6,7];
a = [3,6];
A_a = A(~ismember(A, a));

该方法无需关注数组是否有序，代码简洁易读，能覆盖大多数常规场景。

二、利用有序性的优化方法（更高性能） #

由于两个数组均为有序状态，可给ismember添加sorted参数，让Matlab采用二分查找算法匹配元素，时间复杂度从默认的O(n*m)降至O(n log m)，在数组规模较大时性能提升显著：

A_a = A(~ismember(A, a, 'sorted'));

三、手动双指针实现（极致效率） #

若追求最高性能，可手动实现双指针线性扫描，完全利用数组的有序性，时间复杂度为O(n+m)：

A = [1,3,5,6,7];
a = [3,6];
n = length(A);
m = length(a);
i = 1; j = 1;
A_a = [];

while i <= n && j <= m
    if A(i) < a(j)
        A_a = [A_a, A(i)];
        i = i + 1;
    elseif A(i) == a(j)
        i = i + 1;
        j = j + 1;
    else
        j = j + 1;
    end
end
// 追加A中剩余未处理的元素
A_a = [A_a, A(i:end)];

这种方法彻底避免了重复比较，适合处理超大数组的场景。

为什么不推荐循环朴素方法？ #

你提到的循环执行A_a = A(A~=a(i))的问题在于：每次循环都会重新遍历整个A数组，时间复杂度为O(k*n)（k为子集a的长度），且每次赋值都会创建新数组，内存开销更大，效率远低于上述方法。

内容的提问来源于stack exchange，提问作者Argyll