第2讲 matlab数学运算.ppt

第二讲MATLAB的数值计算,matlab具有出色的数值计算能力，占据世界上数值计算软件的主导地位,数值运算的功能,创建矩阵矩阵运算多项式运算线性方程组数值统计线性插值函数优化微分方程的数值解,一、MATLAB中的变量类型,数值变量字符串变量单元型变量结构型变量,1.数值变量,MATLAB是以矩阵为基本运算单元的，而构成数值矩阵的基本单元就是数值。MATLAB中的变量名必须遵循：变量名区分大小写；变量名的长度不超过63位，超过时给出警告信息；变量名必须以字母开头，其余可包含字母、数字、下划线，但不得使用标点符号。,MATLAB中变量的作用域一般默认为局部变量，仅在当前调用的M文件中有效。如果要定义全局变量，则必须用global来声明。一般情况下，为了和局部变量有所区别，常将全局变量用大写字母表示，但这并不是必需的，只是人为的一种约定而已。,比如以下的例子就是输入了变量x和X，这是两个不同的变量，一个是数字，一个是矩阵。,2.字符串变量,MATLAB中的字符串运算功能非常丰富，特别是符号运算功能的加入，使得字符串函数功能得到了极大增强。在MATLAB中，所有字符串都用英文单引号标识，字符串和字符数组是等价的，字符串中的每个字符(包括空格)都是字符数组的一个元素。,字符串与数值数组的相互转换详情请参阅帮助(help),字符串的操作与执行详情请参阅帮助(help),串检验与进制转换函数详情请参阅帮助(help),3.单元型变量,单元型变量是MATLAB中较为特殊的一种数据类型，本质上也是一种数组，但这种数组和传统数组的区别是：传统数组中所有元素只能是同一种数据类型，而单元型数组可以把不同的数据类型组合在一起，从而形成一种比较复杂的数组。,(1)单元型变量的定义,直接赋值法：单元型变量使用大括号标识，元素之间用逗号分隔。也可以直接对单元型变量的元素直接赋值单元型变量的下标用大括号索引。如A3，B1，5等。由cell函数预先分配存储空间，然后逐个元素进行赋值。例如：Bcell(2,3)，可在内存空间中建立一个单元型空变量B，然后可逐个对其每个元素赋值。,(1)单元型变量的定义,直接赋值法：单元型变量使用大括号标识，元素之间用逗号分隔。也可以直接对单元型变量的元素直接赋值单元型变量的下标用大括号索引。如A3，B1，5等。由cell函数预先分配存储空间，然后逐个元素进行赋值。例如：Bcell(2,3)，可在内存空间中建立一个单元型空变量B，然后可逐个对其每个元素赋值。,值得注意的是：单元型变量的存储并不是以指针的方式来存储的，因此改变元素的值是不会影响到原来所引用的变量值的。单元型变量还可以嵌套，即单元型变量的元素也可以是单元型变量。,(2)单元型变量的相关函数,4.结构型变量,结构型变量是另外一种可以将不同的数据类型组合在一起的特殊数据类型。与单元型数据类型相同它也不是以指针方式传递数据的。不同的是其作用相当于数据库中的记录，可存储一系列相关数据。同一个数据字段（Field）必须具有相同的数据类型，而单元型数据每个元素彼此可以不同。,(1)结构型变量的定义,直接赋值法：结构型变量的使用必须指出结构的属性名，并以操作符“.”来连接结构变量名与属性名。对该属性直接赋值，MATLAB会自动生成该结构变量。如A.b1，B(2,3).a3等。结构型数组的不同元素类型可不同。由struct函数预先分配存储空间方法是：结构型变量struct(元素名1，元素值1，元素名2，元素值2，)例如：c=struct(c1,1,c2,b,c3,abcd),值得注意的是：结构型变量的存储也不是以指针方式存储的。因此改变元素的值就不会影响到所引用变量的值。结构型变量也可以嵌套使用，即结构型变量的元素也可以是结构型变量。,(2)结构型变量的相关函数,二、命令行的基本操作,创建矩阵的方法直接输入法规则：矩阵元素必须用括住矩阵元素必须用逗号或空格分隔在内矩阵的行与行之间必须用分号分隔,矩阵元素可以是任何matlab表达式，可以是实数，也可以是复数，复数可用特殊数i，j输入a=123;456x=2pi/2;sqrt(3)3+5i,矩阵元素,符号的作用,逗号和分号的作用逗号和分号可作为指令间的分隔符，matlab允许多条语句在同一行出现。分号如果出现在指令后，屏幕上将不显示结果。,注意：只要是赋过值的变量，不管是否在屏幕上显示过，都存储在工作空间中，以后可随时显示或调用。变量名尽可能不要重复，否则会覆盖。当一个指令或矩阵太长时，可用续行,冒号的作用用于生成等间隔的向量，默认间隔为1。用于选出矩阵指定行、列及元素。循环语句,2.用matlab函数创建矩阵,空阵matlab允许输入空阵，当一项操作无结果时，返回空阵。rand服从均匀分布的随机矩阵eye单位矩阵(对角元素为1,其他为0)zeros全部元素都为0的矩阵ones全部元素都为1的矩阵,还有伴随矩阵、稀疏矩阵、魔方矩阵、对角矩阵、范德蒙等矩阵的创建，就不一一介绍了。注意：matlab严格区分大小写字母，因此a与A是两个不同的变量。matlab函数名必须小写。,3.矩阵的修改,直接修改可用键找到所要修改的矩阵，用键移动到要修改的矩阵元素上即可修改。指令修改可以用A(,)=来修改。,例如a=120;305;789a=120305789a(3,3)=0a=120305780,还可以用函数subs修改，matlab6.0之后可用find函数配合修改。,例如:单输入参数的情形：设a=980C1=3已经存在于工作空间中.语句y=dsolve(Dy=-a*y)会产生y=exp(-a*t)*C1这时用语句subs(y)会产生ans=3*exp(-980*t)如果是三个参数形如:subs(a+b,a,4)则会返回4+b.,还可以进行多重替换，如:symsabsubs(cos(a)+sin(b),a,b,sym(alpha),2)或者subs(cos(a)+sin(b),a,b,sym(alpha),2)返回cos(alpha)+sin(2)标量表达式的替换:subs(exp(a*t),a,-magic(2)返回exp(-t),exp(-3*t)exp(-4*t),exp(-2*t)多重标量表达式的替换:subs(x*y,x,y,01;-10,1-1;-21)返回0,-12,0,把matlab工作空间中一些有用的数据长久保存下来的方法是生成mat数据文件。save将工作空间中所有的变量存到matlab.mat文件中。,三、数据的保存与获取,默认文件名,savedata将工作空间中所有的变量存到data.mat文件中。savedataab将工作空间中a和b变量存到data.mat文件中。下次运行matlab时即可用load指令调用已生成的mat文件。,loadloaddataloaddataabmat文件是标准的二进制文件，也可用ASCII码形式保存。,即可恢复保存过的所有变量,SAVE.-ASCII以8位ASCII存贮。其他选项如：-DOUBLE以16位ASCII存贮。-TABS以TAB制表符分格ASCII数-APPEND以追加方式存贮变量-MAT以MAT格式存贮而不管扩展名-STRUCT结构性变量存贮专用另外还有V4，V6，REGEXP等选项,矩阵加、减（,）运算规则：相加、减的两矩阵必须有相同的行和列两矩阵对应元素相加减。允许参与运算的两矩阵之一是标量。标量与矩阵的所有元素分别进行加减操作。,四、矩阵运算,2.矩阵乘（）运算规则：A矩阵的列数必须等于B矩阵的行数标量可与任何矩阵相乘。a=123;456;780;b=1;2;3;c=a*bc=143223,d=-1;0;2;f=pi*df=-3.141606.2832矩阵除的运算在线性代数中没有，有矩阵逆的运算，在matlab中有两种矩阵除运算,apa自乘p次幂,方阵,1的整数,3.矩阵乘方an,ap,pa,对于p的其它值,计算将涉及特征值和特征向量，如果p是矩阵，a是标量ap使用特征值和特征向量自乘到p次幂；如a,p都是矩阵，ap则无意义。,a=1,2,3;4,5,6;7,8,9;a2ans=303642668196102126150,当一个方阵有复数特征值或负实特征值时，非整数幂是复数阵。,a0.5ans=0.4498+0.7623i0.5526+0.2068i0.6555-0.3487i1.0185+0.0842i1.2515+0.0228i1.4844-0.0385i1.5873-0.5940i1.9503-0.1611i2.3134+0.2717i,inv矩阵求逆det行列式的值eig矩阵的特征值diag对角矩阵矩阵转置sqrt矩阵开方,4.矩阵的其它运算,5.矩阵的一些特殊操作,矩阵的变维a=1:12;b=reshape(a,3,4)c=zeros(3,4);c(:)=a(:)矩阵的变向rot90:旋转;fliplr:左右翻;flipud:上下翻矩阵的抽取diag:抽取主对角线;tril:抽取主下三角;triu:抽取主上三角矩阵的扩展,关系运算,逻辑运算,关系函数和逻辑函数,数组运算指元素对元素的算术运算，与通常意义上的由符号表示的线性代数矩阵运算不同数组加减(+,-)a+ba-b,5.矩阵的数组运算,对应元素相加减（与矩阵加减等效）,2.数组乘除(，./，.）aba，b两数组必须有相同的行和列两数组相应元素相乘。a=123;456;789;b=246;135;7910;a.*bans=281841530497290,a=123;456;789;b=246;135;7910;a*bans=253746558510985133172,a./b=b.aa.b=b./aa./b=b.a都是a的元素被b的对应元素除a.b=b./a都是a的元素被b的对应元素除例:a=123;b=456;c1=a.b;c2=b./ac1=4.00002.50002.0000c2=4.00002.50002.0000,给出a,b对应元素间的商.,3.数组乘方(.)元素对元素的幂例:a=123;b=456;z=a.2z=1.004.009.00z=a.bz=1.0032.00729.00,常见的基本数学函数,三角函数,指数与对数函数,复数函数,取整函数,矩阵函数,特殊函数,matlab语言把多项式表达成一个行向量，该向量中的元素是按多项式降幂排列的。f(x)=anxn+an-1xn-1+a0可用行向量p=anan-1a1a0表示poly产生特征多项式系数向量特征多项式一定是n+1维的特征多项式第一个元素一定是1,五、多项式运算,例:a=123;456;780;p=poly(a)p=1.00-6.00-72.00-27.00p是多项式p(x)=x3-6x2-72x-27的matlab描述方法，我们可用：p1=poly2str(p,x)函数文件，显示数学多项式的形式p1=x3-6x2-72x-27,2.roots求多项式的根,a=123;456;780;p=poly(a)p=1.00-6.00-72.00-27.00r=roots(p)r=12.12-5.73显然r是矩阵a的特征值-0.39,当然我们可用poly令其返回多项式形式p2=poly(r)p2=1.00-6.00-72.00-27.00matlab规定多项式系数向量用行向量表示，一组根用列向量表示。,3.conv多项式乘运算,例:a(x)=x2+2x+3;b(x)=4x2+5x+6;c=(x2+2x+3)(4x2+5x+6)a=123;b=456;c=conv(a,b)=conv(123,456)c=4.0013.0028.0027.0018.00p=poly2str(c,x)p=4x4+13x3+28x2+27x+18,4.deconv多项式除运算,a=123;c=4.0013.0028.0027.0018.00d=deconv(c,a)d=4.005.006.00,5.多项式微积分,matlab提供了polyder函数多项式的微分。命令格式：polyder(p):求p的微分p=polyder(a,b):求多项式乘积a*b的微分p,q=polyder(a,b):求多项式商a/b的微分例：a=12345;poly2str(a,x)ans=x4+2x3+3x2+4x+5b=polyder(a)b=4664poly2str(b,x)ans=4x3+6x2+6x+4,polyint求多项式函数的不定积分：命令格式：p=polyint(a):求a的不定积分，常数项为0p=polyint(a,k):求a的不定积分，常数项为k例：a=12345;poly2str(a,x)ans=x4+2x3+3x2+4x+5b=polyint(a)b=0.20000.50001.00002.00005.00000poly2str(b,x)ans=0.2x5+0.5x4+x3+2x2+5x,六、代数方程组求解,matlab中有两种除运算左除和右除。对于方程ax=b，a为anm矩阵，有三种情况：当n=m时，此方程成为“恰定”方程当nm时，此方程成为“超定”方程当nm时，此方程成为“欠定”方程matlab定义的除运算可以很方便地解上述三种方程,1.恰定方程组的解,方程ax=b(a为非奇异)x=a-1b矩阵逆两种解:x=inv(a)b采用求逆运算解方程x=ab采用左除运算解方程,方程ax=ba=12;23;b=8;13;x=inv(a)*bx=abx=x=2.002.003.003.00,=,ax=b,例:x1+2x2=82x1+3x2=13,2.超定方程组的解,方程ax=b,mn时此时不存在唯一解。方程解(aa)x=abx=(aa)-1ab求逆法x=abmatlab用最小二乘法找一个准确地基本解。,例:x1+2x2=12x1+3x2=23x1+4x2=3a=12;23;34;b=1;2;3;解1x=ab解2x=inv(aa)abx=x=1.001.0000.00,=,ax=b,3.欠定方程组的解,当方程数少于未知量个数时,即不定情况,有无穷多个解存在。matlab可求出两个解：用除法求的解x是具有最多零元素的解是具有最小长度或范数的解，这个解是基于伪逆pinv求得的。,x1+2x2+3x3=12x1+3x2+4x3=2a=123;234;b=1;2;x=abx=pinv(a)bx=x=1.000.8300.330-0.17,ax=b,七、微分方程求解,微分方程求解的仿真算法有多种，常用的有Euler（欧拉法）、RungeKutta(龙格-库塔法。Euler法称一步法，用于一阶微分方程,当给定仿真步长时：所以yn+1=yn+hf(xn,yn)n=0,1,2y(x0)=y0,RungeKutta法龙格-库塔法：实际上取两点斜率的平均斜率来计算的，其精度高于欧拉算法。龙格-库塔法：ode23ode45,k1=hf(xn,yn)k2=hf(xn+h,yn+k),例：x+(x2-1)x+x=0为方便令x1=x，x2=x分别对x1,x2求一阶导数，整理后写成一阶微分方程组形式x1=x2x2=x2(1-x12)-x1建立m文件解微分方程,建立m文件functionxdot=wf(t,x)xdot=zeros(2,1)xdot(1)=x(2)xdot(2)=x(2)*(1-x(1)2)-x(1)给定区间、初始值;求解微分方程t0=0;tf=20;x0=00.25;t,x=ode23(wf,t0,tf,x0)plot(t,x),figure(2),plot(x(:,1),x(:,