[教育]第三章 MATLAB数值运算.ppt

第三章 MATLAB 数值运算 主要内容 3.1 向量及其运算 向量是矢量运算的基础 行向量 列向量 3.1.1 向量的构造 1．逐个输入 a=[1 3 9 10 15 16] %采用空格和逗号分隔构成行向量 b=[1; 3; 9; 10; 15; 16] %采用分号隔开构成列向量 2．利用冒号表达式“:”生成向量 x=1:2:9 %初值=1，终值=9，步长=2 z=1:5 %初值=1，终值=5，默认步长=1 3．利用函数生成向量 x=linspace(1,10,5) %初值=1，终值=10，元素数目=5 x=logspace(1,3,3) %初值=10，终值=1000，元素数目=3 3.1.2 向量的读取 例： a=linspace(1,10,10) a(3) a(2:7) a(8:-2:3) a(6:end) length(a) size(a) 3.2 矩阵及其运算 MATLAB 以矩阵为基本的运算单元，向量和标量作为特殊的矩阵处理：向量看作只有一行或一列的矩阵；标量看作只有一个元素的矩阵。 3.2.1 矩阵的构造 通过直接输入矩阵的元素构造矩阵： 用中括号[ ]把所有矩阵元素括起来 同一行的不同数据元素之间用空格或逗号间隔 用分号（；）指定一行结束 可分成几行进行输入，用回车符代替分号 数据元素可以是表达式，系统将自动计算结果 通过数组（矩阵）编辑器 通过M文件创建矩阵 利用MATLAB函数创建矩阵 标准矩阵生成函数 3.2.2 矩阵下标与子矩阵提取 A(m, n) 提取第m行，第n列元素 A(:, n) 提取第n列元素 A(m, :) 提取第m行元素 A(m1:m2, n1:n2) 提取第m1行到第m2行和第n1 列到第n2列的所有元素 A(m:end, n) 提取从第m行到最末行和第n列的子块 A(:) 得到一个长列矢量，该矢量的元素按矩 阵的列进行排列 例1： 修改矩阵A中元素的数值 A=[1 2 3 4; 5 6 7 8; 9 10 11 12; 13 14 15 16]; A(1,1)=0;A(2,2)=A(1,2)+A(2,1);A(4,4)=cos(0); A = 0 2 3 4 5 7 7 8 9 10 11 12 13 14 15 1 3.2.3 矩阵的特殊化处理函数 reshape 在总元素数不变的前提下，改变矩阵的行 数和列数 flipud 以矩阵“水平中线”为对称轴，交换上下对 称位置上的矩阵元素 fliplr 以矩阵“垂直中线”为对称轴，交换左右对 称位置上的矩阵元素 rot90 把矩阵逆时针旋转90°的整数倍 用于专门学科的特殊矩阵 魔方矩阵 魔方矩阵有一个有趣的性质，其每行、每列及两条对角线上的元素和都相等。对于n阶魔方阵，其元素由1,2,3,…,n2共n2个整数组成。MATLAB提供了求魔方矩阵的函数magic(n)，其功能是生成一个n阶魔方阵。 用于专门学科的特殊矩阵 例2：将101～125等25个数填入一个5行5列的表格中，使其每行每列及对角线的和均为565。 用于专门学科的特殊矩阵 范得蒙矩阵 范得蒙德(Vandermonde)矩阵最后一列全为1，倒数第二列为一个指定的向量，其他各列是其后列与倒数第二列的点乘积。可以用一个指定向量生成一个范得蒙矩阵。在MATLAB中，函数vander(V)生成以向量V为基础向量的范得蒙矩阵。例如，A=vander([1;2;3;5]) 希尔伯特（Hilbert）矩阵 希尔伯特矩阵的第i行、第j列的元素值为1/(i+j-1)在MATLAB中，生成希尔伯特矩阵的函数是hilb(n)。 使用一般方法求逆会因为原始数据的微小扰动而产生不可靠的计算结果。MATLAB中，有一个专门求希尔伯特矩阵的逆的函数invhilb(n)，其功能是求n阶的希尔伯特矩阵的逆矩阵。 例 求4阶希尔伯特矩阵及其逆矩阵。命令如下：format rat %以有理形式输出H=hilb(4)H=invhilb(4) 3.2.4