数值计算程序的通用性——函数名作为参数传递 在进行数值算法的程序设计时,往往需要实现程序的通用性,即同一个程序可以对不同的问题进行求解,如:计算函数f(x)的数值积分的程序,我们需要写一个通用的数值积分程序——Integr(),使其可以对不同的f(x)函数均可以进行积分计算,而不需要每计算一个函数的积分就重新写一个Integr()程序,仅需要另写一个f(x)的表达式程序即可。
若使用Matlab实现,则显然少不了eval(),和feval()函数的使用,即把需要进行积分计算的函数f(x)写完之后保存为.m文件,然后在积分程序中只需要把f(x)的文件名作为参数传递给Integr()即可。 若是把f(x)函数名作为字符串传入,则在Integr()中只需将f(x)函数名及其参数写为一个字符串,使用eval(字符串)执行即可,如:
function result=integr(a,b,n,fun)
%定义积分函数,a为积分下限,b为积分上限,n为划分区间的个数,fun为需进行计算的函数的文件名
h=(b-a)/n; x=[a:h:b]; y=eval([fun,'(x)']);
%将函数名fun和其参数x合成一个字符串并执行字符串,即执行计算fun(x),注意不能丢了fun的括号。
result=sum(y);
定义需要积分计算的函数fun之后,直接作为integr()的参数引用即可,如任意定义函数如下:
function r=ff(x) r=sin(x).*(1+x.*x);
则,在命令窗口引用即可:将'ff'换为其他函数名字符串即可进行计算其数值积分。
integr(0,1,1000,'ff')
———————————————————————————————————————————————
使用feval()函数则更加方便,integr()函数如下定义:
function result=integr(a,b,n,fun) ;
h=(b-a)/n; x=[a:h:b]; y=fun(x);%fun为传入的函数指针,可以直接作为引用
result=sum(y);
则可以通过如下方式使用:
integr(0,1,1000,@ff)%‘@’即表示ff作为函数指针传入