第三章 Matlab法度模范设计[宝典].docVIP

第三章 Matlab程序设计 前面介绍了M文件的概念和组成，看到了一些用MATLAB语言编写的简单程序．要想实现更强的功能，需要用到循环控制．几乎所有实用的程序都包含循环，熟练使用MATLAB的循环结构的选择结构是编程的基本要求。MATLAB提供四种循环和选择控制结构，它们是:for循环，while循环，if-else-end结构和switch-case-end结构 们经常出现在M文件中． 运算符 MATLAB的运算符可分为三类：算术运算符、关系运算符和逻辑运算符其中算术运算符的优先级最高，其次是关系运算符，再其次是逻辑运算符。算术运算符在前面已经介绍，这里只介绍关系运算符和逻辑运算符．1．关系运算符 关系运算符对于程序的流程控制非常有用．MATLAB共有六个关系运算符，它们分别是 ＜ 小于；＜ = 小于等于；＞ 大于；＞= 大于等于；﹦﹦ 等于；~ =不等于关系运算符可以比较同型矩阵，此时将生成一个0-l矩阵，当相应元素经关系运算为真时，对应位置上生成1，否则为0关系运算符也可以比较标量和矩阵，此时是标量与矩阵的每个元素分别比较，生成一个0－1矩阵． 2．逻辑运算符 MATLAB共有三个逻辑运算符：与（＆）或（|）非（~）。 对于数值矩阵，当元素为0时，逻辑上为假；当元素为非0时。逻辑上为真．同关系运算符一样，逻辑运算符两端的运算数可以是同型矩阵，对两矩阵的相应元素分别运算，结果为一个0－1矩阵．当逻辑表达式的值为真时。赋值1，否则为0．同样，其中一个矩阵也可以是标量。 与（＆）运算： 两个运算数都为真时，结果为真，其他情况下（一真一假或两个都假）结果为假． 或（1）运算： 两个运算数都为假时，结果为假，其他情况下（一真一假或两个都真）为真． 非（～）运算： 只有一个运算数，当该运算数为真时，结果为假，否则，结果为真． for循环 for循环允许一组命令以固定的和预定的次数重复。for循环的一般形式为： for x=表达式1:表达式2：表达式3 语句体end 其中表达式1的值为循环的初值，表达式2的值为步长，表达式3的值为循环的终值．如果省略表达式2，则默认步长为1．该循环体的执行过程如下： l）将表达式1的值赋给x; 2) 对于正的步长，当x的值大于表达式3的值时，结束循环；对于负的步长，当x的值小于表达式3的值时结束循环．否则，执行for和end之间的语句体，然后执行下面的第3）步． 3）x加上一个步长后，返回第2）步继续执行． 例如，程序： for k＝l：4 X（k ）＝l／k； end format rat ％设置输出格式为有理数 x 将输出： X＝ 1 1/2 1/3 1/4 *************************************************************************** 注意： l) for语句可以嵌套使用； 2）当有一个等效的数组方法来解给定的问题时，应避免用for循环．例如，上例可被重写为 n＝1：4； x＝1／n； format rat；x这种方法执行更快，要求较小的输入． while循环 while循环一般用于事先不能确定循环次数的情况．while循环的一般形式为： While 表达式 语句体 end 只要表达式的值为1(真),就执行while与end之间的语句体，直到表达式的值为0(假)时终止该循环．通常，表达式的值为标量，但对数组值也同样有效，此时，数组的所有元素都为真,才执行while与end之间的语句体．例如，程序： n＝0；EPS=1； while(1+EPS) 1 EPS＝EPS／2；n＝n＋l； end n，EPS＝EPS * 2 运行结果： n＝ 53 EPS＝ 2.2204 e-016 这个例子给出了计算MATLAB的特殊变量eps的一种方法，eps是一个可以加到1，在计算机的有限精度下，而使结果大于1的最小数值．这里我们用大写EPS，以便与eps相区别．EPS从1开始，不断被2除，直到EPS＋l不大于1。因为MATLAB用16位数来表示数据，因此，当EPS接近10-16时，它会认为EPS+1不大于 1，于是w