第1章 MTLAB程序设计.pptVIP

第1章 MATLAB程序设计 1.1 M文件 1.2 程序控制结构 1.3 函数文件 Matlab 编程简介 M 文件的建立与打开 命令文件举例 首先建立函数文件fexch.m： function [a,b]=exch(a,b) c=a;a=b;b=c; 然后在MATLAB的命令窗口调用该函数文件： clear;x=1:10;y=[11,12,13,14;15,16,17,18]; [x,y]=fexch(x,y) 一些测试函数 运算优先级 1.2 程序控制结构 顺序结构 例1.2 求一元二次方程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))]); if 条件语句 if 条件语句 例1.4 输入一个字符，若为大写字母，则输出其对应的小写字母；若为小写字母，则输出其对应的大写字母；若为数字字符则输出其对应的数值，若为其他字符则原样输出。 c=input(请输入一个字符,s); if c=A c=Z disp(setstr(abs(c)+abs(a)-abs(A))); elseif c=a c=z disp(setstr(abs(c)- abs(a)+abs(A))); elseif c=0 c=9 disp(abs(c)-abs(0)); else disp(c); end switch 语句 例1.5 某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示)： price200 没有折扣 200≤price500 3%折扣 500≤price1000 5%折扣 1000≤price2500 8%折扣 2500≤price5000 10%折扣 5000≤price 14%折扣 输入所售商品的价格，求其实际销售价格。 程序如下： price=input(请输入商品价格); switch fix(price/100) case {0,1} %价格小于200 rate=0; case {2,3,4} %价格大于等于200但小于500 rate=3/100; case num2cell(5:9) %价格大于等于500但小于1000 rate=5/100; case num2cell(10:24) %价格大于等于1000但小于2500 rate=8/100; case num2cell(25:49) %价格大于等于2500但小于5000 rate=10/100; otherwise %价格大于等于5000 rate=14/100; end price=price*(1-rate) %输出商品实际销售价格 语句格式为： try 语句组1 catch 语句组2 end try语句先试探性执行语句组1，如果语句组1在执行过程中出现错误，则将错误信息赋给保留的lasterr变量，并转去执行语句组2。 例1.6 矩阵乘法运算要求两矩阵的维数相容，否则会出错。先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘。 程序如下： 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 %显示出错原因 循环结构 例1.7 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。 程序如下： for m=100:999