1、1 Text 函数的用法：用法 text(x,y,string)在图形中指定的位置(x,y)上显示字符串stringtext(x,y,z,string) 在三维图形空间中的指定位置(x,y,z)上显示字符串 string2，plot(0,z1,z12,-b,LineWidth,3) 里面表示数组.3,x,y 均为矩阵，plot 命令就是画出 x，y 矩阵对应的二维平面的点形成的曲线。y（：，1）中逗号前是行，逗号后是列，冒号表示从几到几。所以 y（：，1）表示第一列的所有元素。如果是 y(3:5,1)则表示第一列的第 3 到第 5 行对应的元素。只要你的 y 矩阵有 100 列，那你当然可以将
2、 1 改成 100。同理，x 矩阵也可以这样。4 sym 的意思是 symbol，就是后面括号里面是个代数式，要进行符号运算，class（）判断对象是什么类型。5 matlab 控制运算精度用的是 digits 和 vpa 这两个函数xs = vpa(x,n)在 n 位相对精度下，给出 x 的数值型符号结果 xsxs = vpa(x)在 digits 指定的精度下，给出 x 的数值型符号结果 xsdigits 用于规定运算精度，比如：digits(20);这个语句就规定了运算精度是 20 位有效数字。但并不是规定了就可以使用，因为实际编程中，我们可能有些运算需要控制精度，而有些不需要控制。vp
3、a 就用于解决这个问题，凡是用需要控制精度的，我们都对运算表达式使用 vpa 函数。例如：digits(5);a=vpa(sqrt(2);这样 a 的值就是 1.4142，而不是准确的1.4142135623730950488016887242097又如：digits(5);a=vpa(sqrt(2);b=sqrt(2);这样 a 的值是 1.4142，b 没有用 vpa 函数，所以 b 是1.41421356237309504880168872420976 识别变量的指令：findsym(EXPR) 确认表达式 EXPR 中所有自由符号变量findsym(EXPR,N) 从表达式 EXPR
4、中确认出靠 x 最近的 N 个自由符号变量默认自由变量的顺序为：x y w z v.7 matlab 中 size 函数的用法 size：获取数组的行数和列数（1）s=size(A),当只有一个输出参数时，返回一个行向量，该行向量的第一个元素时数组的行数，第二个元素是数组的列数。（2）r,c=size(A),当有两个输出参数时，size 函数将数组的行数返回到第一个输出变量，将数组的列数返回到第二个输出变量。（3）如果在 size 函数的输入参数中再添加一项，并用 1 或 2 为该项赋值，则 size 将返回数组的行数或列数。其中 r=size(A,1)该语句返回的时数组 A 的行数， c=s
5、ize(A,2) 该语句返回的时数组 A 的列数。 8 matlab eig 函数 在 MATLAB 中，计算矩阵 A 的特征值和特征向量的函数是 eig(A)，常用的调用格式有5 种：(1) E=eig(A)：求矩阵 A 的全部特征值，构成向量 E。(2) V,D=eig(A)：求矩阵 A 的全部特征值，构成对角阵 D，并求A 的特征向量构成V 的列向量。(3) V,D=eig(A,nobalance)：与第 2 种格式类似，但第 2 种格式中先对 A 作相似变换后求矩阵 A 的特征值和特征向量，而格式 3 直接求矩阵 A 的特征值和特征向量。(4) E=eig(A,B)：由 eig(A,B
6、)返回 NN 阶方阵 A 和 B 的 N 个广义特征值，构成向量 E(5) V,D=eig(A,B)：由 eig(A,B)返回方阵 A 和 B 的 N 个广义特征值，构成 NN 阶对角阵 D，其对角线上的 N 个元素即为相应的广义特征值，同时将返回相应的特征向量构成 NN 阶满秩矩阵，且满足 AV=BVD。eigFind eigenvalues and eigenvectorsSyntaxd = eig(A)d = eig(A,B)V,D = eig(A)V,D = eig(A,nobalance)V,D = eig(A,B)V,D = eig(A,B,flag)d = eig(A)和 V,D
7、 = eig(A)最为常用 注意，第一列为对应第一个特征值的特征向量。9 函数 subexpr()用法函数 subexpr()将符号表达式中重复出现的字符串用符号变量代替，其具体使用方法如下： Y,SIGMA = subexpr(S,SIGMA)：指定用符号变量 SIGMA 来代替符号表达式中重复出现的字符串；10 函数 subs()用法函数 subs()用指定符号替换符号表达式中的某一特定符号，其具体使用方法如下 : R = subs(S,Old,New)：用新符号变量 New 替代原来符号表达式 S 中的变量 Old。11 limit 函数的用法limit(A,B,C) A:函数 B：待赋
8、值的变量 C：要给变量赋的值12 diff 函数的用法diff(f,v,n) 求 f 对变量 v 的 n 阶导数diff(f,n) 求 f 对默认变量的 n 阶导数diff(f,v) 求 f 对变量 v 的一阶导数diff(f) 求 f 对默认变量的一阶导数13 jacobian 函数的用法这是在命令窗口出来的帮助，其大概意思是：函数 jacobian 就是雅可比矩阵，命令 jacobian（f，v）是计算标量函数或矢量函数 f 关于矢量 v 的雅可比矩阵。其结果就相当于 df/dv。需要注意的是，当 f 是标量时，就相当于是求 f 的梯度。14 taylor 函数的用法r=taylor(f,
9、n,v,a) 把 在 处展开为幂级数)(vfaknkaxf)(!10)(15 legend 函数的基本用法是LEGEND(string1,string2,string3, .)分别将字符串 1、字符串 2、字符串 3标注到图中，每个字符串对应的图标为画图时的图标。17 求和 问题，具体如下bavf)(s=symsum(f,v,a,b) 求通式 f 在指定变量 v 取遍 a,b中所有整数时的和。18 符号表达式的基本操作符号运算中有许多操作指令，如 collect（合并同类项） 、expand（对指定项展开） 、factor（进行因式或因子分解） 、horner（转换成嵌套形式） 、numden
10、（提取公因式） 、simplify（恒等式简化） 、pretty（习惯方式显示）等，其中最常用的是simple(EXPR) 运用包括 simplify 在内的各种指令把 EXPR转换成最简短形式19 ezplot 功 能 简 介 ： 在 matlab 的 命 令 窗 口 中 键 入 help ezplot 命 令 或 者 doc ezplot 即 可 获 得 本 函 数 的 帮 助 信 息 。 EZPLOT 即 ： Easy to use function plotter。 它 是 一 个 易 用 的 一 元 函 数 绘 图 函 数 1 。 特别 是 在 绘 制 含 有 符 号 变 量 的 函
11、 数 的 图 像 时 ， ezplot 要 比 plot 更方 便 。 因 为 plot 绘 制 图 形 时 要 指 定 自 变 量 的 范 围 ， 而 ezplot 无需 数 据 准 备 2， 直 接 绘 出 图 形 。 ezplot 的 调 用 格 式 ： 1、 ezplot(fun) 2、 ezplot(fun,min,max) 3、 ezplot(fun2) 4、 ezplot(fun2,xmin,xmax,ymin,ymax) 5、 ezplot(fun2,min,max) 6、 ezplot(funx,funy) 7、 ezplot(funx,funy,tmin,tmax) 8、
12、ezplot(.,figure_handle) 9、 ezplot(axes_handle,.) 10、 h = ezplot(.)20 dsolve 函数的用法S=dsolve(eq1, eq2, , eqn, cond1, cond2, , condn, v)eq1,eq2为微分方程或微分方程组， cond1为边界条件或初始条件， v为变量，默认的变量为 t. 若应变量为 y ，用“Dny”表示“y 的 n 阶导数” ， Dy 为一阶导数。解在 S.y 中21 如何在 matlab 中调用 maple （转自 onlive 快乐生活）（不用安装 maple 软件就可调用）方法 1： map
13、le(maplestatement)其中 maplestatement 是完整的 maple 语句，由一条或几条命令组成，必须符合 maple 的语法方法 2：maple(function，arg1, arg2,)其中 function 为 maple 中的函数名称，arg1, arg2,是函数function 所用的参数。注：如果方法 1 行不通，可尝试方法 2（个人经验）。基本代数部分如何用 matlab 求阶乘factorial（n） 求 n 的阶乘如何用 matlab 配方 没有发现 matlab 有这一命令，不过我们可以调用 maple 的命令，调用方法如下：首先加载 maple 中
14、的 student 函数库，加载方法为：maple（with（student）然后运行 maple 中的配方命令，格式为：maple（completesquare（f） 把 f 配方，其中 f 为代数表达式或代数方程maple(completesquare（f，x）) 把 f 按指定的变量 x配方，其中 f 同上 maple(completesquare（f，x，y，）) 把 f 按指定的变量 x，y，配方maple(completesquare（f，x，y，）) 把 f 按指定的变量 x，y，配方，如何用 matlab 进行多项式运算（1） 合并同类项 syms 表达式中包含的变量 coll
15、ect(表达式，指定的变量)（2）因式分解 syms 表达式中包含的变量 factor(表达式)（3）展开syms 表达式中包含的变量 expand(表达式)我们也可在 matlab 中调用 maple 的命令进行多项式的运算，调用格式如下：maple（maple 中多项式的运算命令）如何用 matlab 进行分式运算 发现 matlab 只有一条处理分式问题的命令，其使用格式如下：n，d=numden（f） 把符号表达式 f 化简为有理形式，其中分子和分母的系数为整数且分子分母不含公约项，返回结果 n 为分子，d 为分母。注意：f 必须为符号表达式 不过我们可以调用 maple 的命令，调用
16、方法如下：maple（denom（f） 提取分式 f 的分母maple(numer（f）)提取分式 f 的分子maple(normal（f） )把分式 f 的分子与分母约分成最简形式maple(expand（f）)把分式 f 的分子展开,分母不变且被看成单项。maple(factor（f）)把分式 f 的分母和分子因式分解，并进行约分。如何用 Matlab 进行因式分解syms 表达式中包含的变量 factor(表达式)如何用 Matlab 展开syms 表达式中包含的变量 expand(表达式) 如何用 Matlab 进行化简syms 表达式中包含的变量 simplify(表达式) 如何用
17、Matlab 合并同类项syms 表达式中包含的变量 collect(表达式，指定的变量)如何用 Matlab 进行数学式的转换 调用 Maple 中数学式的转换命令，调用格式如下：maple(Maple 的数学式转换命令) 即：maple(convert(表达式，form) 将表达式转换成 form的表示方式 maple(convert(表达式，form, x) 指定变量为 x，将依赖于变量 x 的函数转换成 form 的表示方式（此指令仅对 form 为 exp 与sincos 的转换式有用）如何用 Matlab 进行变量替换syms 表达式和代换式中包含的所有变量 subs(表达式，要替
18、换的变量或式子，代换式) 如何用 matlab 进行复数运算 a+b*i 或 a +b*j 表示复数 a+bi 或 a+bjreal（z） 求复数 z 的实部imag（z） 求复数 z 的虚部abs（z） 求复数 z 的模angle（z） 求复数 z 的辐角，conj（z） 求复数 z 的共轭复数exp（z） 复数的指数函数，表示 ez如何在 matlab 中表示集合a, b, c, 表示由 a, b, c,组成的集合 （注意：元素之间也可用空格隔开）unique(A) 表示集合 A 的最小等效集合（每个元素只出现一次） 也可调用 maple 的命令，格式如下：maple(a, b, c,)
19、表示由 a, b, c,组成的集合 下列命令可以生成特殊的集合：maple(seq(f(i),i=nm) 生成集合f(n), f(n+1), f(n+2), , f(m)如何用 Matlab 求集合的交集、并集、差集和补集 union(A,B) 求集合 A 和 B 的并集intersect(A,B) 求集合 A 和 B 的交集setdiff(A,B) 求集合 A 和 B 的差集 A-Bsetdiff(U,A) 求集合 A 关于全集 U 的补集我们也可以调用 Maple 的相应功能，调用方法如下：maple(用 Maple 求集合的交集、并集、差集和补集的命令 )具体地说，共有下列几个调用命令：
20、maple( A union B ) 求集合 A 和 B 的并集maple( A intersect B ) 求集合 A 和 B 的交集maple(A minus B ) 求差集 A-Bmaple( 全集 I minus A ) 求集合 A 关于全集 I 的补集如何用 matlab 排序sort（v） 将向量 v 的元素从小到大排列（升序排列）22Gca 表示当前图形。Xtick 表示坐标轴刻度。函数 get 返回某些对象属性的当前值。使用函数 get 的最简单语法是 get(handle, PropertyName )。例如：c=get(Hl_a, color )返回具有句柄 Hl_a 对象
21、的颜色。函数 set 改变句柄图形对象属性，使用语法 set(handle, PropertyName ,value)。例如： set(Hf_1, Position ,p_vect)将具有句柄 Hf_1 的图形位置设为向量 p_vect 所指定的值。同样 set(Hl_a, color , r )将具有句柄 Hl_a 的对象的颜色设置成红色。23 plot(1,5,0,0,.r,MarkerSize,20)其中当 x=1 时 y=0,同理当 x=5 时 y=0.24 MATLAB 函 数 heaviside 简 介在 matlab 的 命 令 窗 口 中 键 入 doc heaviside 或
22、 者 help heaviside 可 以 获 得 如 下 帮 助 信 息 ： “heaviside(x) has the value 0 for x 0, and 0.5 for x = 0. heaviside is not a function in the strict sense.” 即 ： 若 令 y=heaviside(x) 则 当 x0 时 ， y 的 值 为 1； 当 x 等 于0 时 ， y=0.5。 这 是 一 个 单 位 阶 跃 函 数 。 从 严 格 意 义 上 来 说 ，heaviside 不 是 一 个 函 数 。 单 位 阶 跃 函 数 的 加 窗 特 性 应
23、用 很 广 泛 。例 一 ： 在 matlab 的 命 令 窗 口 输 入 ： heaviside(0) Enter ans = 0.5000 这 是 Matlab 2011b 中 的 结 果 在 原 来 原 点 处 值 不 存 在 的 情 况 有所 修 改 。 把 下 面 的 代 码 写 到 一 个 m 文 件 中 ， 运 行 ： Y= ; for x=-5:5 y=heaviside(x); Y=Y y; end Y 将 在 matlab 的 命 令 窗 口 中 输 出 ： Y = 0 0 0 0 0 NaN 1 1 1 1 1 例 二 ： 下 面 将 举 一 个 利 用 单 位 阶 跃
24、函 数 加 窗 特 性 的 例 子 ： 阶 跃 函 数 的 加 窗 特 性syms t y; y=cos(t)*(heaviside(t+0.5*pi)-heaviside(t-0.5*pi); ezplot(y); 这 样 就 可 以 获 得 cos(t)在 -pi/2不 用 输 入 ) t=0:0.001:1; y1=sin(10*t); y2=sin(15*t); subplot(211) plot(t,y1) subplot(212) plot(t,y2) 运 行 结 果 见 右 图 。 也 可 以 将 上 述 命 令 写 成 一 个 程 序 段 ： function subplot_
25、sample1() close all figure grid on t=0:0.001:1; y1=sin(10*t); y2=sin(15*t); subplot(211) plot(t,y1) subplot(212) plot(t,y2) 26 最常用的矩阵分析指令如下det(A) 行列式 Adiag(A) 取对角元构成向量，或据向量构成对角阵V, D=eig(A) 特征值分解，使 VDexpm(A) 矩阵指数 Aeinv(A) 矩阵逆 1poly(A) 矩阵的特征多项式rank(A) 矩阵秩27 1函数 collect()函数 collect()将符号表达式中同类项合并，其具体使用方
26、法如下： R=collect(S)：将表达式 S 中的相同次幂的项合并； R=collect(S,v)：将表达式 S 中变量 v 的相同次幂的项合并。2函数 expand()函数 expand()将符号表达式进行展开，其具体使用方法如下： R = expand(S)：将表达式 S 中的各项进行展开。3函数 horner()函数 horner()将符号表达式转换成嵌套形式，其具体使用方法如下： R = horner(S)：将符号多项式矩阵 S 中的每个多项式转换成它们的嵌套形式。4函数 factor()函数 factor()对符号多项式进行因式分解，其具体使用方法如下： R=factor(X)：
27、如果 X 是一个多项式或多项式矩阵，该函数将 X 表示成低阶多项式相乘的形式；如果 X 不能分解成有理多项式乘积的形式，则返回 X 本身。5函数 simplify()函数 simplify()将符号表达式按一定规则简化，其具体使用方法如下： R= simplify(S)：该函数可应用于包含和式、方根、分数的乘方、等符号表达式矩阵 S。6函数 simple()该函数是将符号表达式表示成最简形式，其具体使用方法如下： r = simple(S)：用几种不同的算术简化规则对符号表达式进行简化，并显示中间过程； r,how = simple(S)：不显示中间过程，并附加返回最简形式对应的简化方法 。7
28、. pretty习惯方式显示8. numden获取最小公分母和相应分子28 shading interp 在网格片内采用颜色插值处理，得出表面图显得最光滑！29 material 函数的用法该函数通过设置曲面的材料属性，控制 surface 和 patch 对象的反射系数。函数语法material shiny material dull material metal material(ka kd ks) material(ka kd ks n) material(ka kd ks n sc) material default 参数解析 material shiny 设置反射系数为光泽色调，即对
29、象相对于漫反射和环境光照，具有较高的镜面反射效果。material dull 设置反射系数为阴暗色调，即对象具有更多的漫反射光，没有镜面高亮反射，并且反射光的颜色仅取决于光源颜色。material metal 设置反射系数为金属高亮色调，即对象具有很高的镜面反射效果，很低的环境光线和漫反射。material(ka kd ks) 在一个向量中分别指定对象的环境、漫反射、镜面反射光线强度。material(ka kd ks n) 在一个向量中分别指定对象的环境、漫反射、镜面反射光线强度，以及高光指数。material(ka kd ks n sc) 在一个向量中分别指定对象的环境、漫反射、镜面反射光
30、线强度，以及高光指数和高光颜色。material default 对当前坐标系的 surface 和 patch 对象使用默认材质属性30 xlim,ylim 的用法就是设置当前图的横轴和纵轴的范围，lim 就是 limite 的意思。也可以不用 set 设置，而直接用函数 xlim 和 ylim 设置xlim(gca,-20,20);ylim(gca,-20,20);也可以用函数 axis 设置axis(gca,-20,20,-20,20)31 matlab 的 colormap 函数详解matlab 中，每个 figure 都有（而且仅有）一个 colormap，翻译过来就是色图。COLO
31、RMAP(MAP) 用 MAP 矩阵映射当前图形的色图。COLORMAP(default) 默认的设置是 JET.MAP = COLORMAP 获得当前色图矩阵.COLORMAP(AX,.) 应用色图到 AX 坐标对应的图形，而非当前图形。MAP 实际上是一个 mx3 的矩阵，每一行的 3 个值都为 0-1 之间数，分别代表颜色组成的 rgb 值，1 0 0 代表红色， 0 1 0代表绿色，0 0 1代表蓝色。系统自带了一些 colormap，如：winter、autumn 等。输入 winter，就可以看到它是一个 64x3 的矩阵。用户可以自定义自己的 colormap，而且不一定是 64
32、 维的。0 0 0 is black, 1 1 1 is white,1 0 0 is pure red, .5 .5 .5 is gray, and127/255 1 212/255 is aquamarine.那么颜色在 fill 或 patch，SURFACE 等函数中到底是如何显示的呢？本质上，是把具体的颜色变成 colormap 中的相应 index，也就是行数。这个过程叫做换算映射：将指定的数值颜色向量（矩阵）C，映射到对应的颜色。颜色矩阵 C 的数值范围为Cmin, Cmax, Cmin 和 Cmax 的数值或者为min(min(C) max(max(C)，也可以在 CAXIS
33、中设置。在 matlab 中，图形窗的属性CdataMapping缺省设置值为scaled，也就是线性的映射。Cmin 对应的值映射到 colormap 的第一行，Cmax 对应的值映射到colormap 的最后一行。映射过程如下：首先，需要根据 caxis 取得 Cmin 和 Cmax 两个变量（默认值为 0 和1） ，画图时如果指定了数值颜色向量（矩阵）C，Cmin 和 Cmax 自动设置为 C 中的最大值和最小值。当你想控制时，可以自定义。比如将 Cmax 减小，这样将把所有大于 Cmax 的 C 值，全部都映射到同一个颜色（colormap 中 index 最大的行代表的颜色） 。根据
34、 Cij 在 Cmin 和 Cmax 之间的比例关系，确定对应的颜色的index，默认为线性映射。也就是说，当制定了数值颜色向量（矩阵）C 之后画图，图中颜色的使用范围会自动占满整个颜色范围！实例：clc;clear all;x=0 1 1 0;y=0 0 1 1; %定义四个点 0 0 1 0 1 1 0 1H_F=fill(x,y,0 0.1 0.2 0.6); %定义四个点的 C 值row_cmap = 15; %定义色图矩阵的行数color_map1=zeros(row_cmap,3); %定义色图矩阵color_r = 0:1/(row_cmap-1):1; color_g = 0:
35、1/(row_cmap-1):1;color_b = 0:1/(row_cmap-1):1;color_map1(:,1) = color_r; color_map1(:,2) = color_g;colormap(color_map1);colorbar;%本例中颜色从0 0 0 变化到1 1 0%增加 row_cmap 的值，如变化到 100，则可看到颜色的渐变，而非跳跃型变化。32 S=solve(eq1, eq2, , eqn, v1, v2, , vn)求方程组关于指定变量的解33 数值积分trapz(x,y)%梯形法沿列方向求函数 Y 关于自变量 X 的积分cumtrapz(x,y
36、)%梯形法沿列方向求函数 Y 关于自变量 X 的累计积分quad(fun,a,b,tol)%采用递推自适应 Simpson 法计算积分quad1(fun,a,b,tol)%采用递推自适应 Lobatto 法求数值积分dbquad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol)%二重(闭型)数值积分指令 triplequad(fun,xmin,xmax,ymin,ymax,zmin,zmax,tol)%三重(闭型)数值积分指令34 linspace 函数的用法linspace 是用来生成一个等差数列的 lin 是 linear 的缩写x=linspace(a,b,n
37、)35 MATLAB 函 数 logspace 简 介用 法 ： logspace(a,b,n)， 其 中 a、 b、 n 分 别 表 示 开 始 值 、结 束 值 、 元 素 个 数 。 功 能 ： 生 成 从 10 的 a 次 方 到 10 的 b 次 方 之 间 按 对 数 等 分 的n 个 元 素 的 行 向 量 。 n 如 果 省 略 ， 则 默 认 值 为 50。 在 matlab 的命 令 窗 口 中 输 入 help logspace 或 者 doc logspace 可 以 获 得 更多 帮 助 信 息 。 相 关 函 数 ： linspace 编 辑 本 段 应 用 举 例
38、（ 1） x=logspace(a,b)生 成 有 50 个 元 素 的 对 数 等 分 行 向 量x， 且 x(1)=10 的 a 次 方 ， x(50)=10 的 b 次 方 ； （ 2） x=logspace(a,b,n)生 成 有 n 个 元 素 的 对 数 等 分 行 向 量x， 且 x(1)=10 的 a 次 方 ， x(n)=10 的 b 次 方 ； （ 3） x=logspace(a,pi)生 成 有 50 个 元 素 的 对 数 等 分 行 向 量x， 且 x(1)=10 的 a 次 方 ， x(50)= ；36 isnan 函数来判断变量是否是数值37 ind2sub 函数
39、使用方法 ind2sub 函数可以用来把矩阵元素的 index 转换成对应的下标（determines the equivalent subscript values corresponding to a single index into an array）例如： 一个 4*5 的矩阵 A，第 2 行第 2 个元素的 index 的 6(matlab中 matrix 是按列顺序排列)，可以用 ind2sub 函数来计算这个元素的下标 I,J = ind2sub(size(A),6)38 MATLAB 中，a(2:3,1:2:5)中 1:2:5 是什么意思？2：3 指的是第二行到第三行，步长默
40、认为 1，可以表示为2：1：3。1：2：5 这样看，第一列到第 5 列，但不是以 1 增加，而是已 2 增加，步长是 2。就是表示 1、3、5 共 3 列，比如 1：2：9 就是取1，3，5，7，9。39 B=A(1,2,2,2,1,3,5 )的解释把 A 变成四行，第一行去原 A 第一行的 1，3，5 列，第二行取原 A第二行的 1,3,5 列，第三行取原 A 第二行的 1,3,5 列，第四行取原A 第二行的 1,3,5 列。40 num2str(1 2 3;4 5 6)就是将数值转换为字符串。num2str(1 2 3;4 5 6,1,2)就是将数值转换为字符串的同时将第一列和第二列的字符串合并成一个字符串。41 find 函数的用法find 函数用于返回所需要元素的所在位置 (位置的判定：在矩阵中，第一列开始，自上而下，依次为 1，2，3.,然后再从第二列，第三列依次往后数)find（A）返回矩阵 A 中非零元素所在位置 A = 1 0 4 -3 0 0 0 8 6; X = find(A)X =1 3 4 8 9find（A5）返回矩阵 A 中大于 5 的元素所在位置 find(A5)ans =8 9
