计算方法第2章.pptVIP

第2章 MATLAB语言基础 几个基础概念 数值矩阵 字符串和符号矩阵 基本绘图方法 MATLAB语言编程 几个基础概念 标识符 数据类型 变量名及赋值 几个基础概念 标识符 用来做变量名、常量名、文件名 以英文字母开头 由英文字母、数字、下划线混合而成 大小写区分 几个基础概念 数据类型 数值型（Double Array）：含实数和复数 a=3.2 字符串型（Char Array）：以英文单引号界定 b=ab 符号型（Sym Object）：用sym或syms可以把字符、表达式、方程、矩阵等定义成数学符号，运算结果为数学表达式 c=sym(x) 用class（a）可知a的类型 几个基础概念 变量名及赋值 由赋值符号‘=’来完成赋值 由所赋的值决定变量的类型 三类数据结构的基本单位都是矩阵，单个值和向量是矩阵的特例。 2.1 数值矩阵 永久性数值变量名 数值矩阵的创建 数值矩阵元素的标识和修改 数值矩阵的矩阵算法 数值矩阵的数组算法 永久性数值变量名 pi：圆周率 eps：机器浮点运算误差限 Inf：正无穷大 NaN：不定值 i或j：虚数单位 ans：临时变量，代表未定义名称的变量 数值矩阵的创建 直接输入法 元素置于[]中 一行之间用逗号或空格隔开，每行用分号隔开 续行符为三个连续的英文句点… 元素类型不限 每行元素个数相等 输入完成后加分号，则建立的矩阵不显示，直接使用回车，则立即显示建立的矩阵 数值矩阵的创建 用指令创建特殊矩阵 zeros(n); zeros(m,n); ones(n); ones(m,n) eye(n); eye(m,n); diag(a,k); rand(n); rand(m,n); randn(n); randn(m,n); magic(n); tril(a); 见案例ex2_1_2_2.m 数值矩阵的创建 变换矩阵指令 flipud(a);a上下翻转 fliplr(a);a左右翻转 rot90(a,k);a沿逆时针方向旋转k个90度 rot90(a);k默认为1 reshape(a,m,n);a的元素按照m*n重排 数值矩阵的创建 特殊向量的创建 等差数列 a为初值,b为终值,h为步长默认为1 t=a:h:b; t=(a:h:b); t=[a:h:b]; a为初值,b为终值,n为等分节点数,默认值为100 t=linspace(a,b,n); 数值矩阵的创建 特殊向量的创建 等比数列 a为初值,b为终值,n为等比节点数,默认值为50 t=logspace(log10(a),log10(b),n); 全为1的向量 t=linspace(1,1,n); t=ones(1,n); t=logspace(log10(1),log10(1),n); 练习 创建矩阵 不采用直接赋值法 数值矩阵元素的标识用() a(p); a中序号为p的元素,序号按先列后行计算 a(:); a中所有元素,序号按先列后行排列成的列阵 a(m,n); a中m行n列对应的元素 a(:,n); a中n列对应的元素列阵 a(m,:); a中m行对应的元素行阵 a([p,q,r],n); a中第p,q,r行的第n列对应的元素列阵 a(m,[p,q,r]); a中第p,q,r列的第m行对应的元素行阵 a(p:q,n); a中从p到q行的第n列对应的元素列阵 a(m,p:q); a中从p到q列的第m行对应的元素行阵 a([p,q,r],[w,s]); a中第p,q,r行与第w,s列的交点元素构成的矩阵 a(:,:,n); 三维矩阵a的第n页的平面距阵 数值矩阵的增删修改(标识+赋值) a=magic(3)%创建3阶魔方阵 a(4,:)=[1,2,3]%加一行 a(2,:)=[3,4,6]%修改一行 a(:,1)=[]%删除一列 a(:,3)=[3;5;8;9]%增加一列 练习 已知A=(1:8),用A构造组合出以下矩阵B 数值矩阵的矩阵算法 矩阵维数和转置 format short a=rand(3,4) size(a)%输出a的行数和列数 size(a,1)%输出a的行数 size(a,2)%输出a的列数 b=a‘%求a的共轭转置距阵,也即转置矩阵 c=[a;1,2i,3,4i]%求复数矩阵c d=c e=conj(d)%求d的共轭距阵，也即c的转置矩阵 数值矩阵的矩阵算法 矩阵的加,减,乘法 加减法的操作数维数必须相等 a+b a-b 乘法的操作数的内维数必须相等 a*b mtimes(a,b) a^n mpower(a,n) 数值矩阵的矩阵算法 矩阵的求逆及除法 求方阵的逆矩阵 inv(a) 求矩阵的伪逆矩阵 pinv(b) a左除b inv(a)*b a\b