2010-3循环语句.ppt

1、数学实验三,MATLAB循环语句,MATLAB语言编程介绍,MATLAB中各种命令可以完成许多单一的任务，对于某些较为复杂的问题，仅靠现有的命令或函数来解决，往往是难以达到目的 。为此，要运用MATLAB编程语言编制程序，形成M-文件。 程序是使计算机完成各项运算的命令集，运行一个编制好的程序，计算机会从第一条命令行开始，一行接一行地执行相应的命令，直到终止。 程序编写调试完成后，需要存盘，形成永久性文件，可以随时对它进行调用或修改。 文件名以字母开头，但不能用专用变量名，如pi,ans,eps等。,程序结构,程序一般包含：,a=input(请输入数组 an= ); m,n=size(a);

2、M=a(1); for i=2:n if Ma(i) M=a(i); end end fprintf( M=%.5fn ,M);,数据输入,数据处理,结果输出,M-文件中输入、输出命令,直接赋值输入： 提示对话输入(input命令) x=2; x=input(请输入参数 x=); a=1,2;3,4; a=input(请输入矩阵 a=); s= any string!; s=input(Please input s=); 直接输出： 格式控制输出(fprintf命令) x fprintf(x=%.0f, y=%.5fn,pi,pi); disp(a,b); fprintf(x=%5.0f, y

3、=%10.5fn,pi,pi);,一.For语句,格式： for i=n1:(step):n2 commands; end 作用：i从n1开始，执行命令集commands,遇到end，i=i+step,重复执行，直到i n2. 省略格式： for i=n1:n2这里step=1.,例1：求n个奇数和：s=1+3+5+(2n-1) n=input(please input n=); s=0; for i=1:n s=s+(2*i-1); fprintf(i=%.0f, s=%.0fn,i,s); end 程序wp31,例2：求阶乘：p=12 3 n=n! n=input(请输入 n= ); p=

4、1; for i=1:n p=p*i; fprintf( i=%.0f, p=%.0fn ,i,p); end wp32,例3：求e：e=1+1+1/2!+1/3!+1/n! n=input(请输入 n= ); p=1;e=1; for i=1:n p=p*i; p1=1/p; e=e+p1; fprintf( i=%.0f, p=%.0f, e=%.8f n ,i,p,e); end wp33,二.While语句,格式： while (condition is true) commands; end 作用：当条件成立时，执行命令集commands,直到条件不成立.,例4：求和：s=1+2+n

5、3000 s=0;k=1; while s+k3000 s=s+k; fprintf( k=%.0f, s=%.0fn ,k,s); k=k+1; end wp34,例5：求e：e=1+1+1/2!+1/3!+ (1/n!=1.0e-8 p=p*i; r=1/p; e=e+r; fprintf( i=%.0f, r=%.9f, e=%.9f n ,i,r,e); i=i+1; end wp35,三.选择控制命令(if),单项选择控制 格式： if (condition is true) commands; end 作用：若条件成立，则执行命令集 commands. 否则，不执行。,例6：求n个

6、实数中最大的数M. a=input(请输入数组 an= ); m,n=size(a); M=a(1); for i=2:n if Ma(i) M=a(i); end end fprintf( M=%.5fn ,M); wp36,四.多项选择控制命令(if),多项选择控制 格式： if (condition is true) commands; elseif (condition is true) commands; else commands; end 作用：若条件成立，则执行命令集 commands. 否则，不执行。,例7：求函数 y=x2 sgn(x)的值 x=input(please i

7、nput x=) if x0 z=1; elseif x=0 z=0; else z=-1; end y=x2*z wp37,五.定义函数,格式 function y1,y2=funname(x,y,a,n) 其中： function 为关键字； y1,y2,为输出列表； funname为自定义的函数名； (x,y,a,n,)为函数参数，用来传递相 关数值. 存盘时，要求用funname作为函数名.,例7：建立符号函数sgn(x) function sn=sgn(x) if x0 sn=1; elseif x=0 sn=0; else sn=-1; end 以sgn作为文件名存盘，即建立了函数

8、。 调用： 在命令区执行 : sn=sgn(10)或sn=sgn(-2),例8：建立计算(a+b)n,(a-b)n的函数,function p,q=fun1(a,b,n) p=(a+b).n; q=(a-b).n; fun1.m,例编制程序，从键盘输入a,b,n,计算(a+b)n,(a-b)n wp38,a=input(a=) b=input(b=) n=input(n=) p,q=fun1(a,b,n); fprintf(a+b)n=%.4f,(a-b)n=%.4fn,p,q),本次实验上机任务,教材：李继成等 上机操作实验至实验中所有指令和例题，熟悉Matlab相关命令；,数值计算问题举例

9、,问题：求无理数的近似值 先求 的近似值，再设计通用函数.,例1：用二分法求函数x2-2=0的正实根.,输入初值：,输出结果：STOP,No,No,No,Yes,Yes,Yes,a=1;b=2;K=1000;er=1.0e-8;f=inline(x.2-2);k=1;c=(a+b)/2;while kKr=c;if abs(f(r)1.0e-10fprintf(k=%.0f,r=%.8fn,k,r);break;end,if f(a)*f(r)0 b=r;else a=r;endc=(a+b)/2; fprintf(k=%.0f,r=%.8fn,k,r);if abs(r-c)er break

10、;endk=k+1;end,erfenfa.m,1. 每门课程考试阅卷完毕,任课教师都要对各班的考试成绩进行统计,统计内容包括:全班人数,总得分,平均得分,不及格的人数及90分(包括90分)以上的人数.请编制程序解决这一问题,并自给一组数据验证程序的正确性.要求:使用者在提示下通过键盘输入学生成绩,计算机自动处理后,显示需要的结果.,本次实验自选个题目,任选题;,任选题提交能够正确运行并附有结果的程序,2. 根据我国个人所得税计算方法,编制程序,要求:使用者在系统提示下通过键盘输入月工资薪金收入总数,计算机则在屏幕上显示个人所得税额,界面友好,方便使用. 个人所得税计算方法: 月个人所得税（月

11、工资薪金收入1600）*适用税率速算扣除数 附表：个人所得税税率表（工资、薪金所得适用）,个人所得税税率表（工资、薪金所得适用）,级 数 全月应纳税所得额 税率（）税率 速算扣除数1 不超过500元的 52 超过500元至2000元的部分 10 253 超过2000元至5000元的部分 15 1254 超过5000元至20000元的部分 20 3755 超过20000元至40000元的部分 25 13756 超过40000元至60000元的部分 30 33757 超过60000元至80000元的部分 35 63758 超过80000元至100000元的部分 40 103759 超过100000

12、元的部分 45 15375,、任给一个自然数n，写出1到n之间所有能够表示成连续若干个自然数的和的数及其所有表示方法,、求，中同时满足下列条件的数 （）该数各位数字之和为奇数 （）该数是素数,a=1.;b=2.;r=1.0e-7;ya=fc1(a);yb=fc1(b);k=1;if abs(ya)r root=a; fprintf(k=%.0f, Root x=%.4f, f(x)=%12.7fn,k,root,ya); return;endif abs(yb)r root=b; fprintf(k=%.0f, Root x=%.4f, f(x)=%12.7fn,k,root,yb);,function y=fc1(x) y=1./(x-.3).2+.01)+1./(