matlab2数值计算new1.pptVIP

第二讲 MATLAB的数值计算 —— matlab 具有出色的数值计算能力，占据世界上数值计算软件的主导地位 内容 矩阵及其运算 多项式运算 代数方程组求解 线性插值 数据分析与统计 矩阵及其运算 创建矩阵的方法 直接输入法 规则： 矩阵元素必须用[ ]括住 矩阵元素必须用逗号或空格分隔 在[ ]内矩阵的行与行之间必须用分号分隔 矩阵元素 矩阵元素可以是任何表达式(表达式中不能有未知变量），可以是实数 ，也可以是复数，复数可用特殊函数i，j 输入 a=[1 2 3;4 5 6] x=[2 pi/2;sqrt(3) 3+5i] 2.用matlab函数创建矩阵 空阵 [ ] — matlab允许输入空阵，当一项操作无结果时，返回空阵。 rand —— 随机矩阵 eye —— 单位矩阵 zeros ——全部元素都为0的矩阵 ones ——全部元素都为1的矩阵 还有伴随矩阵、稀疏矩阵、魔方矩阵、对角矩阵、范 德蒙等矩阵的创建，就不一一介绍了 注意： matlab严格区分大小写字母。 matlab函数名必须小写。 3. 矩阵的修改 直接修改 可用?键找到所要修改的矩阵，用?键移动到要修改的矩阵元素上即可修改。 指令修改 可以用A(?,?)= ? 来修改。 例如 a=[1 2 0;3 0 5;7 8 9] a =1 2 0 3 0 5 7 8 9 4.数据的保存与获取 把matlab工作空间中一些有用的数据长久保存下来的方法是生成mat数据文件。 save —将工作空间中所有变量保存matlab.mat中 save data—将工作空间中所有变量保存为data.mat save data a b —将工作空间中变量a和b保存为data.mat save data a b -append—将变量a和b添加到data.mat save data a b -ascii—将变量a和b保存为data.mat的8位ascii文件 save data a b –ascii-double—将变量a和b保存为data.mat的16位ascii文件 用load指令调用已生成的mat文件。 load —— load data —— load data a b —— load data a b -ascii mat文件是标准的二进制文件，还可以ASCII码形式保存。 5.矩阵运算 矩阵加、减（＋,－）运算 规则： ? 相加、减的两矩阵必须有相同的行和列两矩阵对应元素相加减。 ? 允许参与运算的两矩阵之一是标量。标量与矩阵的所有元素分别进行加减操作。 . 矩阵乘（?）运算 规则： A矩阵的列数必须等于B矩阵的行数 标量可与任何矩阵相乘。 a=[1 2 3;4 5 6;7 8 0];b=[1;2;3];c=a*b c =14 32 23 . 矩阵乘方——a^p 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 a^0.5 ans = 0.4498 + 0.7623i 0.5526 + 0.2068i 0.6555 -0.3487i 1.0185 + 0.0842i 1.2515 + 0.0228i 1.4844 - 0.0385i 1.5873 - 0.5940i 1.9503 - 0.1611i 2.3134 + 0.2717i . 矩阵的其它运算 inv —— 矩阵求逆 det —— 行列式的值 eig —— 矩阵的特征值 ‘ —— 矩阵转置 sqrt —— 矩阵开方 .矩阵的一些特殊操作 矩阵的变维 a=[1:12];b=reshape(a,3,4) c=zeros(3,4);c(:)=a(:) 矩阵的变向 rot90:旋转; fliplr:左右翻; flipud:上下翻 矩阵的抽取 diag:抽取主对角线;tril: 抽取主下三角; triu:抽取主上三角 矩阵分析 矩阵行列式det 通常用来判断矩阵是否奇异 A=magic(3) A= 8 1 6 3 5 7 4 9 2 det(A) ans=