有限元教程.doc

第1章 MATLAB基础知识 1.1 MATLAB数据的特点 1．矩阵是MATLAB最基本、最重要的数据对象。单个数据(标量)可以看成是矩阵的特例。 2．MATLAB数据类型 数值数据：双精度型、单精度数、带符号整数和无符号整数。 字符数据。 结构(Structure)和单元(Cell)。 多维矩阵和稀疏矩阵(Sparse)。 1.2 变量和赋值 1.2.1 变量的命名 在MATLAB中，变量名是以字母开头，后接字母、数字或下划线的字符序列，最多19个字符。 在MATLAB中，变量名区分字母的大小写。MATLAB提供的标准函数名以及命令名必须用小写字母。 1.2.2 赋值语句 MATLAB赋值语句有两种格式： (1) 变量=表达式 (2) 表达式 一般地，运算结果在命令窗口中显示出来。如果在语句的最后加分号，那么，MATLAB仅仅执行赋值操作，不再显示运算的结果。 在MATLAB语句后面可以加上注释，注释以%开头，后面是注释的内容。 1.2.3 数据的输出格式 MATLAB用十进制数表示一个常数，具体可采用日常记数法和科学记数法两种表示方法。 数据输出时用户可以用format命令设置或改变数据输出格式。format命令的格式为： format 格式符 注意，format命令只影响数据输出格式，而不影响数据的计算和存储。 1.2.4 预定义变量 在MATLAB工作空间中，还驻留几个由系统本身定义的变量。它们有特定的含义，在使用时，应尽量避免对这些变量重新赋值。 1.2.5 内存变量的管理 1. 内存变量的显示与删除 who和whos这两个命令用于显示在MATLAB工作空间中已经驻留的变量名清单。 clear命令用于删除MATLAB工作空间中的变量。 注意，预定义变量不能被删除。 2. 工作空间浏览器 (1) 工作空间浏览器的启动 (2) 工作空间浏览器的操作 3. 内存变量文件 利用MAT文件(.mat)可以把当前MATLAB工作空间中的一些有用变量长久地保留下来。 MAT文件的生成和装入由save和load命令来完成。常用格式为： save 文件名 [变量名表] [-append][-ascii] load 文件名 [变量名表] [-ascii] 1.3 MATLAB矩阵 1.3.1 矩阵的建立 1. 直接输入法 将矩阵的元素用方括号括起来，按矩阵行的顺序输入各元素，同一行的各元素之间用空格或逗号分隔，不同行的元素之间用分号分隔。例如 A=[1 2 3;4 5 6;7 8 9] 2. 利用M文件建立矩阵 对于比较大且比较复杂的矩阵，可以为它专门建立一个M文件。 3. 利用MATLAB函数建立矩阵 几个产生特殊矩阵的函数：zeros、ones、 eye、rand、randn。 这几个函数的调用格式相似，下面以产生零矩阵的zeros函数为例进行说明。其调用格式是： zeros(m) 产生m×m零矩阵 zeros(m,n) 产生m×n零矩阵。 zeros(size(A)) 产生与矩阵A同样大小的零矩阵 相关的函数有：length(A)给出行数和列数中的较大者，即length(A)=max(size(A))；ndims(A)给出A的维数。 此外，常用的函数还有reshape(A,m,n)，它在矩阵总元素保持不变的前提下，将矩阵A重新排成m×n的二维矩阵。 1.3.2 冒号表达式 冒号表达式的一般格式： e1:e2:e3 还可以用linspace函数产生行向量：linspace(a,b,n) linspace(a,b,n)与a:(b-a)/(n-1):b等价。 13.3 矩阵的拆分 1. 矩阵元素 MATLAB允许用户对一个矩阵的单个元素进行赋值和操作。例如 A(3,2)=200 也可以采用矩阵元素的序号来引用矩阵元素。矩阵元素按列编号，先第一列，再第二列，依次类推。 以m×n矩阵A为例，矩阵元素A(i,j)的序号为(j-1)*m+i。其相互转换关系也可利用sub2ind和ind2sub函数求得。 2. 矩阵拆分 (1)利用冒号表达式获得子矩阵 ①A(:,j)表示取A矩阵的第j列全部元素；A(i,:)表示A矩阵第i行的全部元素；A(i,j)表示取A矩阵第i行、第j列的元素。 ②A(i:i+m,:)表示取A矩阵第i～i+m行的全部元素；A(:,k:k+m)表示取A矩阵第k～k+m列的全部元素，A(i:i+m,k:k+m)表示取A矩阵第i～i+m行内，并在第k～k+m列中的所有元素。 此外，还可利用一