MATLAB程序设计课件.ppt

《MATLAB程序设计课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第5章MATLAB程序设计5.1程序流程控制5.2M文件5.3函数调用和参数传递5.4利用函数句柄执行函数5.5利用泛函命令进行数值分析5.6内联函数5.7M文件性能的优化和加速5.1程序流程控制5.1.1for…end循环结构MATLAB的循环结构有2种：for…end结构和while…end结构。这2种语句结构不完全相同，各有各的特点。语法：for循环变量=array循环体end说明：循环体被循环执行，执行的次数就是array的列数，array可以是向量也可以是矩阵，循环变量依次取array的各列，每取

2、1次循环体执行1次。【例5.1】使用for…end循环编程求出1+3+5…+99的值。sum=0;forn=1:2:100sum=sum+nend计算的结果为：sum=2500。5.1.1for…end循环结构【例5.2】使用for…end循环将单位阵转换为列向量。sum=zeros(6,1);forn=eye(6,6)sum=sum+nend计算结果如下。sum=1111115.1.2while…end循环结构for…end循环的循环次数确定，而while…end循环的循环次数不确定。语法：while表达

3、式循环体end说明：只要表达式为逻辑真，就执行循环体；一旦表达式为逻辑假，就结束循环。表达式可以是向量也可以是矩阵。如果表达式为矩阵则当所有的元素都为真时才执行循环体，如果表达式为nan，MATLAB认为是假，就不执行循环体。5.1.2while…end循环结构【例5.3】根据，求y<3时的最大n值和y值。y=0;n=1;whiley<3y=y+1/(2*n-1);n=n+1;z(n)=y;endmn=n-2%y<3之前的nmy=z(n-1)计算结果如下。mn=56my=2.99445.1.3if…else

4、…end条件转移结构if…else…end结构是最常见的条件转移结构。语法：if条件式1语句段1elseif条件式2语句段2…else语句段n+1end5.1.3if…else…end条件转移结构【例5.4】根据不同的分段表达式，绘制分段函数曲线，曲线如图5.1所示。5.1.3if…else…end条件转移结构x=0:0.5:10;y=zeros(1,length(x));%产生0行向量，y的初始值为0forn=1:length(x)ifx(n)>=8y(n)=1;elseifx(n)>=6y(n)=5-x

5、(n)/2;elseifx(n)>=4y(n)=2;;elsey(n)=sqrt(x(n));endendplot(x,y)axis([01002.5]);5.1.4switch…case开关结构switch…case结构是有多个分支结构的条件转移结构。语法：switch开关表达式case表达式1语句段1case表达式2语句段2…otherwise语句段nend5.1.4switch…case开关结构【例5.5】用switch…case开关结构得出各月份的季节。formonth=1:12;switchmon

6、thcase{3,4,5}season='spring'case{6,7,8}season='summer'case{9,10,11}season='autumn'otherwiseseason='winter'endend5.1.5try…catch…end试探结构MATLAB提供一种试探结构try…catch…end，这种语句结构是其他很多语言所没有的。语法：try语句段1catch语句段2end说明：首先试探性地执行语句段1，如果在此段语句执行过程中出现错误，则将错误信息赋给保留的lasterr变量，

7、并放弃该段语句，转而执行语句段2中的语句；当执行语句段2又出现错误，则终止该结构。5.1.5try…catch…end试探结构【例5.6】用try…catch…end结构进行矩阵相乘运算。n=4;a=magic(n);m=3;b=eye(3);tryc=a*bcatchc=a(1:m,1:m)*bendlasterr计算结果如下。c=162351110976用lasterr函数查看出错原因显示如下。ans=Errorusing==>*Innermatrixdimensionsmustagree.5.1.6流

8、程控制语句1．break命令break命令可以使包含break的最内层的for或while语句强制终止，立即跳出该结构，执行end后面的命令，break命令一般和if结构结合使用。【例5.7】将【例5.3】增加条件，用if与break命令结合，停止while循环。计算值，当y>=3时终止计算。y=0;n=1;whilen<=100ify<3y=y+1/(2*n-1);n=n+1;z(n)=y;elsebreake