《MATLAB程序设计课件》.ppt

例3.6 某商场对顾客所购买的商品实行打折销售， 已知打折标准为： 商品价格小于200元时，折扣率为0; 商品价格大于或等于200元且小于500时，折扣率为0.03; 商品价格大于或等于500且小于1000元时，折扣率为0.05; 商品价格大于或等于1000且小于2500元时，折扣率为0.08; 商品价格大于或等于2500且小于5000元时，折扣率为0.1; 商品价格大于或等于5000元时，折扣率为0.14; 求所售商品的实际销售价格。 目录 ＞＞price=input('请输入商品价格'); switch fix(price/100) case {0,1} rate=0; case {2,3,4} rate=3/100; case num2cell(5:9) rate=5/100; case num2cell(10:24) rate=8/100; case num2cell(25:49) rate=10/100; otherwise rate=14/100; end price=price*(1-rate) 请输入商品价格890 price = 845.5000 3.3.3 try语句 语句格式为： try 语句组1 catch 语句组2 end try语句先试探性执行语句组1，如果语句组1在执行过程中出现错误，则将错误信息赋给保留的lasterr变量，并转去执行语句组2。这种试探性执行语句是其他高级语言所没有的。 目录 例3.7 矩阵乘法运算要求两矩阵的维数相容，否则会出错。先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘。 程序如下： 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 %显示出错原因 目录 C = 7 16 27 40 55 72 ans = Error using ==> * Inner matrix dimensions must agree. 3.4 循环结构 3.4.1 for语句 基本格式 for 循环变量＝起始值：步长：终止值 循环体 end 步长缺省值为1，可以在正实数或负实数范围内任意指定。对于正数，循环变量的值大于终止值时，循环结束；对于负数，循环变量的值小于终止值时，循环结束。循环结构可以嵌套使用。 目录 for i=1:5 for j=1:5 a(i,j)=1/(i+j-1); end end a= 1.0000 0.5000 0.3333 0.2500 0.2000 0.5000 0.3333 0.2500 0.2000 0.1667 0.3333 0.2500 0.2000 0.1667 0.1429 0.2500 0.2000 0.1667 0.1429 0.1250 0.2000 0.1667 0.1429 0.1250 0.1111 例3.8 已知求y的表达式， 　　　　 当n=100时，求y的值。 目录 在实际MATLAB编程中，为提高程序的执行速度，常用向量运算来代替循环操作： n=100; i=1:n; f=1./i.^2; y=sum(f) 程序如下： y=0;n=100; for i=1:n y=y+1/i/i; end y y = 2047/1252 例3.9 求定积分 。 程序如下： 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提供了有关数值积分的标准函数，实际应用中可直接调用这些函数求数值积分。 s = 754/837 for语句更一般的格式： for 循环变量=矩阵表达式 循环体语句 end 执行过程是依次