【讨论】快速删除矩阵中含0的行。

rocwoods

如题！
分两种情况，一种是某行只要含有0就去除，另一种是全部为0才去除。
欢迎大家踊跃发言，看看有多少种方法。

rocwoods

比较容易想到的方法是用循环实现，利用any all等函数性质，但是这样比较麻烦。
有更好的方法。

eight

原帖由 rocwoods 于 2006-12-1 15:13 发表
比较容易想到的方法是用循环实现，利用any all等函数性质，但是这样比较麻烦。
有更好的方法。

第一种：A(any(~logical(A),2),:)=[];
第二种：A(all(~logical(A),2),:)=[];

rocwoods

我想的一种方法没有用all和any
前一种情况:
a1=a(sum(~a,2)==0,:)
后一种情况:
a2=a(sum(~a,2)<n,:) (n为矩阵列数)

eight兄的方法可以把logical去掉。
直接A(any(~A,2),:)=[];
A(all(~A,2),:)=[];
要不会有warning，虽然不是大事。

chengweifeng

不另开新贴，就此向大家请教一个问题吧，问题简化后如下：

1. I=[2 11 10 3 12]；
   
2. Element=[1 2 3 4 5 6 7 8 9];

复制代码

现在我想删除Element中Element和 I 相同的元素2，3 大家有什么好的方法？

[ 本帖最后由 chengweifeng 于 2007-10-24 01:56 编辑 ]

rocwoods

1. intersect(setxor(I, Element),Element)
   
2. 或者
   
3. [a,b,c]=intersect(I,Element);
   
4. Element(c)=[]
   
5. 或者setdiff(Element,intersect(I,Element))
   

复制代码

[ 本帖最后由 rocwoods 于 2007-10-24 10:39 编辑 ]

chengweifeng

That's perfect!
Thank you very much.

huaijuliu

这个能否继续讨论？

lin2009

6# rocwoods

intersect(setxor(I, Element),Element)
或者
[a,b,c]=intersect(I,Element);
Element(c)=[]
或者setdiff(Element,intersect(I,Element))

根据集合差运算规则和setdiff的用法，直接用setdiff即可。
setidff(I, Element)

lin2009

1# rocwoods

一种是某行只要含有0就去除

可以用：
[rows,cols,vals] = find(A==0);
A(rows,:) = []

ajun2100

楼上“rows,cols,vals] = find(A==0);
A(rows,:) = []”
牛~~~

wpaul

楼上“rows,cols,vals] = find(A==0);
A(rows,:) = []”
受益匪浅。