第8章M文件程序设计基础.pptVIP

例：求一元二次方程ax2 +bx+c=0的根。 程序如下： a=input(a=?); b=input(b=?); c=input(c=?); d=b*b-4*a*c; x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)]; disp([x1=,num2str(x(1)),,x2=,num2str(x(2))]); 3．程序的暂停 暂停程序的执行可以使用pause函数，其调用 格式为： pause(延迟秒数) 如果省略延迟时间，直接使用pause，则将暂 停程序，直到用户按任一键后程序继续执行。 若要强行中止程序的运行可使用Ctrl+C命令。 例：一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。 程序如下： for m=100:999 m1=fix(m/100); %求m的百位数字 m2=rem(fix(m/10),10); %求m的十位数字 m3=rem(m,10); %求m的个位数字 if m==m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) end end 例:已知y，当n=100时，求y的值。 程序如下： y=0;n=100; for i=1:n y=y+1/ i^2; end y for语句更一般的格式为： for 循环变量=矩阵表达式 循环体语句 end 执行过程是依次将矩阵的每行中各列元素赋给循 环变量，然后执行循环体语句，直至所有行各列 元素处理完毕。 例:写出下列程序的执行结果。 s=0; a=[12,13,14;15,16,17;18,19,20;21,22,23]; for k=a s=s+k; end disp(s); 2．while语句 while语句的一般格式为： while (条件) 循环体语句 end 其执行过程为：若条件成立，则执行循环体语 句，执行后再判断条件是否成立，如果不成立则 跳出循环。 例:从键盘输入若干个数，当输入0时结束输入，求 这些数的平均值和它们之和。 程序如下： sum=0; n=0; val=input(Enter a number (end in 0):); while (val~=0) sum=sum+val; n=n+1; val=input(Enter a number (end in 0):); end if (n 0) sum mean=sum/n end 与循环结构相关的语句还有break语句和continue 语句。它们一般与if语句配合使用。 break语句用于终止循环的执行。当在循环体内 执行到该语句时，程序将跳出循环，继续执行循 环语句的下一语句。 continue语句控制跳过循环体中的某些语句。当 在循环体内执行到该语句时，程序将跳过循环体 中所有剩下的语句，继续下一次循环。 例: 求[100，200]之间第一个能被21整除的整数。 程序如下： for n=100:200 if rem(n,21)~=0 continue end break end n 4．循环的嵌套 如果一个循环结构的循环体又包括一个循环 结构，就称为循环的嵌套，或称为多重循环结 构。 if语句 在MATLAB中，if语句有3种格式。 (1) 单分支if语句： if 条件 语句组 end 当条件成立时，则执行语句组，执行完之后继续 执行if语句的后继语句，若条件不成立，则直接 执行if语句的后继语句。 (2) 双分支if语句： if 条件 语句组1 else 语句组2 end 当条件成立时，执行语句组1，否则执行语句组2， 语句组1或语句组2执行后，再执行if语句的后继 语句。 例：计算分段函数。 程序如下： x=input(请输入x的值:); if x==10 y=cos(x+1)+sqrt(x*x+1); else y=x*sqrt(x+sqrt(x)); end y (3) 多分支if语句： if 条件1 语句组1 elseif 条件2 语句组2 …… elseif 条件m 语句组m else