《高等应用数学问题的MATLAB求解-第二版》mathchap02

1、7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院第第2章章MATLAB语言程序设计基础语言程序设计基础高等应用数学问题的MATLAB求解清华大学出版社2008CAI课件开发：薛定宇、刘莹莹、董雯彬7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院第第2章章 MATLAB语言程序设计基础语言程序设计基础MATLAB 程序设计语言基础基本数学运算MATLAB语言的流程结构函数的编写与调试二维图形绘制三维图形绘制图像处理简介7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院MATLAB语言有如下的优点MATLAB

2、语言的简洁高效性MATLAB 语言的科学运算功能MATLAB 语言的绘图功能MATLAB 庞大的工具箱与模块集MATLAB 强大的动态系统仿真功能MATLAB 语言是当前国际上自动控制领域的首选计算机语言，也是很多理工科专业最适合的计算机数学语言，本课程以 MATLAB 为主要计算机语言。7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.1 MATLAB 程序设计语言基础程序设计语言基础MATLAB语言的变量与常量数据结构MATLAB的基本语句结构冒号表达式与子矩阵提取7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.1.1 MA

3、TLAB语言的变量与常量语言的变量与常量MATLAB语言的变量名规则由一个字母引导，后面可以为其他字符区分大小写Abc ABc有效 MYvar12, MY_Var12 和 MyVar12_错误的变量名 12MyVar, _MyVar12 MATLAB的保留常量eps, i, j, pi, NaN, Inf, i=sqrt(-1)lastwarn, lasterr7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.1.2 数据结构数据结构数值型数据符号型其他数据结构7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.1.2.1 数值型数

4、据数值型数据双精度数值变量IEEE标准，64位 (8字节)，11指数位，53数值位和一个符号位 至 double() 函数的转换其他数据类型uint8，常用于图像表示和处理，8位int8(),int16(),int32(),uint16(),uint32()7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.1.2.2 符号型符号型符号型，sym(A)，常用于公式推导常用于公式推导和求解析解变量声明显示符号变量的任何精度默认精度：32位小数7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.1求出 p 的 300 位有效数字MA

5、TLAB代码可以任意增加有效数字的位数随着有效数字位数的增加，程序运行的时间变长7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.1.2.3 其他数据结构其他数据结构字符串型数据：用单引号括起来 多维数组：是矩阵的直接扩展，多个下标单元数组：将不同类型数据集成到一个变量名下面，用 表示类与对象：可以定义重载函数7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.1.3 MATLAB的基本语句结构直接赋值语句接赋值语句的基本结构赋值语句的末尾加一个分号可以阻止显示运算结果保留变量： ans存放最近一次无赋值变量语句的运算结果7/5/2

6、021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.2表示矩阵MATLAB代码其他命令7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.3试输入复数矩阵MATLAB代码需要避免的语句7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院函数调用语句函数调用举例函数可以通过不同的方式被调用内在函数， *.m函数匿名函数， 内联inline函数重载函数7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.1.4 冒号表达式与子矩阵提取冒号表达式是定义行向量的有效方法开始于值 ，步长为

7、值 ，终止于值 .默认步长值为17/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.4用不同的步距生 间的向量7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院子矩阵提取的基本语句格式 表示子矩阵要保留的行号构成的向量 表示要保留的列号构成的向量:，表示要提取所有的行或列，取决于它的位置7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.5矩阵A的不同子矩阵提取MATLAB代码7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.2 基本数学运算基本数学运算矩阵的代数

8、运算矩阵的逻辑运算矩阵的比较运算解析结果的化简与变换基本数论运算7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.2.1矩阵的代数运算矩阵的代数运算矩阵转置加减法运算矩阵乘法矩阵的除法矩阵翻转矩阵乘方运算点运算7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院矩阵转置矩阵转置矩阵表示矩阵A，n 行m列，被称作 矩阵厄米转置一般转置7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院加减法运算加减法运算数学表示用C语言编程难于实现，同A*BMATLAB语法注意：任一个变量可以为标量如果矩阵维数不匹配，系统会报错

9、7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院矩阵乘法矩阵乘法数学表示：MATLAB表示：注意：系统自动检测矩阵维数是否匹配7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院矩阵的除法矩阵的除法矩阵左除：求解线性方程组：MATLAB解法：最小二乘解若A为非奇异方阵，则7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院矩阵右除：求解线性方程组：MATLAB解法：最小二乘解若A为非奇异方阵，则更精确的解法7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院矩阵翻转矩阵翻转左右翻转 上下

10、翻转旋转 90o如何旋转180o？7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院矩阵乘方矩阵乘方A为方阵，求矩阵A的x次幂数学描述MATLAB命令7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院点运算点运算矩阵对应元素的直接运算例如： 为 另一个例子： 为7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.6求出矩阵A的全部三次方根，并检验结果MATLAB代码另两个根7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.2.2 矩阵的逻辑运算矩阵的逻辑运算逻辑变量：当前

11、版本有逻辑变量非 0 表示逻辑 1逻辑运算（相应元素间的运算）与运算或运算非运算异或运算7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.2.3 矩阵的比较运算矩阵的比较运算各种允许的比较关系, =, , =, = =,=, find(), all(), any()实例7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.2.4 解析结果的化简与变换解析结果的化简与变换函数 simple() 用于数学公式的化简其他常用化简函数numden(), collect(), expand(), factor()7/5/2021Slide 1

12、(of 11)MATLAB求解东北大学信息学院例例 2.7化简多项式用不同的函数求解7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院变量替换变量替换变量替换语句该函数执行点运算转换成LaTeX表示7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.8试用 对下式进行双线性变换，得出其LaTeX排版格式MATLAB代码7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院结果： 8,frac left( 2,z+1 right) 2z left( 3,z+1 right) left( 5,z+3 right

13、) 3 left( z+1 right) 7LaTeX排版显示7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.2.5 基本数论运算基本数论运算7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.9运用各种函数，对下面的数据进行取整运算-0.2765,0.5772,1.4597,2.1091,1.191,-1.6187MATLAB代码：7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.10假设 的Hilbert矩阵可以由 定义，试对其进行有理数变换。结果：7/5/2021Slide 1 (

14、of 11)MATLAB求解东北大学信息学院例例 2.11试求1856120和1483720的最大公约数与最小公倍数，并求出所得出的最小公倍数的质因数分解。MATLAB求解语句：7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.12试求出11000间的全部质数得到的所有质数7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.3 MATLAB语言的流程结构循环结构转移结构开关结构试探结构7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.3.1 循环结构循环结构for循环结构 为一个向量，循环变量

15、每次从 向量中取一个数值，执行一次循环体的内容，如此下去，直至执行完 向量中所有的分量7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院while循环结构7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.13用循环求解更简单的语句7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.14用循环求解求最小的 m，使下式成立 使用while结构 这里不宜应用for结构 7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.15求如下值使用MATLAB的循环结构和向量

16、化运算方法7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.3.2 转移结构转移结构7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.16用for循环求解最大的 m，使下式成立MATLAB代码远比while结构繁琐7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.3.3 开关结构开关结构7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院7/5/2021Slide 1 (of 11)MATL

17、AB求解东北大学信息学院2.3.4 试探结构试探结构全新结构优点：设置错误陷阱提高算法的执行效率7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.4 函数编写与调试函数编写与调试MATLAB 语言函数的基本结构可变输入输出个数的处理inline函数与匿名函数7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.4.1 MATLAB语言函数的语言函数的基本结构基本结构 nargin, nargout, varargin, varargout7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.17为什么

18、需要函数？问题M-脚本以M-文件的形式被存取若最大值不为 10000， M-文件就需要被修改则需要一种新的文件形式7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.18按例2.17的要求编写一个MATLAB函数M-文件如下举例优点：无需修改程序7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.19编写一个函数生成 Hilbert矩阵要求若只给出一个输入参数，则会自动生成一个方阵在函数中给出合适的帮助信息检测输入和返回变量的个数7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院MATLAB函数

19、如下7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院在线帮助命令生成Hilbert矩阵7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.20利用函数的递归调用，计算阶乘递归调用函数7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院计算11！其他函数7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.21试比较递归算法和循环算法在Fibonacci数列中应用的优劣Fibonacci数列MATLAB代码7/5/2021Slide 1 (of 11)MATLAB求解东北

20、大学信息学院第25项用 for循环结构计算第100项注意：本例不宜采用递归结构7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.4.2 可变输入输出可变输入输出个数的处理个数的处理vararginvarargoutvarargin1, varargin2, vararginn7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.22conv()可以计算两个多项式的积，时使用varargin实现任意多个多项式的积MATLAB代码7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院函数调用7/5/2021S

21、lide 1 (of 11)MATLAB求解东北大学信息学院2.4.3 inline 函数与匿名函数函数与匿名函数inline 函数可以免去编辑M-文件MATLAB7.0以上版本提供的匿名函数在MATLAB工作空间内已经定义的变量可以直接使用7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.5 二维图形绘制二维图形绘制二维图形绘制基本语句其他二维图形绘制语句隐函数绘制及应用图形修饰7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.5.1 二维图形绘制基本语句二维图形绘制基本语句两个序列 和 构造向量:利用这些数据绘图7/5/202

22、1Slide 1 (of 11)MATLAB求解东北大学信息学院其他方法 仍为向量，而 为矩阵，亦即 和 均为矩阵，且 和 的行和列数均相同假设有多对这样的向量或矩阵7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院改变曲线性质选项列表7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.23绘制方程MATLAB代码问题：如何检查曲线？7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院改用不同的步长，观察曲线是否相同换用较小步长在 附近用变步长算法7/5/2021Slide 1 (of 11)MATLAB

23、求解东北大学信息学院例例 2.24绘制饱和函数方程MATLAB绘图语句更简单的命令7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院图像修饰图像修饰每一个窗口、曲线和坐标轴都是一个对象，对象的属性可以通过函数set()来设置，也可以通过函数get()来获取。7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.5.2 其他二维图形绘制语句其他二维图形绘制语句7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.25绘制极坐标函数 和绘制曲线7/5/2021Slide 1 (of 11)MATLAB求解东北

24、大学信息学院例例 2.26以正弦数据为例，试在同一窗口的不同区域用不同的绘图方式绘制出相应的曲线MATLAB求解语句：7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院下述函数可以用于将一个绘图窗口分割成几个不同的区域n和m分别是行数和列数， k则指示区域的序号7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.5.3 隐函数绘制及应用隐函数绘制及应用隐函数绘图语句默认区域是-2p,2p其他语法7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.27试绘制隐函数MATLAB语句扩大范围7/5/2021

25、Slide 1 (of 11)MATLAB求解东北大学信息学院2.5.4 图形修饰图形修饰直接采用工具栏文字修饰 特殊符号表下标分别用 和 _ 表示 优越性MATLAB 7.0 的新功能7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.6 三维图形表示三维图形表示三维曲线绘制三维曲面绘制三维图形视角设置7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.6.1 三维曲线绘制三维曲线绘制三维曲线绘制stem3, fill3, bar3等7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院例例 2.28试绘制参数方程其中，MATLAB绘图7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院使用stem() 函数3D绘图和图形旋转7/5/2021Slide 1 (of 11)MATLAB求解东北大学信息学院2.6.2 三维曲面绘制三维曲面绘制一般曲面绘制其他函数sur