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

《高等应用数学问题的MATLAB求解.ppt》由会员分享，可在线阅读，更多相关《高等应用数学问题的MATLAB求解.ppt（99页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第2章MATLAB语言程序设计基础高等应用数学问题的MATLAB求解清华大学出版社2008CAI课件开发：薛定宇、刘莹莹、董雯彬5/28/2023Slide 1(of 99)东北大学信息学院第2章 MATLAB语言程序设计基础MATLAB 程序设计语言基础基本数学运算MATLAB语言的流程结构函数的编写与调试二维图形绘制三维图形绘制图像处理简介5/28/2023Slide 2(of 99)MATLAB语言有如下的优点MATLAB 语言的简洁高效性MATLAB 语言的科学运算功能MATLAB 语言的绘图功能MATLAB 庞大的工具箱与模块集MATLAB 强大的动态系统仿真功能MATLAB 语言是

2、当前国际上自动控制领域的首选计算机语言，也是很多理工科专业最适合的计算机数学语言，本课程以 MATLAB 为主要计算机语言。5/28/2023Slide 3(of 99)2.1 MATLAB 程序设计语言基础MATLAB语言的变量与常量数据结构MATLAB的基本语句结构冒号表达式与子矩阵提取5/28/2023Slide 4(of 99)2.1.1 MATLAB语言的变量与常量MATLAB语言的变量名规则由一个字母引导，后面可以为其他字符区分大小写Abc ABc有效 MYvar12,MY_Var12 和 MyVar12_错误的变量名 12MyVar,_MyVar12 MATLAB的保留常量eps

3、,i,j,pi,NaN,Inf,i=sqrt(-1)lastwarn,lasterr5/28/2023Slide 5(of 99)2.1.2 数据结构数值型数据符号型其他数据结构5/28/2023Slide 6(of 99)2.1.2.1 数值型数据双精度数值变量IEEE标准，64位(8字节)，11指数位，53数值位和一个符号位 至 double()函数的转换其他数据类型uint8，常用于图像表示和处理，8位int8(),int16(),int32(),uint16(),uint32()5/28/2023Slide 7(of 99)2.1.2.2 符号型符号型，sym(A)，常用于公式推导常用

4、于公式推导和求解析解变量声明显示符号变量的任何精度默认精度：32位小数5/28/2023Slide 8(of 99)例 2.1求出 p 的 300 位有效数字MATLAB代码可以任意增加有效数字的位数随着有效数字位数的增加，程序运行的时间变长5/28/2023Slide 9(of 99)2.1.2.3 其他数据结构字符串型数据：用单引号括起来 多维数组：是矩阵的直接扩展，多个下标单元数组：将不同类型数据集成到一个变量名下面，用 表示类与对象：可以定义重载函数5/28/2023Slide 10(of 99)2.1.3 MATLAB的基本语句结构直接赋值语句接赋值语句的基本结构赋值语句的末尾加一个

5、分号可以阻止显示运算结果保留变量：ans存放最近一次无赋值变量语句的运算结果5/28/2023Slide 11(of 99)例 2.2表示矩阵MATLAB代码其他命令5/28/2023Slide 12(of 99)例 2.3试输入复数矩阵MATLAB代码需要避免的语句5/28/2023Slide 13(of 99)函数调用语句函数调用举例函数可以通过不同的方式被调用内在函数，*.m函数匿名函数，内联inline函数重载函数5/28/2023Slide 14(of 99)2.1.4 冒号表达式与子矩阵提取冒号表达式是定义行向量的有效方法开始于值，步长为值，终止于值.默认步长值为15/28/202

6、3Slide 15(of 99)例 2.4用不同的步距生 间的向量5/28/2023Slide 16(of 99)子矩阵提取的基本语句格式 表示子矩阵要保留的行号构成的向量 表示要保留的列号构成的向量:，表示要提取所有的行或列，取决于它的位置5/28/2023Slide 17(of 99)例 2.5矩阵A的不同子矩阵提取MATLAB代码5/28/2023Slide 18(of 99)2.2 基本数学运算矩阵的代数运算矩阵的逻辑运算矩阵的比较运算解析结果的化简与变换基本数论运算5/28/2023Slide 19(of 99)矩阵的代数运算矩阵转置加减法运算矩阵乘法矩阵的除法矩阵翻转矩阵乘方运算点

7、运算5/28/2023Slide 20(of 99)矩阵转置矩阵表示矩阵A，n 行m列，被称作 矩阵厄米转置一般转置5/28/2023Slide 21(of 99)加减法运算数学表示用C语言编程难于实现，同A*BMATLAB语法注意：任一个变量可以为标量如果矩阵维数不匹配，系统会报错5/28/2023Slide 22(of 99)矩阵乘法数学表示：MATLAB表示：注意：系统自动检测矩阵维数是否匹配5/28/2023Slide 23(of 99)矩阵的除法矩阵左除：求解线性方程组：MATLAB解法：最小二乘解若A为非奇异方阵，则5/28/2023Slide 24(of 99)矩阵右除：求解线性

8、方程组：MATLAB解法：最小二乘解若A为非奇异方阵，则更精确的解法5/28/2023Slide 25(of 99)矩阵翻转左右翻转 上下翻转旋转 90o如何旋转180o？5/28/2023Slide 26(of 99)矩阵乘方A为方阵，求矩阵A的x次幂数学描述MATLAB命令5/28/2023Slide 27(of 99)点运算矩阵对应元素的直接运算例如：为 另一个例子：为5/28/2023Slide 28(of 99)例 2.6求出矩阵A的全部三次方根，并检验结果MATLAB代码另两个根5/28/2023Slide 29(of 99)2.2.2 矩阵的逻辑运算逻辑变量：当前版本有逻辑变量非

9、 0 表示逻辑 1逻辑运算（相应元素间的运算）与运算或运算非运算异或运算5/28/2023Slide 30(of 99)2.2.3 矩阵的比较运算各种允许的比较关系,=,=,=,=,find(),all(),any()实例5/28/2023Slide 31(of 99)2.2.4 解析结果的化简与变换函数 simple()用于数学公式的化简其他常用化简函数numden(),collect(),expand(),factor()5/28/2023Slide 32(of 99)例 2.7化简多项式用不同的函数求解5/28/2023Slide 33(of 99)变量替换变量替换语句该函数执行点运算转

10、换成LaTeX表示5/28/2023Slide 34(of 99)例 2.8试用 对下式进行双线性变换，得出其LaTeX排版格式MATLAB代码5/28/2023Slide 35(of 99)结果：8,frac left(2,z+1 right)2z left(3,z+1 right)left(5,z+3 right)3 left(z+1 right)7LaTeX排版显示5/28/2023Slide 36(of 99)基本数论运算5/28/2023Slide 37(of 99)例 2.9运用各种函数，对下面的数据进行取整运算-0.2765,0.5772,1.4597,2.1091,1.191,

11、-1.6187MATLAB代码：5/28/2023Slide 38(of 99)例 2.10假设 的Hilbert矩阵可以由 定义，试对其进行有理数变换。结果：5/28/2023Slide 39(of 99)例 2.11试求1856120和1483720的最大公约数与最小公倍数，并求出所得出的最小公倍数的质因数分解。MATLAB求解语句：5/28/2023Slide 40(of 99)例 2.12试求出11000间的全部质数得到的所有质数5/28/2023Slide 41(of 99)2.3 MATLAB语言的流程结构循环结构转移结构开关结构试探结构5/28/2023Slide 42(of 9

12、9)2.3.1 循环结构for循环结构 为一个向量，循环变量 每次从 向量中取一个数值，执行一次循环体的内容，如此下去，直至执行完 向量中所有的分量5/28/2023Slide 43(of 99)while循环结构5/28/2023Slide 44(of 99)例 2.13用循环求解更简单的语句5/28/2023Slide 45(of 99)例 2.14用循环求解求最小的 m，使下式成立 使用while结构 这里不宜应用for结构 5/28/2023Slide 46(of 99)例 2.15求如下值使用MATLAB的循环结构和向量化运算方法5/28/2023Slide 47(of 99)2.3

13、.2 转移结构5/28/2023Slide 48(of 99)5/28/2023Slide 49(of 99)例 2.16用for循环求解最大的 m，使下式成立MATLAB代码远比while结构繁琐5/28/2023Slide 50(of 99)2.3.3 开关结构5/28/2023Slide 51(of 99)5/28/2023Slide 52(of 99)2.3.4 试探结构全新结构优点：设置错误陷阱提高算法的执行效率5/28/2023Slide 53(of 99)2.4 函数编写与调试MATLAB 语言函数的基本结构可变输入输出个数的处理inline函数与匿名函数5/28/2023Sli

14、de 54(of 99)2.4.1 MATLAB语言函数的基本结构 nargin,nargout,varargin,varargout5/28/2023Slide 55(of 99)例 2.17为什么需要函数？问题M-脚本以M-文件的形式被存取若最大值不为 10000，M-文件就需要被修改则需要一种新的文件形式5/28/2023Slide 56(of 99)例 2.18按例2.17的要求编写一个MATLAB函数M-文件如下举例优点：无需修改程序5/28/2023Slide 57(of 99)例 2.19编写一个函数生成 Hilbert矩阵要求若只给出一个输入参数，则会自动生成一个方阵在函数中给

15、出合适的帮助信息检测输入和返回变量的个数5/28/2023Slide 58(of 99)MATLAB函数如下5/28/2023Slide 59(of 99)在线帮助命令生成Hilbert矩阵5/28/2023Slide 60(of 99)例 2.20利用函数的递归调用，计算阶乘递归调用函数5/28/2023Slide 61(of 99)计算11！其他函数5/28/2023Slide 62(of 99)例 2.21试比较递归算法和循环算法在Fibonacci数列中应用的优劣Fibonacci数列MATLAB代码5/28/2023Slide 63(of 99)第25项用 for循环结构计算第100

16、项注意：本例不宜采用递归结构5/28/2023Slide 64(of 99)2.4.2 可变输入输出个数的处理vararginvarargoutvarargin1,varargin2,vararginn5/28/2023Slide 65(of 99)例 2.22conv()可以计算两个多项式的积，时使用varargin实现任意多个多项式的积MATLAB代码5/28/2023Slide 66(of 99)函数调用5/28/2023Slide 67(of 99)2.4.3 inline 函数与匿名函数inline 函数可以免去编辑M-文件MATLAB7.0以上版本提供的匿名函数在MATLAB工作空

17、间内已经定义的变量可以直接使用5/28/2023Slide 68(of 99)2.5 二维图形绘制二维图形绘制基本语句其他二维图形绘制语句隐函数绘制及应用图形修饰5/28/2023Slide 69(of 99)2.5.1 二维图形绘制基本语句两个序列 和 构造向量:利用这些数据绘图5/28/2023Slide 70(of 99)其他方法 仍为向量，而 为矩阵，亦即 和 均为矩阵，且 和 的行和列数均相同假设有多对这样的向量或矩阵5/28/2023Slide 71(of 99)改变曲线性质选项列表5/28/2023Slide 72(of 99)例 2.23绘制方程MATLAB代码问题：如何检查曲

18、线？5/28/2023Slide 73(of 99)改用不同的步长，观察曲线是否相同换用较小步长在 附近用变步长算法5/28/2023Slide 74(of 99)例 2.24绘制饱和函数方程MATLAB绘图语句更简单的命令5/28/2023Slide 75(of 99)图像修饰每一个窗口、曲线和坐标轴都是一个对象，对象的属性可以通过函数set()来设置，也可以通过函数get()来获取。5/28/2023Slide 76(of 99)2.5.2 其他二维图形绘制语句5/28/2023Slide 77(of 99)例 2.25绘制极坐标函数 和绘制曲线5/28/2023Slide 78(of 9

19、9)例 2.26以正弦数据为例，试在同一窗口的不同区域用不同的绘图方式绘制出相应的曲线MATLAB求解语句：5/28/2023Slide 79(of 99)下述函数可以用于将一个绘图窗口分割成几个不同的区域n和m分别是行数和列数，k则指示区域的序号5/28/2023Slide 80(of 99)2.5.3 隐函数绘制及应用隐函数绘图语句默认区域是-2 p,2p其他语法5/28/2023Slide 81(of 99)例 2.27试绘制隐函数MATLAB语句扩大范围5/28/2023Slide 82(of 99)2.5.4 图形修饰直接采用工具栏文字修饰 特殊符号表下标分别用 和 _ 表示 优越性

20、MATLAB 7.0 的新功能5/28/2023Slide 83(of 99)2.6 三维图形表示三维曲线绘制三维曲面绘制三维图形视角设置5/28/2023Slide 84(of 99)2.6.1 三维曲线绘制三维曲线绘制stem3,fill3,bar3等5/28/2023Slide 85(of 99)例 2.28试绘制参数方程其中，MATLAB绘图5/28/2023Slide 86(of 99)使用stem()函数3D绘图和图形旋转5/28/2023Slide 87(of 99)2.6.2 三维曲面绘制一般曲面绘制其他函数surfl(),surfc()等高线绘制 contour(),cont

21、ours()5/28/2023Slide 88(of 99)例 2.29给出二元函数如下，绘制3D图像MATLAB绘图表面图5/28/2023Slide 89(of 99)例 2.30试绘制出二元函数绘制3D图像5/28/2023Slide 90(of 99)变步长方法5/28/2023Slide 91(of 99)例 2.31绘制如下分段函数的三维图像5/28/2023Slide 92(of 99)MATLAB绘图5/28/2023Slide 93(of 99)2.6.3 三维图形视角设置两种方法改变图形的视角直接采用工具栏命令语句view()定义为方位角，定义为仰角5/28/2023Sli

22、de 94(of 99)视点定义5/28/2023Slide 95(of 99)例 2.32对函数 从不同的方向上观察其图形MATLAB代码5/28/2023Slide 96(of 99)2.7 图像处理简介图像处理工具箱提供了图像处理的强大功能读图像文件：imread()图像显示：image()、imview()、imshow()和imtool()图像颜色转换：rgb2gray()、rgb2hsv()和hsv2rgb()图像边缘提取：edge()5/28/2023Slide 97(of 99)例 2.33试读取tiantan.jpg文件中的图像并提取其边缘tiantan.jpg文件5/28/2023Slide 98(of 99)MATLAB代码结果5/28/2023Slide 99(of 99)