第3章_Matlab程序设计_part2讲述.ppt

[例3-14] 通过键盘输入百分制成绩，输出对应的等级（A: 90~100, B: 80~89,C:70~79,D:60~69,E: 60) 程序如下： n=input(‘n=‘); if n0 | n100 display(‘Error!’); else i=fix(n/10); switch i case {9,10} display(‘A’) case 8 display(‘B’) case 7 display(‘C’) case 6 display(‘D’) otherwise display(‘E’) end %switch end % if 例：判断用户选择的是哪种插值方法（’bilinear’, ‘cubic’或’nearest’） method = input(Input method: (bilinear,cubic or nearest), s) switch lower(method) case {linear,bilinear} disp(Method is linear) case cubic disp(Method is cubic) case nearest disp(Method is nearest) otherwise disp(Unknown method.) End 练习：利用if-elseif-else结构重写以上两个程序。 3.4.3 try-catch结构 Matlab从5.3版本开始提供了一种新的试探式语句结构。其一般格式为： try 语句组１ catch 语句组２ end 执行时，首先试探性地执行语句组１。如果在执行语句组１的过程中出错，则将错误信息赋给Matlab预定义变量lasterr后转去执行语句组２（如果执行语句组２时也出错，将终止程序的执行）。 [例]矩阵乘法运算要求两矩阵的维数相容，否则会出错。创建一个命令文件，先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘。 A=input(请输入A矩阵：) B=input(请输入B矩阵：) lasterr(); %清除原有的错误信息 try C=A*B; catch C=A.*B; end C disp(lasterr) %显示出错原因 运行程序，观察程序的执行结果。 输入矩阵A：[1 2 3;4 5 6] 输入矩阵B ：[7 8 9;10 11 12] 运行结果：（略） 再运行一次，输入矩阵A：[1 2 3;4 5 6] 输入矩阵B ：[1 2;3 4;5 6] 运行结果：（略） [例]采用单下标方式对n阶魔方阵的元素进行引用。当下标超出元素个数（最大单下标）时，改为引用最后一个元素，并显示出错警告信息。 程序如下： n=input(‘输入魔方矩阵的阶数:’); A=magic(n) index=input(‘输入单下标:’) try a=A(index) catch a=A(end) lasterr end 练习： 用if语句来完成本题。 3.5 循环结构 Matlab中循环语句有两种：for语句和while语句。 3.5.1 for语句 for 语句的一般格式： for 循环变量=数组 循环体语句 end 数组为向量或矩阵。执行时依次将数组的各列赋值给循环变量，然后执行循环体语句，直至各列元素处理完毕。 [例] 建立命令文件,输入以下程序，分析执行过程 A=input(A=) sumA=0; for x=A sumA=sumA+x; end sumA 请分别输入以下数据，观察程序的运行结果，分析程序的功能： 1) [1 2 3;4 5 6] 2) [1 2 3 4 5 6] 3) [1;2;3] 说明：本题可以用sum函数来改写（留作练习）。 for语句最常用的形式 格式： for 循环变量=表达式1：表达式2：表达式3 循环体语句 end 其中，表达式1为循环变量的初值，表达式2的值为步长，表达式3的值为循环变量的终值。当步长为1时，表达式2可以省略。 执行过程 [例3-18]求 ，其中，n=10. 程序如下： y=0; for k=1:10 y=y+1/2^k; end y 说明：若用sum函数，本题可改写为： K=1:10; x=1./2.^K; y=sum(x) 或者： n=1