MATLAB程序设计12261.pptVIP

第三章 MATLAB程序设计 宇航科学技术学院 隋丽 010bit_suili@bit.edu.cn MATLAB的工作模式 指令驱动模式 在MATLAB命令行窗口下用户输入单行指令时， MATLAB立即处理这条指令，并显示结果，这就是MATLAB命令行方式。 命令行方式程序可读性差，而且不能存储，当处理复杂问题和大量数据时很不方便。 M文件模式 将要执行的命令全部写在一个文本文件中，这样能使程序简洁明了，便于对程序的修改和维护。 3.1 M文件 3.1 M文件 用MATLAB语言编写的程序，称为M文件。M文件可以根据调用方式的不同分为两类：脚本文件(Script File)和函数文件(Function File)。 脚本文件和函数文件的扩展名均是.m。 例3-1 分别建立脚本文件和函数文件，将华氏温度f转换为摄氏温度c。 1）首先建立脚本文件并以文件名f2c.m存盘。 clear; %清除工作空间中的变量 f=input(Input Fahrenheit temperature：); c=5*(f-32)/9 2）然后在MATLAB的命令窗口中输入f2c，将会执行该脚本文件，执行情况为： f2c Input Fahrenheit temperature：73 c = 22.7778 首先建立函数文件f2c2.m。 function c=f2c2(f) c=5*(f-32)/9 然后在MATLAB的命令窗口调用该函数文件。 clear; y=input(Input Fahrenheit temperature：); x=f2c2(y) 输出情况为： Input Fahrenheit temperature：73 c = 22.7778 x = 22.7778 3.1.1 脚本M文件 将所有要执行的命令按顺序放到一个扩展名为.m的文本文件中，即可生成一个M文件。 脚本文件就是由一系列的MATLAB指令和命令组成的纯文本格式的M文件。 脚本文件没有输入参数，也没有输出参数。 执行脚本文件时，文件中的指令或者命令按照出现在脚本文件中的顺序依次执行。 脚本文件示例 % 注释行 % M文件示例 % “flower petal” % 以下为代码行 % 计算 theta=-pi:0.01:pi; %设置角度向量，分度为0.01 rho(1,:)=2*sin(5*theta).^2; ％计算频率为5的正弦函数平方 rho(2,:)=cos(10*theta).^3; rho(3,:)=sin(theta).^2; rho(4,:)=5*cos(3.5*theta).^3; for k=1:4 ％设置循环次数 subplot(2,2,k) % 图形输出 polar(theta,rho(k,:)) ％绘制极坐标图 end disp(程序运行结束！) 在MATLAB命令行中运行该脚本文件： script_examp 程序运行结束！ MATLAB会出现相应的图形窗口 3.1.2 函数M文件 函数文件能够接受用户的输入参数，进行计算，并将计算结果作为函数的返回值返回给调用者。 在MATLAB中具有不同类型的函数文件 内装函数文件 系统M函数文件 系统MEX函数文件 用户自定义MEX函数文件 用户自定义的M函数文件 3.1.2 函数M文件 当调用MATLAB的函数时，在命令窗口输入函数名，则该函数的程序即被调用并执行。函数M文件由5部分组成： 函数定义行，由关键字function开头 帮助文本的标题，简称H1行 帮助文本的内容 函数体 附注 例：…… 函数M-文件的结构 function 输出形参表=函数名称（输入形参表） 特定规则： 1） 函数M文件第一行必须以单词function作为引导词； 2） 程序中的变量均为局部变量，不保存在工作空间中，其变量只在函数运行期间有效。 3） 函数名的命名规则与变量的命名规则相同。 4） 函数M文件的名称一般要与函数名一致，否则调用时会出错。 函数调用的一般格式为：在命令窗口中输入 例：下面的程序是一段求n！的函数文件 function f=jch(n) f=1; for i=2:n f=f*i; end 保存为jch.m。 jch(10) ans = 3628800 对于多输入、多输出函数，它的定义行为：function [x,y,z,…]=函数名（a,b,c,…) 设可逆方阵A，编写同时求