[工学]第五章 MATLAB程序设计.pptVIP

第五章 MATLAB程序设计 5.1 M文件介绍 5.2 程序控制结构及语句 5.3 变量作用域 5.4 文件操作 5.5 面向对象的程序设计 5.6 程序设计辅助函数 5.7 Matlab编译 5.8 程序调试 5.9 程序设计优化 5.10 应用程序接口 5.1 M文件介绍 5.1.1 M文件的创建与编辑 5.1.2 M文件的分类 5.2 程序控制结构及语句 5.2.1 常见程序控制流语句 5.2.2 分支结构 5.2.3 循环结构 5.3 变量作用域 5.4 文件操作 5.5 面向对象程序设计 5.6 程序设计辅助函数 5.7 MATLAB编译 5.8 程序调试 5.9 程序设计优化 5.10 应用程序接口 例5.2.2.3-1 矩阵乘法运算要求两矩阵的维数相容，否则会出错。先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘。 程序如下： A=[1,2,3;4,5,6]; B=[7,8,9;10,11,12]; try C=A*B; catch C=A.*B; end C lasterr %显示出错原因 5.2.3.1 for语句 语句格式为： for 循环变量=表达式1:表达式2:表达式3 循环体语句 end 首先计算三个表达式的值，再将表达式1的值赋给循环变量，如果此时循环变量的值介于表达式1和表达式3的值之间，则执行循环体语句，否则结束循环的执行。执行完一次循环之后，循环变量自增一个表达式2的值，然后再判断循环变量的值是否介于表达式1和表达式3之间，如果满足仍然执行循环体，直至不满足为止。这时将结束for语句的执行，而继续执行for语句后面的语句。 注意： for 语句一定要有end结束标志； 循环语句中要用“；”防止中间结果的输出； 书写时采用锯齿结构以增加可读性； 例5.2.3.1-1 已知y的表达式，当n=100时，求y的值。 程序如下： y=0;n=100; for i=1:n y=y+1/i/i; end y 在实际MATLAB编程中，为提高程序的执行速度，常用向量运算来代替循环操作： n=100; i=1:n; f=1./i.^2; y=sum(f) 注意： for 语句一定要有end结束标志； 循环语句中要用“；”防止中间结果的输出； 书写时采用锯齿结构以增加可读性； 例5.2.3.1-1 已知y的表达式 程序如下： y=0;n=100; for i=1:n y=y+1/i/i; end y 当n=100时，求y的值。 例5.2.3.1-2 求定积分 程序如下： a=0;b=3*pi;n=1000; h=(b-a)/n; x=a:h:b; f=exp(-0.5*x).*sin(x+pi/6); for i=1:n s(i)= (f(i)+f(i+1))*h/2; end s=sum(s) 事实上，MATLAB提供了有关数值积分的标准函数，实际应用中可直接调用这些函数求数值积分。 for语句更一般的格式： for 循环变量=矩阵表达式 循环体语句 end 执行过程是依次将矩阵的各整列元素赋给循环变量，然后执行循环体语句，直至各列元素处理完毕。实际上，“表达式1:表达式2:表达式3”是一个仅为一行的矩阵(行向量)，因而列向量是单个数据。 例5.2.3.1-3 已知5个学生4门功课的成绩，求每名学生的总成绩。 程序如下： s=0; a=[65,76,56,78;98,83,74,85;76,67,78,79;98,58,42,73;67,89,76,87]; for k=a s=s+k; end disp(s); 5.2.3.2 while语句 while语句的一般格式为： while (条件) 循环体语句 end 其执行过程为：若条件成立，则执行循环体语句，执行后再判断条件是否成立，如果不成立则跳出循环。 条件为真？ 循环体语句 是 否 结束 循环 y=0; i=1; while 1 f=1/i/i; y=y+f; if i==100 break; end