MATLAB介绍(2)ppt课件

1、MATLABMATLAB工具软件工具软件1;.2科学家和工程技术人员需要的研究工具科学家和工程技术人员需要的研究工具 一种演算纸式的科学工程计算语言一种演算纸式的科学工程计算语言 MATLAB3 1核心部分：几百个核心内部函数。 2大量可选的工具箱。MATLABMATLAB的两个主要部分的两个主要部分4MATLAB的基本组成wMATLAB主程序MATLAB语言、工作环境、句柄图形、数学函数库和应用程序接口wMATLAB的功能性工具箱 如： Simulink：动态系统仿真 wMATLAB的学科性工具箱 如：signal processing toolbox 信号处理工具箱 5例：用一个简单命令求

2、解线性系统 3x1+ x2 - x3 = 3.6 x1+2x2+4x3 = 2.1 -x1+4x2+5x3 = -1.4线性系统求解：（Ax=b）6线性系统求解：（Ax=b）输入指令：A=3 1 -1;1 2 4;-1 4 5;b=3.6;2.1;-1.4;x=Ab得到结果：x = 1.4818 -0.4606 0.384870123456- 1- 0 .8- 0 .6- 0 .4- 0 .200 .20 .40 .60 .81x=linspace(0,6)y1=sin(2*x),y2=sin(x.2),y3=(sin(x).2;plot(x,y1,x, y2,x, y3)例、用简短命令计算并

3、绘制在0 x6范围内的sin(2x)、sinx2、sin2x。8常用专业性工具箱：wControl System 最基本的控制工具箱wSignal Processing 信号处理工具箱wImage Processing 图像处理工具箱wCommunication 通讯工具箱wHigher-Order Spectral Analysis 高阶谱分析wSpline 样条分析工具箱wStatistics 统计分析工具箱w910MATLAB的部分通用命令的部分通用命令clear : 删除内存中的变量。删除内存中的变量。clc： 删除工作窗口的内容。删除工作窗口的内容。clf： 删除图形窗口的内容。删除

4、图形窗口的内容。11M Matlabatlab联机帮助命令联机帮助命令 help 功能提供功能提供matlab大部分主题的在线帮助信息大部分主题的在线帮助信息help 显示显示help 主题一览表主题一览表help plotxyz 显示有关三维做图指令帮助信息显示有关三维做图指令帮助信息12虽然虽然help可以随时提供帮助，但必须知道准确的函数名称。当不能确定函数名称可以随时提供帮助，但必须知道准确的函数名称。当不能确定函数名称时，时，help就无能为力了。就无能为力了。 Lookfor函数函数 它可提供通过一般的关键词，搜索出一组与之相关的命令它可提供通过一般的关键词，搜索出一组与之相关的命

5、令lookfor fourier 寻找含有傅立叶变换的相关指令寻找含有傅立叶变换的相关指令13变量查询函数变量查询函数who与与whos作用都是列出在作用都是列出在matlab工作空间中已经驻留的变量名清单工作空间中已经驻留的变量名清单不同的是不同的是whos在给出驻留变量的同时，还给出他们的维数及性质在给出驻留变量的同时，还给出他们的维数及性质14Matlab的演示功能的演示功能demo 在线演示在线演示15MATLAB的数值运算w数据类型 数据运算：双精度型（Double)w支持的数据类型： 数值型、字符型、单元型和结构型 16复数的表示：w i,j表示虚数单位 z=a+b*i 或 z=a

6、+b*j 当b是一个数字时，可以表示为 z=a+bi 或 z=a+bjMATLAB的矩阵和数组运算matlab 具有出色的矩阵运算能力，占据世界上数值计算软件的重要地位17;.18矩阵的创建:w创建矩阵的方法创建矩阵的方法w直接输入法规则： 矩阵元素必须用 括住 矩阵元素必须用逗号或空格分隔 在 内矩阵的行与行之间必须 用分号分隔19 矩阵元素可以是任何matlab表达式 ，可以是实数 ，也可以是复数，复数可用特殊函数I，j 输入 a=1 2 3;4 5 6 x=2 pi/2;sqrt(3) 3+5i 矩阵元素:20w注意：只要是赋过值的变量，不管是否在屏幕上显示过，都存储在工作空间中，以后可

7、随时显示或调用。变量名尽可能不要重复，否则会覆盖212.用matlab函数创建矩阵w空阵 matlab允许输入空阵，当一项操作无结果时，返回空阵。wrand 随机矩阵weye 单位矩阵wzeros 全部元素都为0的矩阵wones 全部元素都为1的矩阵22向量的生成:w利用冒号生成向量x=i:j 如果IJ，生成X空向量x=i:j:k j是递增步长23矩阵的标识和引用w矩阵的子矩阵通过向量、标量的标识来引用和赋值。 1.向量标识 A（u, v) 2. “01”向量标识 A（L）24例：w找出4阶魔方阵中大于10的元素。 A = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 1

8、5 125规则：规则：w所有字符串都用单引号括起来；所有字符串都用单引号括起来；w字符串中的每个字符都是字符串变量中的一个元素；字符串中的每个字符都是字符串变量中的一个元素；w字符串中的字符以字符串中的字符以ASCIIASCII码形式存储并区分大小码形式存储并区分大小字符串数组字符串数组26w字符串数组的建立。字符串数组的建立。 例：例： a=This is a test. A=中国中国 成都成都 B=China 中國中國 27字符串数组的建立字符串数组的建立w有效字符串连成长字符串有效字符串连成长字符串w多行字符串数组的直接创建多行字符串数组的直接创建w利用字符串操作函数创建字符串数组利用字

9、符串操作函数创建字符串数组w利用转换函数建立字符串数组利用转换函数建立字符串数组w字符串数组的元素标识字符串数组的元素标识28符号的作用: “ ,” 逗号的作用 逗号可作为指令间的分隔符，matlab允许多条语句在同一行出现。 29w把matlab工作空间中一些有用的数据长久保存下来的方法是生成mat数据文件。 save 将工作空间中所有的变量存到matlab.mat文件中。数据的保存与获取默认文件名30save 将工作空间中所有的变量存到.mat文件中。save a b 将工作空间中a和b变量存到文件中。31load load load a b mat文件是标准的二进制文件，还可以ASCII

10、码形式保存。即可恢复保存过的所有变量32w矩阵加、减（,）运算规则： 相加、减的两矩阵必须有相同的行和列两矩阵对应元素相加减。 允许参与运算的两矩阵之一是标量。标量与矩阵的所有元素分别进行加减操作。第二章 矩阵运算及应用矩阵运算及应用332. 矩阵乘（）运算规则：wA矩阵的列数必须等于B矩阵的行数w标量可与任何矩阵相乘。a=1 2 3;4 5 6;7 8 0;b=1;2;3;c=a*bc =14 32 23 34 矩阵除的运算在线性代数中没有，有矩阵逆的运算，在matlab中有两种矩阵除运算。 矩阵左除：AB 定义为：X=AB AX=B 矩阵右除：B/A 定义为：X=B/A XA=B (B/A

11、)=(AB)35 a p a 自乘p次幂 方阵方阵整数整数3. 矩阵乘方当p为正整数时，表示方阵自乘p次；当p为负整数时，表示方阵自乘p次后的逆；当p为零时，为方阵同维的单位阵。36 a=1,2,3;4,5,6;7,8,9;a2 ans =30 36 42 66 81 96 102 126 15037winv 矩阵求逆wdet 行列式的值weig 矩阵的特征值wdiag 对角矩阵 w 矩阵转置wsqrt 矩阵开方4. 矩阵的其它运算 38 数组运算指元素对元素的算术运算，与通常意义上的由符号表示的线性代数矩阵运算不同w 数组加减(.+,.-) a.+b a.- b数组运算数组运算 对应元素相加

12、减（与矩阵加减等效）对应元素相加减（与矩阵加减等效）392. 数组乘除(，./，.）ab a，b两数组必须有相同的行 和列两数组相应元素相乘。a=1 2 3;4 5 6;7 8 9;b=2 4 6;1 3 5;7 9 10;a.*bans = 2 8 18 4 15 30 49 72 90 40a=1 2 3;4 5 6;7 8 9;b=2 4 6;1 3 5;7 9 10;a*bans = 25 37 46 55 85 109 85 133 172 41a./b=b.aa.b=b./aa./b=b.a 都是都是a的元素被的元素被b的对应元的对应元 素除素除a.b=b./a 都是都是b的元素被

13、的元素被a的对应元的对应元 素除素除例例: a=1 2 3;b=4 5 6; c1=a.b; c2=b./ac1 = 4.0000 2.5000 2.0000c2 = 4.0000 2.5000 2.0000 给出a,b对应元素间的商.42例例:a=1 2 3;b=4 5 6;z=a.2z = 1.00 4.00 9.00z=a.bz = 1.00 32.00 729.003. 数组乘方数组乘方(.) 元素对元素的幂元素对元素的幂43关系操作和逻辑操作关系操作和逻辑操作w关系操作符：关系操作符：w、=、 = =(等于等于)、 =(不等于不等于)w关系运算规则：关系运算规则：w两个变量都是标量，

14、则结果为真两个变量都是标量，则结果为真(1)或假或假(0)。w两个变量都是数组，则必须大小相同，结果也是同样大小的数组，数组的元素为两个变量都是数组，则必须大小相同，结果也是同样大小的数组，数组的元素为0或或1。44关系操作关系操作w一个数组和一个标量，则把数组的每个元素分别与标量比较，结果为与数组大小相一个数组和一个标量，则把数组的每个元素分别与标量比较，结果为与数组大小相同的数组，数组的元素为同的数组，数组的元素为0或或1。w、=，仅对参加比较变量的实部进行比较，仅对参加比较变量的实部进行比较，w= = 和和 = ，则同时对实部和虚部进行比较。，则同时对实部和虚部进行比较。45逻辑操作逻辑

15、操作w逻辑操作符：w&(与)、|(或)、(非)和xor(异或)。46a ba & ba | baxor(a,b)00001001011110010111110047运算符优先级运算符优先级w(矩阵转置矩阵转置)、(矩阵幂矩阵幂)和和.(数组转置数组转置)、.(数组幂数组幂)w (逻辑非逻辑非) w *(乘乘)、/(左除左除)、(右除右除)和和.*(点乘点乘)、./(点左除点左除)、.(点右除点右除) w+、-(加减加减)w: (冒号冒号) w、=、= w&(逻辑与逻辑与) 48多项式的表示和创建 在Matlab中，任意的多项式都是用一个行向量表示，将多项式的系数按降幂排

16、列存放在行向量中。多项式p(x)=a0 xn+a1xn-1+an-1x+an对应的行向量为:P=a0 a1 an。49w直接创建多项式例：lp(x)=2x4+5x3-8x+32 5 0 -8 350通过根创建多项式polyw其调用格式如下：lP=poly(A)l如A为向量，则P是多项式的系数行向量，该多项式的根A。51l注意：l如A为nn的矩阵，则P是一个长度为n+1的行向量，它是A的特征多项式的系数；例：A=6 -8 6;1 0 0;0 2 0;p=poly(A),roots(p),D=eig(A)52求根：wroots 调用方式： roots(A)53注意：w求根的精度问题： 例：y=(x

17、-1)6 =x6-x5+15x4-20 x3+15x2-6x+154wa=1 -6 15 -20 15 -6 1wa =w 1 -6 15 -20 15 -6 1w roots(a)wans =w 1.0042 + 0.0025iw 1.0042 - 0.0025iw 1.0000 + 0.0049iw 1.0000 - 0.0049iw 0.9958 + 0.0024iw 0.9958 - 0.0024i55多项式的基本运算w多项式的乘法convw多项式的乘法本质上是多项式系数向量的卷积运算。w卷积的定义为 K c(k)= a(i)b(k+1-i) i=1 w在Matlab，完成此功能的命令

18、是conv：56lC= conv(a,b)。例：求a=x3+3x2+2x1+1和b=4x3+3x2+9x1+10的乘积。a=1 3 2 1;b=4 3 9 10;c=conv(a,b)57w多项式的除法dconvw多项式的除法的运算实质就是多项式系数的解卷积运算。向量a对向量c进行解卷积得到的商向量q和余向量r，满足： knc(k)-r(k)= a(i)q(k+1-i) i=1n在Matlab，完成此功能的命令是dconv：lq,r= dconv(c,a)其中： c=conv(a,q)+r58例:a=1 3 2 1;c=4 15 26 47 51 29 10;q,r=deconv(c,a)59

19、w多项式的加减n对多项式所对应的向量利用数组加减运算规则执行加减运算。注：多项式所对应的向量如果大小不相等，必须用0补齐。例：求a=x3+3x2+2x1+1和b=4x3+3x2+9x1+10的和与差。a=1 3 2 1;b=4 3 9 10;c=a+b,c1=a-b,60w多项式的求导polyder 不仅可以计算单个多项式的导数，还可以计算两个多项式相乘和相除的导数。其调用方式如下：lpolyder(p):返回多项式系数向量p的导数；lpolyder(a,b):返回多项式a*b的导数；lq,d=polyder(b,a):返回b/a的导数，用q/d表示。其中q是结果的分子多项式，d是结果的分母多项式。61例：求a=x3+3x2+2x1+1和b=4x3+3x2+9x1+10的乘和商的导数。a=1 3 2 1;b=4 3 9 10;polyder(a,b)q,d=polyder(b,a)62w多项式的求值polyvaln在Matlab中，多项式求值有两种方式：按数组规则和按矩阵规则计算多项式值。ly=polyval(p,x):按数组规则运算。用来计算多项式p在自变量x处的值，可以是数、向量和矩阵。如果x是向量或矩阵，则该命令将对x的每个元素计算p的值，结果大小和x的大小相同；ly=polyvalm(p,x):按矩阵规则运算。自变量x是一方阵。63例：p=2 3 3 0 6;x