首页 > 其他分享 >Matlab 之 find()函数

Matlab 之 find()函数

时间:2022-10-29 18:01:12浏览次数:77  
标签:函数 元素 matlab ind find Matlab

 当我第一次用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倍。

本文出自:Matlab 之 find()函数 - anzhiwu - 博客园 (cnblogs.com)

招募大量matlab技术人员,有大量matlab需求订单,均为个人短期可以完成,有时间的朋友可以加我微信:xiaoyuer-8988 加好友备注博客园matlab技术即可。有需求也可以联系此微信。

 

标签:函数,元素,matlab,ind,find,Matlab
From: https://www.cnblogs.com/sy01/p/16839277.html

相关文章

  • matlab 中fft的用法
    一.调用方法X=FFT(x);X=FFT(x,N);x=IFFT(X);x=IFFT(X,N)用MATLAB进行谱分析时注意:(1)函数FFT返回值的数据结构具有对称性。例:N=8;n=0:N-1;xn=[43267890];Xk=fft(xn)→Xk......
  • CMake Error at /opt/ros/melodic/share/catkin/cmake/catkinConfig.cmake:83 (find_p
    问题描述ros遇到如下错误。CMakeErrorat/opt/ros/melodic/share/catkin/cmake/catkinConfig.cmake:83(find_package):Couldnotfindapackageconfigurationfilepr......
  • 为什么有时候某函数明明只写了一次,gcc和ld却报错 “multiple definition” ?
    最主要的可能:你在头文件里写了非staticinline的函数定义(注意定义和声明的区别)有函数体的函数->定义,而非声明有初始化的全局变量->定义,而非声明原因:1.有可能是......
  • 入坑MATLAB必会的吐血总结
    本渣想回过头来整理一下MATLAB的一些基本的知识(很多东西比较琐碎,应该系统的梳理梳理),下文中没有提到的,自己用help查即可。此文用来存个档,便于回顾。由于matlab......
  • Matlab GUI界面
    MatlabGUI界面去年做了一些关于MatlabGUI的程序,现在又要做相关的东西,回想起来,当时很多经验没有记录下来,现在回顾起来始终觉得不爽,所以从现在开始,一定要勤写记录。  ......
  • MATLAB命令大全(转载)
    MATLAB命令大全管理命令和函数   help     在线帮助文件   doc     装入超文本说明   what     M、MAT、MEX文件的目录列表   type  ......
  • Python函数
    文章目录​​函数的概念​​​​函数的定义格式​​​​函数的调用格式​​​​函数的调用方式​​​​链式调用​​​​嵌套调用​​​​递归调用​​​​默认参数与变长参......
  • C语言函数指针在形参列表和返回值中的函数声明写法
    tags:C/C++Pointer写在前面最近看网络,发现系统调用中的信号函数的声明有点奇怪,如下:void(*signal(intsig,void(*func)(int)))(int);虽然书中给出了解释,但是奈何......
  • python Tips -----主函数
    if__name__=="__main__":python和C++/Java不一样,没有主函数一说,也就是说python语句执行不是从所谓的主函数main开始的。当运行单个python文件时,如运行a.py,这个时候a......
  • 通过主动调用脱函数抽取类壳
    被动调用以自实现的抽取壳为例,对Test01.proc1函数进行抽空,但这次不同的是此Test01类在apk加载的过程中并不会被加载。通过之前被动调用的方式dumpdex发现此函数还是被抽......