matlab教程(全)ppt课件

MATLAB语言及其应用授课教师:刘晨晨单位：网络工程与信息安全系email:lcczzm@126.comApplicationofMatlabLanguage.课程的作用课程的目的课程的特点课程安排.课程的作用在欧美各高等学校，Matlab成为线性代数、自动控制理论、数字信号处理、时间序列分析、动态系统仿真、图像处理等诸多课程的基本教学工具，成为本科生、硕士生和博士生的必须掌握的基本技能。在 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 研究单位和工业部门，Matlab已被广泛地用于研究和解决各种具体的工程问题。可以预见，Matlab将在我国科学研究和工程应用中发挥越来越大的作用。.讲授MATLAB语言基础入门知识，介绍MATLAB产品的体系、MATLAB桌面工具的使用方法，重点介绍MATLAB的数据可视化、数值计算的基本步骤以及如何使用MATLAB语言编写整洁、高效、 规范 编程规范下载gsp规范下载钢格栅规范下载警徽规范下载建设厅规范下载 的程序。并涉及到一些具体的专业应用工具箱（如：信号处理工具箱、图像处理工具箱等）。通过本课程的学习，了解、熟悉、掌握MATLAB的基本编程方法，并具有初步的利用计算机处理、解决实际问题的能力，为进一步学习后续的专业课程做好准备。本课程的目的（ObjectivesofThisCourse）.课程安排课堂教学：共24学时；（1-12周）上机试验：共24学时。（2-13周，周二7-8节，九实401、402、403）学习成绩：1)上机实验成绩占30%;2)考勤10%;3)考试60%（随堂考试）。主要参考书《精通MATLAB6.5》张志涌等编著，北航出版，2003年《高等应用数学问题的Matlab求解》薛定宇等著，清华大学出版社，2004年《Matlab程序设计与应用》刘卫国主编，高等教育出版社.授课宗旨讲授MATLAB的通用功能。寓教于例，由浅入深。关于科学计算，着重强调理论概念、算法和实际计算三者之间的关系。.第一讲Matlab概述前言Matlab软件概述Matlab的桌面环境及入门知识.1Matlab概述内容Matlab发展历史Matlab产品家族(Matlabfamilyofproducts)体系Matlab语言的特点。目的全面了解Matlab软件包激发对Matlab软件的学习兴趣。.1.1MATLAB的历史及影响70年代中期，CleveMoler博土及其同事在美国国家基金会的帮助下，开发了LINPACK和EISPACK的FORTRAN语言子程序库，这两个程序库代表了当时矩阵运算的最高水平。到了70年代后期，身为美国新墨西哥州大学计算机系系主任的CIeveMoler，在给学生上线性代数课时，为了让学生能使用这两个子程序库，同时又不用在编程上花费过多的时间，开始着手用FORTRAN语言为学生编写使用LINPACK和EISPACK的接口程序，他将这个程序取名为MATLAB，其名称是由MATrix和LABoratory（矩阵实验室）两个单词的前三个字母所合成。在1978年，Malab就面世了。这个程序获得了很大的成功，受到了学生的广泛欢迎。在以后的几年里，Matlab在多所大学里作为教学辅助软件使用，并作为面向大众的免费软件广为流传。CleveMoler.将MATLAB商品化的不是CleveMoler,而是一个名叫JackLittle的人。当免费的MATLAB软件到Stanford大学，JackLittle正在该校主修控制，便接触到了当时MATLAB，直觉告诉他，这是一个具有巨大发展潜力的软件。因此他在毕业沒多久，就开始用C语言重新编写了MATLAB的核心。在Moler的协助下，于1984年成立MathWorks公司，首次推出MATLAB商用版。在其商用版推出的初期，MATLAB就以其优秀的品质（高效的数据计算能力和开放的体系结构）占据了大部分数学计算软件的市场，原来应用于控制领域里的一些封闭式数学计算软件包（如英国的UMIST、瑞典的LUND和SIMNON、德国的KEDDC）就纷纷被淘汰或在MATLAB上重建。JackLittle.CleveMoler至今仍是该公司的首席科学家，他以60多岁的高龄，还常常亲自进行撰写程序的工作，非常令人佩服。如果你有数值运算方面的高水平问题，寄到MathWorks后，大部份还是会由CleveMoler亲自回答。在1994年，Pentium芯片曾发生Fdiv的bug，当时CleveMoler是第一个以软件方式解决此bug的人，曾一时脍炙人口。在公司初创的五年，JackLittle非常辛苦，常常身兼数职（董事长、总经理、推销、程序开发等），但公司一直稳定发展，从当初的一人公司，到1993年的200人，到2000年的500余人，到2005年公司员工达到了1300人，不但打败其他竞争软件，而且前景一片欣欣向荣。根据JackLittle个人说法，MATLAB早期成功的两大因素是：选用了C语言及选定PC为主要平台，这似乎和微软的成功有相互呼应之妙。MathWorks公司，目前仍然是私人企业，并未上市，这和JackLittle个人理念有关，他认为MATLAB的设计方向应该一直是以顾客的需求与软件的完整性为首要目标，而不是以盈利为主要目的，因此MATLAB一直是在稳定中求进步，而不会因为上市而遭受股东左右其发展方向。这也是为什么MATLAB新版本总是姗姗来迟的原因，因为他们不会因为市场的需求而推出不成熟的产品。此外，由于JackLittle保守的个性，也使得MathWorks不曾跨足MATLAB/Simulink以外的行业，当前商场上纷纷扰扰的并购或分家，MathWorks完全是绝缘体。.1992年，支持Windows3.x的MATLAB4.0版本推出，增加了Simulink，Control，NeuralNetwork，SignalProcessing等专用工具箱。1993年11月，MathWorks公司推出了Matlab4.1，其中主要增加了符号运算功能。当升级至Matlab4.2c，这一功能在用户中得到广泛应用。1997年，Matlab5.0版本问世了,实现了真正的32位运算，加快数值计算，图形表现有效。2001年初，MathWorks公司推出了Matlab6.0（R12）。2002年7月，推出了Matlab6.5(R13)，在这一版本中Simulink升级到了5.0，性能有了很大提高，另一大特点是推出了JIT程序加速器，Matlab的计算速度有了明显的提高。2005年9月，推出了MAILAB7.1(Release14SP3)，在这一版本中Simulink升级到了6.3，软件性能有了新的提高，用户界面更加友好。值得 说明 关于失联党员情况说明岗位说明总经理岗位说明书会计岗位说明书行政主管岗位说明书 的是，MatlabV7.1版采用了更先进的数学程序库，即“LAPACK”和“BLAS”。目前，Matlab软件支持多种系统平台，如常见的WindowsNT/XP、UNIX、Linux等。Matlab版本的发展.MATLAB就是这样经过了近30年的专门打造、20多年的千锤百炼，它以高性能的数组运算（包括矩阵运算）为基础，不仅实现了大多数数学算法的高效运行函数和数据可视化，而且提供了非常高效的计算机高级编程语言，在用户可参与的情况下，各种专业领域的工具箱不断开发和完善，MATLAB取得了巨大的成功，已广泛应用于科学研究、工程应用，用于数值计算分析、系统建模与仿真。早在20世纪90年代初，欧美等发达国家的大学就将MATLAB列为一种必须掌握的编程语言。近几年来，国内的很多大学也将MATLAB列为了本科生必修课程。与Maple、Mathematica数学计算软件相比，MATLAB以数值计算见长，而Maple等以符号运算见长，能给出解析解和任意精度解，而处理大量数据的能力远不如MATLAB。MATLAB软件功能之强大、应用之广泛，已成为为21世纪最为重要的科学计算语言。可见学习掌握这一工具的重要性。.MATLAB产品由若干模块组成，不同的模块完成不同的功能，其中有：实际上MATLAB本身就是一个极其丰富的资源库，那么应该从哪一部分开始着手、学习使用MATLAB呢？这就有必要了解这一软件产品的体系结构。MATLABCompiler这种编译器可以将MATLAB程序文件编译生成标准的C/C++语言文件，而生成的标准的C/C++文件可以被任何一种C/C++编译器编译生成函数库或可执行文件，以提高程序的运行效率。1.2MATLAB产品的体系结构MATLABMATLABToolboxesMATLABCompilerSimulinkSimulinkBlocksetsReal-TimeWorkshop(RTW)StateflowStateflowCoder由这些模块产品之间的关系可以图1.1表示。图1.1MATLAB产品体系结构围绕着MATLAB这个计算核心，形成了诸多针对不同应用领域的算法程序包，被称为专用工具箱（Toolbox），这些工具箱的列表以及每个工具箱的使用详见MATLAB在线帮助文档。MATLAB本身所提供的工具箱大概有40多个，另外还有其他公司或研究单位开发提供的工具箱，这些工具箱的总数已有100多个，而且新的工具箱还在不断增加。如果你有特别的应用领域，可以首先到网上查找是否已有相关的工具箱，很可能已有人将你要做的应用程序作成工具箱了。Simulink是窗口图形方式的、专门用于连续时间或离散时间的动态系统建模、分析和仿真的核心。Real-TimeWorkshop是一种实时代码生成工具，它能够根据Simulink模型生成程序源代码，并打包、编译所生成的源代码生成实时应用程序。从现有的Simulink和Stateflow自动生成C语言程序代码的功能、定点运算模块集（Fixed-pointBlockset）与C语言程序代码到VHDL（VeryHighSpeedIntegratedCircuitHardwareDescriptionLanguage，一种标准的硬件 电路 模拟电路李宁答案12数字电路仿真实验电路与电子学第1章单片机复位电路图组合逻辑电路课后答案 设计语言）的自动转换功能，可以看出，高级的系統仿真或低级的芯片算法设计，都可用MATLAB、Simulink、Stateflow及相关的工具箱来完成。Stateflow是基于有限状态机理论针对复杂的事件驱动系统进行建模、仿真的工具。StateflowCoder是基于Stateflow状态图生成高效、优化的程序代码。围绕着Simulink仿真核心所开发的应用程序包，称为模块集（Blocksets）,MATLAB产品提供许多专用模块集，如CommunicationBlockset、DSPBlockset、SimPowerSystemBlockset、SignalProcessingBlockset等，详见MATLAB在线帮助文档。MATLAB是MATLAB产品家族的计算核心与基础，是集高性能数值计算与数据可视化于一体的高效编程语言。.1.3MATLAB编程语言的特点语法 规则 编码规则下载淘宝规则下载天猫规则下载麻将竞赛规则pdf麻将竞赛规则pdf 简单。尤其内定的编程规则，与其他编程语言（如C、Fortran等）相比更接近于常规数学表示。对于数组变量的使用，不需类型声明，无需事先申请内存空间。MATLAB语言主要有以下几个特点：MATLAB基本的语言环境提供了数以千计的计算函数，极大的提高了用户的编程效率。如，一个fft函数即可完成对指定数据的快速傅里叶变换，这一任务如果用C语言来编程实现的话，至少要用几十条C语言才能完成。因此，MATLAB是一个简单易用、功能强大的高效编程语言。平台无关性（可移植性）。MATLAB软件可以运行在很多不同的计算机系统平台上，如WindowsMe/NT/2000/XP、很多不同版本的UNIX以及Linux。无论你在哪一个平台上编写的程序都可以运行在其它平台上，对于MATLAB数据文件也一样，是平台无关的。极大保护了用户的劳动、方便了用户。其绘图功能也是平台无关的。无论任何系统平台，只要MATLAB能够运行，其图形功能命令就能正常运行。MATLAB是一种脚本式（scripted）的解释型语言，无论是命令、函数或变量，只要在命令窗口的提示符下键入，并“回车（Enter）”，MATLAB都予以解释执行。.功能强大数值运算优势符号运算优势(Maple)强大的2D、3D数据可视化功能许多具有算法自适应能力的功能函数.语言简单、内涵丰富语言及其书写形式非常接近于常规数学书写形式；其操作和功能函数指令就是常用的计算机和数学书上的一些简单英文单词表达的，如：help、clear等；完备的帮助系统，易学易用。扩充能力、可开发能力较强MATLAB完全成了一个开放的系统用户可以开发自己的工具箱可以方便地与Fortran、C等语言接口编程易、效率高Matlab以数组为基本计算单元具有大量的算法优化的功能函数.2MATLAB的桌面环境及入门知识启动与退出MATLAB命令窗口及使用数值表示、变量、表达式命令历史窗口工作空间获取在线帮助.2.1启动与退出MATLAB启动MATLAB直接用鼠标双击桌面上MATLAB7.1图标或Windows桌面的“开始”—〉“所有程序”—〉“MATLAB7.1”—〉“MATLAB7.1”。退出MATLAB关闭MATLAB桌面在命令窗口执行quit或exit命令MATLAB缺省桌面（见下页）.菜单栏工具栏文件编辑窗口工作空间历史命令窗口Start菜单命令窗口命令提示符2.1启动与退出MATLAB（续）.2.2命令窗口的使用激活命令窗口。“>>”与闪烁的光标一起表明系统就绪，等待输入。命令窗口脱离MATLAB桌面。简单计算【例2.2-1】计算（1）在MATLAB命令窗口输入以下内容：>>(12+2*(7-4))/3^2（2）按【Enter】键，指令执行。（3）返回的计算结果：ans=2.2.2命令窗口（续）〖说明〗在命令窗口【Enter】键提交命令执行。Matlab所用运算符（如+、-、^等）是各种计算程序中常见的。计算结果中的“ans”是英文“answer”的一种缩写，其含义就是“运算答案”。ans是Matlab的一个预定义变量。.2.2命令窗口（续）简单计算（续）【例2.2-2】计算sin(45ْ)>>sin(45*pi/180)Matalb中正弦函数sin就是常见的正弦函数。它的参数值是以“弧度”为单位的。pi也是Matalb的预定义变量。pi=3.14159…Matlab对字母大小写是敏感的。>>sqrt(2*exp(4.92+0.5)+1)Matalb中开平方—sqrt(x)，是英文squareroot的缩写。Matalb中指数函数exp(x)，常见的表达方式。ans=0.7071ans=21.2781.2.2命令窗口（续）“clc”清除窗口显示内容的命令。>>y=2*sin(0.3*pi)/(1+sqrt(5))y=0.5000>>y=2*cos(0.3*pi)/(1+sqrt(5))y=0.3633命令行编辑“↑”键调回已输入过命令。修改。.2.2命令窗口（续）【例2.2-5】计算半径为5.2m的圆的周长和面积。>>radius=5.2;%圆的半径>>area=pi*5.2^2,circle_len=2*pi*5.2area=84.9487circle_len=32.6726以上两例，命令行中用到了等号“=”。计算结果不再赋给“ans”，而是赋给用户指定的变量y、area、circle_len。无论是预定义变量还是用户自定义变量都被存储在系统的工作空间内，即系统定义的一个存储窗口变量的内存空间。Who、whos命令用来显示工作空间的变量clear命令用来清除工作空间的变量。.2.2命令窗口（续）>>who>>whosYourvariablesare:anscircle_lenyarearadiusNameSizeBytesClassans1x18doublearrayarea1x18doublearraycircle_len1x18doublearrayradius1x18doublearrayy1x18doublearrayGrandtotalis5elementsusing40bytes.2.2命令窗口（续）>>cleary>>whoYourvariablesare:anscircle_lenarearadius>>clearansarea>>whosYourvariablesare:NameSizeBytesClasscircle_len1x18doublearrayradius1x18doublearrayGrandtotalis2elementsusing16bytes.2.2命令窗口（续）数值显示格式设置缺省显示格式：简洁的短（shortg）格式窗口命令及语法格式：format显示格式关键字如：formatlong%15位数字显示常见通用命令命令含义clc清除命令窗口的显示内容clear清除Matlab工作空间中保存的变量who或whos显示Matlab工作空间中的变量信息dir显示当前工作目录的文件和子目录清单cd显示或设置当前工作目录type显示指定m文件的内容help或doc获取在线帮助quit或exit关闭/推出MATALB.2.3工作空间查看工作空间内存变量，可以由who、whos。命名新变量。修改变量名删除变量绘图保存变量数据装入数据.2.4历史窗口历史窗口：首先记录每次启动时间并记录在命令窗口输入命令，此次运行期间，输入的所有命令被记录为一组，并以此次启动时间为标志。使用历史窗口：可以查看命令窗口输入过的命令或语句可以选择一条或多条命令执行拷贝、执行、创建M文件等。要清除历史记录，可以选择Edit菜单中的ClearCommandHistory命令.2.5当前目录窗口和搜索路径当前目录窗口：指Matlab运行时的工作目录。只有在当前目录和搜索路径下的文件、函数才可以被运行和调用。如果没有特殊指明，数据文件也将存放在当前目录下；用户可以将自己的工作目录设置成当前目录，从而使得所有操作都在当前目录中进行。.搜索路径：指Matlab执行过程中对变量、函数和文件进行搜索的路径。在File菜单中选择SetPath命令或在命令窗口输入pathtool命令，出现搜索路径设置对话框：！修改完搜索路径后，需要进行保存。.2.6获取在线帮助MATLAB提供的帮助信息有两类简单纯文本帮助信息helplookfor（条件比较宽松）例：inverse窗口式综合帮助信息（文字、公式、图形）dochelpwin.【功能演示-1】求方程的全部根。p=[2,0,-3,71,-9,13];％建立多项式系数向量x=roots(p)；求根x=-3.49141.6863+2.6947i1.6863-2.6947i0.0594+0.4251i0.0594-0.4251i.【功能演示-2】求解线性方程组a=[2,3,-1;8,2,3;45,3,9];％建立系数矩阵ab=[2;4;23];％建立列向量bx=inv(a)*bx=0.55310.2051-0.2784.符号计算symsxyz％建立符号变量[x,y,z]=solve(2*x+3*y-z-2,8*x+2*y+3*z-4,45*x+3*y+9*z-23)x=151/273y=8/39z=-76/273.【功能演示-3】求解定积分求解I=quad('x.*log(1+x)',0,1)ans=0.250或symsxnt(x*log(1+x),0,1)ans=1/4.【功能演示-4】多项式曲线拟合考虑如下x-y一组实验数据：x=[1,2,3,4,5,6,7,8,9,10]y=[1.2,3,4,4,5,4.7,5,5.2,6,7.2]注：一次多项式拟合：p1=polyfit(x,y,1)三次多项式拟合：p3=polyfit(x,y,3)plot原始数据、一次拟合曲线和三次拟合曲线x2=1:0.1:10;y1=polyval(p1,x2)y3=polyval(p3,x2)plot(x,y,’*’,x2,y1,‘:’,x2,y3).拟合曲线图由图可见，三次拟合结果较好。.2.3数值表示、变量及表达式数值的记述Matlab的数只采用习惯的十进制表示，可以带小数点和负号;其缺省的数据类型为双精度浮点型（double）。例如：3-100.0011.3e101.256e-6变量命令规则变量名、函数名对字母的大小写是敏感的。如myVar与myvar表示两个不同的变量。变量名第一个字母必须是英文字母。变量名可以包含英文字母、下划线和数字。变量名不能包含空格、标点。变量名最多可包含63个字符（6.5及以后的版本）。.2.3数值表示、变量及表达式（续）Matlab预定义的变量〖说明〗每当MATLAB启动完成，这些变量就被产生。MATLAB中，被0除不会引起程序中断，给出报警的同时用inf或NaN给出结果。用户只能临时覆盖这些预定义变量的值，Clear或重启MATLAB可恢复其值。变量名意义ans最近的计算结果的变量名epsMATLAB定义的正的极小值=2.2204e-16pi圆周率πinf∞值，无限大i或j虚数单元，sqrt(-1)NaN非数，0/0、∞/∞.2.3数值表示、变量及表达式（续）运算符和表达式〖说明〗Matlab用“\”和”/”分别表示“左除”和“右除”。对标量而言，两者没有区别。对矩阵产生不同影响。MATLAB表达式的书写规则与“手写方式”几乎完全相同。表达式按与常规相同的优先级自左至右执行运算。优先级：指数运算级别最高，乘除次之，加减最低。括号改变运算的次序。.2.3数值表示、变量及表达式（续）复数及其运算MATLAB中复数的表达：z=a+bi，其中a、b为实数。MATLAB把复数作为一个整体，象计算实数一样计算复数。>>z1=3+4*i,z2=1+2*i,z3=2*exp(i*pi/6),z=z1*z2/z3>>z_real=real(z),z_image=imag(z),>>z_angle=angle(z),z_length=abs(z),.课堂总结对MATALB整个软件产品进行了概述历史发展软件产品家族体系构成M语言的特点MATLAB的入门知识MATLAB桌面MATLAB数值的记述和数据显示格式变量及其命名规则运算符及表达式常用的MATLAB命令在线帮助的使用.思考题简述MATALB软件的功能及特点。简述MATALB变量的命名规则。熟悉课件中的例子。.Matlab语言及其应用第2讲.实验安排.工作空间Matlab系统界面.2.2命令窗口（续）【例2.2-5】计算半径为5.2m的圆的周长和面积。>>radius=5.2;%圆的半径>>area=pi*5.2^2,circle_len=2*pi*5.2area=84.9487circle_len=32.6726以上两例，命令行中用到了等号“=”。计算结果不再赋给“ans”，而是赋给用户指定的变量y、area、circle_len。无论是预定义变量还是用户自定义变量都被存储在系统的工作空间内，即系统定义的一个存储窗口变量的内存空间。Who、whos命令用来显示工作空间的变量clear命令用来清除工作空间的变量。.2.2命令窗口（续）>>who>>whosYourvariablesare:circle_lenarearadiusNameSizeBytesClassarea1x18doublearraycircle_len1x18doublearrayradius1x18doublearrayGrandtotalis3elementsusing24bytes.2.2命令窗口（续）>>clearradius>>whoYourvariablesare:circle_lenarea>>cleararea>>whosYourvariablesare:NameSizeBytesClasscircle_len1x18doublearrayGrandtotalis1elementsusing8bytes.2.2命令窗口（续）数值显示格式设置缺省显示格式：简洁的短（shortg，7位有效数字）格式窗口命令及语法格式：format显示格式关键字如：formatlong%15位数字显示常见通用命令命令含义clc清除命令窗口的显示内容clear清除Matlab工作空间中保存的变量who或whos显示Matlab工作空间中的变量信息dir显示当前工作目录的文件和子目录清单cd显示或设置当前工作目录type显示指定m文件的内容help或doc获取在线帮助quit或exit关闭/推出MATALB.2.3工作空间查看工作空间内存变量，可以由who、whos。命名新变量修改变量名删除变量绘图保存变量数据装入数据.2.4历史窗口历史窗口：首先记录每次启动时间并记录在命令窗口输入命令，此次运行期间，输入的所有命令被记录为一组，并以此次启动时间为标志。使用历史窗口：可以查看命令窗口输入过的命令或语句可以选择一条或多条命令执行拷贝、执行、创建M文件等。要清除历史记录，可以选择Edit菜单中的ClearCommandHistory命令.2.5当前目录窗口和搜索路径当前目录窗口：指Matlab运行时的工作目录。只有在当前目录和搜索路径下的文件、函数才可以被运行和调用。如果没有特殊指明，数据文件也将存放在当前目录下；用户可以将自己的工作目录设置成当前目录，从而使得所有操作都在当前目录中进行。.搜索路径：指Matlab执行过程中对变量、函数和文件进行搜索的路径。在File菜单中选择SetPath命令或在命令窗口输入pathtool命令，出现搜索路径设置对话框：！修改完搜索路径后，需要进行保存。.2.6获取在线帮助MATLAB提供的帮助信息有两类简单纯文本帮助信息helplookfor（条件比较宽松）例：inverse窗口式综合帮助信息（文字、公式、图形）dochelpwin.【功能演示-1】求方程的根p=[2,0,-3,71,-9,13];％建立多项式系数向量x=roots(p)；求根x=-3.49141.6863+2.6947i1.6863-2.6947i0.0594+0.4251i0.0594-0.4251i.【功能演示-2】求解线性方程组a=[2,3,-1;8,2,3;45,3,9];％建立系数矩阵ab=[2;4;23];％建立列向量bz=inv(a)*bz=0.55310.2051-0.2784.符号计算symsxyz％建立符号变量[x,y,z]=solve(2*x+3*y-z-2,8*x+2*y+3*z-4,45*x+3*y+9*z-23)x=151/273y=8/39z=-76/273.【功能演示-3】求解定积分求解I=quad('x.*log(1+x)',0,1)ans=0.250或symsxint(x*log(1+x),0,1)ans=1/4.【功能演示-4】多项式曲线拟合考虑如下x-y一组实验数据：x=[1,2,3,4,5,6,7,8,9,10]y=[1.2,3,4,4,5,4.7,5,5.2,6,7.2]注：一次多项式拟合：p1=polyfit(x,y,1)三次多项式拟合：p3=polyfit(x,y,3)plot原始数据、一次拟合曲线和三次拟合曲线x2=1:0.1:10;y1=polyval(p1,x2)y3=polyval(p3,x2)plot(x,y,’*’,x2,y1,‘:’,x2,y3).拟合曲线图由图可见，三次拟合结果较好。.第2章Matlab数据及其运算Matlab数据类型非常丰富，除数值型、字符型等基本数据类型外，还有结构体、单元等更为复杂的数据类型。各种数据类型都以矩阵形式存在，矩阵是Matlab最基本的数据对象，并且矩阵的运算是定义在复数域上的。.2.1数值表示、变量及表达式数值的记述Matlab的数只采用习惯的十进制表示，可以带小数点和负号;其缺省的数据类型为双精度浮点型（double）。例如：3-100.0011.3e101.256e-6变量命令规则变量名、函数名对字母的大小写是敏感的。如myVar与myvar表示两个不同的变量。变量名第一个字母必须是英文字母。变量名可以包含英文字母、下划线和数字。变量名不能包含空格、标点。变量名最多可包含63个字符（6.5及以后的版本）。.2.1数值表示、变量及表达式（续）Matlab预定义的变量〖说明〗每当MATLAB启动完成，这些变量就被产生。MATLAB中，被0除不会引起程序中断，给出报警的同时用inf或NaN给出结果。用户只能临时覆盖这些预定义变量的值，Clear或重启MATLAB可恢复其值。变量名意义ans最近的计算结果的变量名epsMATLAB定义的正的极小值=2.2204e-16pi圆周率πinf∞值，无限大i或j虚数单元，sqrt(-1)NaN非数，0/0、∞/∞.2.1数值表示、变量及表达式（续）运算符和表达式〖说明〗Matlab用“\”和”/”分别表示“左除”和“右除”。对标量而言，两者没有区别。对矩阵产生不同影响。MATLAB表达式的书写规则与“手写方式”几乎完全相同。表达式按与常规相同的优先级自左至右执行运算。优先级：指数运算级别最高，乘除次之，加减最低。括号改变运算的次序。.2.1数值表示、变量及表达式（续）复数及其运算MATLAB中复数的表达：z=a+bi，其中a、b为实数。MATLAB把复数作为一个整体，象计算实数一样计算复数。【例2.3-1】复数z1=3+4i，z2=1+2i,z3=计算z1=3+4*i,z2=1+2*i,z3=exp(i*pi/6),z=z1*z2/z3z_real=real(z),z_image=imag(z),z_angle=angle(z),z_length=abs(z),.2.2Matlab矩阵(数组)的表示数组的概念一维数组变量的创建二维数组变量的创建数组元素的标识与寻访数组运算多维数组.2.2.1.数组(array)的概念数组定义：按行(row)和列(column)顺序排列的实数或复数的有序集，被称为数组。数组中的任何一个数都被称为这个数组的元素，由其所在的行和列标识，这个标识也称为数组元素的下标或索引。Matlab将标量视为1×1的数组。对m行、n列的2维数组a：计为m×n的数组a；*行标识、列标识均从1开始；行标识从上到下递增；列标识从左到右递增。a(3,4)=34rowisfirsta=.数组的分类一维数组，也称为向量(vector)。行向量(rowvector)、列向量(columnvector)。二维数组(矩阵matrix)。多维数组。有效矩阵：每行元素的个数必须相同，每列元素的个数也必须相同。2.2.1.数组(array)的概念.2.2.1、数组(array)的概念（续）行向量列向量a(2,1)=3a(1,2)=2b(3)=3c(2)=2.2.2.2、创建一维数组变量第一种方法：使用方括号“[]”操作符【例2-1】创建数组(行向量)a=[13pi3+5i]>>a=[13pi3+5*i]%ora=[1,3,pi,3+5*i]a=1.00003.00003.14163.0000+5.0000i所有的向量元素必须在操作符“[]”之内；向量元素间用空格或英文的逗点“,”分开。第二种方法：使用冒号“:”操作符【例2-2】创建以1~10顺序排列整数为元素的行向量b。>>b=1:10b=12345678910.2.2.2、创建一维数组变量（续）【例2-3】键入并执行c=1:2:10和d=1:2:9>>c=1:2:10c=13579>>d=1:2:9d=13579利用冒号“:”操作符创建行向量的基本语法格式：x=Start:Increment:EndStart表示新向量x的第一个元素；新向量x的最后一个元素不能大于End；Increment可正可负，若负，则必须Start>End；若正，则必须Start>a=[123];b=[234];>>c=[a;b];>>c1=[ab];第二种方法：函数方法函数ones(生成全1矩阵)、zeros(生成全0矩阵)、reshape“helpelmat”获得基本的矩阵生成和操作函数列表【例3-4】创建全1的3x3数组。>>ones(3)【例3-5】创建全1的3x4数组。>>ones(3,4).2.2.3、创建二维数组变量（续）【例3-5】reshape的使用演示>>a=-4:4a=-4-3-2-101234>>b=reshape(a,3,3)b=-4-12-303-214数组元素的排列顺序，从上到下按列排列，先排第一列，然后第二列，…要求数组的元素总数不变。.4、数组元素的标识与寻访数组元素的标识“全下标（index）”标识经典数学教科书采用“全下标”标识法：每一维对应一个下标。如对于二维数组，用“行下标和列下标”标识数组的元素，a(2,3)就表示二维数组a的“第2行第3列”的元素。对于一维数组，用一个下标即可，b(2)表示一维数组b的第2个元素，无论b是行向量还是列向量。“单下标”（linearindex）标识所谓“单下标”标识就是用一个下标来表明元素在数组的位置。对于二维数组，“单下标”编号：设想把二维数组的所有列，按先后顺序首尾相接排成“一维长列”，然后自上往下对元素位置执行编号。两种“下标”标识的变换：sub2ind、ind2sub.4、数组元素的标识与寻访（续）【例4-1】单下标的使用>>a=zeros(2,5);>>a(:)=-4:5a=-4-2024-3-1135注意数组的排列顺序。.4、数组元素的标识与寻访（续）元素与子数组的寻访与赋值【例4-3】一维数组元素与子数组的寻访与赋值>>a=linspace(1,10,5)a=1.00003.25005.50007.750010.0000>>a(3)%寻访a的第3个元素ans=5.5000>>a([125])%寻访a的第1、2、5个元素组成的子数组ans=1.00003.250010.0000.4、数组元素的标识与寻访（续）>>a(1:3)%寻访前3个元素组成的子数组ans=1.00003.25005.5000>>a(3:-1:1)%由前3个元素倒序构成的子数组ans=5.50003.25001.0000>>a(3:end)ans=5.50007.750010.0000>>a(3:end-1)ans=5.50007.7500%第3个及其后所有元素构成的子数组函数end作为参数使用，返回最后一个元素的下标.>>a([12355321])ans=1.00003.25005.500010.000010.00005.50003.25001.0000数组元素可以被任意重复访问，构成长度大于原数组的新数组。>>a(6)???Indexexceedsmatrixdimensions.下标值超出了数组的维数，导致错误>>a(2.1)???Subscriptindicesmusteitherberealpositiveintegersorlogicals.下标值只能取正整数或逻辑值4、数组元素的标识与寻访（续）.4、数组元素的标识与寻访（续）>>a(3)=0%修改数组a的第3元素值为0a=1.00003.250007.750010.0000>>a([25])=[11]a=1.00001.000007.75001.0000可以修改指定数组元素的值一次可以修改多个数组元素的值要修改的数组元素的个数应与送入数组的元素个数相同.4、数组元素的标识与寻访（续）【例4-3】二维数组元素与子数组的寻访与赋值>>a_2=zeros(2,4)%创建2x4的全0数组a_2=00000000>>a_2(:)=1:8a_2=13572468>>a_2([258])%单下标方式寻访多个元素ans=258注意元素的排列顺序.4、数组元素的标识与寻访（续）>>a_2([258])=[102030]a_2=13207104630>>a_2(:,[23])=ones(2)%双下标方式寻访并修改a_2=1117101130二维数组可以“单下标”方式或“全下标”方式访问、赋值；“单下标”方式赋值时，等号两边涉及的元素个数必须相等；“全下标”方式赋值时，等号右边数组的大小必须等于原数组中涉及元素构成的子数组的大小。.4、数组元素的标识与寻访（续）>>a_2(:,end)ans=730>>a_2(:,end-1)ans=11>>a_2(:,end:-1:3)ans=71301>>a_2(end,:)ans=101130>>a_2(end,[2:4])ans=1130>>a_2([46])=6:7a_2=1117106730>>a_2(end,[2:end-1])Whatistheresult?.4、数组元素的标识与寻访（续）【例4-4】size、length函数>>a=ones(4,6)*6>>m=size(a)>>len=length(a)>>b=1:5;>>length(b)>>c=b’>>length(c)size函数返回变量的大小，即变量数组的行列数length函数返回变量数组的最大维数.4、数组元素的标识与寻访（续）双下标到单下标的转换【例4-5】sub2ind函数-双下标转换为单下标>>A=[172418;222714;461320];>>A(:,:,2)=A-10>>A(2,1,2)>>sub2ind(size(A),2,1,2)>>A(14).4、数组元素的标识与寻访（续）单下标到双下标的转换【例4-6】ind2sub函数-双下标转换为单下标>>b=zeros(3);>>b(:)=1:9>>IND=[3456]>>[I,J]=ind2sub(size(b),IND).Matlab语言及其应用第3讲.2.2.2、创建一维数组变量第一种方法：使用方括号“[]”操作符【例2-1】创建数组(行向量)a=[13pi3+5i]>>a=[13pi3+5*i]%ora=[1,3,pi,3+5*i]a=1.00003.00003.14163.0000+5.0000i所有的向量元素必须在操作符“[]”之内；向量元素间用空格或英文的逗点“,”分开。第二种方法：使用冒号“:”操作符【例2-2】创建以1~10顺序排列整数为元素的行向量b。>>b=1:10b=12345678910.2.2.2、创建一维数组变量（续）【例2-3】键入并执行c=1:2:10和d=1:2:9>>c=1:2:10c=13579>>d=1:2:9d=13579利用冒号“:”操作符创建行向量的基本语法格式：x=Start:Increment:EndStart表示新向量x的第一个元素；新向量x的最后一个元素不能大于End；Increment可正可负，若负，则必须Start>End；若正，则必须Start>a=[123];b=[234];>>c=[a;b];>>c1=[ab];第二种方法：函数方法函数ones(生成全1矩阵)、zeros(生成全0矩阵)、reshape“helpelmat”获得基本的矩阵生成和操作函数列表【例3-4】创建全1的3x3数组。>>ones(3)【例3-5】创建全1的3x4数组。>>ones(3,4).2.2.3、创建二维数组变量（续）【例3-5】reshape的使用演示>>a=-4:4a=-4-3-2-101234>>b=reshape(a,3,3)b=-4-12-303-214数组元素的排列顺序，从上到下按列排列，先排第一列，然后第二列，…要求数组的元素总数不变。.4、数组元素的标识与寻访数组元素的标识“全下标（index）”标识经典数学教科书采用“全下标”标识法：每一维对应一个下标。如对于二维数组，用“行下标和列下标”标识数组的元素，a(2,3)就表示二维数组a的“第2行第3列”的元素。对于一维数组，用一个下标即可，b(2)表示一维数组b的第2个元素，无论b是行向量还是列向量。“单下标”（linearindex）标识所谓“单下标”标识就是用一个下标来表明元素在数组的位置。对于二维数组，“单下标”编号：设想把二维数组的所有列，按先后顺序首尾相接排成“一维长列”，然后自上往下对元素位置执行编号。两种“下标”标识的变换：sub2ind、ind2sub.4、数组元素的标识与寻访（续）【例4-1】单下标的使用>>a=zeros(2,5);>>a(:)=-4:5a=-4-2024-3-1135注意数组的排列顺序。.4、数组元素的标识与寻访（续）元素与子数组的寻访与赋值【例4-3】一维数组元素与子数组的寻访与赋值>>a=linspace(1,10,5)a=1.00003.25005.50007.750010.0000>>a(3)%寻访a的第3个元素ans=5.5000>>a([125])%寻访a的第1、2、5个元素组成的子数组ans=1.00003.250010.0000.4、数组元素的标识与寻访（续）>>a(1:3)%寻访前3个元素组成的子数组ans=1.00003.25005.5000>>a(3:-1:1)%由前3个元素倒序构成的子数组ans=5.50003.25001.0000>>a(3:end)ans=5.50007.750010.0000>>a(3:end-1)ans=5.50007.7500%第3个及其后所有元素构成的子数组函数end作为参数使用，返回最后一个元素的下标.>>a([12355321])ans=1.00003.25005.500010.000010.00005.50003.25001.0000数组元素可以被任意重复访问，构成长度大于原数组的新数组。>>a(6)???Indexexceedsmatrixdimensions.下标值超出了数组的维数，导致错误>>a(2.1)???Subscriptindicesmusteitherberealpositiveintegersorlogicals.下标值只能取正整数或逻辑值4、数组元素的标识与寻访（续）.4、数组元素的标识与寻访（续）>>a(3)=0%修改数组a的第3元素值为0a=1.00003.250007.750010.0000>>a([25])=[11]a=1.00001.000007.75001.0000可以修改指定数组元素的值一次可以修改多个数组元素的值要修改的数组元素的个数应与送入数组的元素个数相同.4、数组元素的标识与寻访（续）【例4-3】二维数组元素与子数组的寻访与赋值>>a_2=zeros(2,4)%创建2x4的全0数组a_2=00000000>>a_2(:)=1:8a_2=13572468>>a_2([258])%单下标方式寻访多个元素ans=258注意元素的排列顺序.4、数组元素的标识与寻访（续）>>a_2([258])=[102030]a_2=13207104630>>a_2(:,[23])=ones(2)%双下标方式寻访并修改a_2=1117101130二维数组可以“单下标”方式或“全下标”方式访问、赋值；“单下标”方式赋值时，等号两边涉及的元素个数必须相等；“全下标”方式赋值时，等号右边数组的大小必须等于原数组中涉及元素构成的子数组的大小。.4、数组元素的标识与寻访（续）>>a_2(:,end)ans=730>>a_2(:,end-1)ans=11>>a_2(:,end:-1:3)ans=71301>>a_2(end,:)ans=101130>>a_2(end,[2:4])ans=1130>>a_2([46])=6:7a_2=1117106730>>a_2(end,[2:end-1])Whatistheresult?.4、数组元素的标识与寻访（续）【例4-4】size、length函数>>a=ones(4,6)*6>>m=size(a)>>len=length(a)>>b=1:5;>>length(b)>>c=b’>>length(c)size函数返回变量的大小，即变量数组的行列数length函数返回变量数组的最大维数.4