软件实习-第一章matlab编程基础资料教程文件

软件实习--第一章MATLAB编程基础资料考核成绩分为：“字母制”第一章MATLAB编程基础第二章MATLAB数值计算第三章MATLAB绘图第四章MATLAB在电子科学与技术专业中的应用教学内容教学目的：

通过本课程的学习，使同学们初步掌握MATLAB的基本知识，学会运用MATLAB语言进行数值计算和计算机仿真的基本技能，为后续的课程、毕业设计和从事生产及科学研究打下较好的基础。

第一章MATLAB编程基础简介1、

MATLAB(MATrixLABoratory)是当前最流行、功能最强大的应用软件和编程语言之一，它集科学计算、可视化和可编程等功能集于一体，而且具有易学易用的特点，可利用计算机的运算能力来处理设计者手工处理所遇到的不便与繁琐。

MATLAB语言是由美国的CleverMoler博士于1980年开发的,设计者的初衷是为解决“线性代数”课程的矩阵运算问题.

MATLAB已经发展成为多学科、多种工作平台的功能强大的大型软件。在各高校，MATLAB已经成为光电子科学与技术等专业的基本教学工具，是应该掌握的基本技能之一。2MATLAB的主要特点是：（1）有高性能数值计算的高级算法，特别适合矩阵代数领域；（2）有大量事先定义的数学函数，并且有很强的用户自定义函数的能力；（3）有强大的绘图功能以及具有教育、科学和艺术学的图解和可视化的二维、三维图；（4）基于HTML的完整的帮助功能；（5）适合个人应用的强有力的面向矩阵(向量)的高级程序设计语言；（6）与其它语言编写的程序结合和输入输出格式化数据的能力；（7）有在多个应用领域解决难题的工具箱。

3、计算机仿真的概念仿真（Simulation）的定义：所谓仿真就是模型实验，其基本思想是利用物理的或数学的模型来类比模仿现实过程，以寻求对真实过程的认识。它所遵循的基本原则是相似性原理。计算机仿真的定义：计算机仿真是基于所建立的系统仿真模型，利用计算机对系统进行分析与研究的方法。计算机仿真的三要素1、三要素（1）系统：研究的对象（2）模型：系统的抽象（3）计算机：工具与手段

2、基本步骤包括三个基本的内容：建模数学模型计算机模型仿真实验结果分析3应用仿真技术的意义1、经济大型、复杂系统直接实验是十分昂贵的，如：空间飞行器的一次飞行实验的成本约在1亿美元左右，而采用仿真实验仅需其成本的1/10~1/5，而且设备可以重复使用。2、安全某些系统（如载人飞行器、核电装置等），直接实验往往会有很大的危险，甚至是不允许的，而采用仿真实验可以有效降低危险程度，对系统的研究起到保障作用。3、快捷提高设计效率：比如电路设计，服装设计等等。4、具有优化设计和预测的特殊功能对一些真实系统进行结构和参数的优化设计是非常困难的，这时仿真可以发挥它特殊的优化设计功能。在非工程系统中（如社会、管理、经济等系统），由于其规模及复杂程度巨大，直接实验几乎不可能，这时通过仿真技术的应用可以获得对系统的某种超前认识。

基于仿真的方法的优点：可进行任意精确地建模，能够比前述两种方法在更大的范围内对系统性能进行评价。可以容易地合并数学和实验模型，可以将器件的实测特性和实际信号用于分析和设计。仿真得到的波形还可用作测试信号对硬件的性能进行验证。成本低，速度快，灵活性高。

缺点是：在某些情况下，计算复杂度很高（合理地选择建模和仿真方法可降低复杂度）了解MATLAB的基本知识熟悉MATLAB的上机环境掌握利用MATLAB进行基本运算的方法初步具备将一般数学问题转化成对应的计算机模型并进行处理的能力第一章MATLAB语言基础工作区（命令）窗口命令历史窗口发射台窗口当前工作目录当前目录浏览器工作空间浏览器MATLAB窗口环境当MATLAB启动时，展现在屏幕上的界面为MATLAB的默认界面，如图所示。

MATLAB主窗口是MATLAB的主要工作界面。主窗口除了嵌入一些子窗口外，还主要包括菜单栏和工具栏。

1．菜单栏

在MATLAB6.5主窗口的菜单栏，共包含File、Edit、View、Web、Window和Help6个菜单项。(1)File菜单项：File菜单项实现有关文件的操作。

(2)Edit菜单项：Edit菜单项用于命令窗口的编辑操作。

(3)View菜单项：View菜单项用于设置MATLAB集成环境的显示方式。(4)Web菜单项：Web菜单项用于设置MATLAB的Web操作。

(5)Window菜单项：主窗口菜单栏上的Window菜单，只包含一个子菜单Closeall，用于关闭所有打开的编辑器窗口，包括M-file、Figure、Model和GUI窗口。

(6)Help菜单项：Help菜单项用于提供帮助信息。2．工具栏

MATLAB主窗口的工具栏共提供了命令按钮。这些命令按钮均有对应的菜单命令，但比菜单命令使用起来更快捷、方便。

子窗口命令窗口是MATLAB的主要交互窗口，用于输入命令并显示除图形以外的所有执行结果。

MATLAB命令窗口中的“>>”为命令提示符，表示MATLAB正在处于准备状态。在命令提示符后键入命令并按下回车键后，MATLAB就会解释执行所输入的命令，并在命令后面给出计算结果。命令历史记录窗口

在默认设置下，历史记录窗口中会自动保留自安装起所有用过的命令的历史记录，并且还标明了使用时间，从而方便用户查询。而且，通过双击命令可进行历史命令的再运行。如果要清除这些历史记录，可以选择Edit菜单中的ClearCommandHistory命令。帮助窗口

进入帮助窗口可以通过以下3种方法：

(1)单击MATLAB主窗口工具栏中的Help按钮。

(2)在命令窗口中输入helpwin、helpdesk或doc。

(3)选择Help菜单中的“MATLABHelp”选项。二、MATLAB命令窗口1、启动MATLAB命令窗口计算机安装好MATLAB之后，双击MATLAB图标，就可以进入命令窗口，此时意味着系统处于准备接受命令的状态，可以在命令窗口中直接输入命令语句。MATLAB语句形式》变量＝表达式； 通过“=”符号将表达式的值赋予变量。当键入回车键时，该语句被执行。语句执行之后，窗口自动显示出语句执行的结果。如果希望结果不被显示，则只要在语句之后加上一个分号（；）即可。此时尽管结果没有显示，但它依然被赋值并在MATLAB工作空间中分配了内存。三、变量和数值显示格式

1、变量（1）变量的命名：变量的名字必须以字母开头（不能超过19个字符），之后可以是任意字母、数字或下划线；变量名称区分字母的大小写；变量中不能包含有标点符号。（2）一些特殊的变量ans：用于结果的缺省变量名i、j：虚数单位pi：圆周率nargin：函数的输入变量个数eps：计算机的最小数nargout：函数的输出变量个数inf：无穷大realmin：最小正实数realmax：最大正实数nan：不定量flops：浮点运算数（3）变量操作在命令窗口中，同时存储着输入的命令和创建的所有变量值，它们可以在任何需要的时候被调用。如要察看变量a的值，只需要在命令窗口中输入变量的名称即可：》a2、数值显示格式任何MATLAB的语句的执行结果都可以在屏幕上显示，同时赋值给指定的变量，没有指定变量时，赋值给一个特殊的变量ans，数据的显示格式由format命令控制。format只是影响结果的显示，不影响其计算与存储；MATLAB总是以双字长浮点数（双精度）来执行所有的运算。如果结果为整数，则显示没有小数；如果结果不是整数，则输出形式有：formatshort：短格式（5位定点数）99.1253formatlong：长格式（15位定点数99.12345678900000formatshorte：短格式e方式9.9123e+001formatlonge：长格式e方式9.912345678900000e+001formatbank：2位十进制99.12formathex：十六进制格式例:formatshort,pians=3.1416MATLAB的基本运算可分为三类：算术运算关系运算逻辑运算四、简单的数学运算1、常用的数学运算符＋，—，*（乘），/（左除），\（右除），^（幂）在运算式中，MATLAB通常不需要考虑空格；多条命令可以放在一行中，它们之间需要用分号隔开；逗号告诉MATLAB显示结果，而分号则禁止结果显示。2、常用数学函数abs,sin,cos,tan,asin,acos,atan,sqrt,exp,imag,real,sign,log,log10,conj（共扼复数）等五、MATLAB的工作空间1、MATLAB的工作空间包含了一组可以在命令窗口中调整（调用）的参数who：显示当前工作空间中所有变量的一个简单列表 whos：则列出变量的大小、数据格式等详细信息clear：清除工作空间中所有的变量clear变量名：清除指定的变量2、保存和载入workspace（1）savefilenamevariables将变量列表variables所列出的变量保存到磁盘文件filename中Variables所表示的变量列表中，不能用逗号，各个不同的变量之间只能用空格来分隔。未列出variables时，表示将当前工作空间中所有变量都保持到磁盘文件中。（2）loadfilenamevariables将以前用save命令保存的变量variables从磁盘文件中调入MATLAB工作空间。用load命令调入的变量，其名称为用save命令保存时的名称，取值也一样。Variables所表示的变量列表中，不能用逗号，各个不同的变量之间只能用空格来分隔。未列出variables时，表示将磁盘文件中的所有变量都调入工作空间。3、退出工作空间quit或exit六、文件管理文件管理的命令，包括列文件名、显示或删除文件、显示或改变当前目录等。（what、dir、type、delete、cd、which）what：显示当前目录下所有与matlab相关的文件及它们的路径。dir：显示当前目录下所有的文件which：显示某个文件的路径cdpath：由当前目录进入path目录cd..：返回上一级目录cd：显示当前目录typefilename：在命令窗口中显示文件filenamedeletefilename：删除文件filenameMATLAB语言编程简单，使用方便。例如：clc%清屏clear%从内存中清除变量和函数moreon％moreoff：不允许分页moreon：允许分页。echoon％echoon/off：控制是否在屏幕上显示程序内容。%求矩阵与矩阵的乘积，矩阵与向量的乘积A=[567;946;436]B=[345;579;731]X=[578]'C=A*BY=A*Xmoreoffechooff数值矩阵的生成不管是任何矩阵（向量），我们可以直接按行方式输入每个元素：同一行中的元素用逗号（，）或者用空格符来分隔，且空格个数不限；不同的行用分号（；）分隔。所有元素处于一方括号（[]）内；矩阵生成不但可以使用纯数字（含复数），也可以使用变量（或者说采用一个表达式）。B=[123；234；345]》y=[2,4,5；368]y=245368》a=1;b=2;c=3;》x=[5bc;a*ba+cc/b]x=5.0002.0003.0002.0004.0001.500◆直接输入法创建简单矩阵。MATLAB提供了丰富的绘图功能helpgraph2d可得到所有画二维图形的命令helpgraph3d可得到所有画三维图形的命令加、减运算运算符：“＋”和“－”分别为加、减运算符。运算规则：对应元素相加、减，即按线性代数中矩阵的“十”，“一”运算进行。例>>A=[1,1,1;1,2,3;1,3,6]>>B=[8,1,6;3,5,7;4,9,2]>>C=A+B>>D=A-B结果显示：C=92747105128D=-70-5-2-3-4-3-64>>A=[1,2;3,4;5,6];B=[1,2,1;3,4,2];A*Bans=7105152211233417乘法运算符：*运算规则：按线性代数中矩阵乘法运算进行，即放在前面的矩阵的各行元素，分别与放在后面的矩阵的各列元素对应相乘并相加。1．两个矩阵相乘>>A=3;>>B=[2,3,4,5,6];>>A.*Bans=69121518点运算是两个维数相同矩阵对应元素之间的运算，在有的教材中也定义为数组运算。数组乘法：A.*B表示A与B对应元素相乘。A=[1,2;3,4;5,6];B=[1,2;1,3;4,2];c=A.*Bc=1431220123．向量点积函数dot格式C=dot(A,B)%若A、B为向量，则返回向量A与B的点积，A与B长度相同；若为矩阵，则A与B有相同的维数。C=dot(A,B,dim)%在dim维数中给出A与B的点积例 >>X=[-102]；>>Y=[-2-11];>>Z=dot(X,Y)则显示：Z=44．向量叉乘在数学上，两向量的叉乘是一个过两相交向量的交点且垂直于两向量所在平面的向量。在Matlab中，用函数cross实现。函数cross格式C=cross(A,B)%若A、B为向量，则返回A与B的叉乘，即C=A×B，A、B必须是3个元素的向量；若A、B为矩阵，则返回一个3×n矩阵，其中的列是A与B对应列的叉积，A、B都是3×n矩阵。C=cross(A,B,dim)%在dim维数中给出向量A与B的叉积。A和B必须具有相同的维数，size(A,dim)和size(B,dim)必须是3。例1-24计算垂直于向量(1,2,3)和(4,5,6)的向量。

>>a=[123];>>b=[456];>>c=cross(a,b)结果显示：

c=-36-35．混合积混合积由以上两函数实现：例1-25计算向量a=(1,2,3)、b=(4,5,6)和c=(-3,6,-3)的混合积解：>>a=[123];b=[456];c=[-36-3];>>x=dot(a,cross(b,c))结果显示：x=54注意：先叉乘后点乘，顺序不可颠倒。a\b运算等效于求a*x=b的解如：a=[12;34]；b=[35;59]》a\b=[-1-1;23]Matlab提供了两种除法运算：左除（\）和右除（/）。一般情况下，x=a\b是方程a*x=b的解，而x=b/a是方程x*a=b的解。而a/b等效于求x*b=a的解。只有方阵才可以求幂。如：a=[12;34]；b=[35;59]a/b=[-0.500.50;3.50–1.50]》a.\b=[3.002.50;1.672.25]如：a=[12;34]；b=[35;59]》a./b=[0.330.40;0.600.44]如：a=[12;34]；b=[35;59]矩阵乘方运算符：^运算规则：（1）当A为方阵，P为大于0的整数时，A^P表示A的P次方，即A自乘P次；P为小于0的整数时，A^P表示A-1的P次方。(2)标量的数组乘方P.^A，标量的数组乘方定义为数组乘方：A.^P：表示A的每个元素的P次乘方。>>A=[568;7910];B=[1247;4155];A<Bans=100010>>A=[568;7910];>>B=[1247;4155];>>A>Bans=0111017．矩阵逻辑运算设矩阵A和B都是m×n矩阵或其中之一为标量，在MATLAB中定义了如下的逻辑运算：（1）矩阵的与运算格式A&B或and(A,B)

说明A与B对应元素进行与运算，若两个数均非0，则结果元素的值为1，否则为0。（2）或运算格式A|B或or(A,B)

说明A与B对应元素进行或运算，若两个数均为0，则结果元素的值为0，否则为1。（3）非运算格式~A或not(A)

说明若A的元素为0，则结果元素为1，否则为0。（4）异或运算格式xor(A,B)

说明A与B对应元素进行异或运算，若相应的两个数中一个为0，一个非0，则结果为0，否则为1。例>>A=[0234;1350],B=[1053;1505]A=02341350B=10531505>>C1=A&B,C2=A|B,C3=~A,C4=xor(A,B)C1=00111100C2=11111111C3=10000001C4=11000011A=[816;357;492]矩阵转置运算符：′运算规则：若矩阵A的元素为实数，则与线性代数中矩阵的转置相同。若A为复数矩阵，则A转置后的元素由A对应元素的共轭复数构成。若仅希望转置，则用如下命令：A.′。》a=[123;456]'a=142536》b=[1+2i2-7i]'b=1.0000-2.0000i2.0000+7.0000i》b=[1+2i2-7i].'b=1.0000+2.0000i2.0000-7.0000i》a=[123;456].'a=142536方阵的行列式函数det格式d=det(X)%返回方阵X的多项式的值例>>A=[123;456;789]A=123456789>>D=det(A)D=0矩阵与矩阵操作基础x=linspace(0,4*pi,100)%将[0,4π]区间100等分，产生了一个100维向量

x=Columns1through600.12690.25390.38080.50770.6347Columns7through120.76160.88851.01551.14241.26931.3963Columns13through181.52321.65011.77711.90402.03092.1579Columns19through242.28482.41172.53872.66562.79252.9195Columns25through303.04643.17333.30033.42723.55413.6811Columns31through363.80803.93494.06194.18884.31574.4427Columns37through424.56964.69654.82354.95045.07735.2043Columns43through485.33125.45815.58515.71205.83895.9659Columns49through546.09286.21976.34676.47366.60056.7275Columns55through606.85446.98137.10837.23527.36217.4890Columns61through667.61607.74297.86987.99688.12378.2506Columns67through728.37768.50458.63148.75848.88539.0122

Columns73through789.13929.26619.39309.52009.64699.7738Columns79through849.900810.027710.154610.281610.408510.5354Columns85through9010.662410.789310.916211.043211.170111.2970Columns91through9611.424011.550911.677811.804811.931712.0586Columns97through10012.185612.312512.439412.5664

logspace(1,3,10)ans=1.0e+003*Columns1through60.01000.01670.02780.04640.07740.1292Columns7through100.21540.35940.59951.0000例产生一个3×4随机矩阵>>R=rand(3,4)R=0.95010.48600.45650.44470.23110.89130.01850.61540.60680.76210.82140.7919例产生一个在区间[10,20]内均匀分布的4阶随机矩阵>>a=10;b=20;>>x=a+(b-a)*rand(4)x=19.218119.354710.578911.388917.382119.169013.528712.027711.762714.102718.131711.987214.057118.936510.098616.0379命令逆函数inv格式Y=inv(X)%求方阵X的逆矩阵。若X为奇异阵或近似奇异阵，将给出警告信息。例1-43求的逆矩阵方法一>>A=[123;221;343];>>Y=inv(A)或Y=A^(-1)则结果显示为

Y=1.00003.0000-2.0000-1.5000-3.00002.50001.00001.0000-1.0000逆矩阵与行列式计算求逆：inv(A)；求行列式：det(A)要求矩阵必须为方阵》a=[123;456;235];》b=inv(a)b=-2.33330.33331.00002.66670.3333-2.0000-0.6667-0.33331.0000》det(a)ans=-3矩阵的迹函数trace格式b=trace(A)%返回矩阵A的迹，即A的对角线元素之和。矩阵的秩函数rank格式k=rank(A)%求矩阵A的秩k=rank(A,tol)%tol为给定误差矩阵对角线元素的抽取函数diag格式X=diag(v,k)%以向量v的元素作为矩阵X的第k条对角线元素，当k=0时，v为X的主对角线；当k>0时，v为上方第k条对角线；当k<0时，v为下方第k条对角线。X=diag(v)%以v为主对角线元素，其余元素为0构成X。v=diag(X,k)%抽取X的第k条对角线元素构成向量v。k=0：抽取主对角线元素；k>0：抽取上方第k条对角线元素；k<0抽取下方第k条对角线元素。v=diag(X)%抽取主对角线元素构成向量v。例>>v=[123];>>x=diag(v,-1)x=0000100002000030>>A=[123;456;789]A=123456789>>v=diag(A,1)v=262．上三角阵和下三角阵的抽取函数tril%取下三角部分格式L=tril(X)%抽取X的主对角线的下三角部分构成矩阵LL=tril(X,k)%抽取X的第k条对角线的下三角部分；k=0为主对角线；k>0为主对角线以上；k<0为主对角线以下。函数triu%取上三角部分格式U=triu(X)%抽取X的主对角线的上三角部分构成矩阵UU=triu(X,k)%抽取X的第k条对角线的上三角部分；k=0为主对角线；k>0为主对角线以上；k<0为主对角线以下。例>>A=ones(4)%产生4阶全1阵A=1111111111111111>>L=tril(A,1)%取下三角部分L=1100111011111111>>U=triu(A,-1)%取上三角部分U=1111111101110011对于小数构成的矩阵A来说，如果我们想对它取整数，有以下几种方法：（1）按-∞方向取整函数floor格式floor(A)%将A中元素按-∞方向取整，即取不足整数。（2）按+∞方向取整函数ceil格式ceil(A)%将A中元素按+∞方向取整，即取过剩整数。（3）四舍五入取整函数

round格式round(A)%将A中元素按最近的整数取整，即四舍五入取整。（4）按离0近的方向取整函数

fix格式fix(A)%将A中元素按离0近的方向取整例1-55>>A=-1.5+4*rand(3)A=2.30050.44390.3259-0.57542.0652-1.42600.92741.54841.7856>>B1=floor(A),B2=ceil(A),B3=round(A),B4=fix(A)B1=200-12-2011B2=31103-1122B3=200-12-1122B4=20002-10114、了解矩阵超越函数在MATLAB中exp、sqrt等命令也可以作用到矩阵上，但这种运算是定义在矩阵的单个元素上的，即分别对矩阵的每一个元素进行计算。超越数学函数可以在函数后加上m而成为矩阵的超越函数，例如：expm,sqrtm。矩阵的超越函数要求运算矩阵为方阵。函数

sin、sinh功能正弦函数与双曲正弦函数格式Y=sin(X)%计算参量X（可以是向量、矩阵，元素可以是复数）中每一个角度分量的正弦值Y，所有分量的角度单位为弧度。Y=sinh(X)%计算参量X的双曲正弦值Y注意：sin(pi)并不是零，而是与浮点精度有关的无穷小量eps，因为pi仅仅是精确值π浮点近似的表示值而已；对于复数Z=x+iy，函数的定义为：sin(x+iy)=sin(x)*cos(y)+i*cos(x)*sin(y)，，，

例2-1x=-pi:0.01:pi;plot(x,sin(x))x=-5:0.01:5;plot(x,sinh(x))图形结果为图2-1。图2-1正弦函数与双曲正弦函数图

图2-1正弦函数与双曲正弦函数图函数asin、asinh功能反正弦函数与反双曲正弦函数格式Y=asin(X)%返回参量X（可以是向量、矩阵）中每一个元素的反正弦函数值Y。若X中有的分量处于[-1,1]之间，则Y=asin(X)对应的分量处于[-π/2,π/2]之间，若X中有分量在区间[-1,1]之外，则Y=asin(X)对应的分量为复数。Y=asinh(X)%返回参量X中每一个元素的反双曲正弦函数值Y例2-2x=-1:.01:1;plot(x,asin(x))x=-5:.01:5;plot(x,asinh(x))图形结果为图2-2。

图2-2反正弦函数与反双曲正弦函数图函数

cos、cosh功能余弦函数与双曲余弦函数格式Y=cos(X)%计算参量X（可以是向量、矩阵，元素可以是复数）中每一个角度分量的余弦值Y，所有角度分量的单位为弧度。我们要指出的是，cos(pi/2)并不是精确的零，而是与浮点精度有关的无穷小量eps，因为pi仅仅是精确值π浮点近似的表示值而已。Y=conh(X)%计算参量X的双曲余弦值Y说明若X为复数z=x+iy，则函数定义为：cos(x+iy)=cos(x)*cos(y)+i*sin(x)*sin(y)，，，

例2-3x=-pi:0.01:pi;plot(x,cos(x))x=-5:0.01:5;plot(x,cosh(x))图形结果为图2-3。图2-3余弦函数与双曲余弦函数图

图2-3余弦函数与双曲余弦函数图函数acos、acosh功能反余弦函数与反双曲余弦函数格式Y=acos(X)%返回参量X（可以是向量、矩阵）中每一个元素的反余弦函数值Y。若X中有的分量处于[-1,1]之间，则Y=acos(X)对应的分量处于[0,π]之间，若X中有分量在区间[-1,1]之外，则Y=acos(X)对应的分量为复数。Y=asinh(X)%返回参量X中每一个元素的反双曲余弦函数Y例2-4x=-1:.01:1;plot(x,acos(x))x=-5:.01:5;plot(x,acosh(x))图形结果为图2-4。

图2-4反余弦函数与反双曲余弦函数图函数

tan、tanh功能正切函数与双曲正切函数格式Y=tan(X)%计算参量X（可以是向量、矩阵，元素可以是复数）中每一个角度分量的正切值Y，所有角度分量的单位为弧度。我们要指出的是，tan(pi/2)并不是精确的零，而是与浮点精度有关的无穷小量eps，因为pi仅仅是精确值π浮点近似的表示值而已。Y=tanh(X)%返回参量X中每一个元素的双曲正切函数值Y例2-5x=(-pi/2)+0.01:0.01:(pi/2)-0.01;%稍微缩小定义域plot(x,tan(x))x=-5:0.01:5;plot(x,tanh(x))图形结果为图2-5。

图2-5正切函数与双曲正切函数图

函数atan、atanh功能反正切函数与反双曲正切函数格式Y=atan(X)%返回参量X（可以是向量、矩阵）中每一个元素的反正切函数值Y。若X中有的分量为实数，则Y=atan(X)对应的分量处于[-π/2,π/2]之间。Y=atanh(X)%返回参量X中每一个元素的反双曲正切函数值Y。

图2-6反正切函数与反双曲正切函数图例2-6x=-20:0.01:20;plot(x,atan(x))x=-0.99:0.01:0.99;plot(x,atanh(x))图形结果为图2-6。函数

cot、coth功能余切函数与双曲余切函数格式Y=cot(X)%计算参量X（可以是向量、矩阵，元素可以是复数）中每一个角度分量的余切值Y，所有角度分量的单位为弧度。Y=coth(X)%返回参量X中每一个元素的双曲余切函数值Y例2-7x1=-pi+0.01:0.01:-0.01;%去掉奇点x=0x2=0.01:0.01:pi-0.01;%做法同上plot(x1,cot(x1),x2,cot(x2))plot(x1,coth(x1),x2,coth(x2))图形结果为图2-7。

图2-7余切函数与双曲余切函数图函数

acot、acoth功能反余切函数与反双曲余切函数格式Y=acot(X)%返回参量X（可以是向量、矩阵）中每一个元素的反余切函数YY=acoth(X)%返回参量X中每一个元素的反双曲余切函数值Y例2-8x1=-2*pi:pi/30:-0.1;x2=0.1:pi/30:2*pi;%去掉奇异点x=0plot(x1,acot(x1),x2,acot(x2))x1=-30:0.1:-1.1;x2=1.1:0.1:30;plot(x1,acoth(x1),x2,acoth(x2))图形结果为图2-8。