MATLAB基本语法详细说明

1、Matlab语法第一节 基本数值计算1. 变量：分为数值变量和字符变量2. 常量：计算机中不变的量。如i、j、pi、NaN(不确定)、Inf(无穷大)3. 字符变量：将字符串作为变量。有三种方法表示： (1) 用单引号' ' (2) 用函数sym(' ') (3) 用命令symbs4. 举例 x=2 % 将2赋给变量x y=3; % 有;表示在命令窗口不显示y的值 z=x2 -y % 数值计算。输出结果为1 f='sin(x)' % 用单引号定义一个字符变量 g=sym('cos(y)') % 用函数sym(' '

2、)定义一个字符变量 syms a b % 用命令syms定义字符变量。一般用于多符号变量的定义 u=2*a % 字符计算。输出结果为2*a w=b2-1 % 字符计算。输出结果为b2-1 fg=f+g % 字符计算。输出结果为sin(x)+cos(y) uw=u*w % 字符计算。输出结果为2*a*(b2-1) u/w % 字符计算。输出结果为2*a/(b2-1)第二节 矩阵构造及运算 Matlab中数据的结构形式就是一个矩阵。如x=2是一个1×1的矩阵 1. 矩阵的建立 (1) 直接输入法。 (2) 冒号法(1×N)。 (3) 函数法(特殊矩阵)。 (4) 矩阵的编辑(A

3、rray Editor)。 2. 向量 向量是1×N的特殊矩阵，即只有一行或者一列，称为N维向量。 3. 向量的点积与叉积 点积：dot(A,B) 叉积：cross(A,B) 4. 举例： x=1 2 3; 4 5 6; 7 8 9 % .，元素之间用空格，行之间用分号。 y=0: 0.1 :1 % 不用.，初值、步长、终值间用冒号。 w=eye(3) % 3阶单位矩阵。 u=rand(3) % 3阶随机矩阵。元素在0-1之间。 u1=rand(2, 3) % 2*3阶随机矩阵。元素在0-1之间。 q=randn(3) % 3阶随机矩阵。元素在0-1之间。元素符合正态分布。 q1=r

4、andn(2, 3) % 2*3阶随机矩阵。元素在0-1之间。元素符合正态分布。 s=magic(3) % 魔方阵。各行各列以及对角元素的和相等。 ss=zeros(3) % 3阶全零阵。 uu=ones(3) % 3阶全1阵。 5. 矩阵的加、减、乘、除 (1) 参与加、减运算的矩阵必须同维 (2) A/B时，A、B列数必须相同。AB时，A、B行数必须相同。 (3) 矩阵左乘与右乘不同，左除与右除不同。(A*B不等于B*A，AB不等于A/B)。 (4) 标量或函数与矩阵的运算等于该标量或函数与矩阵的每一个元素的运算。 6. 对矩阵的操作 6.1 对矩阵的元素操作： (1) 提取矩阵A的第n行

5、第m列的元素，表示为A(n, m) (2) 提取矩阵A的第n行的所有元素，表示为A(n, :) (3) 提取矩阵A的第m列的所有元素，表示为A(:, m) (4) 将矩阵A的第n行第m列的元素重新赋值b，表示为A(n, m)=b (5) 将矩阵A的第n行的所有元素重新赋值b，表示为A(n, :)=b (6) 将矩阵A的第m列的所有元素重新赋值b，表示为A(:, m)=b (7) 将矩阵A的第n行第m列的元素删除，表示为A(n, m)= (8) 将矩阵A的第n行的所有元素删除，表示为A(n, :)= (9) 将矩阵A的第m列的所有元素删除，表示为A(:, m)= 6.2 矩阵的部分操作： Fli

6、plr(A) % 矩阵左右翻转 Flipud(A) % 矩阵上下翻转 Flipdim(A, m) % 矩阵沿特定维(m)翻转 Rot90(A, k) % 矩阵逆时针旋转k*90度 Tiag(A, k) % 取矩阵对角线元素 Tril(A, k) % 取矩阵的下三角部分 Triu(A, k) % 取矩阵的上三角部分 注：k可以是正、负整数，但绝对值一般不大于矩阵的维数。 6.3 矩阵分析操作： (1) 方阵A的行列式值，可表达为 det(A) (2) 矩阵A的秩，可表达为 rank(A) (3) 行列式值不为堆的方阵A，求逆矩阵可表达为 inv(A) (4) 矩阵A的转置阵，可表达为 A'

7、; (5) 矩阵A的特征向量与特征值：V, D=eig(A)。其中V和D分别为A的特征向量和特征值 6.4 矩阵的数组运算： 矩阵乘：A*B。A的列数与B的行数要相等。 矩阵的数组乘：A.*B。表示为A、B矩阵的对应元素一一相乘，即Aij * Bij。A与B的维度要相同。 矩阵的数组除：A./B或者A.B 举例：x=-10 : 2 :10 y=sin(x) % correct y=sin(x2) % incorrect - Matrix x must be square y=sin(x.2) % correct 6.5 矩阵元素的关系运算与逻辑运算(与、或、非) < 小于 > 大于

8、 <= 小于或等于 >= 大于或等于 = 等于 = 不等于 6.6 矩阵的多维数组形式 (1) 函数cat的使用 (2) B=cat(dim, a1, a2,.) (3) 意义：将多个同维数组a1, a2,.构成一个高维数组B。dim是高维数组B的维数。 它必须等于或大于a1, a2,.的阶次 举例：a=1 2 3 ; 2 3 4 ; 3 5 6 b=a+10 % 新构成一个矩阵 a(:, :, 2)=b % 将新矩阵赋给a的第二层 a(:, :, 1) % 显示a的第一层 c=cat(3, a, b) % 用函数来构成高维数组 6.7 多项式及其运算式 多项式的创建： (1) 由

9、1×N的N维向量 P=a0 a1 a2 an 表示 Pn=a0xn + a1xn-1 + a2xn-2+ . +an多项式。 (2) 由函数poly(A)定义 如果A为二维或以上的矩阵，poly(A)表示由A的特征根确定的多项式。如果A为一维矩阵 poly(A)表示由A的元素为多项式的根确定的多项式。 举例： %方法一 p=1 -2 3 %直接给出多项式p poly2sym(p) %给出p多项式的表达式 %方法二 a=1 2; -2 4 ps=poly(a) % 计算a的特征根确定的多项式 poly2sym(ps) % 给出ps多项式的表达式 %方法三 x=-1 2 px=poly(

10、x) % 以x的元素为多项式的根确定的多项式 poly2sym(px) % 给出ps多项式的表达式1.变量赋值标量：a=10向量：a=1,3,5,6;a=1:2:9;a=1 7 6 2矩阵：a=1 2 3;4 5 6;7 8 92.程序控制语句（1）判断语句if  表达式1                      a=4;b=6;   命令 &#

11、160;                 if a>b  elseif  表达式2             'a大于b'   命令          

12、0;          elseif a<b。                        'a小于b' else 命令         

13、0;         elseif  a=bend                          a等于b            

14、60;               else                               error('不可能的事情发生了')&

15、#160;                           end（2）分支语句switch和case  switch(a)      case 0      case 1   

16、0;  case 2      .      otherwise  enda=5;b=mod(5,3);switch(b)   case 0      模3为0   case 1      模3为1   case 2      模3为2&#

17、160;  otherwise      不可能end(3)循环语句for n=3:32    r(n)=n;end while 条件    命令或函数end当条件不满足时，循环语句执行中间的命令或函数。当条件满足时，跳出循环执行下面的命令。n=1;while n<30r(n)=n;n=n+1;end(4) 跳出循环语句有时可能需要在循环体中的某个点跳出，这里通过continue和break实现continue语句循环后，直接从该点跳到循环体的开始；break

18、语句执行后，直接从该点跳出循环体.n=0;while n<30n=n+1;if n=15 continue;endr(n)=n;endn n=30 n=0;while n<30n=n+1;if n=15 break;endr(n)=n;endn n=151.循环语句forfor  i=s1:s3：s2    循环语句组end解释：首先给i赋值s1；然后，判断i是否介于s1与s2之间；如果是，则执行循环语句组,i=i+s3（否则，退出循环.）；执行完毕后，继续下一次循环。例：

19、求1到100的和，可以编程如下：   sum=0   for i=1:1:100     sum=sum+i   end   这个程序也可以用while语句编程。      注：for循环可以通过break语句结束整个for循环.2.循环语句while  例：sum=0;i=1;      while(i<=100)       sum=sum+i;i=i+1;      end      3.if语句  if(条件)