中南大学matlab总复习.ppt

科学计算与MATLAB 主讲：唐建国 中南大学材料科学与工程学院 2013.11 第一章 绪 论 1 科学计算 数值计算的定义 数值计算的研究过程 数值计算的误差分析 数值计算的稳定性与收敛性 2 MATLAB简介 MATLAB的类别 MATLAB的概述 3 课程的要求与学习方法 课程的内容 课程的要求 课程的学习与考核 教材与主要参考书 小结 如果存在一适当小的正数εr ，使得 则称εr为相对误差限。 例：x 15, ε x 2; εr x 2/15 13.33% y 1000,ε y 5 , εr y 5/1000 0.5% 2.有效数字： 定义：若近似值x的误差限是某一位的半个单位，该位到x的第一位非零数字共有n位，就说x有n位有效数字，x可表示为 其中，a1≠0.且 a1, a2, …, an 都是 0 ～ 9 中的任一整数。 其绝对误差限满足： 例. 测量一物体的长度为954cm，问测量数据的相对误差限多大？ 解:因实际问题所截取的近似数，其绝对误差限一般不超过最小刻度的半个单位，故当x 954cm时，有ε x ＝0.5cm， 而x的相对误差 er x ≤0.5/954 0.0005241…. 0.00053 0.053 % 故εr x 0.053 %. 第二讲 ＭＡＴＬＡＢ基础知识（数值、符号计算） 1 MATLAB的启动与运行 Matlab的启动和退出 Matlab的工作窗口 Matlab的指令行的操作 Matlab的帮助系统 2 MATLAB的矩阵与数值计算 Matlab数据类型 Matlab的变量及其命名 Matlab矩阵及其运算 Matlab的数组关系/逻辑运算 Matlab的多项式运算 3 MATLAB的符号计算 小结 数值显示格式 任何MATLAB的语句的执行结果都可以在屏幕上显示，同时赋值给指定的变量，没有指定变量时，赋值给一个特殊的变量ans，数据的显示格式由format命令控制。 format只是影响结果的显示，不影响其计算与存储；MATLAB总是以双字长浮点数（双精度）来执行所有的运算。 如果结果为整数，则显示没有小数；如果结果不是整数，则输出形式有： format short ：短格式（5位定点数）99.1253 format long：长格式（15位定点数 ） 99.1座机电话号码00000 format short e：短格式e方式 9.9123e+001 format long e：长格式e方式 9.座机电话号码e+001 format bank：2位十进制 99.12 format hex：十六进制格式 如：a [1 2;3 4]；b [ 3 5; 5 9] 》c a+b d a-b 》c d 4 7 -2 -3 8 13 -2 -5 》a*b [13 23; 29 51] 》a/b [-0.50 0.50;3.50 –1.50] 》a\b [-1 -1;2 3] 》a^3 [37 54; 81 118] 》a.*b [3 10;15 36] 》a./b [0.33 0.40;0.60 0.44] 》a.\b [3.00 2.50;1.67 2.25] 》a.^3 [1 8; 27 64] 转置：对于实矩阵用（’）符号或（.’）求转置结果是一样的；然而对于含复数的矩阵，则（’）将同时对复数进行共轭处理，而 （.’）则只是将其排列形式进行转置。 逆矩阵与行列式计算 求逆：inv A ； 求行列式：det A 要求矩阵必须为方阵 5、矩阵分解 （1）奇异值分解 [U,S,V] svd A 例：a 9 8 6 8 可以验证： u*u’ I v*v’ I u*s*v’ a （2）特征值分解 [V,D] eig A 例： a 9 8 6 8 [v,d] eig a v 0.7787 -0.7320 0.6274 0.6813 d 15.4462 0 0 1.5538 （3）正交分解 [Q,R] qr A 例： a 9 8 6 8 [q,r] qr a q -0.8321 -0.5547 -0.5547 0.8321 r -10.8167 -11.0940 0 2.2188 （4）三角分解 [L,U] lu A 将A做对角线分解，使得A L*U,其中L为下三角矩阵，U为上三角矩阵。 注意：L实际上是一个“心理上”的下三角矩阵，它事实上是一个置换矩阵P的逆矩阵与一个真正下三角矩阵L1（其对角线元素为1）的乘积。 [L1,U1,P] lu A 例：a [1 2 3;4 5 6;7 8 9] 比较： [l1,u1,p] lu a [l,u] lu a