MATLAB7.X程序设计 第3章 MATLAB 程序设计(2)

《MATLAB7.X程序设计 第3章 MATLAB 程序设计(2)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第3章MATLAB程序设计(2)3.5循环语句3.5.1while循环while循环的一般形式：while条件表达式语句组；end2MATLAB7.X程序设计例3-16从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。3MATLAB7.X程序设计例3-17已知（1）n为任意给定的整数，求y的值。（2）y≥1.5时，求n的值。4MATLAB7.X程序设计3.5.2for循环for循环变量=表达式1：表达式2：表达式3循环体语句endfor语句的格式为：5MATLAB7.X程序设计for语句的更一般的格

2、式为：for循环变量=矩阵表达式循环体语句end执行时依次将矩阵的各列元素赋给循环变量，然后执行循环体语句，直至各列元素处理完毕。6MATLAB7.X程序设计例3-18已知当m=100时，求y的值。7MATLAB7.X程序设计例3-19设求8MATLAB7.X程序设计例3-20求[100，200]之间第一个能被13整除的整数。例3-21使用for循环语句逐一显示构架数组的域值。例3-22输出全部三位的水仙花数。例3-23编程求出Fibonacci数列中第一个大于10000的元素及序号。9MATLAB7.X程序设计例3

3、-24编程列出n阶Hilbert矩阵。例3-25求某正整数范围内的全部素数。10MATLAB7.X程序设计3.5.3循环程序的等价形式1数组运算的应用MATLAB的基本数据单元为数组，定义了矩阵运算和数组运算，循环程序在一些特殊情况下可以用数组运算替代。11MATLAB7.X程序设计例3-26求1-100的平方、平方根和立方根。方法1：square=zeros(1,100);square_root=zeros(1,100);cube_root=zeros(1,100);forii=1:100square(ii)=ii

4、^2;square_root(ii)=ii^(1/2);cube_root(ii)=ii^(1/3);end方法2：ii=1:100;square(ii)=ii.^2;square_root(ii)=ii.^(1/2);cube_root(ii)=ii.^(1/3);12MATLAB7.X程序设计2逻辑数组的应用下面介绍逻辑数组和循环程序等价方面的应用。13MATLAB7.X程序设计例3-27编程求3阶魔方矩阵大于5的元素的平方根。方法1：a=magic(3);forii=1:size(a,1)forjj=1:siz

5、e(a,2)ifa(ii,jj)>5a(ii,jj)=sqrt(a(ii,jj));endendenda14MATLAB7.X程序设计方法2：a=magic(3);b=a>5;%b是逻辑数组a(b)=sqrt(a(b))分析：方法1和方法2程序的执行结果相同，b=a>5和a(b)=sqrt(a(b))语句等价于二重循环。15MATLAB7.X程序设计如果要求同时对小于等于5的元素求平方呢？a=magic(3);b=a>5;a(b)=sqrt(a(b));a(~b)=a(~b).^2;a16MATLAB7.X程序设计例

6、3-28创建向量1，2，...,10000，求其中大于5000的所有元素的平方根。要求证明用逻辑数组方法执行时间短。17MATLAB7.X程序设计3.6函数文件3.6.1函数文件的格式函数文件是一种特殊类型的运行在自己独立的工作空间的M-文件，通过在输入形参表接收数据，通过输出形参表返回结果到调用文件。函数文件的一般格式为：function输出形参表=函数名（输入形参表）注释行…函数体语句…(return)18MATLAB7.X程序设计说明：1）函数说明行2）H1行3）helptext4）编写和修改记录5）函数体注意

7、：1)、2）、3）和4）是注释行19MATLAB7.X程序设计例3-29编写一个函数文件，求Fibonacci数组中第一个大于整数n的元素及序号。定义函数文件f0：functionf0(n)a(1)=1;a(2)=1;i=2;whilea(i)<=na(i+1)=a(i-1)+a(i);i=i+1;endi,a(i)调用时在命令窗口输入：f0(10000)20MATLAB7.X程序设计例3-30求n！,n=1,2,…,10。定义函数文件ffactor：functionf=ffactor(n)ifn<=1f=1;els

8、ef=ffactor(n-1)*n;endreturn;%返回21MATLAB7.X程序设计编写下述的命令文件中调用函数文件ffactor.mfori=1:10fac(i)=ffactor(i);endfac22MATLAB7.X程序设计3.6.2全局变量和临时变量1全局变量globalXYZ2局部变量presistentXYZ23MATLAB