基础篇matlab2010超强教程

第1MATLAB正常运行所必须具备的基础条件；二是简明地介绍MATLAB及其操作桌面Desktop的基本使用方法；第1MATLAB正常运行所必须具备的基础条件；二是简明地介绍MATLAB及其操作桌面Desktop的基本使用方法；三是全面介绍MATLAB的帮助系统。本章的前两节讲述：MATLABMATLAB环境的启动。因为指令窗性。本章专设第1.9MATLAB的帮助体系和求助方法。作者建议：不管读者此前是否使用MATLAB，都不要忽略本章的安装和工具包选MATLAB运行良好的先决条件。MATLAB本身可适应于许多机种和系PCUnixPC机用户来说，常常需要自己安装MATLABMATLABR2010a（MATLAB7.10）版要求WinXP或WindowsVista平台。下面介绍从光盘上安装MATLAB的方法。现，那么可以在<我的电脑>或<资源管理器>setup.exe应用程序，使“安装向导”启Typical⚫MATLAB软件所在的目◼◼该目录包含MATLAB运作所需的所有文件，如启动文件、各种工具包等。⚫ 自动生成的供用户使用的◼MATLABC:\DocumentsandSettings\acer\MyDocuments文件夹下。（注意：这文件夹名中的acer会随电脑不同而变。）该工作目录C:\DocumentsandSettings\acer\MyDocuments\MATLAB被自动记录在MATLAB的搜索路径中。因此，在这目录上的M文件、MAT文件、MDL文件等都能被MATLAB搜索到。◼◼ 操作桌面的启MATLAB当MATLAB安装到硬盘上以后，一般会在Windows桌面上自动生成MATLAB程序图标。在这种情况下，只要直接点击那图标即可启动MATLAB，打开如图1.2-1的MATLAB1假如WindowsMATLAB图标，那么点击matlab\假如WindowsMATLAB图标，那么点击matlab\。DesktopMATLABR2010aDesktopMATLAB工作界面。其默认形式，如图1.2-1所示。该桌面的上层铺放着三个最常用的界面：指令窗（CommandWindow）、当前目录（CurrentDirectory）浏览器、MATLAB工作内存空间（Workspace）浏览器、历史指令（CommandHistory）窗。1.2-1Desktop⚫该窗是进行各种MATLAB操作的最主要窗口。在该窗内，可键入各种送给MATLAB运作的指令、函数、表达式；显示除图形外的所有运算结果；运行错误时，在该浏览器中，展示着子目录、M文件、MATMDL文件等。对该界面上的M文件，可直接进行复制、编辑和运行；界面上的MAT数据文件，可直接送入MATLAB工作内存。此外，对该界面上的子目录，可进行Windows平台的各种标准操概况信息。比如该窗会展示：M函数文件的H1行内容，最基本的函数格式；所包含的内该浏览器默认地位于当前目录浏览器的后台。该窗口罗列出MATLAB工作空间⚫⚫2所有指令、文字都允许复制、重运行及用于产生M文件。⚫⚫所有指令、文字都允许复制、重运行及用于产生M文件。⚫⚫MATLAB指令窗（CommandWindow）的基本表现形态和操作方式。本书作者相信，通过本节的文字解释，读者将对MATLAB使用方法有一个良好的初始感受。CommandWindow指令窗简指令窗默认地位于MATLAB桌面的右方（见图1.2-1）。假如，用户希望1.3-1⚫⚫或选中指令窗菜单{Desktop：DockCommandWindow}便可。1.3-1[12274)]32的算术运算结果。本例演示：最初步的指令输入形式（1）用键盘在MATLAB指令窗中输入以下内>>(12+2*(7-3ans2⚫⚫>>”是“指令输入提示符”，它是自动生成的。本书在此后的输入指令前将不再带提>>理由是：（A）为使本ans2⚫⚫>>”是“指令输入提示符”，它是自动生成的。本书在此后的输入指令前将不再带提>>理由是：（A）为使本书简洁；（B）本MATLAB的M-book写成，而在M-book中运行的指令前是没有提示符的。一条指令输入结束后，必须按[Enter]键，那指令才被执行。⚫⚫⚫S=⚫⚫⚫MATLAB用3个或3个以上的连续黑点表示“续行”，即表示下一行是上一行的继续。指令执行后，变量S被保存在MATLAB的工作空间（Workspace）中，以备后用。如直到本MATLAB指令窗被关闭为止。学习MATLAB，有必要系统介绍一些基本规定。本节先介绍关于变量的若干规定。1-的数值采用习惯的十进制表示，可以带小数点或负号。以下记述3其相对精度是eps（MATLAB的一个预定义变量），大约保持有效数字16位。数值范103082⚫变量名、函数名是对字母大小写敏感的。如变量myvar和MyVar表示两个不同的变量sinMATLAB定义的正弦函数名，但SIN,Sin等都不是变量名的第一个字符必须是英文字母，最多可包含63个字符（英文、数字和下连符）。如myvar201是合法的变量名。变量名中不得包含空格、标点、运算符，但可以包含下连符。如变量名my_va_201是合法的，且读起来更方便。而my,va201MATLAB默认MATLAB为一些数学（MathContants）预定义了变量1.3-1⚫⚫3尽可能不对表1.3-1所列预定义变量名重新赋值，以免产生混淆。41.3-1MATLAB⚫“临时”覆盖。所谓“临时”是指：假如使用clearMATLAB内存中的变量，MATLAB指令窗被关闭后重新启动，那么所有的预定义变量将被重置为默认值，不在遵循IEEE0除是允许的。它不会导致程序执行的中断，只是在给出警告信息的同时，用一个特殊名称（Inf，NaN）记述。这个特殊名称将在以关于它们的更详细的帮助信息，可在MATLAB帮助浏览器左侧Contents页的⚫⚫1.3-1MATLAB⚫“临时”覆盖。所谓“临时”是指：假如使用clearMATLAB内存中的变量，MATLAB指令窗被关闭后重新启动，那么所有的预定义变量将被重置为默认值，不在遵循IEEE0除是允许的。它不会导致程序执行的中断，只是在给出警告信息的同时，用一个特殊名称（Inf，NaN）记述。这个特殊名称将在以关于它们的更详细的帮助信息，可在MATLAB帮助浏览器左侧Contents页的⚫⚫Contants>formatshortRMAd=RMAs=%%%%%int64整数类型时最大正（默认）整数类型时int16整数类型时最大正IMA64=IMA32=IMA32=e1=e2%%表达2时的绝对精度ans5预定义变 预定义变 浮点数相对精度NaN如0/0/i或虚单元ij Inf4运算符和表达（1）经典教科书上的算术运算符（ArithmeticOperations）MATLAB中的表达方式，见1.3-2。1.3-2MATLAB表达式的基本运算⚫只返还一个“主解”。要得复数的全部方根，必须专门编写程序（1.3-6）。因为MATLAB面向矩阵/数组设计，标量被看作(14运算符和表达（1）经典教科书上的算术运算符（ArithmeticOperations）MATLAB中的表达方式，见1.3-2。1.3-2MATLAB表达式的基本运算⚫只返还一个“主解”。要得复数的全部方根，必须专门编写程序（1.3-6）。因为MATLAB面向矩阵/数组设计，标量被看作(11)的矩阵/数组点”。（参见例1.3-9，例1.3-10。更详细说明请看第3章）关于它们的更详细的帮助信息，可在MATLAB帮助浏览器左侧Contents页的<MATLAB/UserGuide/ProgrammingFundamentals/BasicProgramComponents/Operators/ArithmeticOperations>节点找到。MATLABArithmetic⚫⚫⚫⚫⚫书写表达式的规则与“手写算式”几乎完全相⚫⚫⚫⚫⚫5面向复数设计的运算——zabireiMATLAB指令如下。z的实部arcos。zbrsinz的模a2b2zarctanbaz1.3-4z143iz212iz32e6z1z1=4+%合法，但建议少用或6数学表达矩阵运算数组运算加aa+a+减aa-a-乘aa*a.*除aa/bb\abb幂a^a.^(((z14⚫在这种书写格式中，4i4和i⚫z2=1+2*iz2z14⚫在这种书写格式中，4i4和i⚫z2=1+2*iz2%运算符构成的直角坐标%运算符构成的极坐标表1+z3+z+real_z=image_zmagnitude_zangle_z_radianangle_z_degree%弧度单%度数单算在复数域上进行；指令后“分号”的作用；复数加法的几何意义；展示MATLAB的可视能力（让读者感受，但不要求理解）%在一个物理行中，允许输入多条指%但各指令间要用“分号”或“逗号%指令后采用“分号”，使运算结果不clf,holdon%clf清空图形窗。逗号用来分隔两holdoff,gridon,axisequalz125.0000+71.3-21.3-6MATLAB计算38–2吗（1.3-3）？本例演示：MATLAB运1.3-21.3-6MATLAB计算38–2吗（1.3-3）？本例演示：MATLAB运 %3r_a1.0000+（2）38先构造一个多项式p(r)r3p=[1,0,0,-a]pp(r的系数向%指令末尾的“英文状态分号”使该指令运行后，不显示结果R%求多项式1.0000+1.0000-%计算复根%产生参变量在0到2*pi间的一组采样%画一个半径为R%注意“英文状态逗号”在不同位置的hold%画第一象plot(R([1,3]),'o','MarkerSize',15,'Color','b')%画另两个方86543210 %保证屏幕显示呈1.3-3（-8）⚫⚫6本例有助于MATLAB的计算特%保证屏幕显示呈1.3-3（-8）⚫⚫6本例有助于MATLAB的计算特对复数进行方根运算时只给出处于“第一象限”的那个根面向数组设计的运算——在MATLAB中，标量数据被看作(11)的数组（Array）数据。所有的数据都被存放在适当大小的数组中。为加快计算速度（运算的向量化处理），MATLAB对以数组形式存储算例展示MATLAB的计算特点，更详细的叙述请见第3章。 AR=（3）在指令执行AR指令窗中将显示以下结1234在MATLAB中，不必事先对数组维数及大小做任何说明，内存将自动配二维数组输入的三大要素：数组标识符“[]”；元素分隔符空格或逗号“,”；数组行间分⚫MATLAB对字母大小写是敏感比如本例中的数组赋给了变量AR，而不是Ar,aR,ar93210 AI5678⚫⚫1 3A2示：复数数组的生成 指令对数组元素“并行操作”的实质A1.0000-2.0000-%形成复数3.0000-4.0000-A_real1234A_image forforAa1(m,n)=angle(A(m,n))*180/pi以度为单位计算幅Am1= Aa1-78.6901--71.5651-Am2=Aa2-78.6901--71.5651-⚫函数real,imag,abs,angle是同时、并行地作用于数组的每个元素。对4个元素运算所Aa2-78.6901--71.5651-⚫函数real,imag,abs,angle是同时、并行地作用于数组的每个元素。对4个元素运算所需MATLAB来说，应该尽量摒弃“循环”处理，而采用“向量化”处理方式。⚫t1.3-10ye3sin3tt的取值范围是[0,4]（1.3-4）。本例演示：展示数组运算的优点；展示MATLAB的可视化能力。%定义自变量t的取值数%计算与自变量相应的y数组。注意：乘法符前面的小黑%绘制曲1.3-⚫耗费机时的“循环计算”。关于数组运算的详细叙述请见第3章。⚫3 2B【例1.3-1151.3-9）。本例演示矩阵运算指令的简捷性y10 t%复数数组的又一种输入%注意标点符号的%矩阵乘B+ 2.0000+5.0000%复数数组的又一种输入%注意标点符号的%矩阵乘B+ 2.0000+5.0000+ 4.0000-C-39.0000i30.0000-62.0000-42.0000i40.0000-⚫AB的行数相等，所以可以进行A乘B。那是由于MATLAB的设计者采用了“面向对象”编程技术。⚫⚫前一节借助算例，使读者对MATLAB指令窗的使用方法有了一个直观的感受。本MATLAB，更方便地使用MATLAB1从MATLAB7.0起，指令窗中的字符、数值等采用更为醒目的分类显对于输入指令中的if,for,endMATLAB关键词自动地采用蓝色字体显示。2运算结果的显5位有效数字。实际上，MATLAB的数值了比较简洁、紧凑地显示数值输出，才默认地采用formatshortg格式显示出5位有效数中进行选择，都可获得所需的数值计算结果显示格式。MATLAB数值计算结果显示格式的类型见表1.4-1。1.4-1指令 7位；对于大于1000的实数，用5位314.159314.1590；format3{File：Preferences}下拉菜单项，引出一个参数设置对话框；在此弹中的地位极其重要。为此，把各标点的作用归纳成表1.4-2MATLAB常用3{File：Preferences}下拉菜单项，引出一个参数设置对话框；在此弹中的地位极其重要。为此，把各标点的作用归纳成表1.4-2MATLAB常用标点的功 ,.;:%'formatlong从formatshort和formatshorte中自动选formatlong从formatlong和formatlonge中自动选formatformatformat空格表示+formatformatformatformatshort显示格式是默认的显示格式该表中实现的所有格式设置仅在MATLAB的当前执行过程中有效1.4-3⚫⚫表1.4-3所列的指令是基本的，它们对MATLAB各版都适用种“等价”仅对“人机交互”过程而言。至于这些指令在M文件中的作用仍是不可替cd及dirMATLAB操作桌面上或当前目录浏览器中的“浏览键（Browser）”替代。关于当前目录浏览器的使用，请看第1.6节。clear清除内存变量的操作，可以等价地在“工作空间浏览器”交互界面中实现。关于工作空间浏览器的使用，请看第1.7节。⚫⚫⚫1.4-3⚫⚫表1.4-3所列的指令是基本的，它们对MATLAB各版都适用种“等价”仅对“人机交互”过程而言。至于这些指令在M文件中的作用仍是不可替cd及dirMATLAB操作桌面上或当前目录浏览器中的“浏览键（Browser）”替代。关于当前目录浏览器的使用，请看第1.6节。clear清除内存变量的操作，可以等价地在“工作空间浏览器”交互界面中实现。关于工作空间浏览器的使用，请看第1.7节。⚫⚫⚫ 图标。关于M文件编辑器的使用，请看第指 指 打开M文件编辑列出指定目录下的文件和子目显示指定M文件的([{_“At”@MATLABContents<MATLAB/UserGuideProgrammingFundamentals/BasicProgramComponents/SymbolReference>节点找到。为了操作方便，MATLAB不但允许用户在指令窗中对输入的指令行进行各种编辑和运1.4-4MATLAB1.4-4所列的操作对MATLAB各版均适用事实上，MATLAB把指令窗中输入的所有指令都记录在内存中专门开辟的“为了操作方便，MATLAB不但允许用户在指令窗中对输入的指令行进行各种编辑和运1.4-4MATLAB1.4-4所列的操作对MATLAB各版均适用事实上，MATLAB把指令窗中输入的所有指令都记录在内存中专门开辟的“指令历史空间（CommandHistory）”中，只要用户对它们不进行专门的删除操作，它们既不会因为用户对指令窗进行“清屏”操作（即运行clc指令）而消失，也不会因用户对“工作空间”进行“清除内存变量”（即运行clear指令）而消失。和调用。关于“历史指令浏览器”的使用请看第1.5节。2sin(0.3y11 y1=2cos(0.3y21 []键调回已输入过的指令y1=2*sin(0.3*pi)/(1+sqrt(5))；然后移动光y1y2sin改成cos；再按[Enter]键，就可得到结果。即y2=CommandHistory历史指时蹦出的思想“火花”可即刻通过计算加以验证。历史指令窗（CommandHistory）就是为键 键 →CommandHistory历史指令窗简记录的指令行都能被复制，或再运行。关于历史指令窗的功能详见表1.5-1。1.5-1CommandHistory历史指令窗简记录的指令行都能被复制，或再运行。关于历史指令窗的功能详见表1.5-1。1.5-11.5-1⚫⚫Directory、路径设置器和文件管简捷操作方单行或多行指令点亮单行指令；按鼠标右键引出现场菜单；选中鼠标左键双击点亮多行指令；按鼠标右键引出现场菜单；选中见到相应结果。（详见例1.5-1）把多行指令写成M文点亮多行指令键引出现{CreateM-File}菜单项引出书写着这些指令的M文件编辑调试器；再进行相应操作，即可得所需M文件。当在指令窗中运行一条指令时，MATLAB是怎样从庞大的函数和数据库中，找到所需怎样能与MATLAB原有环境融为一体呢？这就是本节要介绍的内容。CurrentDirectory当前目录当在指令窗中运行一条指令时，MATLAB是怎样从庞大的函数和数据库中，找到所需怎样能与MATLAB原有环境融为一体呢？这就是本节要介绍的内容。CurrentDirectory当前目录件、文件夹列表及文件描述区等。此外，MATLAB还为当前目录窗设计了一个专门的操作菜单。借助该菜单可方便地打开或运行M文件、装载MAT文件数据等。详见表1.6-1。1.6-11.6-1简捷操作方运行M文点亮待运行文件；按鼠标右键引出现场菜单；选中{Run}菜单项，即可使该M文件运编辑M文点亮待运行文件；按鼠标右键引出现场菜单；选中{Open}菜单项，此M文件就出现在编辑/鼠标左键双击M把MAT文件全部{Open}菜单项，此文件的数据就全部装入工作内存鼠标左键双击MAT把MAT文件部分{ImportData}菜单项，引出数据预览选择对话框⚫MATLAB启动后的默认当前目录通常是：C:\DocumentsandSettings\acer\My可靠的。MathWorks公司之所以设计这样一个目录，就是供用户使用的。⚫MATLAB启动后的默认当前目录通常是：C:\DocumentsandSettings\acer\My可靠的。MathWorks公司之所以设计这样一个目录，就是供用户使用的。⚫（1）用户MATLABR2010aC:\DocumentsandSettings\acer\MyDocuments\MATLABMATLAB文件。假若用户想另建一个工作目录，采用Windows规范操作就可实现（2）开设的“C:\DocumentsandSettings\acer\MyDocuments\MATLAB”设置成当前目录。MATLAB操作桌面右上方，或当前目录浏览器左上方，都有一个当前目录设置区。M文件中。假设c:\mydircdc:\mydir。注意：以上方法设置的当前目录，只是在当前开启的MATLAB环境中有效。一重新启动，以上设置操作必须重新进MATLAB的搜索把这些目录按优先次序设计为“搜索路径”上的各个节点。此后，MATLAB工作时，就沿当用户从指令窗送入一个名为cont的指令后，MATLAB的基本搜索过程大抵如下。检查MATLAB内存，看cont是不是变量；假如不是变量，则进行下一步。⚫⚫⚫检查cont是不是内建函数（Built-inFunction）在当前目录上，检查是否有名为cont的M文件存在；假如不是，再往下执行。在MATLAB搜索路径的其他目录中，检查是否有名为cont的M文件存在。凡不在搜索路径上的内容，不可能被搜索。（B）指令existwhichload执行时，也都遵循MATLAB搜索点击[Finish]键，就完成操作。1假如用户有多个目录需要同时与 1假如用户有多个目录需要同时与 21.6-2⚫假如在设置后，点击了对话框下方[Save]按键。⚫3所谓永久有效修改是指：所进行的修改不因MATLAB的关闭而消path设置路利用path指令设置路径的方法对任何版本的MATLAB都适用假设待纳入搜索路径把c:\my_dir设置在搜索路径的尾端把c:\my_dir⚫path指令扩展的搜索路径仅在当前MATLAB环境下有效。也就是说：若用户退出当前MATLAB后，再重新启动MATLAB，那么在前一环境下用path所定义的扩展搜索路径⚫工作空间浏览器和变量编辑该浏览器的功用，详见表1.7-1。1.7-11.7-1该浏览器的功用，详见表1.7-1。1.7-11.7-1点击图标，在工作空间中生成一个“unnamed”的新变量；双击该VariableEditor变量编辑器（1.7-3）；在变量 则变量内含的数据就显示在“VariableEditor”变量编辑器中。（参见第1.7-3节）选择待保存到文件的（一个或多个）菜单中的{SaveWorkspaceAs}项，便可把那些变量保存到MAT数据文件。（详见第1.7-4节之一）在“工作空间浏览器”中，用鼠标点亮所需图示的变量y和t在“工作空间浏览器”中，用鼠标点亮所需图示的变量y和t标点中“填色面图”栏（参见图1.7-1），即可。ytarea(y,t)。显然，指令area的两个输入量位置不正确。为纠正这种错误，用鼠标点击[输入量位置交换键]，参见图1.7-2数组(t,y)⚫观特征的卓越能力。正是出于这种考虑提供了“图示数据”的多种途径本节要介绍管理工作空间的三个常用指令：who,clear,packwho1【例1.7-2】在指令窗中运用who,whos查阅MATLAB内存who,whos在指令窗中运行后的显示结果如下Yourvariables点 点 ansy8s⚫⚫⚫,ansy8s⚫⚫⚫,2clearall清除工作空间中的var1和var2清除工作空间中所有的变量、全局变量、编译过的MMEX清除工作空间中名为fun1和fun2clearfun13整理工作空间内存碎片的指令在MATLAB运行期它会自动地为产生的变量分配内存，也会为使用到的M函数配内存。有时对于容量较大的变量，会出现“Outofmemory”的错误。此时，可能使用clear指令清除若干内存中的变量也无济于事。产生这种问题的一个原因是：MATLAB存放一个变pack指令也许VariableEditor变量1.7-3双击工作空间浏览器中的变量图标，将引出如图1.7-3所示的变量编辑器Variablebrush）”与图形窗中的图形相关联（参见第5.5节）。brush）”与图形窗中的图形相关联（参见第5.5节）。1⚫⚫借助工作空间浏览器产生保存变量的MAT点击工作空间浏览器工具条图标，或选中弹出菜单中的{SaveAs}项，就弹出Windows标准的目录和文件名输入对话窗动生成的用户工作目录，即C:\DocumentsandSettings\acer\MyDocuments\MATLAB目MAT数据文件的产生借助输入向导ImportWizard向工点击工作空间浏览器上的图标，或MATLABDesktop下拉菜单{File>ImportData}，引出Windows标准的目录和文件选择对话窗。在用户希望的目录上，选中MAT数据文件（如光盘For2010a\mfiles子目录上的的t和y）；再点击[Finish]按键，变量t和y就被装载到工作空间。⚫⚫2⚫⚫1.7-4save利用saveloadMATLAB3savesaveFileNamev1把全部内存变量保存为FileName.mat文件把变量v1,v2保存为FileName.mat文件--把变量v1,v2添加到FileName.mat把变量v1v2保存为FileName8ASCII文件把变量v1,v2保存为--把变量v1,v2添加到FileName.mat把变量v1v2保存为FileName8ASCII文件把变量v1,v2保存为FileName16ASCII文件把FileName.mat文件中的全部变量装入内存把FileName.mat文件中的v1,v2变量装入内存把FileNameASCII文件中的v1,v2变量装入内---⚫⚫文件名可以带路径，也可以带扩展-asciiASCII格式处理。生成的（不带扩展名的）ASCII文件可以在任何“文字处理器”中被修改。如果数据较多的变量需要进行修改，那么ASCII格式的数如果指令后没有-ascii⚫⚫cdc:\my_dir%在C盘上创建目录%c:\my_dir成为当前目%选择内存中的X,Y,Z变量保存为saf.mat%显示目录.safYZ（2）清空内存，从saf.mat向内存装载变量Z%清除内存中的全部%把saf.mat文件中的Z变量装入内%检查内存中有什么saf⚫saveload重新获取这组数据。这种处理模式常和脚本编写初的方法就显得烦琐和笨拙。M脚本文件就是设计来解决这个矛盾的。Editor/DebuggerM文件编辑器简1.8-1MM文件时才启动。MMM1.8-1MM文件时才启动。MMM文件进行交互式调试；M文件编辑器不仅可处理带.m扩展名的文件，而且可以阅读和编辑其他ASCII码文M文件编辑器的启动方法有以下几种⚫MATLAB桌面窗口输入指令edit，都可以打开空白的M文件编辑器⚫MATLAB桌面图标，或选中菜单项{File>Open}，可引出“Open”文件[Open]键，就可引出展示相应文件的M文件编用鼠标左键双击当前目录窗中的所需MM文件编辑⚫MMATLAB中去执行。（3）与在指令窗中直接运行指令一MATLAB基本工作空间中。（4）文件扩展名是1.8-1】编写解算例1.3-10题目M脚本文件，并运行之。⚫⚫在历史指令窗中，找到算例1.3-10的运行指令，并把它们选中点亮，参见图1.8-2。WindowsC:\DocumentsandSettings\acer\MyDocuments\MATLABexm010801为文件名进行保存，于是就得到了exm010801.m文件。因为C:\DocumentsandSettings\acer\MyDocuments\MATLAB是MATLAB安装时自动⚫⚫成“在搜索路径上的”目录，所以直接点击M文件编辑器上的工具图标可以得到如图1.3-4的曲线。1.8-2利用历史指令创建M⚫⚫1.8-2利用历史指令创建M⚫⚫帮助系统及其使读者接触、学习MATLAB的起因可能不同，借助MATLAB所想解决的问题也可能MATLAB老用户很想知道的是：MATLAB新版本有什么新特点、新功能。再如对科怎么用；或想解某个具体问题，不知道MATLAB有哪些指令可用。MATLAB作为一个优秀的科学计算软件，其帮助系统考虑了不同用户的不同需求，构MATLABMATLAB的1.9-1MATLAB的帮助体帮助形 直接从指令窗中，通过help指令获得；所有包含在M文件之中的帮助注释内容。HTML形式；系统叙述位于matla\help目录下；通过帮助浏览器获得；HTMLXMLM文件，是Web网mathworks公司网站；MATLAB操作界面下拉菜单{Help>WebResources>}；本书不作介绍。1helpFunNamehelpwinFunNamedocFunName⚫1helpFunNamehelpwinFunNamedocFunName⚫helpM文件帮助注释区的内容。这部分资源用纯文本形式写成。它简扼docHTML文件构成的帮助子系统中进行的。HTMLM文件资源M文件帮助注释详细。该子系统，由于采用“超链接”机理，因此检⚫⚫⚫2lookforKeyWord对MH1行进行单词条检索对HTML子系统进行多词条检⚫⚫lookfor搜索的资源是docsearch指令使用格文件帮助注释区中的第一行（简称H1行）◼格式：docsearch('Word1Word2')docsearch('Word1OR格式：docsearch('"Word1Word2"')将对由Word1Woed2◼◼⚫搜索到的内容也比较详细。该词条搜索指令的功能与帮助导航器中的“Search搜索窗”Help1帮助浏览器 Brower）搜索的资源是Mathworks专门创建的HTML随“机”帮（1）MATLABDesktop⚫⚫（1）MATLABDesktop⚫⚫ 方法二：选中下拉菜单项{Help>MATLAB⚫在指令窗中运行 或helpdesk1.9-12⚫◼◼◼Contents分类目录活页SearchResults检索结果活页⚫◼◼⚫MSIMULINK类◼M码类目录的说“一级目录”分别是：ReleaseNotes目录，发布新产品、⚫◼◼⚫MSIMULINK类◼M码类目录的说“一级目录”分别是：ReleaseNotes目录，发布新产品、版本升级老版本的修订、兼容性等信息；InstallationMATLAB及各工具包的安装和激活信息；MATLABMATLAB功能、规则、和基本函数指令等关于M码工具包的分类目录。各工具包目录按英文字母表排序；通用性较强的工具包有：OptimizationToolbox优化工具包，StatisticsToolbox统计工具包，SymbolicMathToolbox符号计算工具包。◼SIMULINK的分类目◼⚫◼◼◼◼◼◼⚫在搜索框里，既可以输入函数指令名（如inv），又可以输入各种专业词条（如inverse⚫◼格式：Word1Word2，或Word1ORWord2格式："Word1Word2"将对由Word1Woed2◼◼laplacethranform，按[Enter]键进行搜索，那么帮助浏览器⚫◼◼◼Product产品类型排列方如果用鼠标◼Product产品类型排列方如果用鼠标点击帮助浏览器左侧的[Product]按键1.9-23，或选中下拉菜单项当点击MATLAB工作台的工具条图标简捷的“通往各种帮助文件的超链接通道”（参看图1.9-1的右侧）。具体如下：⚫⚫⚫⚫⚫PDF文件超链接通道⚫MathWorks网站资源超链接通供用户查阅⚫MathWorks网站资源超链接通供用户查阅MATLAB各种资源或问题解（1）1.9-3⚫ByCategoryAlphabeticalList按字母排列的函数指令帮助通ObjectProperties图形对象属性及舌质指令帮助MATLAB通道特点：层次清楚；属性分列清晰；超链接交互查阅方便。⚫⚫（2）1.9-4⚫适用场合：MATLAB新手不必看；需要了解MATLAB历史演进较多、较重要的MATLAB历史版本编写的文件资源启用通道特点：对各版变化、补丁修正描述清晰。GeneralReleaseNotesforR2010a使用场合：MATLAB新手可浏览新版本对环境的要求；MATLAB设计较大型软件库的用户； 历史版本写成文件运行产生不明原因错误的场⚫（3）1.9-5⚫GettingStarted快速入通道特点：所包含的网络接口通道，可直接观看Mathworks1.9-5⚫GettingStarted快速入通道特点：所包含的网络接口通道，可直接观看Mathworks制作的入门视频；⚫⚫⚫ExamplesinDocumentationTips编程技⚫（4）1.9-6⚫MATLABDemosMATLAB界面各图标、菜单的功能，以及交互操作手法；了解、学习典型的M码文件的指令运用和编写技巧；⚫MATLABDemosMATLAB界面各图标、菜单的功能，以及交互操作手法；了解、学习典型的M码文件的指令运用和编写技巧；通道特点：Video视频演示需网络支持；GUI资源、及部分M码资源所提供的帮助具有独特性研人员适宜），还有“编程技巧”型（MATLAB较熟悉者适宜）（5）PDF1.9-7PDF⚫PrintableversionsPDF格式帮助文适用场合：需要系统阅读、学习UserGuider用户指南。通道特点：需要得到网络支持；（6）MathWorks1.9-8MathWorksMATLAB产品、第三方相关产品信息；向MathWorks公司进行技术咨询。通道特点：需网络支持。1.5e2，1.5e31500相同？ MATLAB环境中，比1大的最小数是多少a8指令clear,clf,clc指令clear,clf,clc同。”（2）MATLAB指令窗中显示的数值有效位数不超过7位。”258S=[123;456;789]%整个指令在中文状态下1.3-5编写一个解题用的M第2章符号对象和符号表达12【例第2章符号对象和符号表达12【例2.1-1】asapi+5^(1/2)Ca=Csa=ans=34symsuvwza5ans=[a5,u,v,w,ans[w,z,v,u,result_1=1a5/3-(v*z)/3-result_2=-(v-(v^2+4*a5*u*sin(3)-(va5/3-(v*z)/3-result_2=-(v-(v^2+4*a5*u*sin(3)-(v+(v^2+4*a5*u*sin(3)-symsabxXYans=[X,Y,a,b,c,d,t,x,ans=[x,y,t,d,c,b,a,X,ans=xans=[x,y,???Errorusing==>atError:argumentmustbeof'Type::Arithmetical'sym.sym>sym.sym/char2symat337Scell=scalarsym(x);sym.sym>sym.sym/symcharat182Scell=char2sym(x);sym.sym>sym.symatS=E4E52symsabtuvxyA[a+b*x,[x/exp(t),ans=xu+v+symsabtuvxyA[a+b*x,[x/exp(t),ans=xu+v+符号对象的识Mn1324McMs=[a,[c,SizeMn=22SizeMc1SizeMs292CMn=CMc=3ans1ans1ans1whosans1ans1ans1whosBytes183260123symsxclearrf1/4-(79^(1/2)*i)/4+symsxrealrfr=clearxsymsxWarning:Explicitsolution>Insolveat98rg=[emptysymnotbe4symsxclearrg=-(19^(1/2)*i)/2(19^(1/2)*i)/2clearallsymssaswsxsysymsAsymssymsxclearrg=-(19^(1/2)*i)/2(19^(1/2)*i)/2clearallsymssaswsxsysymsAsyms88ans={A,B,ans{A,B,symsBclear5ans{A,ans{A,2.1-62.1-⚫2.1-⚫72.1-符号数字及表达式的操2.1-符号数字及表达式的操12sa2^(1/2)+8Digits=formatlongasa_Plus_asa_Minus_aDigits=formatlongasa_Plus_asa_Minus_a%%%%sa32sa5=sa48=a=a2=ans0ans0symsxf(12/x+6/x^2+1/x^3+8)^(1/3)g1=((2*x+91A=sym('[ab;cA[a,[c,d]V=[D1A=sym('[ab;cA[a,[c,d]V=[D[sigma(a^2-2*a*d+ans=[(a/2+d/2%%d^2+sigma/2)/c-d/c,(a/2d/2+sigma/2)/c-[[[a/2+d/2-variablesa/2d/2+ADVw=(a^2-2*a*d+d^2Dw=[a/2+d/2-+a/2+d/2+[RVD[[w=+++d/2-w/2)/c-d/c,a/2+d/2-2*a*d+d^2+-+a/2+2symsabx;f=b+f1b+a*log(y)ans=symsabx;f=b+f1b+a*log(y)ans=f2b+(311*sin(x))/100ans=f3=3^(1/2)+5formatcompactf4=ans=f5[b,a+b,b]ans=f6Columns1through Columns9through Columns17through 2.2-ans62.2-321 2.2-ans62.2-321 54321 符号微积symstxk符号微积symstxk1Lf=Lf1=symsatx;df[[-t*sin(x),1/x]dfdt2=[0,[ dfdxdt[0,[-sin(x),symsx1x2;Jf[[[-sin(x1)*sin(x2),symsdpositive%df_p=%%df_ndf_n0dfdx=dfdx0=0%%holdongridon%%%%df_ndf_n0dfdx=dfdx0=0%%holdongridon%%2.3-10 xgcos(x+sin(y(x)))=%%dgdxdgdx1=subs(dgdx,'diff(y(x),x)','dydx')gcos(x+sin(y(x)))=%%dgdxdgdx1=subs(dgdx,'diff(y(x),x)','dydx')dgdx1-sin(x+sin(y(x)))*(dydx*cos(y(x))+1)=dydx=-sin(x++cos(y(x))*sin(x+symsxr%%+++x^5/24+x^4/6+x^3/2+x^2+8x7x6x+--5x4x3x2+++--+x+xR=%%Rx+x^2+x^3/2+x^4/6+x^5/24+x^6/120+x^7/720+x^8/5040+83 x+x+-24x5x6x+--7x+--x9+----+O(x)【例2.3-7】TL1=(x^6*y^2)/12-x^6/6+(x^4*y^3)/12-(x^4*y)/2-(x^2*y^6)/720+(x^2*y^4)/24-(x^2*y^2)/2+x^2-y^7/5040+y^5/120-y^3/6+y序列/级数的符号求symsnks1=1-序列/级数的符号求symsnks1=1-1/(n+s2=piecewise([abs(x)<1,s3[pi^2/8,-symsabxf1=s1(x^2*(log(x)-1/2))/2s1=x^2*(log(x)/2-%【例2.3-10】disp('')disp('Theintegraloffis')f2[a*x,b*x^2][1/x,Thef|||||| b 2a||||symsxyzWarning:Explicitintegralcouldnotbefound.F2=(14912*2^(1/4))/4641-(6072064*2^(1/2))/348075+symsxyzWarning:Explicitintegralcouldnotbefound.F2=(14912*2^(1/4))/4641-(6072064*2^(1/2))/348075+VF2=(64*2^(3/4))/225symsarthetaphiwarningoff%%L+ ++L_2pilog(2*pi+(4*pi^2+1)^(1/2))/2+pi*(4*pi^2+1)^(1/2)L_2pi_vpa=gridonholdon%title('')holdoff%%2.3-2阿基米德螺线（粗红）和螺线长度函数（细蓝微分方程的符号解clear%S=dsolve('Dx=y,Dy=-x')%disp('')Sx:y:[x+yC1*cos(t)-cleary%%2.3-2阿基米德螺线（粗红）和螺线长度函数（细蓝微分方程的符号解clear%S=dsolve('Dx=y,Dy=-x')%disp('')Sx:y:[x+yC1*cos(t)-cleary%%y螺线长度-420 xC3*x-clf,hold%%fork=-2:0.5:2boxon%%%C3*x-clf,hold%%fork=-2:0.5:2boxon%%%–+y=0','的解2.4-1y(31*x^4)/468-x^3/3+holdontitle(['x*D2y-3*Dy=x^2',',y(1)=0,y(5)=0'])holdoffyny微分方程(y2–xyy0的解8奇解6420 x02.4-2符号变换和符号卷Fourier变换及其反变symstwUT=pi*dirac(-w)-02.4-2符号变换和符号卷Fourier变换及其反变symstwUT=pi*dirac(-w)-Ut(pi+pi*(2*heaviside(t)SUt=%-holdon%%x*D2y-3*Dy=x2,5432 0 x%gridon2.5-1Heaviside(t)symsAtw%gridon2.5-1Heaviside(t)symsAtw%%%%YwA*((1/exp((tao*w*i)/2))*(-pi*dirac(-w)+i/w)-exp((tao*w*i)/2)*(-pi*dirac(-w)+i/w))Yw_fy=Yw_fy_e=%%Yt(A*transform::fourier(sin((tao*w)/2)/w,w,t))/piYt_e=-A*(heaviside(t-tao/2)-heaviside(t+ut10 tkk=find(tn==-t3/2|tn==t3/2);%holdon%holdoffgridyt13=+--2.5-2Heaviside(t)kk=find(tn==-t3/2|tn==t3/2);%holdon%holdoffgridyt13=+--2.5-2Heaviside(t)构造的矩形10 2.5-3symstx%%F1=(1/exp(w*x*i))/(1+w*i)G1=%%2.5-3symstx%%F1=(1/exp(w*x*i))/(1+w*i)G1=%%t),t,-F2=-exp(t^2*i)/(-1+t*i)F3=-(1/exp(t*w*i))/(-1+%%Laplace变换及其反变symstsa(2sin((33210 w(2abs(sin((33210 wf1=F1b/((a+s)^2+%symaclearans=heaviside(t-a)F2=%f1=F1b/((a+s)^2+%symaclearans=heaviside(t-a)F2=%%-t,symsapositiveF3=%%F4=piecewise([b<0,F5=ft_F5=dirac(a-t)%%[0%%%F6=piecewise([-1<%%++%%n=gamma(n+1)/s^(n+Zsymsnzclear%%%%gn6-6z223z+-symswTz%%%pretty(FW),disp('')%z2zgn6-6z223z+-symswTz%%%pretty(FW),disp('')%z2z-w)z+inv_FW=symsnzclear%zz-1symsnzclearinv_KD=%%symsnzclear%%FDpiecewise([kinZ_,inv_FD=piecewise([kinZ_,f(k)/z^k],[Otherwise,f(k)*kroneckerDelta(k-n,0)],symsaznclear%%gn【例2.5-6】symsTt%gn【例2.5-6】symsTttaoyt=-(1/exp(t)-1/exp(t/T))/(T-%%%%symssyt-(1/exp(t)-1/exp(t/T))/(T-符号矩阵分析和代数方程symsa11a12a21a22A[a11,[a21,a22]DA=a11*a22-a12*a21IA=[a22/(a11*a22-a12*a21),-a12/(a11*a22[-a21/(a11*a22-a12*a21),--D=(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)EA=a11/2+a22/2-D/2a11/2+a22/2+symst A[GA=1/2,[(3^(1/2)*cos(t))/2-sin(t)/2,cos(t)/2-[cos(t)/2+(3^(1/2)*sin(t))/2,sin(t)/2+%%%%A[GA=1/2,[(3^(1/2)*cos(t))/2-sin(t)/2,cos(t)/2-[cos(t)/2+(3^(1/2)*sin(t))/2,sin(t)/2+%%%%%holdonLstr=['旋转前的v1';'旋转前的v2';'旋转后的holdu1';'旋转后%%grid5.6-1GivensA=sym([11/21/2-1;11-11;1-1/4-11;-8-1%Givens10 旋转前的旋转前的 旋转旋转 889%%%%disp(['d','n','p','%%dnp889%%%%disp(['d','n','p','%%dnpq[1,8,8,【例2.6-4】S=y:[2x1sym]z:[2x1%(v+2*u*w+(v^2+4*u*w*v-4*u*w)^(1/2))/(2*u)(v+2*u*w-(v^2+4*u*w*v-4*u*w)^(1/2))/(2*u)--(v+2*u*w+(v^2+4*u*w*v--(v+2*u*w-(v^2+4*u*w*v-symsdnp%% [z/3-2,8,(4*z)/3-4,clearall,symsx;s代数状态方程求符号传递函symsG1G2G3G4H1H2 0,symsG1G2G3G4H1H2 0, b=[ disp([blanks(5),'传递函数Y2Ua传递函数Y2UaG1G2G3G2G3H1-G3G4H2+G2G3G4+1此G3G4H1 H20,-0,-0, Y2Ub=c*((eye(size(A))-A)\b);disp([blanks(5),'传递函数Y2Ub为'])传递函数Y2UbG1G2G3G2G3G4H1-G3G4H2+G1G2G3G4+1syms%%参数具体化的传递函数Y2Uc2100s300+%%参数具体化的传递函数Y2Uc2100s300+421s+15735s2s+++s+⚫符号计算结果的可视1symstgridongridontitle('s=\inty=%%%%%%s2%%%%shadinginterp%%%%(2cos((31/2t)/2))/(32%%%%shadinginterp%%%%(2cos((31/2t)/2))/(30 ts= t2.8-ezsurf在参变量格式下绘制的2.8-ezsurf在参变量格式下绘制的symsxyrealdisp('')2%%%1exp(x)+fxint=-+%%%%legend('f(x)','\int^x_0f(x)2.8-3gy=subs(finverse(fx),x,y)gylog(-(y+1)/(y-%gyintpiecewise([y<1,log(1-y^2)2.8-3gy=subs(finverse(fx),x,y)gylog(-(y+1)/(y-%gyintpiecewise([y<1,log(1-y^2)+y*log(y+1)-y*log(1-y)],[1log(y^2-1)+y*log(-(y+1)/(y-1))-pi*i])gf=x%%%%holdonholdofflegend('直接法计算反函数积分'互补法求反函数积分函数及其积分函 xf(x)00 x2.8-4TL1(x^6*y^2)/12-x^6/6+(x^4*y^3)/12-(x^4*y)/2-(x^2*y^6)/7202.8-4TL1(x^6*y^2)/12-x^6/6+(x^4*y^3)/12-(x^4*y)/2-(x^2*y^6)/720+(x^2*y^4)/24-(x^2*y^2)/2+x^2-y^7/5040+y^5/120-y^3/6+yFxy=sin(x^2+y)%Fxy_TL1sin(x^2+y)-+(x^2*y^2)/2-(x^2*y^4)/24-(x^4*y^3)/12+shadinginterp%%%%积分0 y2.8-52.8-52.8-6Taylor2.8-6Taylorshadinginterp%2.8-7shadinginterp符号计算资源深入利符号表达式、串操作及数值计算M符号计算资源深入利符号表达式、串操作及数值计算M符号工具包资源表达式转换成M1把符号包资源转换成M码函数的示s(a^4*b^3+3*yR*a^4-a^3*b^4-3*yL*b^4)/(3*a^4-3*b^4)-x^3/3+(x^4*(a^3-b^3+3*yL-3*yR))/(3*a^4-3*b^4)Hs=%%%ynColumnsthroughynColumnsthroughColumn5.9-2functiony=%%%%%%%%%t7=a.^2;t8=b.^2;t9=t7.^2;t11=t12=3.*t10;t13=t12-3.*t9;t14=1./t13;y=t14.*(3.*t10.*yL-3.*t9.*yR+a.*t10.*t7-b.*t8.*t9)-(x.*t11)./3-t14.*(3.*yL-3.*yRyxy-3y=x2,y(1)=0,5 0 xa.*t7-ym=Columns1through% Columna.*t7-ym=Columns1through% Column◼〖functiony=%%%%%%%%iferror('输入量数目7)if%y=y0y=exm050902_ZZY(de2,x,3,6,-5.9-3用exm050902_ZZYmfun调MuPAD5.9-3用exm050902_ZZYmfun调MuPAD特殊函symstgt_0=limit(gt,%gridon%yx*D2y-50 x5.10-1fx=piecewise([x<1,%integralcouldnotbeLi(x)],[Otherwise,t=%%%fxMfunColumns18 5.10-1fx=piecewise([x<1,%integralcouldnotbeLi(x)],[Otherwise,t=%%%fxMfunColumns18 Columns9through Columns17through holdholdoff% t5.10-2被积函数曲线gt和积分曲线%fx_matlabColumns1through Columns9through Columns17through 5.10-2被积函数曲线gt和积分曲线%fx_matlabColumns1through Columns9through Columns17through %3/7+0.1,sym(3/7+0.1),sym('sin(w*t)'),sym('a*exp(-X)'),（提示：关于符号变量的假设要注意 44.50正实根的程序。注意：只要正实根，不要出现x2axa20a0a=@,b=sym(@),c=sym(@,'d'),d=sym('@')在此，@分别代表具体数值7/3,pi/3,pi*3^(1/3)vpa(abs(b-d)),vpa(abs(c-d))等来观察。；而异同通过vpa(abs(a-d)),0 ta13 5Aa23a33f(k)kk（提示：symsumk实际上，这就是根据定义求Z变换问题。x12k2 。（提示：理论结果为lnx；注意限定性假设2k1xa13 5Aa23a33f(k)kk（提示：symsumk实际上，这就是根据定义求Z变换问题。x12k2 。（提示：理论结果为lnx；注意限定性假设2k1x1ksintdy。（2）8．（1）y(t)和t0。2tesinxdx64位有效数字的积分值。（提示：intvpa,52(xy)dydx 11sinttx11．在[0,2yx)0y(1)y(n)sinn212n00313.有序列x(k)akh(k)bk，（k0ab），求这两个序列的卷积ky(k)h(n)x(kn。（提示：symsum,h(te3t，求该系统在输入u(tcostt0（提示：直接卷积法，变换法均可f(t)Aet,0的Fourier变换。（提示：注意限定ttt16f(t0s的Laplaces33s26s18．利用符号运算证明Laplace变换的时域求导性质：Ldf(tsLf(t)f(0)。（ f(k)kek的Zx2y21,xy2的解。（提示：正确使用0） 22．采用代数状态方程法求图p2-2所示结构框图的传递函 。（提示：列出正 x0） 22．采用代数状态方程法求图p2-2所示结构框图的传递函 。（提示：列出正 xAxbUYcxdUc(IA)1bdYc(IA)1fgNYU图p2-01p2-3x 24xat224xat2btx(02 34ff(0)g(0)1第3章数值数组及向量化运数值计算的特点和地symstx%%%%ft=sx=-2*sin(x)+第3章数值数组及向量化运数值计算的特点和地symstx%%%%ft=sx=-2*sin(x)+%%%%%%%%%%%ans ans-19.9833-19.7907-3.1-1在区间[0,5]150 xsymstxft=Warning:Explicitintegralcouldsx=nott=%%%%%%%%%holdonholdoffans=%3.1-2在区间[0,4]数symstxft=Warning:Explicitintegralcouldsx=nott=%%%%%%%%%holdonholdoffans=%3.1-2在区间[0,4]数值数组的创建和寻23210 x一维数组的创1%%%==1234560=ColumnsColumnsthroughthrough0b1= b21%%一维数组的创1%%%==1234560=ColumnsColumnsthroughthrough0b1= b21%%c1=[2c1=Column3.0000%through%%c2a=2.7358;%%C3.5000+5.4716+23.2-3】试用变量编辑器，把如下(36)的数组输入MATLAB内存，并命名为A1833.2-13M文件3.2-2利用M4利函数创3.2-13M文件3.2-2利用M4利函数创建数3.2-5】ans=11%111111%%ansD100%0100014ans111%%ans100010001%%ans131 ans111%%ans100010001%%ans131 13through 1311313111333311131313【例3.2-6】%%%%A001200340056007800900Aans8ans8%%%%ans12125656ans%5%ans789%%%A1345789)%B159%%ans789%%%A1345789)%B159%%%%L11A0001000100345789a=%%%12345678A123412567834A5678b=14B10%%04D110041004100410046100410]0410041004D1=044004400440044【例3.2-8】A=123456789B100410]0410041004D1=044004400440044【例3.2-8】A=123456789B321%654987C789%456123%D987654321数组运vr=[0.89,1.20,%9.70,ir=[0.028,0.040,%fork=1:L%fork=1:Lrm=%0.118,0.258,0.299,0.257,%%%%%%vr=[0.89,1.20,4.27,3.62,7.71,8.99,7.92,9.70,%7r=ColumnsColumnsrm=0.145,%0.258,0.299,0.257,r=ColumnsColumnsrm=0.145,%0.258,0.299,0.257,%%throughthroughforii=1:Mfor%%%%%%%%%ans=0%83.3-1“非数”和“空3.3-1“非数”和“空”数a=bc=ans=ans=ans=ans=1%%%%9R%LR00%00010010% disp('非数在二维数组R中的位置R%LR00%00010010% disp('非数在二维数组R中的位置%disp(['单下标时的第',int2str(si(1)),'和第',int2str(si(2)),'个元素si=9=21=352135==非数在二维数组R中的位置ab==cd==f=ans=ans1ans%%%1whichaaisavariable.ans=2ans%%%00A012345%A0145关系操作和逻辑1whichaaisavariable.ans=2ans%%%00A012345%A0145关系操作和逻辑操【例3.5-1】A%%123456789B987654321r0111000000r1000010000%%%%3.5-1A%%00123L1L211111000000100AB0C1%0012310103.5-1A%%00123L1L211111000000100AB0C1%00123101010000%%%%%y残缺图 正确图 3.5-2要求在闭区间[0,210个等距采样点的一维数组。试用两种不同的指令实的位置，分别求出它们的“全下标”和“单下标”。（概率双位（即取-1，+1）取值的随机码”程序指令，并给出-1码的数目。（提示：rand, 4A3.5-2要求在闭区间[0,210个等距采样点的一维数组。试用两种不同的指令实的位置，分别求出它们的“全下标”和“单下标”。（概率双位（即取-1，+1）取值的随机码”程序指令，并给出-1码的数目。（提示：rand, 4AB1=A.^(0.5),B2=A^(0.5),34 得结果不同。（1）请分别写出根据B1,B2恢复原矩阵A的程序。（2）用指令检验所先运行clear,formatlong,rand('twister',1),A=rand(3,3)，然后根据A对角阵B，其相应元素由A的对角元素构成；另一个矩阵C0，而其余元素与对应的A阵元素相同。（提示：diag）7．先运行指令y10 10 10 10 tZ=sin(X).*sin(Y)./X./Y;产生矩阵Z。（1）请问Z中有多少个“非数”数据？（Z=sin(X).*sin(Y)./X./Y;产生矩阵Z。（1）请问Z中有多少个“非数”数据？（2）用指令surf(X,Y,Z);shadinginterp观察所绘的图形。（3）请写出绘制相应的“无裂缝”图形的全部指令。（提示：isnan,sum,eps）9k9kkk 456 Sak行。例如k3时，A阵为 的和是一个(110)行数组615 87，并把它保存为Sa的第3行。问题：该段（sum对行数组的作用而设计。仔细观察以下程序运行后所得到Sa正确吗？fork=10:-1:1;A=reshape(1:10*k,k,10);Sa(k,:)=sum(A);end;Sa）4数值计4数值计L1=L2=%%symsholdonholdoff%%1holdonholdoff%4.1-2【例4.1-3】%21 0holdonholdoff%4.1-2【例4.1-3】%21 0 t10 thold%%title('[0,2\pi]')xlabel('t'),boxoffholdoffholdontitle('[end-10,end]')hold%%title('[0,2\pi]')xlabel('t'),boxoffholdoffholdontitle('[end-10,end]')4.1-3diff和gradient【例4.1-4】%%%%3[0, [end-10,11 dxdtdif 6.