Matlab 第10讲函数文件
10.4 一些特殊函数的使用
例5 匿名函数使用示例。 >> A=[2,3,4]; >> B=[5,6,7]; >> sum1=@(x,y)(A*x+B*y) sum1 =
@(x,y)(A*x+B*y) >> sum1(5,7) ans =
45 57 69
10.4 一些特殊函数的使用
多项式的表示
函数M文件：基本内容包含以下几个方面。 1. 函数声明行 只有函数M文件才有函数声明行，函数声明行位于M文件的第一行，通过
function关键字指定该文件是什么函数文件，并指定函数名和输入、输 出参数。函数声明行的语法表达式如下：
function [out1,out2,……] = myfun (in1,in2,……) 参数说明
如果在若干个M文件中，都把某个变量定义为全局变量，则这些函数将公共 使用这一变量。所有函数都可以对它进行存取和修改操作。定义全局变 量是M文件间传递信息的一种手段，全局变量给函数间的数据带来了方 便，但却破坏了函数对变量的封装，降低了程序的可读性，因而在结构 化程序设计中，全局变量是不受欢迎的。特别是当程序较大，子程序较 多时，全局变量将个程序调试和维护带来不便，故不提倡使用全局变量 。
10.3 函数句柄
主函数和子函数
例4 %ex_subfun.m% function[avg,med] =ex_subfun(x) %主函数 n =length(x); avg=mean(x,n); med=median(x,n); function a=mean(x,n) %子函数，计算平均值 a =sum(x)/n; function m=median(x,n) %子函数，计算中值 x =sort(x); if rem(n,2) ==1 m=x((n+1)/2); else m=(x(n/2)+x(n/2+1))/2; end
10.4 一些特殊函数的使用
多项式的表示
函数 roots polyval polyvalm polyder poly
Polyfit Residue
功能 多项式求根 多项式求值 矩阵多项式求值 多项式求导 求矩阵的特征多项式；或者求一个多项 式，其根为指定的数值 多项式曲线拟合 求解余项
10.4 一些特殊函数的使用
当函数有一个以上输出变量时，输出变量包含在括号内。例如 ，[V,D] = eig(A)。不要把这个句法与等号右边的[V,D] 相混淆。右边的[V,D] 是由数组V和D 所组成。
10.2 变量及函数以及函数的调用
变量及函数
MATLAB中的变量分为局部变量和全局变量。函数文件中的变量都是局部 的，即一个函数文件中定义的变量不能被另一个函数文件或其它M文件 使用。当函数调用完毕后，该函数文件中定义的所有局部变量都将被释 放，即全部被清除。
10.1 m文件和函数文件
2. H1行 H1表示帮助文件的第一行，紧跟函数定义行。H1行的内容应该是有关该文
件的重要信息，因此H1行应尽量言简意赅。当使用lookfor查找与某个 单词相关的函数时，lookfor便在H1行进行搜索。 3. 帮助文档 帮助文档是为了对M文件更加详细的说明，它的内容应包括解释m文件实现 的功能、M文件中各参数的含义以及创作版权等信息。帮助文档紧跟着 H1行，以“%”为开头。当在命令窗口中输入“help filename”时，帮助 文档会在H1行后显示。 4. 注释 注释是为了便于程序阅读和修改对程序内容的解释说明。注释以“%”为开头 ，可以在程序的任何位置。
x.^2. arglist:为参数列表，它指定函数的输入参数列表，对于多个输入参数的情
况，通常要用逗号分隔各个参数。 符号@：是MATLAB中创建函数句柄的操作符，表示创建由输入参数列表
arglist和expr确定的函数句柄，并把这个函数句柄返回给变量 fhandle,这样以后就可以通过fhandle来调用定义好的这个函数。
Function y=Fab(n) if (n==1) y=0; elseif (n==2) y=1; else y=Fab(n-1)+2*Fab(n-2);
end
10.3 函数句柄
主函数和子函数
主函数一般是指与保存文件同名的函数/在指令窗口或其它函数 中可直接调用的函数。MATLAB中的主函数是针对其内部 嵌套函数和子函数而言的。每一个函数M文件的第一行定义 的函数就是主函数，除了主函数外，还可以编写多个嵌套 函数或子函数，以便在主函数功能实现中进行调用。一个M 文件只能包含一个主函数，并且一般将M文件名和主函数名 设为一致。
使用MATLAB函数时，例如inv, abs, angle和sqrt，MATLAB获取传 递给它的变量，利用所给的输入，计算所要求的结果。然后，把这些结果 返回。由函数执行的命令，以及由这些命令所创建的中间变量，都是隐含 的。所有可见的东西是输入和输出，也就是说函数是一个黑箱。
MATLAB提供了一个创建用户函数的结构，并以M文件的文本形式存储在 计算机上。
Matlab 程序设计第10讲 函数文件
主讲人：王佐才 博士 研究员/博导
合肥工业大学土木与水利工程学院
提纲
m文件和函数文件 变量及函数以及函数调用 函数句柄 一些特殊函数的使用
10.1 m文件和函数文件
MATLAB提供了各种功能强大的函数库，用户可以通过这些 函数库可以实现各种复杂的运算。然而对于构建更为复杂的系 统而言，用户需要精心设计程序，构建出完成某一功能的高级 程序。通过本章的学习，能够更好的认识MATLAB编写程序 ，运用自如。
fliplr.m文件中。
函数可以有零个或更多个输入参量。函数可以有零个或更多个 输出参量。
函数可以按少于函数M文件中所规定的输入和输出变量进行调 用，但不能用多于函数M文件中所规定的输入和输出变量数 目。如果输入和输出变量数目多于函数M文件中function语 句一开始所规定的数目，则调用时自动返回一个错误。
q=quad(fun,a,b),指定允许的误差，该命令运行更快，但结果精确度低。 MATLAB中二重积分由dblquad函数来实现。二重积分调用格式:
q=dblquad(fun,xmin,xmax,ymin,ymax),函数的参数uad函数积分应用。 >> F=@(x,y)y*sin(x)+x*cos(y); >> Q=dblquad(F,pi,2*pi,0,pi)
10.3 函数句柄
函数句柄：递归调用 例2 利用函数的递归调用计算n！
% 函数文件myfactor.m function y=myfactor(n) if (n<=1) y=1; else y=n*myfactor(n-1): end
10.3 函数句柄
例3 编写Fab函数，求Fab(n)=Fab(n-1)+2Fab(n-2) (n>2)，Fab(1)=0, Fab(2)=1。
10.2 变量及函数以及函数的调用
MATLAB对变量名称有3个规定：
变量名称可以包含英文字母、数字、下划线，但变量的开头必须是一个英文 字母。
MATLAB对变量名称区分大小写。
变量名称不能超过31字符。 全局变量之间要以空格分隔，而不能以逗号分隔。通过全局变量可以实现
MATLAB工作区变量与多个函数空间的共享，任何一处对全局变量的修 改都会直接改变此全局变量的数值。 。
在MATLAB中，多项式用一个行向量表示，向量中的元素为该 多项式的系数，按照降序排列。如多项式9+7+4x+3可以 表示为向量p=[9 7 4 3]。用户可用创建 向量的方式创建 多项式，再将其显示为多项式，如：
>> p=[9 ,7,4,3]; >> y=poly2sym(p) y= 9*x^3 + 7*x^2 + 4*x + 3
0.6435 y1 =
0 >> [x2,y2]=fzero(h_f,[pi/2,pi]) x2 =
2.4981 y2 =
0 >> [x3,y3]=fzero(h_f,0) x3 =
0.6435 y3 =
0
10.4 一些特殊函数的使用
数值积分函数命令
MATLAB中一元函数对的积分用quad调用格式: q=quad(fun,a,b),采用递归自适应方法计算fun函数在区间[a,b]上的积分。
Q= -9.8696
10.4 一些特殊函数的使用
自编函数文件 例8 编写函数文件FEM求解如下图所示的弹簧系统在静力作用
下的各质量块的位移（忽略摩擦力的影响）。其输入变量 为：k=[k1,k2,k3,…,kn+1], p=[p1,p2,p3,…,pn], 输 出变量为每一个质量块的位移量x=[x1,x2,x3,…,xn]
5. M文件正文 这是M文件代码的主要部分，通常包括程序流程控制、变量赋值、函数调用
、数学运算等。
10.1 m文件和函数文件
例 1：分别创建命令文件和函数文件，将华氏温度f 转化成摄氏温度c。c=5/9(f-32)
10.1 m文件和函数文件
几点注意 函数名和文件名必须相同。例如，函数fliplr存储在名为
10.4 一些特殊函数的使用
匿名函数通常是很简单的函数，不像一般的M函数需要通过M文件编写。匿 名函数是面向命令行代码的函数形式，它通常只由一句很简单的声明语 句组成。与一般M文件函数一样的是，匿名函数也可以接受多个输入和 输出参数。
匿名函数创建的语法形式如下：
functionhandle=@(arglist) expression 参数说明： expression:通常是一个简单的MATLAB变量表达式，实现函数的功能，如
function ： 关键字，表示定义一个函数。 out1,out2,…… ： 指定输出参数列表，多个参数之间用逗号隔开。 in1,in2,…… : 输出参数列表。 myfun ： 指定该函数的函数名。函数名必须以字母开头，可以包含 任意的字母和下划线。 M文件保存时一般采用函数名作为文件名，这样在使用时有利于理解和 调用。
如何规范地编写一个MATLAB函数文件
如何规范地编写⼀个MATLAB函数⽂件在matlab中，M⽂件分为脚本⽂件和函数⽂件。
如果M⽂件的第⼀个可执⾏语句以function开头，那这个M⽂件就是函数⽂件。
函数⽂件内定义的变量为局部变量，只在函数⽂件内部起作⽤，当函数⽂件执⾏完后，这些内部变量将被清除。
本⽂介绍如何规范地编写⼀个函数⽂件。
通常，函数⽂件由函数声明⾏、H1⾏、在线帮助⽂本区、编写和修改记录、函数主体等⼏个部分组成。
格式如下：function 输出形参表 = 函数名(输⼊形参表)在线帮助⽂本区，其中第⼀⾏为H1⾏编写和修改记录函数主体函数声明⾏1. 由关键字function引导，指明这是⼀个函数⽂件，并定义函数名、输⼊参数和输出参数。
2. 函数名应当与⽂件名⼀致（保存函数⽂件时会默认以函数名作为⽂件名予以保存），如果两者不⼀致，MATLAB将以⽂件名为准，以后调⽤时使⽤⽂件名调⽤该函数。
3. 在脚本⽂件中编写函数时，需以end结尾；独⽴的函数⽂件不需以end结尾。
4. MATLAB中的函数⽂件名必须以字母开头，可以是字母、下划线、数字的任意组合，但不可以超过31个字符。
H1⾏1. 紧随函数声明⾏之后的以“%”开头的第⼀注释⾏。
2. H1⾏包括⼤写的函数名和函数功能简要描述，采⽤lookfor命令可在命令⾏窗⼝显⽰H1⾏的信息。
3. 建议在编写H1注释⾏时，尽量采⽤英⽂表述，这是为了之后的使⽤过程中关键词检索的⽅便。
在线帮助⽂本区1. 包括H1⾏以及H1⾏之后的连续的以“%”开头的注释⾏。
2. 通常包括函数输⼊变量和输出变量的含义以及调⽤说明。
3. 采⽤help命令可在命令⾏窗⼝显⽰在线帮助⽂本区的信息。
编写和修改记录1. 与在线帮助⽂本区以⼀个空⾏相隔。
2. 该⾏以“%”开头，记录了编写及修改函数⽂件的作者、⽇期、版本等内容，以⽅便后来的使⽤查询或修改。
函数主体1. 应与编写和修改记录以⼀个空⾏相隔。
2. 这部分内容包括了所有实现该函数⽂件功能的MATLAB指令。
MATLAB教程第9、10讲
函数文件：tran.m: function [rho,theta] = tran(x,y) rho = sqrt(x*x+y*y); theta = atan(y/x);
2014-7-29
y = input(‘please input y=:’); [rho,the] = tran(x,y); rho the
2014-7-29 Application of Matlab Language 5
5.3.2 函数调用
函数调用的一般格式是：
[输出实参表] = 函数名(输入实参表） 注意：函数调用时，各实参出现的顺序、个数，应与函数定 义时相同。 例5.11 利用函数文件，实现直角坐标(x,y)与极坐标(ρ,θ)之 间的转换。 调用tran.m的命令文件main1.m:
9
5.3.4 全局变量与局部变量
Matlab中，函数文件中的变量是局部变量。 如在若干函数中，都把某一变量定义为全局变量，那么这些函数将 共用这个变量。 全局变量的作用域是整个Matlab的工作空间，所有函数都可以对它 进行存取和修改。 全局变量用global命令定义，格式为： global 变量名 例5.13 全局变量应用示例。 先建立函数文件wadd.m，该函数将输入的参数加权相加： function f = wadd(x,y) BETA = 2; global ALPHA BETA s = wadd(1,2) f = ALPHA*x + BETA*y; 输出为： 在命令窗口中输入： s= global ALPHA BETA 5 ALPHA = 1; 2014-7-29 Application of Matlab Language
matlab培训讲义
MATLAB 解方程讲师：施六五邮箱：shiliuwu@主要内容：●界面设置●编辑和调试●解方程●拟合●积分●优化一、界面的设置⏹Perference设置字体，属性等⏹Set path 添加路径,增加具体的子函数⏹Desktop→Desktop Layout→Default ⏹Command window(tab type)⏹Command history⏹帮助窗口（doc docsearch look for ）二、编辑与调试⏹M文件：函数文件和脚本文件⏹文件的命名：不能与matlab内部函数同名⏹多行注释：1.Ctrl+R,Ctrl+T 2.if o end,3.%{ %}⏹自动缩进：Ctrl+I⏹双击变量查看该程序的所有的该变量⏹做标记Ctrl+F2更多右键信息⏹设置断点，F5,F10,F11,Go until Cursor ⏹分块执行%%⏹程序中断：命令窗口输入Ctrl+C⏹cell模式%% 按住Ctrl+Enter解方程●线性方程●非线性方程●隐函数方程●带参数的方程求解●微分方程●积分方程●优化线性方程组Ax=b x=A\b矩阵的分解求解线性方程组矩阵分解是指根据一定的原理用某种算法将一个矩阵分解成若干个矩阵的乘积。
常见的矩阵分解有LU分解、QR 分解、Cholesky分解，以及Schur分解、Hessenberg分解、奇异分解等。
（具体可以查看帮助文件）非线性方程求解函数：fzero 求解单变量非线性方程的根基本用法：z=fzero('fname',x0,tol,trace)例子：初值的选取Data cursor 获取坐标ar-1-0.8-0.6-0.4-0.200.20.40.60.81texp(-3 t) sin(4 t+2)+4 exp(-0.5 t) cos(2 t)-0.5函数：fsolve 基本用法：X=fsolve(@fun,x0,options)clc;clear;close allf{1}=inline('x-0.6*sin(x)-0.3*cos(y)'); f{2}=inline('y-0.6*cos(x)+0.3*sin(y)'); hold onh=arrayfun(@(i)ezplot(f{i},[-1,1,-1,1]),[1,2]); set(h(1),'color','r')-1-0.8-0.6-0.4-0.200.20.40.60.81-1-0.8-0.6-0.4-0.200.20.40.60.81xyy-0.6 cos(x)+0.3 sin(y) = 0function y=myfun(x)y(1)=x(1)-0.6*sin(x(1))-0.3*cos(x(2)); y(2)=x(2)-0.6*cos(x(1))+0.3*sin(x(2)); end隐函数的求解给定一个隐函数的自变量，通过数值方法求解得到因变量。
MATLAB7.0使用详解-第10章 SIMULINK高级仿真技术
10.3.1
使能子系统（Enable Subsystem）
• 使能子系统（Enable Subsystem）将控制信号分 为允许（enable）和禁止（disenable）两种。当 控制信号为正时，即控制信号为允许状态，系统 可以执行使能子系统中的模块；否则就禁止模块 功能。
10.3.2 触发子系统（Triggered Subsystem）
10.2 子系统的封装
• 使用子系统技术可以很好的优化系统模型得界面，是系统 模型的可读性更强。在对系统进行仿真时，首先要对系统 模块参数进行设置，对子系统也是如此，需要对子系统所 有模块进行合适的参数设置。 • 在第8章的例子中，子系统中模块的参数是逐一设置的，这 会给用户带来极大的不便。子系统的封装技术则是对子系 统的包装，使其成为一个真正意义上的SIMULINK模块，通 过封装，可以简化模型，用户使用一个动态设置的对话框 代替多个静态对话框，省去了对子系统内部结构的了解， 为用户的直接调用提供了方便。 • 一个已经创建的子系统的封装方法主要是通过使用下图所 示的封装编辑器来实现的，在封装编辑器中用户可以设置 封装参数设置、图标、初始以及文本。归纳起来，即设定 子系统初始值及特性、生成模块图标以及创建模块帮助文 档和描述信息。
10.4 SIMULINK的模型调试
• 如同诸多系统设置平台一样，SIMULINK也具有界 面友好、功能强大的调试功能。SIMULINK提供了 一个图形化的调试界面，简化了调试操作。 simulink调试器是一个定位和诊断模型错误的工 具．通过SIMULINK调试工具，用户可以采用多种 方式对模型进行调试，发现其中可能存在的问题 并加以修改，从而使得用户的模型设计、仿真、 分析更加快速、便捷。
ቤተ መጻሕፍቲ ባይዱ
《matlab函数库》课件
数值计算
支持矩阵和数组运算，统计分 析，数值优化等。
图形可视化
绘制二维和三维图形，支持多 种图形格式和交互操作。
工具箱扩展
提供多种专业领域的工具箱， 如信号处理、控制系统、图像 处理等。
Matlab的应用领域
科学研究
物理、化学、生物、材料科学等领域的研究 工作。
工程设计
机械、电子、航空航天、汽车等领域的建模 、仿真和分析。
04
信号处理函数
01
信号处理函数用于分析和处理各 种信号。
03
使用信号处理函数，您可以对音 频、图像、视频和其他类型的数
据进行信号处理和分析。
02
这些函数包括用于信号滤波、频 谱分析、特征提取和信号生成等
的函数。
201 4
04
这些函数可以帮助您提取信号中 的有用信息，如频率、幅度和相
位等。
图像处理函数
应用示例
数值优化函数在许多领域都有广泛应用，如机器 学习、数据挖掘、控制系统等。例如，在机器学 习中，可以使用数值优化函数来训练神经网络， 以最小化预测误差。
详细描述
数值优化函数包括fminunc、fmincon等，它们 采用不同的算法和约束条件，适用于不同的最优 化问题。这些函数通过不断迭代和调整变量的值 ，逐渐逼近最优解。
03
CATALOGUE
Matlab函数库进阶
控制系统函数
01
控制系统函数用于分析和设计控 制系统。
02
这些函数包括用于系统建模、分 析和控制的函数。
控制系统函数可以帮助您模拟和 控制各种类型的系统，如线性时 不变系统、非线性系统等。
03
使用这些函数，您可以进行系统 稳定性分析、控制设计、信号处
matlab的函数文件
matlab的函数文件
Matlab的函数文件是一种用于存储Matlab代码的文件类型。
它们包含一个或多个函数定义，可以被其他程序或用户调用。
这些函数可以用于执行各种任务，包括数学计算、数据分析和图形绘制等。
Matlab函数文件通常具有以下格式：
function [output1,output2,...] =
myFunction(input1,input2,...)
% Comments about the function
% More comments
% Code for the function
output1 = someOperation(input1,input2);
output2 = anotherOperation(input1,input2);
end
在这个示例中，“myFunction”是函数的名称，“input1”和“input2”是输入参数，“output1”和“output2”是输出参数。
函数体内包含
了代码，执行了一些操作并将结果输出。
在函数文件中，可以包含多个函数定义，但只能有一个主函数。
函数文件通常以“.m”为扩展名，例如“myFunction.m”。
要使用函数文件，可以通过调用函数名称来调用函数，例如“myFunction(input1,input2)”将执行函数并返回输出参数。
函数
文件也可以被其他程序或脚本文件调用。
为了保持代码的可读性，建议在函数文件中添加注释并使用有意义的变量名。
总之，Matlab的函数文件是一种有用的工具，可以帮助用户轻松地组织和重复使用代码，从而提高效率和准确性。
(2021年整理)MatLab常用函数大全
（完整）MatLab常用函数大全编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（（完整）MatLab常用函数大全）的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。
本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快业绩进步，以下为（完整）MatLab常用函数大全的全部内容。
1、求组合数求kC，则输入：nnchoosek（n，k）例：nchoosek(4，2） = 6.2、求阶乘求n！。
则输入：Factorial（n).例：factorial(5） = 120.3、求全排列perms（x).例：求x = ［1，2,3］; Perms(x),输出结果为：ans =3 2 13 1 22 3 12 1 31 2 31 3 24、求指数求a^b:Power（a，b) ；例:求2^3 ；Ans = pow(2，3) ；5、求行列式求矩阵A的行列式：det（A）;例：A=［1 2;3 4］；则det(A） = —2 ;6、求矩阵的转置求矩阵A的转置矩阵：A’转置符号为单引号.7、求向量的指数求向量p=[1 2 3 4]’的三次方：p.^3例：p=[1 2 3 4］'A=[p,p.^2,p。
^3，p。
^4]结果为：注意:在p与符号”^"之间的"。
”不可少.8、求自然对数求ln(x）：Log(x)例:log(2) = 0.69319、求矩阵的逆矩阵求矩阵A的逆矩阵:inv(A)例：a= ［1 2；3 4];则10、多项式的乘法运算函数conv（p1,p2）用于求多项式p1和p2的乘积。
这里，p1、p2是两个多项式系数向量。
例2—2 求多项式43x x23-+的乘积.810+-和2x x命令如下：p1=[1,8,0，0，-10］；p2=［2，—1，3]；c=conv（p1,p2）11、多项式除法函数［q，r］=deconv（p1，p2）用于多项式p1和p2作除法运算,其中q返回多项式p1除以p2的商式,r返回p1除以p2的余式.这里，q和r仍是多项式系数向量。
(完整版)Matlab第十讲教案
西南科技大学本科生课程备课教案计算机技术在安全工程中的应用——Matlab入门及应用授课教师：徐中慧班级：专业：安全技术及工程第十章 MATLAB自定义函数课型：新授课教具：多媒体教学设备，matlab教学软件一、目标与要求✧通过解说与实例练习，掌握matlab创建函数M文件的方法✧掌握matlab中全局变量与局部变量的定义与用法✧通过解说与实例练习，掌握在matlab主函数M文件中创建子函数✧在实例练习过程中，回顾利用伪码编写简单程序的方法✧掌握通过创建matlab函数M文件解决生活中的计算问题二、教学重点与难点本堂课教学的重点在于引导学生掌握matlab中函数M文件的创建及应用。
本堂课的难点在于理解matlab中函数M文件主函数与子函数的区别及调用，局部变量与全局变量的定义与应用范围的区别。
三、教学方法本课程主要通过讲授法、演示法、练习法等相结合的方法来引导学生掌控本堂课的学习内容。
1)通过讲授法向学生讲述创建matlab函数M文件的基本方法、全局变量与局部变量的定义及用法等。
2)通过运用多媒体设备现场演示matlab创建函数M文件的应用实例。
3)在掌握创建matlab函数M文件基本方法的基础上，采用练习法引导学生创建函数M文件解决实际问题。
四、教学内容课后习题五（1）拉力测试装置在测试过程中，被测样本受均匀外力的作用产生形变。
下图中显示的是一组拉力测试数据。
根据以下公式计算应力与形变：00l l F A l σε-=和＝ 其中，σ是产生的应力，单位为lbf/in 2(psi);F 为施加的外力，单位为lbf;A 为样本的截面积，单位为in 2；ε为产生的形变，单位为in/in ；l 为样本的长度；0l 为样本的原始长度。
（a ）测试样本是直径为0.505in 的金属杆，根据直径可以计算出金属杆的截面积，进一步利用所提供的数据计算金属杆的应力和形变。
（b ）以形变为x 轴，应力为y 轴，作x-y 线图。
matlab函数文件的编写和调用
matlab函数文件的编写和调用Matlab是一种广泛应用于科学计算和工程领域的高级编程语言和环境。
它提供了丰富的函数库和工具箱，可以用于各种数学计算、数据分析和可视化等任务。
在Matlab中，函数文件的编写和调用是非常重要的，本文将介绍如何编写和调用Matlab函数文件。
我们需要了解什么是Matlab函数文件。
Matlab函数文件是一种独立的文件，以.m为扩展名，包含了一系列的Matlab命令。
这些命令可以用来实现特定的功能，如计算某个数学函数的值、进行数据处理、实现算法等。
通过编写函数文件，我们可以将复杂的计算任务分解为多个简单的步骤，提高代码的可维护性和重用性。
编写Matlab函数文件的基本步骤如下：1. 打开Matlab编辑器。
在Matlab主界面的“主页”选项卡中，点击“新建脚本”按钮，即可打开Matlab编辑器。
2. 编写函数代码。
在编辑器中，我们可以编写函数的代码。
函数的代码由一系列的语句组成，用于实现特定的功能。
在编写代码时，我们可以使用Matlab提供的函数和操作符，也可以根据需要自定义函数和变量。
在函数代码中，我们可以使用注释来说明代码的目的和实现方法，提高代码的可读性。
3. 保存函数文件。
在编辑器中，点击“文件”菜单中的“保存”选项，选择保存函数文件的路径和名称，并设置文件的扩展名为.m。
保存函数文件后，我们可以在Matlab命令窗口或其他函数文件中调用该函数。
调用Matlab函数文件的方法如下：1. 在命令窗口中直接调用。
在Matlab命令窗口中，我们可以直接输入函数名和参数，即可调用对应的函数。
例如，如果我们有一个名为“myfunc”的函数文件，可以在命令窗口中输入“myfunc(arg1, arg2, ...)”来调用该函数。
2. 在其他函数文件中调用。
在Matlab中，我们可以在一个函数文件中调用其他函数文件。
在调用其他函数文件之前，我们需要确认被调用的函数文件已经保存在Matlab的当前工作目录下或在Matlab 的搜索路径中。
MATLAB如何定义函数文件PPT课件
X function dy = myfun2(t,x,y,z) dy = zeros(3,1); dy(1) = y * z; dy(2) = -x * z; dy(3) = -0.51 * x * y;
-
6
说明
function dy = myfun2(t,y)
odefun
变量属性必须 一一对应！
2、主文件 main3.m
clear;
[T,Y]=ode45(@myfun3,[0,10],[1,1]);
-
8
高阶常微分方程
高阶常微分方程
变量替换 化为
一阶常微分方程组
例： Van der Pol 初值问题
令
x1 y, x2
dy，则原方程可化为 dt
-
参数怎么处理？
用全局变量传递
9
参数传递
1、函数文件 verderpol.m
function dy = myfun2(t,y)
ddyy == z[eryo(s2(3),*1y)(;3)%; dy-myu(s1t)*bey(a3)c;olu.m.n.vector!
dy(1) dy(2)
==-0y-.(y52()11)**y*(y1(y)3(*)3y;)(;2)];
dy(3) = -0.51 * y(1) * y(2);
如果是常微分方程组，y 就是列向量！
dy 必须是列向量，长度为方程组的个数，通常与y的长度相同！
函数中的输入参数和输出参数是形参，名字可以任意取，但必须 满足上述条件。即输入参数有两个，第一个表示自变量，第二个 是由因变量组成的列向量，输出参数必须是列向量。
-
7
例
例：解初值问题：