《2.Matlab 数值计算.pptVIP

第二讲 MATLAB的数值计算 —— matlab 具有出色的数值计算能力，占据世界上数值计算软件的主导地位 数值运算的功能 创建矩阵 矩阵运算 多项式运算 线性方程组 数值统计 线性插值 函数优化 微分方程的数值解 一、命令行的基本操作 创建矩阵的方法 直接输入法 2.用matlab函数创建矩阵 空阵 [ ] — matlab允许输入空阵，当一项操作无结果时，返回空阵。 rand —— 随机矩阵 eye —— 单位矩阵 zeros ——全部元素都为0的矩阵 ones ——全部元素都为1的矩阵 diag ——产生对角矩阵 也可用linspace函数产生行向量。其调用格式为： linspace(a, b, n) 其中a和b是生成向量的第一个和最后一个元素，n是元素总数。 例 》a=linspace(1 , 10 , 10) a= 　 1 2 3 4 5 6 7 8 9 10 还有伴随矩阵、稀疏矩阵、魔方矩阵(magic)、对角矩阵、范德蒙等矩阵的创建，就不一一介绍了。 注意：matlab严格区分大小写字母，因此a与A是两个不同的变量。 matlab函数名必须小写。 二、数据的保存与获取 把Matlab工作空间中一些有用的数据长久保存下来的方法是生成mat数据文件。 ?save —— 将工作空间中所有的变量存到matlab.mat文件中。 ?save data——将工作空间中所有的变量存到data.mat文件中。 ?save data a b ——将工作空间中a和b变量存到data.mat文件中。 下次运行Matlab时即可用load指令调用已生成的mat文件。 ?load —— ?load data —— ?load data a b —— mat文件是标准的二进制文件，还可以ASCII码形式保存。 三、矩阵运算 矩阵加、减（＋,－）运算 规则： ? 相加、减的两矩阵必须有相同的行和列两矩阵对应元素相加减。 ? 允许参与运算的两矩阵之一是标量。标量与矩阵的所有元素分别进行加减操作。 2. 矩阵乘（?）运算 规则： A矩阵的列数必须等于B矩阵的行数 标量可与任何矩阵相乘 a=[1 2 3;4 5 6;7 8 0];b=[1;2;3];c=a*b c =14 32 23 d=[-1;0;2]; f=pi*d f = -3.1416 0 6.2832 矩阵除的运算在线性代数中没有，有矩阵逆的运算，在matlab中有两种矩阵除运算。 两种除法：\和/，分别表示左除和右除。如果A矩阵是非奇异方阵，则A\B和B/A运算可以实现。 A\B等效于A的逆左乘B矩阵，而B/A等效于A矩阵的逆右乘B矩阵。 对于矩阵来说，左除和右除表示两种不同的除数矩阵和被除数矩阵的关系。对于矩阵运算，一般A\B≠B/A。 3. 矩阵乘方—— a^n,a^p,p^a a ^ p —— a 自乘p次幂 a=[1,2,3;4,5,6;7,8,9];a^2 ans =30 36 42 66 81 96 102 126 150 4. 矩阵的其它运算 inv —— 矩阵求逆 det —— 行列式的值 eig —— 矩阵的特征值 diag —— 对角矩阵 ’ —— 矩阵转置 sqrt —— 矩阵开方 7. 矩阵的数组运算 数组运算指元素对元素的算术运算， 与通常意义上的由符号表示的线性代数 矩阵运算不同。 数组加减(.+,.-) a.+b a.- b 2. 数组乘除(??，./，.\） a??b —— a，b两数组必须有相同的行 和列两数组相应元素相乘。 a=[1 2 3;4 5 6;7 8 9]; b=[2 4 6;1 3 5;7 9 10]; a.*b ans = 2 8 18 4 15 30 49 72 90 a=[1 2 3;4 5 6;7 8 9]; b=[2 4 6;1 3 5;7 9 10]; a*b ans = 25 37 46 55 85