(精品)第二章matlab语言基础.ppt
《(精品)第二章matlab语言基础.ppt》由会员分享,可在线阅读,更多相关《(精品)第二章matlab语言基础.ppt(51页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、CH2、MATLAB语言基础第一节 使用MATLAB的窗口环境一、MATLAB语言的显著特点1、具有强大的矩阵运算能力:Matrix Laboratory(矩阵实验室),使得矩阵运算非常简单。2、是一种演算式语言MATLAB的基本数据单元是既不需要指定维数,也不需要说明数据类型的矩阵(向量和标量为矩阵的特例),而且数学表达式和运算规则与通常的习惯相同。因此MATLAB语言编程简单,使用方便。例 exp2_1.m二、MATLAB命令窗口1、启动MATLAB命令窗口计算机安装好MATLAB之后,双击MATLAB图标,就可以进入命令窗口,此时意味着系统处于准备接受命令的状态,可以在命令窗口中直接输入
2、命令语句。MATLAB语句形式变量表达式;通过等于符号将表达式的值赋予变量。当键入回车键时,该语句被执行。语句执行之后,窗口自动显示出语句执行的结果。如果希望结果不被显示,则只要在语句之后加上一个分号(;)即可。此时尽管结果没有显示,但它依然被赋值并在MATLAB工作空间中分配了内存。2、命令行编辑器(1)方向键和控制键可以编辑修改已输入的命令 :回调上一行命令 :回调下一行命令(2)命令窗口的分页输出more off:不允许分页 more on:允许分页more(n):指定每页输出的行数回车前进一行,空格键显示下一页,q结束当前显示。(3)多行命令()如果命令语句超过一行或者太长希望分行输入
3、,则可以使用多行命令继续输入。S=1-12+13+4+9-4-18;三、变量和数值显示格式1、变量(1)变量的命名:变量的名字必须以字母开头(不能超过19个字符),之后可以是任意字母、数字或下划线;变量名称区分字母的大小写;变量中不能包含有标点符号。(2)一些特殊的变量ans:用于结果的缺省变量名 i、j:虚数单位pi:圆周率 nargin:函数的输入变量个数eps:计算机的最小数 nargout:函数的输出变量个数inf:无穷大 realmin:最小正实数realmax:最大正实数 nan:不定量flops:浮点运算数(3)变量操作在命令窗口中,同时存储着输入的命令和创建的所有变量值,它们可
4、以在任何需要的时候被调用。如要察看变量a的值,只需要在命令窗口中输入变量的名称即可:a2、数值显示格式任何MATLAB的语句的执行结果都可以在屏幕上显示,同时赋值给指定的变量,没有指定变量时,赋值给一个特殊的变量ans,数据的显示格式由format命令控制。format只是影响结果的显示,不影响其计算与存储;MATLAB总是以双字长浮点数(双精度)来执行所有的运算。如果结果为整数,则显示没有小数;如果结果不是整数,则输出形式有:format(short):短格式(5位定点数)99.1253format long:长格式(15位定点数 99.12345678900000format short
5、e:短格式e方式 9.9123e+001format long e:长格式e方式 9.912345678900000e+001format bank:2位十进制 99.12format hex:十六进制格式四、简单的数学运算(例exp2_2.m)1、常用的数学运算符,*(乘),/(左除),(右除),(幂)在运算式中,MATLAB通常不需要考虑空格;多条命令可以放在一行中,它们之间需要用分号隔开;逗号告诉MATLAB显示结果,而分号则禁止结果显示。2、常用数学函数abs,sin,cos,tan,asin,acos,atan,sqrt,exp,imag,real,sign,log,log10,co
6、nj(共扼复数)等 五、MATLAB的工作空间1、MATLAB的工作空间包含了一组可以在命令窗口中调整(调用)的参数who:显示当前工作空间中所有变量的一个简单列表whos:则列出变量的大小、数据格式等详细信息clear:清除工作空间中所有的变量clear 变量名:清除指定的变量2、保存和载入workspace(1)save filename variables将变量列表variables所列出的变量保存到磁盘文件filename中Variables所表示的变量列表中,不能用逗号,各个不同的变量之间只能用空格来分隔。未列出variables时,表示将当前工作空间中所有变量都保持到磁盘文件中。缺
7、省的磁盘文件扩展名为“.mat”,可以使用“-”定义不同的存储格式(ASCII、V4等)(2)load filename variables将以前用save命令保存的变量variables从磁盘文件中调入MATLAB工作空间。用load 命令调入的变量,其名称为用save命令保存时的名称,取值也一样。Variables所表示的变量列表中,不能用逗号,各个不同的变量之间只能用空格来分隔。未列出variables时,表示将磁盘文件中的所有变量都调入工作空间。3、退出工作空间quit 或 exit六、文件管理文件管理的命令,包括列文件名、显示或删除文件、显示或改变当前目录等。(what、dir、ty
8、pe、delete、cd、which)what:显示当前目录下所有与matlab相关的文件及它们的路径。dir:显示当前目录下所有的文件which:显示某个文件的路径cd path:由当前目录进入path目录cd.:返回上一级目录cd:显示当前目录type filename:在命令窗口中显示文件filenamedelete filename:删除文件filename七、使用帮助1、help命令,在命令窗口中显示MATLBA的所有函数都是以逻辑群组方式进行组织的,而MATLAB的目录结构就是以这些群组方式来编排的。help matfun:矩阵函数数值线性代数help general:通用命令he
9、lp graphics:通用图形函数help elfun:基本的数学函数help elmat:基本矩阵和矩阵操作help datafun:数据分析和傅立叶变换函数help ops:操作符和特殊字符help polyfun:多项式和内插函数help lang:语言结构和调试help strfun:字符串函数help control:控制系统工具箱函数2、helpwin:帮助窗口3、helpdesk:帮助桌面,浏览器模式4、lookfor命令:返回包含指定关键词的那些项5、demo:打开示例窗口矩阵生成不但可以使用纯数字(含复数),也可以使用变量(或者说采用一个表达式)。矩阵的元素直接排列在方括号
10、内,行与行之间用分号隔开,每行内的元素使用空格或逗号隔开。大的矩阵可以用分行输入,回车键代表分号。y=2,4,53 6 8y=2 4 5 3 6 8一、矩阵的输入第二节 MATLAB矩阵运算及多项式处理a=1;b=2;c=3;x=5 b c;a*b a+c c/bx=5.000 2.000 3.000 2.000 4.000 1.5001、在命令窗口中输入(1)用线性等间距生成向量矩阵(start:step:end)a=1:2:10a=1 3 5 7 9其中start为起始值,step为步长,end为终止值。当步长为1时可省略step参数;另外step也可以取负数。2、语句生成(2)a=lin
11、space(n1,n2,n)在线性空间上,行矢量的值从n1到n2,数据个数为n,缺省n为100。a=linspace(1,10,10)a=1 2 3 4 5 6 7 8 9 10(4)一些常用的特殊矩阵单位矩阵:eye(m,n);eye(m)零矩阵:zeros(m,n);zeros(m)一矩阵:ones(m,n);ones(m)对角矩阵:对角元素向量 V=a1,a2,an A=diag(V)随机矩阵:rand(m,n)产生一个mn的均匀分别的随机矩阵(3)a=logspace(n1,n2,n)在对数空间上,行矢量的值从10n1到10n2,数据个数为n,缺省n为50。这个指令为建立对数频域轴坐标
12、提供了方便。a=logspace(1,3,3)a=10 100 1000eye(2,3)ans=1 0 0 0 1 0zeros(2,3)ans=0 0 0 0 0 0ones(2,3)ans=1 1 1 1 1 1V=5 7 2;A=diag(V)A=5 0 0 0 7 0 0 0 2eye(2)ans=1 0 0 1zeros(2)ans=0 0 0 0ones(2)ans=1 1 1 1如果已知A为方阵,则V=diag(A)可以提取A的对角元素构成向量V。1、转置:对于实矩阵用()符号或(.)求转置结果是一样的;然而对于含复数的矩阵,则()将同时对复数进行共轭处理,而(.)则只是将其排列
13、形式进行转置。a=1 2 3;4 5 6a=1 4 2 5 3 6a=1 2 3;4 5 6.a=1 4 2 5 3 6b=1+2i 2-7ib=1.0000-2.0000i 2.0000+7.0000ib=1+2i 2-7i.b=1.0000+2.0000i 2.0000-7.0000i二、矩阵的运算+;-;*;和/;.*;.;./;.如:a=1 2;3 4;b=3 5;5 9c=a+b d=a-bc=d=4 7 -2 -3 8 13 -2 -5a*b=13 23;29 51a/b=-0.50 0.50;3.50 1.50ab=-1-1;2 3a3=37 54;81 118a.*b=3 10
14、;15 36a./b=0.33 0.40;0.60 0.44a.b=3.00 2.50;1.67 2.25a.3=1 8;27 642、四则运算与幂运算只有维数相同的矩阵才能进行加减运算。注意只有当两个矩阵中前一个矩阵的列数和后一个矩阵的行数相同时,才可以进行乘法运算。ab运算等效于求a*x=b的解;而a/b等效于求x*b=a的解。只有方阵才可以求幂。点运算是两个维数相同矩阵对应元素之间的运算,在有的教材中也定义为数组运算。3、逆矩阵与行列式计算求逆:inv(A);求行列式:det(A)要求矩阵必须为方阵4、了解矩阵超越函数在MATLAB中exp、sqrt等命令也可以作用到矩阵上,但这种运算是
15、定义在矩阵的单个元素上的,即分别对矩阵的每一个元素进行计算。超越数学函数可以在函数后加上m而成为矩阵的超越函数,例如:expm,sqrtm。矩阵的超越函数要求运算矩阵为方阵。a=1 2 3;4 5 6;2 3 5;b=inv(a)b=-2.3333 0.3333 1.0000 2.6667 0.3333 -2.0000 -0.6667 -0.3333 1.0000det(a)ans=-3A(m,n):提取第m行,第n列元素A(:,n):提取第n列元素A(m,:):提取第m行元素A(m1:m2,n1:n2):提取第m1行到第m2行和第n1列到 第n2列的所有元素(提取子块)。A(:):得到一个长
16、列矢量,该矢量的元素按矩阵的列进行排列。矩阵扩展:如果在原矩阵中一个不存在的地址位置上设定一个数(赋值),则该矩阵会自动扩展行列数,并在该位置上添加这个数,而且在其他没有指定的位置补零。消除子块:如果将矩阵的子块赋值为空矩阵,则相当于消除了相应的矩阵子块。三、矩阵的操作1、矩阵下标MATLAB通过确认矩阵下标,可以对矩阵进行插入子块,提取子块和重排子块的操作。2、矩阵的大小m,n=size(A,x):返回矩阵的行列数m与n,当x=1,则只返回行数m,当x=2,则只返回列数n。length(A)=max(size(A):返回行数或列数的最大值。rank(A):求矩阵的秩a=1 2 3;3 4 5
17、;m,n=size(a)m=2n=33、了解矩阵操作函数:flipud;fliplr;rot90length(a)ans=3max(size(a)ans=3rank(a)ans=2在MATLAB中,多项式使用降幂系数的行向量表示,如:多项式p=poly(r)p=1 -12 -0 25 116四、多项式处理(1)多项式的建立与表示方法r=roots(p)r=11.7473 2.7028 -1.2251+1.4672i -1.2251-1.4672i表示为:p=1-12 0 25 116,使用函数roots可以求出多项式等于0的根,根用列向量表示。若已知多项式等于0的根,函数poly可以求出相应多
18、项式。(2)多项式的运算相乘conva=1 2 3;b=1 2 c=conv(a,b)=1 4 7 6conv指令可以嵌套使用,如conv(conv(a,b),c)相除deconvq,r=deconv(c,b)q=1 2 3 商多项式r=0 0 0 余多项式求多项式的微分多项式polyderpolyder(a)=2 2求多项式函数值polyval(p,n):将值n代入多项式求解。polyval(a,2)=11(3)*多项式的拟合多项式拟合又称为曲线拟合,其目的就是在众多的样本点中进行拟合,找出满足样本点分布的多项式。这在分析实验数据,将实验数据做解析描述时非常有用。命令格式:p=polyfit
19、(x,y,n),其中x和y为样本点向量,n为所求多项式的阶数,p为求出的多项式。例exp2_15.m(4)*多项式插值多项式插值是指根据给定的有限个样本点,产生另外的估计点以达到数据更为平滑的效果。该技巧在信号处理与图像处理上应用广泛。所用指令有一维的interp1、二维的interp2、三维的interp3。这些指令分别有不同的方法(method),设计者可以根据需要选择适当的方法,以满足系统属性的要求。Help polyfun可以得到更详细的内容。y=interp1(xs,ys,x,method)在有限样本点向量xs与ys中,插值产生向量x和y,所用方法定义在method中,有4种选择:n
20、earest:执行速度最快,输出结果为直角转折linear:默认值,在样本点上斜率变化很大spline:最花时间,但输出结果也最平滑cubic:最占内存,输出结果与spline差不多例exp2_16.m五*、MATLAB数据处理1、矩阵分解(1)奇异值分解U,S,V=svd(A)例:a=9 8 6 8可以验证:u*u=Iv*v=Iu*s*v=a求矩阵A的奇异值及分解矩阵,满足U*S*V=A,其中U、V矩阵为正交矩阵(U*U=I),S矩阵为对角矩阵,它的对角元素即A矩阵的奇异值。u,s,v=svd(a)u=0.7705 -0.6375 0.6375 0.7705s=15.5765 0 0 1.5
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 精品 第二 matlab 语言 基础
限制150内