《matlab用法详解.ppt》由会员分享，可在线阅读，更多相关《matlab用法详解.ppt（46页珍藏版）》请在三一办公上搜索。
1、第二讲 MATLAB的数值计算,matlab 具有出色的数值计算能力，占据世界上数值计算软件的主导地位,数值运算的功能,创建矩阵矩阵运算多项式运算线性方程组微分方程的数值解,一、命令行的基本操作,创建矩阵的方法直接输入法规则：矩阵元素必须用 括住 矩阵元素必须用逗号或空格分隔 在 内矩阵的行与行之间必须 用分号分隔,矩阵元素可以是任何matlab表达式，可以是实数，也可以是复数，复数可用特殊函数I，j 输入 a=1 2 3;4 5 6 x=2 pi/2;sqrt(3)3+5i,矩阵元素,符号的作用,逗号和分号的作用 逗号和分号可作为指令间的分隔符，matlab允许多条语句在同一行出现。分号如果
2、出现在指令后，屏幕上将不显示结果。,注意：只要是赋过值的变量，不管是否在屏幕上显示过，都存储在工作空间中，以后可随时显示或调用。变量名尽可能不要重复，否则会覆盖。当一个指令或矩阵太长时，可用续行,冒号的作用 用于生成等间隔的向量，默认间隔为1。用于选出矩阵指定行、列及元素。循环语句,2.用matlab函数创建矩阵,空阵 matlab允许输入空阵，当一项操作无结果时，返回空阵。rand 随机矩阵eye 单位矩阵zeros 全部元素都为0的矩阵ones 全部元素都为1的矩阵,还有伴随矩阵、稀疏矩阵、魔方矩阵、对角矩阵等矩阵的创建，就不一一介绍了。注意：matlab严格区分大小写字母，因此a与A是两
3、个不同的变量。matlab函数名必须小写。,3.矩阵的修改,直接修改 可用键找到所要修改的矩阵，用键移动到要修改的矩阵元素上即可修改。指令修改 可以用A(,)=来修改。,例如a=1 2 0;3 0 5;7 8 9a=1 2 0 3 0 5 7 8 9a(3,3)=0a=1 2 0 3 0 5 7 8 0,把matlab工作空间中一些有用的数据长久保存下来的方法是生成mat数据文件。save 将工作空间中所有的变量存到matlab.mat文件中。,二、数据的保存与获取,默认文件名,save data将工作空间中所有的变量存到data.mat文件中。save data a b 将工作空间中a和b变
4、量存到data.mat文件中。下次运行matlab时即可用load指令调用已生成的mat文件。,load load data load data a b mat文件是标准的二进制文件，还可以ASCII码形式保存。,即可恢复保存过的所有变量,矩阵加、减（,）运算规则：相加、减的两矩阵必须有相同的行和列两矩阵对应元素相加减。允许参与运算的两矩阵之一是标量。标量与矩阵的所有元素分别进行加减操作。,三、矩阵运算,2.矩阵乘（）运算规则：A矩阵的列数必须等于B矩阵的行数标量可与任何矩阵相乘。a=1 2 3;4 5 6;7 8 0;b=1;2;3;c=a*bc=14 32 23,d=-1;0;2;f=pi
5、*df=-3.1416 0 6.2832 矩阵除的运算在线性代数中没有，有矩阵逆的运算，在matlab中有两种矩阵除运算x=ab 采用左除运算解方程,x=inv(a)b 采用求逆运算解方程,inv 矩阵求逆det 行列式的值eig 矩阵的特征值diag 对角矩阵；diag（）矩阵转置;注意输入时在英文状态下sqrt 矩阵开方,3.矩阵的其它运算,关系运算,数组运算指元素对元素的算术运算，与通常意义上的由符号表示的线性代数矩阵运算不同 数组加减(.+,.-)a.+b a.-b,4.矩阵的数组运算,对应元素相加减（与矩阵加减等效）,2.数组乘除(，./，.）ab a，b两数组必须有相同的行 和列两
6、数组相应元素相乘。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,a=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,a./b=b.aa.b=b./aa./b=b.a 都是a的元素被b的对应元 素除a.b=b./a 都是a的元素被b的对应元 素除例: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.00
7、00,给出a,b对应元素间的商.,3.数组乘方(.)元素对元素的幂例:a=1 2 3;b=4 5 6;z=a.2z=1.00 4.00 9.00z=a.bz=1.00 32.00 729.00,四、代数方程组求解,matlab中有两种除运算左除和右除。对于方程ax+b，a 为anm矩阵，有三种情况：当n=m时，此方程成为“恰定”方程 当nm时，此方程成为“超定”方程 当nm时，此方程成为“欠定”方程 matlab定义的除运算可以很方便地解上述三种方程,1.恰定方程组的解,方程ax+b(a为非奇异)x=a-1 b 矩阵逆两种解:x=inv(a)b 采用求逆运算解方程 x=ab 采用左除运算解方程
8、,方程ax=ba=1 2;2 3;b=8;13;x=inv(a)*b x=ab x=x=2.00 2.00 3.00 3.00,=,a x=b,例:x1+2x2=8 2x1+3x2=13,练习,求解方程组,matlab语言把多项式表达成一个行向量，该向量中的元素是按多项式降幂排列的。f(x)=anxn+an-1xn-1+loa0 可用行向量 p=an an-1 a1+a0表示poly 产生特征多项式系数向量特征多项式一定是n+1维的特征多项式第一个元素一定是1a 个准确地基本解。,五、多项式运算,例:a=1 2 3;4 5 6;7 8 0;p=poly(a)p=1.00-6.00-72.00-
9、27.00 p是多项式p(x)=x3-6x2-72x-27的matlab描述方法，我们可用：p1=poly2str(p,x)函数文件，显示数学多项式的形式p1=x3-6 x2-72 x 27,2.roots 求多项式的根,a=1 2 3;4 5 6;7 8 0;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,con
10、vs多项式乘运算,例:a(x)=x2+2x+3;b(x)=4x2+5x+6;c=(x2+2x+3)(4x2+5x+6)a=1 2 3;b=4 5 6;c=conv(a,b)=conv(1 2 3,4 5 6)c=4.00 13.00 28.00 27.00 18.00p=poly2str(c,x)p=4 x4+13 x3+28 x2+27 x+18,4.deconv多项式除运算,a=1 2 3;c=4.00 13.00 28.00 27.00 18.00d=deconv(c,a)d=4.00 5.00 6.00,5.多项式微分,matlab提供了polyder函数多项式的微分。命令格式：pol
11、yder(p):求p的微分polyder(a,b):求多项式a,b乘积的微分p,q=polyder(a,b):求多项式a,b商的微分例：a=1 2 3 4 5;poly2str(a,x)ans=x4+2 x3+3 x2+4 x+5b=polyder(a)b=4 6 6 4poly2str(b,x)ans=4 x3+6 x2+6 x+4,六、微分方程求解,微分方程求解的仿真算法有多种，常用的有Euler（欧拉法）、Runge Kutta(龙 格-库塔法。Euler法称一步法，用于一阶微分方程,当给定仿真步长时：所以 yn+1=yn+hf(xn,yn)n=0,1,2 y(x0)=y0,Runge
12、Kutta法龙格-库塔法：实际上取两点斜率的平均 斜率来计算的，其精度高于欧拉算法。龙格-库塔法：ode23 ode45,k1=hf(xn,yn)k2=hf(xn+h,yn+k),例：x+(x2-1)x+x=0为方便令x1=x，x2=x分别对x1,x2求一阶导数，整理后写成一阶微分方程组形式 x1=x2 x2=x2(1-x12)-x1建立m文件解微分方程,建立m文件function xdot=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=0 0.25;t,x=o
13、de23(wf,t0,tf,x0)plot(t,x),figure(2),plot(x(:,1),x(:,2),命令格式:T,Y=ODE23(ODEFUN,TSPAN,Y0)建立m文件function dxdt=wf(t,x)dxdt=x(2);x(2)*(1-x(1)2)-x(1);求解微分方程t,x=ode23(wf,0 30,0 0.25);plot(t,x);figure(2)plot(x(:,1),x(:,2),八、数据分析与插值函数,max 各列最大值 mean 各列平均值sum 各列求和std 各列标准差var 各列方差sort 各列递增排序,九、拟合与插值,1.多项式拟合x0=
14、0:0.1:1；y0=-.447 1.978 3.11 5.25 5.02 4.66 4.01 4.58 3.45 5.35 9.22;p=polyfit(x0,y0,3)p=56.6915-87.1174 40.0070-0.9043xx=0:0.01:1;yy=polyval(p,xx);plot(xx,yy,-b,x0,y0,or),2.插值插值的定义是对某些集合给定的数据点之间函数的估值方法。当不能很快地求出所需中间点的函数时，插值是一个非常有价值的工具。Matlab提供了一维、二维、三次样条等许多插值选择,table1 table2 intep1 interp2 spline 利用已知点确定未知点粗糙 精确集合大的 简化的,插值函数,小 结 本节介绍了matlab语言的数值运算功能，通过学习应该掌握：如何创建矩阵、修改矩阵符号的用法矩阵及数组运算多项式运算线性方程组与微分运算,
