MATLAB计算与仿真I.pptVIP

例3：求阶乘：p=1×2 × 3 × … × n=n! n=input(请输入 n= ); p=1; for i=1:n p=p*i; fprintf( i=%.0f, p=%.0f\n ,i,p); end 例4：求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 例5：求pi：pi/4=1-1/3+1/5-1/7+1/9-1/11+…… n=input(请输入 n= ); pi=0;k=-1; for i=1:2:n k=(-1)*k; pi=pi+k*1/i; fprintf( i=%.0f, pi=%.8f \n ,i,4*pi); end M-文件中条件循环命令(while命令) 格式： while (condition is true) commands; end 作用：重复执行命令集commands. 例6：求和：s=1+2+…+n3000 s=0;k=1; while s3000 s=s+k; fprintf( k=%.0f, s=%.0f\n ,k,s); k=k+1; end 例7：求e：e=1+1+1/2!+1/3!+… (1/n!1.0e-8) p=1;e=1;r=1;i=1; while r=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 M-文件中选择控制命令(if命令) 单项选择控制 格式： if (condition is true) commands; end 作用：若条件成立，则执行命令集 commands. 否则，不执行。 例8：求n个实数中最大的数M. a=input(请输入数组 a[n]= ); [m,n]=size(a); M=a(1); for i=2:n if Ma(i) M=a(i); end fprintf( M=%.5f, a(%.0f)=%.5f \n ,M,i,a(i)); end M-文件中选择控制命令(if命令) 多项选择控制 格式： if (condition is true) commands; elseif (condition is true) commands; else commands; end 作用：若条件成立，则执行命令集 commands. 否则，不执行。 例9：建立符号函数sign(x) x=input(x=); if x0 sn=1; elseif x==0 sn=0; else sn=-1; end fprintf(x=%.5f, sn=%.0f\n,x,sn); Matlab编程 (续) 程序设计技术 循环向量化计算模式 目的: 向量计算替代循环模式 实例: i初始化