matlab_ch2_数值计算-104()分析.ppt

2 Matlab的数值计算 Matlab 具有出色的数值计算能力，占据世界上数值计算软件（Maple、Mathmatica、MathCAD）的主导地位。 数值运算的功能 创建矩阵 矩阵运算 多项式运算 线性方程组 一、矩阵的创建 Matlab中矩阵是运算的基本单元，且定义在复数域上的 矩阵元素可以是任何matlab表达式 ，可以是实数 ，也可以是复数，复数可用特殊函数i或j 表示 Matlab中所有的矩阵事先都不必定义维数大小，系统会根据用户的输入自动配置，在运行中自动调整矩阵的维数。 一、矩阵的创建（cont.） 矩阵可以通过以下几种形式创建： 以直接列出元素形式输入； 通过语句和函数产生； 从外部文件装入； 建立在M文件中。 1.直接输入 ① 当需要的矩阵维数比较小。 ② 矩阵元素必须用[ ]括住。 ③ 矩阵元素必须用逗号或空格分隔。 ④ 在[ ]内矩阵的行与行之间必须用分号或用enter键分隔。 1.直接输入（Cont.） a=[1 2 3;4 5 6] b=[1,2,3;4,5,6] c=[1,2,3 4,5,6] 矩阵中的元素可以数字或表达式，但表达式中不可含未知变量。x=[2 pi/2;sqrt(3) 3+5i] 符号的作用 逗号和分号的作用 ? matlab允许多条语句在同一行出现，逗号和分号可作为指令间的分隔符。 ?分号如果出现在指令后，屏幕上将不显示结果。 符号的作用（cont.） 符号的作用（cont.） 冒号的作用（n:s:m） ?产生从n到m步长（公差）为s的一系列的值。s可以是正数、负数、小数，默认值是1。 ?用于选出矩阵指定行、列的元素。 ?循环语句。 A=[1:3;1:2:5;0.3:0.1:0.5] A = 1.0000 2.0000 3.0000 1.0000 3.0000 5.0000 0.3000 0.4000 0.5000 符号的作用（cont.） A=[1 2 3;4 5 6;7 8 9] 例：利用for循环求1！+2！+3！+ ??? +5！的值 sum=0; for i=1:5 pdr=1; for k=1:i pdr=pdr*k; end sum=sum+pdr; end 2.矩阵编辑器 当输入的矩阵很大，不适合手工输入。Matlab提供了一个矩阵编辑器（matrix editor）。 需要预先定义一个变量。 打开workspace窗口，选中需要编辑的变量，双击，就可以打开矩阵编辑器。用户可以改变矩阵的大小和元素的值。 打开 desktop\workspace窗口 双击，打开矩阵编辑器 3.由函数创建矩阵 Matlab提供了大量的函数用于创建一些常用的矩阵。 对角阵、单位阵、零矩阵等在数值计算中常用的矩阵，以及由它们经变换或截取产生的新矩阵。 zeros(m,n) 零矩阵。 ones(m,n) 元素全部为1的矩阵。 eye(m,n) 单位阵。 rand(m,n) 正态分布的随机矩阵。 diag(M) ,M是一个矩阵，diag()函数取矩阵的对角元素产生一个列向量；M是一个列向量则产生一个对角阵。 4.从外部文件装入 *.mat文件， save ， load命令。 save mydata.mat 把工作空间中所有变量存入文件。 save yourdata.mat a b c 把工作空间中a b c变量存入文件。 load mydata.mat 把文件中的变量读入到工作空间中。 load mydata .mat a b c 把文件中的变量a b c读入到工作空间中。 *.xls文件，首先创建一个excel文档，然后再matlab命令窗口中用xlsread(‘filename’)调用。 在office excel环境下编辑保存book.xls文件。 在matlab环境下读取数据xlsread(‘book.xls’) 5.复数矩阵的创建 在matlab中复数单位为i=sqrt(-1)，在工作空间中显示为0+1.0000i 。 复数可由下面语句生成：z=a+b*i（或z=a+bi）。 当复数作为矩阵元素时，复数内不能留有空格，否则会被认为是两个分开的数。 A= [1+5*i,2+6*i;3+7*i,4+8*i] B= [1,2;3,4]+i*[5,6;7,8] C=[1+5i,2+6i;3+7i,4+8i] 6.矩阵的子阵 矩阵的子阵可以通过标