matlab find函数

Posted 2023-04-18

x1=-2:2;
k1=-2:2;
x2=[1,-1,1];
k2=-1:1;
k=min([k1,k2]):max([k1,k2]);
f1=zeros(1,length(k));
f2=zeros(1,length(k));
f1(find((k>=min(k1))&(k<=max(k1))==1))=x1;
f2(find((k>=min(k2))&(k<=max(k2))==1))=x2;
f=f1+f2;
stem(k,f,'filled');
axis([min(min(k1),min(k2))-1,max(max(k1),max(k2))+1,min(f)-0.5,max(f)+0.5]);
其中f1(find((k>=min(k1))&(k<=max(k1))==1))=x1;是什么意思？

举个例子讲下find函数的用法：

例如A=[1 2 3;4 5 6;7 8 9]

    1 2 3

    4 5 6

    7 8 9

[m n]=find(A<4)

将m和n合并看[m n],为：

     1     1

     1     2

     1     3

即(1,1),(1,2),(1,3)

也就是说A矩阵中(1,1),(1,2),(1,3)处元素小于4.

find函数是找出矩阵A中满足一定条件的元素，返回的为这些元素的行列标号，分别将这些满足要求的元素的行标号存于m向量中，满足要求的元素的列标号存于向量n中，且m和n存储的元素标号是对应的。

参考技术A 功能：
找到非零元素的索引和值
语法：
1. ind = find(X)
2. ind = find(X, k)
3. ind = find(X, k, 'first')
4. ind = find(X, k, 'last')
5. [row,col] = find(X, ...)
6. [row,col,v] = find(X, ...)
说明：
1. ind = find(X)
找出矩阵X中的所有非零元素，并将这些元素的线性索引值（linear indices：按列）返回
到向量ind中。如果X是一个行向量，则ind是一个行向量；否则，ind是一个列向量。如果
X不含非零元素或是一个空矩阵，则ind是一个空矩阵。
2. ind = find(X, k) 或 ind = find(X, k, 'first') 返回第一个非零元素k的索引值。
k必须是一个正数，但是它可以是任何数字数值类型。 参考技术B 建议你在命令窗口输入 doc find，出来的帮助文件介绍得很清楚。
[C,index]=find(a>=max(a)),对于这句指令，首先得清楚max(a),max(a)好像是求矩阵每一列的最大值，返回一个行向量，你给的a是单行矩阵，max(a)就返回其中的最大值。然后是find（），即寻找a中大于等于最大值的元素，C是满足条件的元素，index是元素在矩阵中的行列位置……不知道这样解答你清不清楚，总之很多指令的用法都可以在帮助文件里面找到的 参考技术C 从你的需要查找的数据来看，语句u=find（a>=max(a)）所得到的应该只是得到a向量中最大值所在原始向量中的序号。
对于行向量或列向量，u=find（a>=max(a)）得到的都只是一个值，a(u)都代表的是这个值，我刚试了下没出现楼主出现的问题。

Matlab 之 find()函数

　　当我第一次用matlab语言编写一个工程项目时，发现自己编写的脚本里循环特别多，导致编程效率很低，这让我特别苦恼。有一次导师让我阅读他编写的一个Matlab脚本，并按照新要求对其进行更改。我发现脚本里多次用到find()函数，猛然间豁然开朗，原来有大部分循环可以用find(）函数解决！对我而言，find()函数是我从Matlab小白开始进阶的标志，仅以此文，献给find()函数！

　　[1] find()函数基本功能

　　find()函数的基本功能是返回向量或者矩阵中不为0的元素的位置索引。

1 >> X = [1 0 4 -3 0 0 0 8 6];
2 >> ind = find(X)
3 
4 ind =
5 
6      1     3     4     8     9

　　其有多种用法，比如返回前k的不为0的元素：

1 >> ind = find(X,2)
2 
3 ind =
4 
5      1     3

　　也可以写成：

1 >> ind = find(X,2,‘first‘)
2 
3 ind =
4 
5      1     3

　　返回后k个不为0的元素：

1 >> ind = find(X,2,‘last‘)
2 
3 ind =
4 
5      8     9

　　若X是一个矩阵，索引该如何返回呢？

 1 >> X = [1 -3 0;0 0 8;4 0 6]
 2 
 3 X =
 4 
 5      1    -3     0
 6      0     0     8
 7      4     0     6
 8 
 9 >> ind = find(X)
10 
11 ind =
12 
13      1
14      3
15      4
16      8
17      9

　　这是因为在Matlab在存储矩阵时，是一列一列存储的，我们可以做一下验证:

1 >> X(4)
2 
3 ans =
4 
5     -3

　　假如你需要按照行列的信息输出该怎么办呢？

>> [r,c] = find(X)
r =
     1
     3
     1
     2
     3

c =
     1
     1
     2
     3
     3

　　如果你还需要输出具体的元素值：

 1 >> [r,c,v] = find(X)
 2 r =
 3      1
 4      3
 5      1
 6      2
 7      3
 8 
 9 c =
10      1
11      1
12      2
13      3
14      3
15 
16 v =
17      1
18      4
19     -3
20      8
21      6

　　[2] 进阶技巧

　　（1）find()函数的功能是找到向量或者矩阵中不为0的元素，那如果需要找到其中满足一定条件的元素，比如，等于4的元素该怎么办呢？

1 >> X = [1 0 4 -3 0 0 0 8 6];
2 >> ind = find(X == 4)
3 
4 ind =
5 
6      3

　　在Matlab中，有一个logical数据类型，和C++中的bool值相同。find()函数将logical值0也视为数值0，logical值1视为数值1，因此可以用上面的方法找到矩阵中满足一定条件的元素。

　　（2）在程序中，我们经常会以向量中是否包含某一元素为判断条件，比如X中是否有等于9的元素，这时也可以应用find()函数：

>> if isempty(find(X == 9))
        log = 0
    else
        log = 1
    end

log =
     0        

　　这里需要补充说明下，当输入中没有非0元素时，findf()返回一个空数组。

　　[3] 高级技巧

　　在进阶技巧（2）中，假如向量X的阶数很高，比如1000万维，如果我们用上面的方法，运行速度会很慢。仔细分析这个语句，我们其实并不希望找到X中等于9的全部元素，我们只想找到其中是否有等于9的元素，所以，我们只需作如下更改：

1 isempty(find(X == 9,1))

　　我们做一个小实验：

1 X = randi(100,1,10000000);
2 tic
3 for i = 1:100
4     isempty(find(X == 9));
5 end
6 toc

　　输出是：

1 Elapsed time is 1.954537 seconds.

　　更改后：

1 X = randi(100,1,10000000);
2 tic
3 for i = 1:100
4     isempty(find(X == 9,1));
5 end
6 toc

　　输出为：

1 Elapsed time is 0.757994 seconds.

　　运行效率约是原来的2.5倍。