第2章_MATLAB数据及其运算.ppt

MATLAB语言及应用 电气学院 叶满园 myye@4 第二章 MATLAB数据及其运算 2.1矩阵及其赋值 2.1.1 矩阵 MATLAB中最基本的数据结构是矩阵(matrix)。 1*1的矩阵----标量(scalar): [5] 只有一行或一列的矩阵-----向量(vector): [1 3 5 7] MATLAB中一个有趣的例子 魔方阵 2.1.2 矩阵元素赋值 矩阵不需维数说明和类型定义，存储单元完全由计算机自动分配。每个变量代表一个矩阵，矩阵的每个元素都看作复数，默认的数据类型为双精度型。 （1）输入矩阵最简单的方法是输入矩阵的元素表， 每行的元素之间用空格或逗号隔开， 用“；”号作为元素表中每一行的结束符， 并用[ ]将元素表括起来。 例如，在命令窗口输入语句： a=[1 2 3；4 5 6；7 8 9] 则结果为：a = 1 2 3 4 5 6 7 8 9 ? 内存变量文件 Save 文件名[变量名表] [-append] [-ascii] load 文件名[变量名表] [-ascii] 将工作空间的变量am和D存储 Save mydata am D Load mydata 常用预定义变量及其含义 重点注意Inf在电容参数赋值时的应用 矩阵赋值的其它方式 （2）向量的生成 向量的生成归纳为：from:step:to结构。 如：　y=[0：2：10] 得：　y=[0,2,4,6,8,10] 如果step=1，则可省略： 　　如：　 x=1:5 则生成一个行向量， x=[1 2 3 4 5 ] （3）向量还可以采用linspace函数产生 　　linspace(a,b,n) 　　★注意上述两种方法的区别！ （4）从外部数据文件装载矩阵 （5）使用函数或程序生成矩阵:如ones(3,2), magic(3)等 2.1.3 矩阵元素的使用 x=a(1,3) +a(2,2) +a(3,1) 矩阵变量名(下标或序号) a（3，2）=8 （或a（6）=8 ）； 下标(Subscript): (行号,列号),行号从上到下，列号从左到右。 序号(Index)：矩阵元素的存储是按列进行的，顺序是从上到下，从左到右。 2.1.4 矩阵拆分 （1）取子矩阵 在一个下标处单独一个冒号，表示取出该下标处所对应的全部行或列。 a( : , k ) 表示取矩阵a的第k列全部元素。 a( m , : ) 表示取矩阵a的第m行全部元素。 a( m : m+k , :) 表示取a矩阵的第m~m+k行元素的全部元素 。 a( : ) 将矩阵a 的所有元素按列堆叠起来，成为一个列向量。 例如： b=a（：，3） %取出A矩阵的第三列。 　　　　　　 （注：%为注释符） c= a（1：2，：） %取出A矩阵的1到2行。 d= a（1：2，1：2）%取出A矩阵1到2行的1至2列。 （2）利用空矩阵删除矩阵的元素 a=[ ] a的维数为0。 例：a( 2 , : )= [ ]; 得: a = 8 1 6 4 9 2 2.1.5 复数（Complex Number）操作 复数由实部（real part）和虚部（imaginary part）组成。 虚数单位用i或j来表示。 6+5i = 6+5j 例： R = [1,2;4,5]; I = [6,7;8,9]; ri=R+i*I ri = 1.0000 + 6.0000i 2.0000 + 7.0000i 4.0000 + 8.0000i 5.0000 + 9.0000i 或：ri=[1+6i,2+7i;4+8i,5+9i] 2.2 表达式（Expressions） 表达式表达式由变量、常量、运算符和函数组成。 MATLAB对键入的表达式进行翻译和计算，语句的形式通常为： 变量=表达式； 或简单地写作：表达式； ；分号是可选的，如果该语句的最后一个字符是分号，则在执行词句时不在屏幕上显示结果。 rho = (sqrt(5)-1)/2 得：rho =0.6180 a = abs(3+4i) 得：a =5 若一个表达式太长，可用三个句号将其延伸到下一行： z = 10*sin(pi/3)* ... sin(pi/3); 2.2.1 变量（Variables） 变