MATLAB软件设计上机指南.doc

MATLAB软件设计 应用技术学院 2008-6 第一章 MATLAB程序设计 一、M文件 用户如想灵活应用matlab去解决实际问题，充分调用matlab的科学技术资源，就需要编辑m文件 包含matlab语言代码的文件称为m文件，其扩展名为m。 编辑m文件可使用各种文本编辑器。 m文件的类型是普通的文本文件，我们可以使用系统认可的文本文件编辑器来建立m文件。如dos下的edit，windows的记事本和word等。 具体的创建方法： 在matlab命令窗口点击file 菜单 new ( m-file 2．. m文件 m文件的语法类似于c语言，但又有其自身特点。它只是一个简单的ASCII码文本文件，执行程序时逐行解释运行程序，matlab是解释性的编程语言。 m文件有两类(独立的m文件 — 称命令文件、可调用m文件 — 称函数文件 (1). 命令文件 — 简单的m文件 命令文件实际上是一串指令的集合，与在命令窗口逐行执行文件中的所有指令，其结果是一样的。没有输入输出参数 命令文件包括两部分：注释文件和程序文件 %RANK Number of linearly independent rows or columns. % K = RANK(X) is the number of singular values of X % that are larger than MAX(SIZE(X)) * NORM(X) * EPS. % K = RANK(X,tol) is the number of singular values of X that % are larger than tol. % Copyright (c) 1984-94 by The MathWorks, Inc. s = svd(x); if (nargin == 1) tol = max(size(x)) * max(s) * eps; end r = sum(s tol); (2).matlab内置函数文件 matlab自定义的函数文件称内置函数文件 调用内置函数的方法：使用函数名并给出相应的入口、出口参数即可。 例如：sin.m函数——用type sin查不到。 调用格式：y=sin(2*x) 实际应用中： x=0:2*pi/180:2*pi; y=sin(2*x) plot(x,y) 函数m文件的格式： function 返回变量=函数名（输入变量） 注释说明语句段 程序语句段 特定规则： 函数m文件第一行必须以单词function作为引导词，必须遵循如下形式： function 因变量=函数名(自变量) m文件的文件名必须是函数名 .m。 程序中的变量均为局部变量，不保存在工作空间中。 其变量只在函数运行期间有效 二、微分方程 一般微分方程式描述系统内部变量的变化率如何受系统内部变量和外部激励，如输入，的影响。当常微分方程式能够解析求解时，可用MATLAB的符号工具箱中的功能找到精确解。在本书的后面将介绍该工具箱的一些特点。 在微分方程难以获得解析解的情况下，可以方便地在数值上求解。为了说明起见，考虑描述振荡器的经典的范得波（Var der Pol）微分方程。 与所有的数值求解微分方程组的方法一样，高阶微分方程式必须等价地变换成一阶微分方程组。对于上述微分方程，通过重新定义两个新的变量，来实现这种变换。 令y1=x 且 y2=dy/dx 则 dy1/dt=y2 根据这个微分方程组，可用MATLAB的函数ode23和ode45求出系统随时间变化的运动情况。调用这些函数时，需要编写一个函数M文件，给定当前时间及y1和y2的当前值，该函数返回上述导数值。MATLAB中，这些导数由一个列向量给出。在本例中，这个列向量为yprime。同样，y1和y2合并写成列向量y。所得函数M文件是： function yprime=vdpol(t , y); % VDPOL(t , y) returns derivatives of the Van der Pol equation: % % x ‘‘-mu *(1-x ^2)*x ‘+x=0 (‘ = d/dx , ‘‘ = d^2/dx^2) % % let y(1)=x and y(2)=x‘ % % then y(1) ‘ = y(2) % y(2) ‘ = MU*(1-y(1)^2)*y(2)-y(1) gl