2024年4月28日发(作者:)

MATLA自定义函数及局部变量

2009-11-20 09:17

在开始学习MATLA的时候并没有发现这个软件有着这么强大的功能,随着课题 的

不断深入, 也在逼迫着自己不断的去应用新的公式并开发新的算法, 这就牵涉 到

了如何在MATLA中自定义函数的问题,随之而来自然就是所有编程语言所面 临的

问题,函数调用、局部变量等等。下面就我自己整理的一些心得与大家交流。 希望

对你也有所帮助。

1、编写自定义函数时尽量分以下四部分:

(1) 函数定义行: function[out1,out2,..]=filename(in1,in2,..)

和输出参数个数分别由nargin和nargout两个MATLA保留的变量来给出。

(2) 第一行帮助行,以 %开头,作为 lookfor 指令搜索的行

(3) 函数体说明及有关注解:以( %)开头,用以说明函数的作用及有关内

容。如果不希望显示某段信息,可在它的前面加空行

(4) 函数体:函数体内使用的除返回和输入变量这些在 function 语句中直 接引

用的变量以外的所有变量都是局部变量, 即在该函数返回之后, 这些变量会 自动

在MATLA的工作空间中清除掉。如果希望这些中间变量成为在整个程序中 都起作

用的变量,则可以将它们设置为全局变量。

例如下面就是一个标准的自字义函数

function A=myhilb(n, m)

% MYHILB 是一个示范性的 M-function.

% A=MYHILB(N, M)会生成一个 NX

M

的 Hilbert 矩阵 A.

% A=MYHILB(N会生成一个 NX

N

的 Hilbert 矩阵.

% MYHILB(N,M) 仅仅显示一个 Hilbert 矩阵,而不会返回任何矩阵

%这些内容在用 help 时不会显示

if nargout>1, error('Too many output arguments.'); end

if nargin==1, m=n;

elseif nargin==0 | nargin>2

error('Wrong number of iutput arguments.');

,输入

1 / 4

end

A1=zeros(n,m);

for i=1: n

for j=1:m

A1(i,j)=1/(i+j-1);

end

end

if nargout==1, A=A1;

elseif nargout==0, disp(A1);

end

把该函数保存在MATLAB勺某个搜索路径中

(

比如work目录中

)

,在命令窗口 输

入:

>> help myhilb

MYHILB是一个示范性的 M-function.

A=MYHILB(N, M)会生成一个 NX

M

的 Hilbert 矩阵 A.

A=MYHILB(N会生成一个 NX

N

的 Hilbert 矩阵.

MYHILB(N,M) 仅仅显示一个 Hilbert 矩阵,而不会返回任何矩阵。

>> lookfor myhilb

MYHILB是一个示范性的 M-function.

2、MATLAB^全局变量与局部变量有什么区别?

答:如果一个函数内的变量没有特别声明, 那么这个变量只在函数内部使用,

即为局部变量。 如果两个或多个函数共用一个变量 (或者说在子程序中也要用到 主

程序中的变量, 注意不是参数) ,那么可以用 global 来将它声明为全局变量。

全局变量的使用可以减少参数传递, 合理利用全局变量可以提高程序执行的 效

率。

2 / 4

如果需要用到其他函数的变量。就要利用在主程序与子程序中分别都声明全 局

变量的方式实现变量的传递。否则函数体内使用的都为局部变量。比如下例:

文件名:

exam1.m

global a

x=1:100;

a=2;

文件名:prods.m

fun ctio n

result=prods(x)

global a

result=a*sum(x);

c=prods(x)

>> examl

10100

3、怎样编写递归函数?

答:在建立数学模型时,经常会遇到利用递归表达的式子, 这种数学表达式 简

单明了,在MATLAB^实现该模型,对应的需要建立递归函数。写一个好的递 归函

数,关键是找到递归的边界条件,和在中间某步骤中会发生的各种情况, 然 后再对

应到程序中即可。试看Hanoi塔例子:要将塔座A上的n个圆盘全部转移 到塔座C

上,如n=1,直接把这一片移到目标柱子即可;若 n>1,可以先把塔座 A上的个圆

盘想法转移到塔座 B上,然后把塔座A上的最后一个大圆盘转移到塔 座C上,最后

再把塔座B上的个圆盘转移到塔座 C上。把上面思想对应到递归程 序中,就得到如

下代码:

fun ctio n

result=ha no i(disk nu m,begi npillar,midpillar,e ndpillar,sch)

if disknum==1%添加一行移动方式,递归的边界

result=[sch;1,begi npillar,e ndpillar];

else%下面一句相当于把上面n-1片移到中间柱子

result=ha no i(disk nu m-1,beg in pillar,e ndpillar,midpillar,sch);

%然后把最后一片移到目标柱子上

result=[result;disk nu m,beg in pillar,e ndpillar];

3 / 4

%把中间当作第一根 , 原来第一根当作中间柱子 , 继续移动

result=hanoi(disknum-1,midpillar,beginpillar,endpillar,result);

end

4 / 4