第3章--MATLAB程序设计.pptVIP

3.3.4 全局变量与局部变量 全局变量用global命令定义，格式为： global 变量名 例3-18 全局变量应用示例。 先建立函数文件wadd.m，该函数将输入的参数加权相加。 function f=wadd(x,y) global ALPHA BETA f=ALPHA*x+BETA*y; 在命令窗口中输入： global ALPHA BETA ALPHA=1; BETA=2; s=wadd(1,2) 3.4 程序举例 例3-19 猜数游戏。首先由计算机产生[1,100]之间的随机整数，然后由用户猜测所产生的随机数。根据用户猜测的情况给出不同提示，如猜测的数大于产生的数，则显示“High”，小于则显示“Low”，等于则显示“You won”，同时退出游戏。用户最多可以猜7次。 例3-20 用筛选法求某自然数范围内的全部素数。 素数是大于1，且除了1和它本身以外，不能被其他任何整数所整除的整数。用筛选法求素数的基本思想是：要找出2~m之间的全部素数，首先在2~m中划去2的倍数(不包括2)，然后划去3的倍数(不包括3)，由于4已被划去，再找5的倍数 (不包括5)，…，直到再划去不超过的数的倍数，剩下的数都是素数。 例3-21 设，求s=。 求函数f(x)在[a,b]上的定积分，其几何意义就是求曲线y=f(x)与直线x=a，x=b，y=0所围成的曲边梯形的面积。为了求得曲边梯形面积，先将积分区间[a,b]分成n等分，每个区间的宽度为h=(b-a)/n，对应地将曲边梯形分成n等分，每个小部分即是一个小曲边梯形。近似求出每个小曲边梯形面积，然后将n个小曲边梯形的面积加起来，就得到总面积，即定积分的近似值。近似地求每个小曲边梯形的面积，常用的方法有：矩形法、梯形法以及辛普生法等。 例3-22 Fibonacci数列定义如下： f1=1 f2=1 fn=fn-1+fn-2 (n2) 求Fibonacci数列的第20项。 例3-23 根据矩阵指数的幂级数展开式求矩阵指数。 五、matlab程序流的控制 matlab的控制语句同c语言有相似之处, 但没有c 语言复杂、灵活和多变。因而语法比较简单，容易掌握。 echo 指令— 用来控制m文件在执行过程中是否 显示 echo on — 打开所有命令文件的显示方式 echo off — 关闭所有命令文件的显示方式 echo — 在以上两者间切换 echo对于命令文件函数文件有所不同，命令文件用法简单，函数对所有命令起作用；函数文件用法较复杂，具体请help查询。 echo file on — 打开file函数文件的显示方式 echo file off — 关闭file函数文件的显示方式 echo file — 切换file函数文件的显示方式 echo on all — 打开所有函数文件的显示方式 echo off all —关闭所有函数文件的显示方式 input —— 提示用户从键盘输入数值、字符串、表达式。 例：n=input(How many apples= ) How many apples= 5 n = 5 pause —— 暂停，等待用户响应 可令m文件执行暂停，并在继续执行前按任一个键 。 keyboard —— 与input类似，功能更强 keyboard 调用键盘与调用m文件一样，当keyboard在m文件中出现时，它停止文件的执行并控制键盘，用于修改程序。 break —— 中断 break 中断for, while循环语句的执行，在嵌套循环结构中， break进从最里层循环推出。 menu —— 产生一个拥护可选择输入的菜单 格式：k=menu( ‘标题’，‘菜单1’，‘菜单2’ …） 六、文件操作 matlab文件输入、输出（I/O）函数允许把其它格式的数据文件直接读进matlab，它是以c语言文件I/O函数为基础的。 打开和关闭文件（fopen, fclose) 格式：fid=fopen(‘数据文件’，‘r’) ‘r’ —— 读 ‘w’ —— 写 ‘a’ —— 附加 ‘rt’ —— 读写 文件一旦打开，即可进行读写操作。 sta=fclose(fid) —— 关闭上述打开文件 sta=fclose(all) —— 关闭所有打开的文件