假如我們編寫了一個積分通用程序,想使它更具有通用性,那么可以把被積函數(shù)也作為一個參數(shù)。在c/c++中,可以使用函數(shù)指針來實現(xiàn)上邊的功能,在matlab中如何實現(xiàn)呢?使用函數(shù)句柄——這時類似于函數(shù)指針的一項功能。
假如我們現(xiàn)在已經(jīng)編寫了一個辛普森方法積分的函數(shù)如下,我們要把函數(shù)作為參數(shù)傳給它(當然啦,matlab早已經(jīng)提供給我們非常牛x的積分函數(shù)了,這里只是個例子而已):
function g=myintegrate(f,N,a,b)
%integration using simpson method
n=N/2;
h=(b-a)/N;
x=a:h:b;
fv=f(x);
s=sum(fv(1:2:2*n-1))*4;
s=s+sum(fv(2:2:2*n))*2;
s=s+f(a)+f(b);
s=s*h/3;
g=s;
現(xiàn)在要做的工作就是如何把上面的f函數(shù)換成我們想要的隨便什么函數(shù)。
定義函數(shù)句柄:
handle=@ functionname
之后就可以把handle作為參數(shù)傳遞給上面的函數(shù)(或者其他需要使用函數(shù)作為參數(shù)的函數(shù))了。
例如在0~pi/2上對sin(x)進行積分,可以如下輸入:
>>fh=@sin;
>>myintegrate(fh,100,0,pi)
ans =
0.9947
函數(shù)句柄可以指定為任意的matlab函數(shù)——當然也包括我們自定義的函數(shù),比如我們在myfun.m文件中定義了如下函數(shù):
function f=myfun(x)
f=x;
那么,
>>fh2=@myfun;
>>myintegrate(fh2,100,0,1)
ans =
0.4967