matlab学习6-程序设计1.pptVIP

* 第五章 MATLAB程序设计 第一节 M文件 1.M文件 大量的MATLAB代码（指令）罗列成文件，送MATLAB系统中解释运行。M文件的后缀名为M 2.M文件的分类 （1）命令式文件（脚本script）。 （2）函数式文件（function） 3.两文件的区别 （1）命令式文件没有输入输出参量，函数式文件有。 （2）命令式文件运行中可以调用工作空间的数据，运行中产生的所有变量为全局变量。函数式文件不能调用工作空间的数据，运行中产生的所有变量为局部变量。 （3）命令式文件可以直接运行，函数文件不能，只能调用。 4.函数式文件（function）的构成 （1）函数定义行： Function[输出参量]=函数名（输入参量） （2）函数体： 完成函数的功能。 （3）函数说明。%提示 （4）函数行注。%提示 function y=average(x) y=sum(x)/length(x); z=1:99; y=average(z) 5.例子 注：存入当前子目录 average funave2 第二节 程序的流程控制 1.关系运算符 ==等于，=小于等于，~=不等于，=大于等于， 小于，大于 2.条件分支语句（if-else-end） 如果-那么-否则 （1）结构： if...elseif ...else...end 注意这一结构的条件优先问题。 （2）条件的嵌套。 exlin function f=exlin(x) if x0 f=0; elseif x1 f=x; elseif x2 f=2-x; else f=0; end %exlin(-1) %exlin(1) %exlin(0.5) %exlin(1.2) %exlin(3) 3.硬循环语句（for-end） （1）结构：for（循环变量）end （2）例 example2 n=5; %每个for和if都必须对应一个end for i=1:n for j=1:n if i==j a(i,j)=2; else a(i,j)=0; end end end a 4.条件循环语句（while-end） 当-就- （1）结构：while（表达式）end （2）例：求阶乘大于或等于99^99的最小整数。 example3 n=1; while prod(1:n)99^99; %prod()向量元素的乘积 n=n+1; end n 5.分支循环语句（switch-case-otherwise-end） （1）结构：switch（开关量） （2）switch 语句是一个控制语句，它通过将控制传递给其体内的一个 case 语句来处理多个选择 example4 a=input(a=?) switch a case 1 disp(It is raning) case 0 disp(It do not know) case -1 disp(It is not raning) otherwise disp(It is raning?) end 第三节 程序的调试 1.语法和逻辑两种错型 2.根据出错信息调试 命令窗的提示，注意一般情况开始时不加；号调试 减少矩阵的维数来检验程序的维数 3.设置断点来判断 Breakpoints 4.利用keyboard命令来判断 m文件中调试比较方便 run 当出现k 时 return 5.变量的观测法 x=1:0.1:4; y=1:0.1:4 keyboard [X,Y]=meshgrid(x,y); R=(X+Y).^2; Z=sin(R)./R keyboard meshc(Z); axis vis3d off; 6.代码运行的计时方法 （1）整段程序代码的计时 tic... toc...表示计算tic和toc之间的时间 example5 tic a=rand(300); inv(a); %计算逆矩阵 toc %注：程序中不需要显示结果的就不显示 tic n=100000; a=1:n; total=sum(1./a) toc （2）也可以用 etime（t1，t2）来计算t1,t2之间的时间差来完成上述功能 example5 t0=clock %调用时钟函数 a=rand(300); inv(a); time=etime(clock,t0) （3）也可以用cputime变量来完成上述功能 * * * 电子12