matlab软件设计上机指南

《matlab软件设计上机指南》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、MATLAB软件设计应用技术学院2008-6第一章MATLAB程序设计一、M文件v用户如想灵活应用matlab去解决实际问题，充分调用matlab的科学技术资源，就需要编辑m文件v包含matlab语言代码的文件称为m文件，其扩展名为m。v编辑m文件可使用各种文本编辑器。vm文件的类型是普通的文本文件，我们可以使用系统认可的文本文件编辑器来建立m文件。如dos下的edit，windows的记事本和word等。v具体的创建方法：1.在matlab命令窗口点击file菜单new®m-file2．.m文件m文件的语法类似于c语言，但又有其自身特点。它只是一个简单的ASCII码文本

2、文件，执行程序时逐行解释运行程序，matlab是解释性的编程语言。m文件有两类独立的m文件—称命令文件、可调用m文件—称函数文件(1).命令文件—简单的m文件命令文件实际上是一串指令的集合，与在命令窗口逐行执行文件中的所有指令，其结果是一样的。没有输入输出参数v命令文件包括两部分：注释文件和程序文件%RANKNumberoflinearlyindependentrowsorcolumns.%K=RANK(X)isthenumberofsingularvaluesofX%thatarelargerthanMAX(SIZE(X))*NORM(X)*EPS.%K=RANK(X,

3、tol)isthenumberofsingularvaluesofXthat%arelargerthantol.%Copyright(c)1984-94byTheMathWorks,Inc.s=svd(x);if(nargin==1)tol=max(size(x))*max(s)*eps;endr=sum(s>tol);(2).matlab内置函数文件matlab自定义的函数文件称内置函数文件调用内置函数的方法：使用函数名并给出相应的入口、出口参数即可。例如：sin.m函数——用typesin查不到。调用格式：y=sin(2*x)实际应用中：x=0:2*pi/180:2*

4、pi;y=sin(2*x)plot(x,y)函数m文件的格式：function返回变量=函数名（输入变量）注释说明语句段程序语句段特定规则：1.函数m文件第一行必须以单词function作为引导词，必须遵循如下形式：function<因变量>=<函数名>(<自变量>)2.m文件的文件名必须是<函数名>.m。3.程序中的变量均为局部变量，不保存在工作空间中。其变量只在函数运行期间有效二、微分方程一般微分方程式描述系统内部变量的变化率如何受系统内部变量和外部激励，如输入，的影响。当常微分方程式能够解析求解时，可用MATLAB的符号工具箱中的功能找到精确解。在本书的后面将介绍该

5、工具箱的一些特点。在微分方程难以获得解析解的情况下，可以方便地在数值上求解。为了说明起见，考虑描述振荡器的经典的范得波（VarderPol）微分方程。与所有的数值求解微分方程组的方法一样，高阶微分方程式必须等价地变换成一阶微分方程组。对于上述微分方程，通过重新定义两个新的变量，来实现这种变换。令y1=x且y2=dy/dx则dy1/dt=y2根据这个微分方程组，可用MATLAB的函数ode23和ode45求出系统随时间变化的运动情况。调用这些函数时，需要编写一个函数M文件，给定当前时间及y1和y2的当前值，该函数返回上述导数值。MATLAB中，这些导数由一个列向量给出。在本

6、例中，这个列向量为yprime。同样，y1和y2合并写成列向量y。所得函数M文件是：functionyprime=vdpol(t,y);%VDPOL(t,y)returnsderivativesoftheVanderPolequation:%%x‘‘-mu*(1-x^2)*x‘+x=0(‘=d/dx,‘‘=d^2/dx^2)%%lety(1)=xandy(2)=x‘%%theny(1)‘=y(2)%y(2)‘=MU*(1-y(1)^2)*y(2)-y(1)globalMU%choose0

7、(1)^2)*y(2)-y(1)];%outputmustbeacolumn给定这个完整地描述微分方程的函数，计算结果如下：>>globalMU%defineMUasaglobalvariableintheCommandWorkspace>>MU=2;%setglobalparametertodesiredvalue>>[t,y]=ode23(‘vdpol‘,0,30,[1;0]);%to=0,tf=30,yo=[1;0]>>y1=y(:,1);%firstcolumnisy(1)versustimepointsint>>y2