MATLAB及其应用-第4讲.pptVIP

第4讲 M文件和面向对象编程 4.1 M文件 4.2 MATLAB控制流 4.3 脚本文件和函数文件 4.1 M文件 一、M文件概述 用MATLAB语言编写的程序，称为M文件。 二、M文件的建立和编辑 1、建立新的M文件 (1)菜单操作：File → New→ M-file。 (2)命令操作：在Command Window中输入edit。 (3)命令按钮操作：单击工具栏上的“新建” 按钮。 2、编辑已有的M文件 例1：分别建立命令文件和函数文件，将华氏温度F转换为摄氏温度C。 命令文件： F2C1.m clear F=input(Input Fahrenheit temperature：F=); C=5*(F-32)/9 函数文件：F2C2.m function C=F2C2(F) C=5*(F-32)/9 4.2 MATLAB控制流 一、程序结构 二、控制程序流的其他常用指令 一、 程序结构 1. 顺序结构 例： %用符号计算求微分 syms a t x; f=[a, t^3; t*cos(x), log(x)]; df=diff(f)； dfdt2=diff(f, t, 2)； dfdxdt=diff(diff(f, x), t)； 注意： （1）循环结构可嵌套使用； （2）为得到高效代码，应尽量提高代码的向量化程度，避免使用循环结构； （3）循环指令前最好对数组进行预定义，以提高运行效率； （4）使用while循环时，一定要注意避免进入死循环。 例3：若数组a中a1=a2=1,且ak+2=ak+ak+1，试求该数组中第一个大于10000的元素 程序如下： clear a(1)=1; a(2)=1; i=2; while a(i)=10000 a(i+1)=a(i-1)+a(i); i=i+1; end i, a(i) 例4. 计算分段函数值。 程序如下： clear x=input(请输入x的值:); if x=0 y=cos(x+1)+sqrt(x*x+1); else y=x*sqrt(x+sqrt(x)); end y 例5 输入一个字符，若为大写字母，则输出其后继字符，若为小写字母，则输出其前导字符，若为数字字符则输出其对应的数值，若为其他字符则原样输出。 程序如下： clear c=input(请输入一个字符,s); if c=A c=Z disp(char(abs(c)+1)); elseif c=a c=z disp(char(abs(c)-1)); elseif c=0 c=9 disp(abs(c)-abs(0)); else disp(c); end 例6 某商场对顾客所购买的商品实行打折销售，已知打折标准，求所售商品的实际销售价格。 程序如下： 例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 %显示出错原因 二、 控制程序流的其他常用指令 1. return 指令 强制MATLAB结束执行当前函数并把控制转出 2. keyboard指令： 当程序遇到keyboard时，MATLAB将控制权交给键盘，用户可以从键盘输入各种合法的MATLAB指令。只有当用户使用回车结束输入后，控制权才交还给程序。 3. input指令和disp指令 input指令： 把MATLAB的控制权暂时交给用户。此后，用户通过键盘键入数值、字符串或表达式，并经“回车”把键入内容输入工作空间，同时把“控制权”交还给MATLAB。 常用的格式： v=input(‘message’) %将用户键入的内容赋给变量v v=input(‘message’,’s’) %将用户键入的内容作为字符串赋给变量v disp指令： 格式：disp(输出项) 输出项可以为字符串，也可为矩阵 例8：求一元二次方程ax2+bx+c=0的根。 程序如下： clear a=input(‘a=?’); b=input(‘b=?’); c=input(‘c=?’); d=