实验1熟悉matlab环境和基本操作

1、实验目的：1.熟悉Matlab环境，掌握Matlab的主要窗口及功能；2. 学会Matlab的帮助使用；3. 掌握向量、矩阵的定义、生成方法和基本运算；4. 掌握Matlab的基本符号运算；5. 掌握Matlab中的二维图形的绘制和控制。实验内容：1. 启动Matlab,说明主窗口、命令窗口、当前目录窗口、工作空间窗口、历史窗口、图形窗口、M文件编辑器窗口的功能。2. 实例操作Matlab的帮助使用。3. 实例操作向量、矩阵的定义、生成方法和基本运算。4. 实例操作Matlab的基本符号运算。5. 实例操作Matlab中的二维图形绘制和控制。实验仪器与软件：1.CPS频在2GH以上，内存在51

2、2Mb以上的PC;2. Matlab7及以上版本。实验讲评：实验成绩:评阅教师：年月曰、Matlab环境及主要窗口的功能运行Matlab安装目录下的matlab.exe文件可启动Matlab环境，其默认布局如下图其中，1. 主窗口的功能是：主窗口不能进行任何计算任务操作，只用来进行一些整体的环境参数设置，它主要对6个下拉菜单的各项和10个按钮逐一解脱。命令窗口的功能是：对MATLAB搜索路径中的每一个M文件的注释区的第一行进行扫描，一旦发现此行中含有所查询的字符申，则将该函数名及第一行注释全部显小在屏幕上。2. 历史窗口的功能是：历史窗口显示命令窗口中的所有执行过的命令，一方面可以查看曾经执行

3、过的命令，另一方面也可以重复利用原来输入的命令行，可以从命令窗口中直接通过双击某个命令行来执行该命令，3. 当前目录窗口的功能是：显示当前目录下所有文件的文件名、文件类型、和最后修改的时间，同时还提供搜索功能，在该窗口下，可以改变当前目录4. M文件编辑器窗口的功能是：MATLAN提供了一个内置既有编辑和调试功能好的程序编辑器；编辑器窗口也有菜单和工具栏，是编辑和调试程序非常方便6.发行说明窗口功能是：该窗口显示MATLAB总包和已安装的工具箱的帮助、演示、GUI工具和产品主贞等4个内容。工作空间窗口的功能是：该窗口显示所有目前内存中MATLAB变量的变量名、数字结构、字节数以及类型，不同的变

4、量类型分别对应不同的变量名图标7. 图形窗口的功能是：利用图形窗口和工具栏中的选项，可以对图形进行线性、颜色、标记三维视图、光照和坐标轴等的设置8. GUI(GraphacalUserInterface)®口功能是：二、Matlab的帮助使用Matlab提供的联机帮助系统使用户在没有任何资料的情况下就能掌握它的使用和基本操作，作为Matlab的用户应熟练掌握其联机帮助系统的使用，下面是Matlab联机帮助系统的使用方法。FileEdit如ie网QqFarvaritesQesktopWindowHelpContentsAdgSe«rchD«mcsroReleaseN

5、otes以Installation耕MATLABExcelLinkMATLABBuilderYorCOMMATLABBuilderforExcel:(MATLABCompiler&MATLABDistributedComputin金MATLABReportGeneratorMATS日WehServer BioinfbrmaticsToolbox.，CommunicationsToolbox ControlSystemToolbox；CiurveFittingToolbox3DataAcquisitionToolbox偻DatabaseToolbox DatsfeedToolboxDis

6、tributedCarmputingToolboxBeginHereRelease14withServicePack3BeginHereTitie:Relea.5«14.ReleaseNotesHighlightsnewfeatures,installalionnotes,bugfiices,ardcompatibilityissuesIfYouAreUpgradingfromaPreviousRelease.IfYouAreUsingMATLABforthef.4s.m通过命令窗口中直接输入help命令将会显示当前帮助系统只能怪包含的所有项目help三、向量的定义、生成和基本运算1

7、:向量的生成a：逐个元素直接输入，向量元素需要用“”括起来，元素之间可以用空格、逗号或分号分隔。用空格和逗号分隔生成的行向量用分号分隔生成列向量例如：h=345678f=3;4;5;6;7;8b:利用冒号表达式创建通过设定“步长（step）”生成一维行向量，通过格式为：x=x0：step:xn。x0表达向量的首元素值，xn表示尾元素数值限，step表示从第二个元素开始，每一个元素与前一个元素的差值。step=1时，可以省略此项的输入，直接写成x=xo:xno例：y=0:10:100c：定数线性采样生成设定总点数n下，均匀采样生成一维行向量。通用格式为x=linspace(a,b,n)。a,b分

8、别是生成向量的第一个和最后一个元素，n是采样总点数。该指令生成的数组相当于由a:(a-b)/(n-1):b生成的数组。缺省n时，生成100维的行向量。clear%清除工空间中的所有变量x=linspace(6,66,8)y=6:60/7:66z=linspace(6,66)d：定数对数采样生成向量设定总点数n下，经“常用对数”均匀采样生成一维行向量。通用格式为x=logspace(a,b,n)。生成数组的第一个元素值为10a,最后一个元素值为10b,n为采样总点数，缺省时，生成50维的行向量。例如：clear%清除工作空间的所有变量x=logspace(1,8,8)y=1:7/7:8xx=10

9、.Ayz=logspace(1,8)2:向虽元素的引用格式为：向量名(下标范围或元素所满足的条件)。例:clearrand('state',0)%把均匀分布伪随机发生器置为初始状态x=rand(1,8)%产生(1X8)的均匀分布随机数组x(7)%引用数组x的第7个元素y=x(125)%引用数组x的第一、一、五个兀素z=x(1:3)%引用数组x的前二个兀素w=x(3:end)%引用数组x的从第三个元素以后的元素v=x(3:-1:1)%由数组x的前3个元素倒排构成的了数组u=x(find(x>0.5)i%数组x中大于0.5的兀素构成的子数组t=x(12344321)%重复引用

10、数组 3:向量与标量、向量与向量的运算四则运算符号有(+*/)a:标量a与向量x进行四则运算是a分别与x中的每个元素进行四则运算并生一个与x等长的向量。例如clearx=4567891011y=3*x+3b:等长的两个向量才能进行四则运算，向量x与y进行四则运算是这两个向量的对应元素分别进行四则运算并生成一个与它们等长的向量。例如clearx=123456y=x*2z=x+yw=x.*yn=x./yd=x.y 幕运算(.八)a:向量x与标量a的籍运算是对x的每一个元素施行籍运算，例如clearx=123456y=x.A3z=3.Ayb:向量x与向量y的籍运算是元素对元素的籍运算。例如：clea

11、rx=123456y=x*2z=x.Ayb=y.Ax 指数运算、对数运算与开方运算等在MATLAB,数组的运算实质上是数组内部每个元素的运算，因此，数组的指数运算、对数运算与开方运算等与标量运算完全一样，运算函数分别为“exp”、“log”、“sqrt”等。例如：clearx=35791113y=exp(x)z=log(x)t=sqrt(x)四、矩阵的定义、生成和基本运算1. 矩阵的创建a:逐个元素直接输入把矩阵元素需用门”括起来，同行元素之间用空格或逗号分隔，行与行之间用分号或回车符分隔矩阵元素可为运算表达式，无任何元素的矩阵称为空矩阵。例如x=123;456;789y=1,2,3;4,5,

12、6;7,8,9k=sin(pi/7),cos(pi/4)T=编写M文件创建大矩阵对于大型矩阵，可通过编写脚本式M文件，然后运行该文件来创建。例如：编写一名为Example10.m的M文件,内容如下。%Example10.m%编写一M文件创建矩阵的示例文件。emn=4563442245636;0976658645;2950512436;143854259178;45954145245233ans=1/2*kA2-1/2*kans=1/2柘八2恢八4+1/2柘恢八2ans=1/6*piA2 通过函数创建特殊矩阵%Example.m咻写一M文件，通过函数创建特殊矩阵的示例文件。%由函数zeros创建

13、全0矩阵。N=4;M=3;A=12345;23561;44425;B1=zeros(M,N)%生成MXN阶全0阵。C1=zeros(size(A)%生成与A同阶的全0阵。A2=ones(N)%生成NXN阶全1阵。B2=ones(M,N)%生成MXN阶全1阵。C2=ones(size(A)%生成与A同阶的全1阵。%曲函数eye创建单位矩阵。A2=eye(N)%生成NXN阶单位矩阵C2=eye(size(A)%生成与A同阶单位矩阵。%曲函数rand或randn创建随机矩阵。A3=rand(N)%生成NXN阶均匀分布的随机阵，元素值在(0.0,1.0)区间内B3=rand(M,N)%生成MXN阶均匀

14、分布的随机阵。C3=rand(size(A)%生成与A同阶阶均匀分布的随机阵。H=hilb(N)%生成NXN阶Hilbert矩阵。 矩阵元素的引用相对位置引用格式：变虽名(行标，列标)绝对位置引用格式：变虽名(绝对位置索引)clearrand('state',0)A=rand(5,3)A(4)%引用距阵A的第四个元素A(2,3)%引用矩阵A的第二行第三列元素2. 矩阵元素的抽取 抽取行clearrand('state',0)A=rand(5,6)抽取矩阵A的第四行抽取矩阵A的第二行和第三行抽取矩阵A的第三行和第三行赋值给B抽取矩阵A的第三行至最后一行赋值给BA(

15、4,:)%A(23,:)%B=A(33,:)%C=A(3:end,:)% 抽取列clearrand('state',0)A=rand(5,6)A(:,3)%抽取矩阵A的第三列A(:,13)%抽取矩阵A的第一列和第三列B=A(:,31)%抽取矩阵A的第三列和第一列赋值给B 抽取块clearrand('state',0)A=rand(3,4)B=A(12,23)%抽取矩阵A的第一、二行与第二、三列交叉的元素赋值给B 抽取矩阵对角线上的元素clearrand('state',0)A=rand(6)%产生(6X6)的均匀分布随机数组V=diag(A)%抽

16、取矩阵A的主对角线上的元素赋值给向量VD=diag(V)%以向量V为对角线元素生成对角矩阵D1=diag(V,2)D2=diag(V,-2)U=diag(A,1)%抽取矩阵A的主对角线上方第一条对角线的元素赋值给向量UL=diag(A,-1)%抽取矩阵A的主对角线下方第一条对角线的元素赋值给向量L 抽取矩阵上三角部分和下三角部分clearrand('state',0)A=rand(5)%U=triu(A,1)%产生(5X5)的均匀分布随机数组从矩阵A的主对角线上方第一条对角线开始抽取A的上三角部分U=triu(A,-1)%L1=tril(A,1)从矩阵A的主对角线下方第一条对角

17、线开始抽取A的上三角部分L2=tril(A,-1)3. 矩阵的基本数学运算矩阵的四则运算(+-*/)与线性代数理论一致，其中,AB=inv(A)*B=AA-1*BclearA=3230;21-46B=-13;41;60;79C=A+B'D=A*BE=B/DF=DA矩阵与常数间的运算(+-*八A)同线性代数理论一致，需注注的是，当进行数除时，常数通常只能做除数。clearA=1230;21-46C=A+2D=A*2E=A/2F=2AG=A(12,12)八2矩阵的数组运算(.+.-.*./.A)是指同维数组间对应元素之间的加、减、乘、除和窑运算，其中“.+"和“-”分别与“+&q

18、uot;和“-"相同，所以，“.+"和“.-"一般不用。clearA=7280;51-46B=A+2C=A.*BD=A./BE=B.AF=A.A2矩阵的基本初等运算clearA=1230;21-46;-1341;6079A(2,:)=A(2,:)*2%2乘A的第二行A(1,:)=A(1,:)+A(2,:)%2乘A的第二行，加到A的第一行 A(23,:)=A(32,:)%交换A的第二行和第三行矩阵的逆运算clear A=1230;21-46;-1341;6079B=inv(A)矩阵的行列式运算clear A=1230;21-46;-1341;6079B=det(A)

19、矩阵的指数运算clearA=1230;21-46;-1341;6079 B=expm(A)矩阵的对数运算clearA=4232;2146;1341;6479B=expm(A)C=logm(B)D1=logm(A) D2=log(A)矩阵的开方运算clearA=9232;2146;1341;6479B=AA2C=sqrtm(B)B1=sqrtm(A)B2=sqrt(A)4. 矩阵的一些特殊操作 变维方法：“：”和函数reshape”。reshape(A,M,N)%将已知矩阵变维成M<N阶矩阵reshape(A,M,N,p,)%将已知矩阵变维成将NXPX阶矩阵cleara=1:12;A=re

20、shape(a,3,4)c=zeros(2,6);c(:)=a(:) 矩阵的变向。cleara=1:12;A=reshape(a,3,4)k=3;dim=1;A1=rot90(A);A2=rot90(A,k);A3=fliplr(A);A4=flipud(A);A5=flipdim(A,dim);矩阵的扩展与收缩。cleara=1:12;A=reshape(a,3,4)B=eye(3,2)C=ones(2,6)D=AB;C%利用小矩阵的组合来生成大矩阵D(6:10,9:10)=4%利用对矩阵标识块的赋值命令生成大矩阵D(:,3:end)=%将矩阵标识块置为空以收缩矩阵五、多项式的定义、生成和基

21、本运算1.多项式的表小对于多项式an1X翥P/-n.n1n(x)a°xax用行向虽Pa0,ai,ani,an表示，把多项式问题转化为向虽问题。 多项式的创建直接输入系数向虽由于在MATLAB中的多项式是以向虽形式储存的，因此，直接输入多项式对应的向虽，MATLAB会自动将向虽元素按降幕顺序分配给各项系数值，向虽可以为行向虽，也可以是列向虽。例如：输入多项式：x35x26x33p=1-56-33;poly2sym(p)%poly2sym通过特征多项式创建也就是从矩阵求其特征多项式获得。例如：A=123;456;789p=poly(A)poly2sym(p) 由多项式的根创建多项式roo

22、t=7-3+6i-3-4i;p=poly(root)poly2sym(p)2. 多项式运算 求多项式的值一般调用函数polyval进行计算，例如：p=11155125;b=62;0-1;polyval(p,b) 求多项式的根求多项式的根可以有两种方法，一种是直接调用函数roots求解多项式的所有根；另一种是通过建立多项式的伴随矩阵再求其特征值的方法得到多项式的所有根。两种方法求得的根是相等的。例如：p=2-56-19;roots(p)P=compan(p)eig(p)3. 多项式的乘除法运算多项式的乘法由函数conv来实现，除法则由函数deconv来实现，例如：p=2-56-19;poly2s

23、ym(p)d=3-90-18;poly2sym(d)pd=conv(p,d)poly2sym(pd)p1=deconv(pd,d)多项式的微分多项式的微分由函数polyder来实现，例如：p=2-56-19;poly2sym(p)Dp=polyder(p)poly2sym(Dp)多项式拟合多项式拟合的实现，一面可以由矩阵的除法求解超定方程来进行；另一方面可调用函数polyfit来实现，调用方法如下：p,s=polyfit(X,Y,n)其中，X、Y为拟合数据，n为拟合多项式的次，p为拟合多项式的系数向虽，s为拟合多项式系数向虽的结构信息，例如：x=0:pi/20:pi/2;y=sin(x);p=

24、polyfit(x,y,5)x1=0:pi/30:pi*2;y1=sin(x1);y2=polyval(p,x1);plot(x1,y1,'b-',x1,y2,'r+')legend('原曲线',拟合曲线')axis(07-1.29)五、Matlab的基本符号运算1：符号表达式的生成用引号来生成符号表达式，例如仁'exp(x)'仁'a*xA2+b*x+c=0'f='Dy-y=x'用sym来生成符号表达式，例如f=sym('exp(x)');仁sym('a*xA2+b*

25、x+c=0')用函数syms来生成符号函数，例如symsyu;p=exp(y/u)2符号表达式的运算提取分子、分母，例如f=sym('a*xA2/(b-x)')n,d=numden(f)3符号表达式的基本运算4符号表达式的高级运算a:符号表达式的复合函数运算通过compose来实现；例如symsxyt;f=1/xA3;g=tan(y);compose(g,f)compose(g,f,t)b:符号表达式的反函数运算通过函数finverse来实现；例如f=sym(1/sin(x);g=finverse(f)c:符号表达式的符号和运算通过函数symsum来实现；例如k=sym

26、('k');symsum(k)symsum(k,0,n-1)symsum(1/k人2,1,inf)3符号与数值间的转换及符号的可变精度运算a：符号表达式转换成数值表达式p='1+sqrt/2'eval(p)b:数值表达式转换成符号表达式；例如p=1.7071;n=sym(p)n=17071/100004:符号表达式的简化a:见符号表达式类似于数学课本中的形式显示；例如symx;f=taylor(exp(-x)f=1-x+1/2*xA2-1/6*xA3+1/24*xA4-1/120*xA5pretty(f)b:合并符号表达式中的同类项；例如symsxy;f=sym(xA2*y+y*x-xA2-2*x);f=collect(f)f=(y-1)*xA2+(y-2)*xf=collect(f,y)仁(xA2+x)*y-xA2-2*xc:对符号表达式进行因式分解；例如symsx;f=sym(xA7-1);f=factor(f)f=(x-1)*(xA2+x+1)*(