高等应用数学问题的MATLAB求解2

4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,第2章MATLAB语言程序设计基础,高等应用数学问题的MATLAB求解,清华大学出版社2008,CAI课件开发：薛定宇、刘莹莹、董雯彬,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,第2章MATLAB语言程序设计基础,MATLAB程序设计语言基础基本数学运算MATLAB语言的流程结构函数的编写与调试二维图形绘制三维图形绘制图像处理简介,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,MATLAB语言有如下的优点MATLAB语言的简洁高效性MATLAB语言的科学运算功能MATLAB语言的绘图功能MATLAB庞大的工具箱与模块集MATLAB强大的动态系统仿真功能MATLAB语言是当前国际上自动控制领域的首选计算机语言，也是很多理工科专业最适合的计算机数学语言，本课程以MATLAB为主要计算机语言。,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.1MATLAB程序设计语言基础,MATLAB语言的变量与常量数据结构MATLAB的基本语句结构冒号表达式与子矩阵提取,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.1.1MATLAB语言的变量与常量,MATLAB语言的变量名规则由一个字母引导，后面可以为其他字符区分大小写AbcABc有效MYvar12,MY_Var12和MyVar12_错误的变量名12MyVar,_MyVar12MATLAB的保留常量eps,i,j,pi,NaN,Inf,i=sqrt(-1)lastwarn,lasterr,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.1.2数据结构,数值型数据符号型其他数据结构,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.1.2.1数值型数据,双精度数值变量IEEE标准，64位(8字节)，11指数位，53数值位和一个符号位至double()函数的转换其他数据类型uint8，常用于图像表示和处理，8位int8(),int16(),int32(),uint16(),uint32(),4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.1.2.2符号型,符号型，sym(A)，常用于公式推导常用于公式推导和求解析解变量声明显示符号变量的任何精度默认精度：32位小数,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.1,求出p的300位有效数字MATLAB代码可以任意增加有效数字的位数随着有效数字位数的增加，程序运行的时间变长,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.1.2.3其他数据结构,字符串型数据：用单引号括起来多维数组：是矩阵的直接扩展，多个下标单元数组：将不同类型数据集成到一个变量名下面，用表示类与对象：可以定义重载函数,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.1.3MATLAB的基本语句结构,直接赋值语句接赋值语句的基本结构赋值语句的末尾加一个分号可以阻止显示运算结果保留变量：ans存放最近一次无赋值变量语句的运算结果,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.2,表示矩阵MATLAB代码其他命令,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.3,试输入复数矩阵MATLAB代码需要避免的语句,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,函数调用语句函数调用举例函数可以通过不同的方式被调用内在函数，*.m函数匿名函数，内联inline函数重载函数,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.1.4冒号表达式与子矩阵提取,冒号表达式是定义行向量的有效方法开始于值，步长为值，终止于值.默认步长值为1,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.4,用不同的步距生间的向量,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,子矩阵提取的基本语句格式表示子矩阵要保留的行号构成的向量表示要保留的列号构成的向量:，表示要提取所有的行或列，取决于它的位置,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.5,矩阵A的不同子矩阵提取MATLAB代码,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.2基本数学运算,矩阵的代数运算矩阵的逻辑运算矩阵的比较运算解析结果的化简与变换基本数论运算,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.2.1矩阵的代数运算,矩阵转置加减法运算矩阵乘法矩阵的除法矩阵翻转矩阵乘方运算点运算,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,矩阵转置,矩阵表示矩阵A，n行m列，被称作矩阵厄米转置一般转置,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,加减法运算,数学表示用C语言编程难于实现，同A*BMATLAB语法注意：任一个变量可以为标量如果矩阵维数不匹配，系统会报错,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,矩阵乘法,数学表示：MATLAB表示：注意：系统自动检测矩阵维数是否匹配,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,矩阵的除法,矩阵左除：求解线性方程组：MATLAB解法：最小二乘解若A为非奇异方阵，则,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,矩阵右除：求解线性方程组：MATLAB解法：最小二乘解若A为非奇异方阵，则更精确的解法,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,矩阵翻转,左右翻转上下翻转旋转90o如何旋转180o？,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,矩阵乘方,A为方阵，求矩阵A的x次幂数学描述MATLAB命令,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,点运算,矩阵对应元素的直接运算例如：为另一个例子：为,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.6,求出矩阵A的全部三次方根，并检验结果MATLAB代码另两个根,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.2.2矩阵的逻辑运算,逻辑变量：当前版本有逻辑变量非0表示逻辑1逻辑运算（相应元素间的运算）与运算或运算非运算异或运算,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.2.3矩阵的比较运算,各种允许的比较关系,=,=,=,=,find(),all(),any()实例,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.2.4解析结果的化简与变换,函数simple()用于数学公式的化简其他常用化简函数numden(),collect(),expand(),factor(),4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.7,化简多项式用不同的函数求解,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,变量替换,变量替换语句该函数执行点运算转换成LaTeX表示,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.8,试用对下式进行双线性变换，得出其LaTeX排版格式MATLAB代码,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,结果：8,fracleft(2,z+1right)2zleft(3,z+1right)left(5,z+3right)3left(z+1right)7LaTeX排版显示,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.2.5基本数论运算,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.9,运用各种函数，对下面的数据进行取整运算-0.2765,0.5772,1.4597,2.1091,1.191,-1.6187MATLAB代码：,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.10,假设的Hilbert矩阵可以由定义，试对其进行有理数变换。结果：,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.11,试求1856120和1483720的最大公约数与最小公倍数，并求出所得出的最小公倍数的质因数分解。MATLAB求解语句：,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.12,试求出11000间的全部质数得到的所有质数,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.3MATLAB语言的流程结构,循环结构转移结构开关结构试探结构,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.3.1循环结构,for循环结构为一个向量，循环变量每次从向量中取一个数值，执行一次循环体的内容，如此下去，直至执行完向量中所有的分量,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,while循环结构,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.13,用循环求解更简单的语句,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.14,用循环求解求最小的m，使下式成立使用while结构这里不宜应用for结构,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.15,求如下值使用MATLAB的循环结构和向量化运算方法,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.3.2转移结构,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.16,用for循环求解最大的m，使下式成立MATLAB代码远比while结构繁琐,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.3.3开关结构,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.3.4试探结构,全新结构优点：设置错误陷阱提高算法的执行效率,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.4函数编写与调试,MATLAB语言函数的基本结构可变输入输出个数的处理inline函数与匿名函数,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.4.1MATLAB语言函数的基本结构,nargin,nargout,varargin,varargout,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.17,为什么需要函数？问题M-脚本以M-文件的形式被存取若最大值不为10000，M-文件就需要被修改则需要一种新的文件形式,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.18,按例2.17的要求编写一个MATLAB函数M-文件如下举例优点：无需修改程序,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.19,编写一个函数生成Hilbert矩阵要求若只给出一个输入参数，则会自动生成一个方阵在函数中给出合适的帮助信息检测输入和返回变量的个数,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,MATLAB函数如下,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,在线帮助命令生成Hilbert矩阵,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.20,利用函数的递归调用，计算阶乘递归调用函数,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,计算11！其他函数,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.21,试比较递归算法和循环算法在Fibonacci数列中应用的优劣Fibonacci数列MATLAB代码,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,第25项用for循环结构计算第100项注意：本例不宜采用递归结构,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.4.2可变输入输出个数的处理,vararginvarargoutvarargin1,varargin2,vararginn,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.22,conv()可以计算两个多项式的积，时使用varargin实现任意多个多项式的积MATLAB代码,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,函数调用,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.4.3inline函数与匿名函数,inline函数可以免去编辑M-文件MATLAB7.0以上版本提供的匿名函数在MATLAB工作空间内已经定义的变量可以直接使用,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.5二维图形绘制,二维图形绘制基本语句其他二维图形绘制语句隐函数绘制及应用图形修饰,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.5.1二维图形绘制基本语句,两个序列和构造向量:利用这些数据绘图,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,其他方法仍为向量，而为矩阵，亦即和均为矩阵，且和的行和列数均相同假设有多对这样的向量或矩阵,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,改变曲线性质选项列表,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.23,绘制方程MATLAB代码问题：如何检查曲线？,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,改用不同的步长，观察曲线是否相同换用较小步长在附近用变步长算法,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.24,绘制饱和函数方程MATLAB绘图语句更简单的命令,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,图像修饰,每一个窗口、曲线和坐标轴都是一个对象，对象的属性可以通过函数set()来设置，也可以通过函数get()来获取。,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.5.2其他二维图形绘制语句,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.25,绘制极坐标函数和绘制曲线,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.26,以正弦数据为例，试在同一窗口的不同区域用不同的绘图方式绘制出相应的曲线MATLAB求解语句：,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,下述函数可以用于将一个绘图窗口分割成几个不同的区域n和m分别是行数和列数，k则指示区域的序号,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.5.3隐函数绘制及应用,隐函数绘图语句默认区域是-2p,2p其他语法,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.27,试绘制隐函数MATLAB语句扩大范围,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.5.4图形修饰,直接采用工具栏文字修饰特殊符号表下标分别用和_表示优越性MATLAB7.0的新功能,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.6三维图形表示,三维曲线绘制三维曲面绘制三维图形视角设置,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.6.1三维曲线绘制,三维曲线绘制stem3,fill3,bar3等,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,例2.28,试绘制参数方程其中，MATLAB绘图,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,使用stem()函数3D绘图和图形旋转,4/30/2020星期五,2008-8-15,12:22:05,高等应用数学问题的MATLAB求解东北大学信息学院,2.6.2三维曲面绘制,一般曲面绘制其他函数surfl(),surfc()等高线绘制contour(),contours(),4/