1、程序设计补充材料第一节 M 文件1、M 文件大量的 Matlab 代码（指令）罗列成文件，送 MATLAB 系统中解释运行。2、M 文件的分类（1）命令式文件（脚本 script）（2）函数式文件（function）3、两类文件的区别（1）命令式文件没有输出输入参量，函数式文件有（2）命令式文件运行中可以调用工作空间的数据，运行中产生的所有变量为全局变量。函数式文件不能调用工作空间的数据，运行中产生的所有变量为局部变量（3）命令式文件可以直接运行，函数式文件不能，只能调用4、函数式文件（function）的构成（1）函数定义行Function输出参量=gauss(输入参量)（2）函数体完成函数
2、的功能（3）函数说明（4）函数行注5、例子（funave1） （ funave2）funave1：functionaverage=funave(vectory);average=sum(sum(vectory)/length(vectory);funave2:(定义二维数组时输出的平均值)functionave1,ave2=funave2(vectory1,vectory2)if nargin=1ave1=sum(sum(vectory1)/length(vectory1);endif nargin=2ave1=sum(sum(vectory1)/length(vectory1);ave2=s
3、um(sum(vectory2)/length(vectory2);end 第二节 程序的流程控制1.关系运算符= 等于 =大于等于 大于2.条件分支语句（if-else-end）如果-那么-否则（1）结构If.else if.else.end注意这一结构的条件优先问题（2）条件的嵌套。(exlin.m)exlin:function f=exlin(x)if x 时 return5.变量的鼠标观测法6.代码运行的计时方法（1）整段程序代码的计时tic.toc 表示计算 tic 与 toc 之间的时间（exno57）（2）也可以用etime(t1,t2)来计算 t1,t2 之间的时间差来完成上述
4、功能(exno58)（3）也可以用 cputtime 变量来完成上述功能(exno59)exno57:tica=rand(300);inv(a); %计算逆矩阵toc%注：程序中不需要显示结果的就不显示，可以节省时间exno58:t0=clocka=rand(300);inv(a);%计算逆矩阵elapsed_time=etime(clock,t0)7.程序代码的优化（1）程序代码的向量化 （exno61）（2）用矩阵结构进行运算(exno62)一般情况下，完全采用矩阵运行的方式，MATLAB 的程序与 C 语言基本相同。这必须对矩阵非常熟悉，例如x=1 2 3;1 2 1a=4 5 6希望将
5、 a 中的每一个元素乘以 x 的每一列，用 diag(x)（3）矩阵的预先配置(exno63)(enxo64)exno61:t0=cputimen=200000;total=0;for i=1:ntotal=total 1/iendtotalt1=cputime-t0exno62:ticn=100000;a=1:n;total=sum(1./a)tocexno63:tica=1 2 3;4 5 6;7 8 9for i=1:100y(i)=det(ai)endtocexno64:tica=1 2 3;4 5 6;7 8 9;y=zeros(1,100);%分配内存for i=1:100y(i)=det(ai);endtoc
