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
发布评论