matlab大学义PPT课件.ppt

.,1,MATLAB语言,哈尔滨工程大学理学院,MATLAB语言,.,2,MATLAB语言,哈尔滨工程大学理学院,研究偶极子(Dipole)的电势（Electricpotential）和电场强度（Electricfielddensity）。设在处有电荷，在处有电荷。那么在电荷所在平面上任何一点的电势和场强分别为其中又设电荷,.,3,MATLAB语言,哈尔滨工程大学理学院,一些命令：1.“showworkspace”命令是打开MATLAB的工作空间浏览器的命令。工作空间浏览器能够显示可进行操作的变量的集合以及各变量的各种属性。也可以用键入命令的方式Who仅显示变量名Whos显示同工作空间一样的信息2.ShowGraphicsPropertyEditor显示图形属性的编辑器,.,4,MATLAB语言,哈尔滨工程大学理学院,3.“setpath”命令该命令允许用户编辑m文件的目录。也可以用MATLAB提供的命令来完成。,.,5,MATLAB语言,哈尔滨工程大学理学院,数据格式表,.,6,MATLAB语言,哈尔滨工程大学理学院,以上的数据格式是在“Preference”对话框的“CommandWindow选项卡中设定。也可用命令直接给数据定义格式命令形式为：format类型此格式命令只改变显示的数据，并不改变MATLAB的存储数据时的格式。2.1MATLAB的运行方式有两种运行方式：一、直接交互命令行操作方式二、.m文件编辑运行操作方式,.,7,MATLAB语言,哈尔滨工程大学理学院,在MATLAB中.m文件分为两种：一、类似DOS的批处理文件，包含了一连串的MATLAB的指令，既不接受输入参数，也不返回输出参数，称为脚本（script）。二、函数（Function），函数接受用户的输入参数，执行函数的具体运算，最后输出结果。1.脚本文件：包括注解和指令两部分。注解部分必须以“”开头。脚本文件使用与存放的变量和数据都在工作空间中，指向工作空间来操作。2.函数文件：该文件主要包含以下五部分,.,8,MATLAB语言,哈尔滨工程大学理学院,（1）函数定义行：位于函数的第一行，格式为：functionA1,A2,A3,=函数名（B1,B2,B3,)A1,A2,A3,为输出参数，B1,B2,B3,为输入参数如无输入或输出参数，可以不写相应的参数。例如：functiontide(x1)functiony1,y2=tide()注：函数文件的名称必须与函数名称或函数文件中的主函数名称相同。（2）H1行：H1行紧接着函数定义行简要说明函数名和功能，是帮助文本的第一行，以“”开始。lookforfunction-name命令找到相关函数，但只显示H1行,.,9,MATLAB语言,哈尔滨工程大学理学院,（3）函数说明：函数说明位于H1行后，用来具体说明函数的功能和用法，同样以“”开头。Helpdirection-name就可显示该目录下所有函数的H1行（4）函数主体：函数主体在函数说明后，是所有进行具体运算和赋值程序的代码。（5）注解：在MATLAB中，还可以对某一行的具体内容进行注解。注：MATLAB通过判断.m文件的第一行是否有关键字“function”来识别该文件是否为函数文件。如果一个m文件是脚本文件，那么在该文件中不允许再出现函数定义语句。但可以调用在同一工作空间的函数文件中定义的函数,.,10,MATLAB语言,哈尔滨工程大学理学院,MATLAB数据类型MATLAB定义了六种数据类型：字符char双精度数值double稀疏数据sparse存储型storage单元数组cell结构数组structstorage又包括：int8（8位整形）unit8（无符号8位）int16unit16int32unit32该类型只用于内存的有效存储，可以对这种类型的数据进行基本操作，但不能对它们执行任何数学运算，如进行要转换数据类型。一般常用的为：chardouble,.,11,MATLAB语言,哈尔滨工程大学理学院,变量与常量变量MATLAB语言并不要求对所使用变量进行事先声明，也不需要指定变量类型，它会自动根据所赋予变量的值或对变量所进行的操作来确定变量的类型；在赋值过程中，如果变量已经存在，MATLAB语言将使用新值代替旧值，并以新的变量类型代替旧的变量类型。命名原则*变量名区分大小写；*变量名长度不超过31位；*变量名以字母开头，变量名中可以包含字母、数字、下划线，但不能使用标点。,.,12,MATLAB语言,哈尔滨工程大学理学院,在没有特殊声明的情况下，MATLAB语言将所识别的一切变量视为局部变量。若要定义全局变量，在给该变量前加关键字global。一般全局变量常用大写的英文字符表示。,MATLAB语言中的常量,.,13,MATLAB语言,哈尔滨工程大学理学院,字符、字符串1.在MATLAB中，所有的字符串都用单引号设定后输入或赋值。2.字符串的每个字符（包括空格）都是字符数组的一个元素。3.在MATLAB中，字符串和字符数组（或矩阵）基本上是等价的。4.字符数组的生成：函数char可以用来生成字符数组（或矩阵）。5.字符串和数组之间的转换,.,14,MATLAB语言,哈尔滨工程大学理学院,字符串转换为数值代码，此功能可由函数double来实现；字符数组转换为字符串，此功能可由函数cellatr来实现；数值数组和字符串之间的转换函数表,.,15,MATLAB语言,哈尔滨工程大学理学院,5.执行字符串执行字符串的功能在MATLAB中由函数eval来实现向量向量的生成直接输入向量格式：a=s1,s2,s3注意：用空格和逗号分隔生成行向量；用分号分隔生成列向量,.,16,MATLAB语言,哈尔滨工程大学理学院,利用冒号表达式生成向量格式：a=s1:s2:s3S1为起始值，s2为步距，s3为终止值如果省略s2的值，则步距取默认值1线性等分向量的生成格式：d=linspace(x1,x2,n)作用：生成通过维数定义的等间隔为（x2-x1）/n的向量对数等分向量的生成格式：d=logspace(x1,x2,n)作用：生成对数等分向量。d(1)=10 x1,d(n)=10 x2,.,17,MATLAB语言,哈尔滨工程大学理学院,2.向量的运算点积：由函数dot来实现格式：dot(a,b)或dot(a,b,dim)叉积:由函数cross实现格式：d=cross(a,b,dim)作用：返回a和b的dim维向量的叉积。a和b必须同维。如果dim不写，则默认dim3混合积：两函数套用例：dot(a,cross(b,c),.,18,MATLAB语言,哈尔滨工程大学理学院,2.4矩阵与数组MATLAB中最重要的功能就是对矩阵的运算。此时矩阵可以是数学意义上的矩阵，也可以是标量或向量。标量为11矩阵，向量为一个一维的矩阵。MATLAB可以进行实、复数矩阵的数学计算。矩阵的输入规则：矩阵的元素要写在方括号中；每行内的元素用逗号或空格分开；行与行之间用分号分开或用回车符分隔。例如：x=1,2,3;4,5,6;7,8,9,.,19,MATLAB语言,哈尔滨工程大学理学院,4.矩阵大小可不预先定义；5.矩阵元素可为运算表达式；6.无任何元素的空矩阵亦合法；7.若不想获得中间结果，可以“;”结束。,矩阵的运算矩阵的四则运算矩阵的加减“”、“”：要求加减的两矩阵同阶对应元素相加；标量与矩阵的加减法是矩阵的每个元素与该标量相加减。,.,20,MATLAB语言,哈尔滨工程大学理学院,矩阵的乘法，“*”矩阵的除法左除x=AB是方程A*X=B的解右除x=B/A是方程X*A=B的解其中B/A=(AB)2.矩阵与常数间的运算是同此矩阵的各元素之间进行运算。3.矩阵的逆运算格式：D=inv(a),.,21,MATLAB语言,哈尔滨工程大学理学院,4.矩阵的乘方运算格式：AB作用:(1)如果A和B都是矩阵，则不能运算；(2)如果B是个正整数，结果是A自乘B次；(3)如果B是个负整数，则首先对A取逆，然后它自乘B次(4)如果B是个标量，则对A中每个元素B次方。矩阵的转置格式：A矩阵为复矩阵，则转置为共轭转置。,.,22,MATLAB语言,哈尔滨工程大学理学院,数组的运算数组的加减：对应元素相加减，符号数组的乘法：对应元素相乘，要求A，B必须同维，除非一个为标量。格式：A.*B3.数组的除法：同乘法格式：右除A./B左除A.B4.数组的乘方：格式：A.B说明：如果A和B都是数组，则表示且A和B长度相同；如果B为标量，则,.,23,MATLAB语言,哈尔滨工程大学理学院,在MATLAB中矩阵的输入有下列方式：显示列表输入；外部文件加载；在.m文件中创建；通过MATLAB函数库生成。注：输入时，元素可以为常数，也可以为变量，表达式或函数。其维数可以扩大或缩小。,.,24,MATLAB语言,哈尔滨工程大学理学院,显示列表输入：按矩阵的格式规定输入，对于简单的矩阵非常方便。对复杂的矩阵，可将矩阵元素分行输入，此时用回车键代替分号。外部数据文件加载：用户可以使用load命令加载外部数据文件来创建矩阵。例如将生成的矩阵保存为file.mat文件，然后在命令窗口中用loadfile.mat加载在.m文件中创建：如果用户将一个数据文件保存为.m文件，则用户在命令窗口中输入文件名，将显示矩阵。通过MATLAB函数库生成：MATLAB为用户提供了四个产生基本矩阵的函数和一些能够产生特殊矩阵的函数。,.,25,MATLAB语言,哈尔滨工程大学理学院,四个产生基本矩阵的函数：zeros函数ones函数rand函数randn函数产生特殊矩阵的函数有：magic函数diag函数triu函数eye函数size函数（1）zeros函数作用：能够生成元素全为零的矩阵格式：矩阵名zeros（m，n）（2）ones函数作用：能够生成元素全为1的矩阵,.,26,MATLAB语言,哈尔滨工程大学理学院,格式：矩阵名ones（m，n）（3）rand函数作用：能够生成统一分布随机元素的矩阵，元素数值范围为（0.0,1.0）格式：矩阵名rand（m，n）（4）randn函数作用：能够生成正态分布随机元素的矩阵，元素数值范围为（0.0,1.0）格式：矩阵名rand（m，n）,.,27,MATLAB语言,哈尔滨工程大学理学院,（5）magic函数作用：能够自动生成一种特殊的n阶矩阵，这些矩阵的行、列和对角线上的元素之和都相等。格式：矩阵名magic（n）（6）diag函数作用：自动生成指定矩阵的对角矩阵。格式：Y=diag（D）（7）triu函数作用：自动生成指定矩阵的上三角矩阵。格式：Y=triu（D）,.,28,MATLAB语言,哈尔滨工程大学理学院,（8）tril函数作用：自动生成指定矩阵的下三角矩阵。格式：Y=tril（D）（9）eye函数和size函数作用：size函数用来返回指定矩阵的行数和列数。eye函数产生指定行数和列数的单位矩阵。两个函数共同使用可产生与指定矩阵同阶的单位矩阵。格式：Yeye（size（D）（10）repmat函数,.,29,MATLAB语言,哈尔滨工程大学理学院,作用：创建mn维块矩阵，每一块为已知矩阵。格式：Yrepmat（D，m，n）或Yrepmat（D,m,n）矩阵组合一个矩阵中的元素也可以时一个小矩阵。但此时应该注意维数的问题。（11）reshape函数格式：Yreshape（D,m,n）作用：用来对矩阵D进行变维。m，n为变维后的矩阵的行数和列数，mn必须与原矩阵的元素个数相等。,.,30,MATLAB语言,哈尔滨工程大学理学院,冒号表达式给行向量赋值格式：a=s1:s2:s3s1为起始值，s2为步距，s3为终止值。如果s2为负值，则要求s1大于s3的值，否则产生一个空向量。如果省略s2的值，则步距取默认值1注：冒号主要用于以下几个方面（1）输入行向量；（2）从向量、矩阵和数组中挑选指定的行、列或元素；（3）产生简单的表格。,.,31,MATLAB语言,哈尔滨工程大学理学院,矩阵的基本函数运算特征值函数eig和eigs格式：Beig(A)返回矩阵A的特征值，并依次构成B的列向量。格式：F,D=eig(A)返回两个矩阵F和D，其中F为特征向量矩阵，D为特征值矩阵。格式：F,D=eig(A,nobalance)同上，只是这种方式不采用预先平衡的方式。,.,32,MATLAB语言,哈尔滨工程大学理学院,格式：Deig(A,B)若B是一个与A同阶的方阵，则该格式将返回一个表示广义特征值的向量。eigs函数是使用迭代法求解特征值和特征向量的函数。2.秩函数：格式：rank(D)3.迹函数：矩阵所有对角线上元素的和称为矩阵的迹。格式：trace(D),.,33,MATLAB语言,哈尔滨工程大学理学院,特殊矩阵空阵在MATLAB中定义为空阵。一个被赋予空阵的变量具有一下性质：在MATLAB工作内存中确实存在被赋空阵的变量；空阵中不包括任何元素，它的阶数是0 x0;空阵可以在MATLAB的运算中传递；可以用clear从内存中清除空阵变量。注意：空阵不是“0”，也不是“不存在”。它可以用来使矩阵按要求进行缩维。例如：a(:,2,5)=,.,34,MATLAB语言,哈尔滨工程大学理学院,多项式运算多项式的表示方法：,用一下的行向量表示：P=a0,a1,，an用函数poly2sym可以将多项式向量表示成为符号形式。系数向量的直接输入法。特征多项式输入法由函数poly实现由根创建多项式由函数poly实现,.,35,MATLAB语言,哈尔滨工程大学理学院,2.求多项式的值两种形式：一种在输入变量值代入多项式计算时是以数组为单元，此时的计算函数为polyval；另一种是以矩阵为计算单元，进行矩阵式运算，此时的计算函数为polyvalm。这两种计算在数值上有很大的差别，这主要源于矩阵计算和数组计算的差别。3.求多项式的根一种是直接调用MATLAB的函数roots，求多项式的所有根；另一种是通过建立多项式的伴随矩阵再求其特征值的方法得到多项式的所有根。4.多项式的乘除法运算多项式的乘法由函数conv实现；,.,36,MATLAB语言,哈尔滨工程大学理学院,多项式的除法由函数deconv来实现。5.多项式的微分函数polyder6.多项式的拟合函数polyfit调用格式：Polyfit(X,Y,n)其中X，Y为拟合数据，n为拟合多项式的阶数。p,s=polyfit(X,Y,n)其中p为拟合多项式系数向量，s为拟合多项式系数向量的结构信息。,.,37,MATLAB语言,哈尔滨工程大学理学院,.,38,MATLAB语言,哈尔滨工程大学理学院,符号运算功能,符号表达式的生成与字符串的生成方式一样。但这种方式对空格很敏感。因此，不要在字符间乱加空格符。用sym命令生成格式：f=sym(表达式）3.用syms命令生成格式：syms变量f表达式,.,39,MATLAB语言,哈尔滨工程大学理学院,符号函数的运算复合函数运算格式：compose(f,g)返回当f=f(x)和g=g(y)时的复合函数f(g(y)格式：compose（f,g,z)返回的复合函数以z为自变量。2.反函数的运算格式：g=finverse(f),.,40,MATLAB语言,哈尔滨工程大学理学院,符号矩阵的运算此运算和矩阵的运算基本一致符号矩阵的简化因式分解格式：factor（s）2.符号矩阵的展开格式：expand(s)3.同类式合并格式：collect(S,v),.,41,MATLAB语言,哈尔滨工程大学理学院,将符号矩阵S中的各元素的v的同幂项系数合并。格式：collect(S)对默认变量进行同类项合并4.符号简化格式：simple（S)寻找符号矩阵或符号表达式的最简型格式：simplify(S)简化符号矩阵的每一个元素,.,42,MATLAB语言,哈尔滨工程大学理学院,5.分式通分格式：n,d=numden(A)把A的各元素转换为分子和分母都是整系数的最佳多项式型。n为分母，d为分子6.符号极限格式：limit(F,x,a)计算符号表达式F在xa条件下的极限值7.符号积分格式：int(S),int(S,v),