第6章 MATAB数学实验.docVIP

第6章 MATLAB数学实验 没有大胆的猜测，就做不出伟大的发现。 ——牛顿 MATLAB是Matrix Laboratory的缩写，是一个集数值计算、图形处理、符号运算、文字处理、数学建模、实时控制、动态仿真和信号处理等功能为一体的数学应用软件，而且该系统的基本数据结构是矩阵，又具有数量巨大的内部函数库和多个工具箱 图6.1 命令窗口是和MATLAB编译器连接的主要窗口，“”为运算提示符，表示MTALAB处于准备状态，当在提示符后输入一段正确的运算式时，只需按Enter键，命令窗口中就会直接显示运算结果。如图6.2所示。 图6.2 数据的默认格式为五位有效数字，可用format命令改变输出格式。Help是获取帮助的命令，在它之后应该跟一个主题词。例如，help format，系统就会对format的用法提供说明，因此它对初学者是非常有用的。 当需要编写比较复杂的程序时，就需要用到M文件，而只需要将所有命令按顺序放到一个扩展名.m的文本文件下，每次运行只需输入该M文件的文件名即可。下面介绍MATLAB程序设计中常用的程序控制语句和命令。 ●顺序结构是最简单的程序结构，在编写程序后，系统将按照程序的物理位置顺序执行。 ●选择语句。在编写程序时，往往需要根据一定的条件来执行不同的语句，此时，需要使用分支语句来控制程序的进程，通常使用if-else-end结构来实现这种控制。if-else-end结构是： if 表达式　　1 　else 　　　2 end 此时，如果表达式为真，则系统执行语句1；如果表达式为假，则系统执行语句2. 例如比较a和b的大小，其中a=40，b=10. 程序设计： clear a=40; b=10; if ab disp(‘ba’) else disp(‘ab’) end 运行结果： ab ●分支语句。另外，MATLAB语句中还提供了switch-case-otherwise-end分支语句，其使用格式如下： switch 开关语句 　 case 　 case …） 执行语句 ，...，执行语句 ……. otherwise 执行语句 ，...，执行语句 　　　nd 在上面的分支语句中，当某个条件语句的内容与开关语句的内容相匹配时，系统将执行其后的语句，如果所有的条件语句与与开关条件都不相符合时，系统将执行otherwise后的语句。 ●循环语句。当遇到许多有规律的重复运算时，可以方便地使用以下两种循环语句。 ◆for循环基本格式是： for i= 表达式， 执行语句，…，执行语句 end 上述结构是对循环次数的控制。 例 求1+2+ … +100的值。 程序设计： 》sum=0; 》for i=1:100 Sum=sum+i; end 》sum sum= 5050 for循环可以重复使用，即可以多次嵌套。 ◆while循环的判断控制可以是逻辑判断语句，因此，它的循环次数可以是一个不定数，这样就赋予了它较for循环更广泛的用途。其使用格式如下： while 表达式　　　　　　　　　　　End。 解 在命令窗口中输入 》A=[1 1 2；-1 0 3；4 -5 6] A= 6.1.2 矩阵的结构操作 输入矩阵后，可以对矩阵进行的主要操作包括矩阵的扩充，矩阵元素的提取和变换，矩阵元素的部分删除等。下面对其作简单介绍。 1.矩阵的扩充 例如，用下述命令可以在上述矩阵A下面再加上一个行向量： 》A（4，：）=[1 3 2] A= 下述命令可以在上述矩阵A下面再加上一个列向量： 》A（：，4）=[-1 0 3 2] A= 2．矩阵元素的提取和变换 可以用下述命令提取上述矩阵A的第3行第1列的元素： 》A（3,1） ans= 4 可以用下述命令提取上述矩阵A的第1行第3列的元素： 》A（：，[1，3]） ans= 可以用下述命令提取矩阵的上三角和下三角部分及对角线元素： triu（A）——提取矩阵A的上三角部分； tril（A）——提取矩阵A的下三角部分； diag（A）——提取矩阵A的对角线元素。 可以用下列命令对矩阵进行翻转和旋转： fliplr（A）——矩阵A左右翻转； flipud（A）——矩阵A上下翻转； rot90（A）——矩阵A整体逆时针旋转90度。 3．矩阵元素的删除