Automatically define anonymous functions.

Question

丰源 2022-10-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1830578-automatically-define-anonymous-functions

评论： Bjorn Gustavsson 2022-10-19

采纳的回答： Bjorn Gustavsson

Here are my code:

clear

C=[-1,1,1+1/4*i,-1+1/4*i];

syms x;

f=@(x) x.^2 + 1/4;

df=diff(f,x);

fun = df/f;

%fun=@(x) (df/f);

fun =@(x) ((2.*x)./(x.^2 + 1/4));

a=abs(integral( fun,1,1,'Waypoints',C))

I need to compute complex line integrals by the function integral(fun,xmin,xmax,Name,Value). And this function requires that the fun parameter must have an function handler. But my input parameter fun is defined by df/f, which will be changed through variable f.

If I use the annotated code fun=@(x) (df/f), the MATLAB will report an error that "The input function must return a 'double' or 'single' value. Find 'sym'."

Therefore I have to define the anonymous function manually, which is very annoying. Is there any method for automatically defining the anonymous functions? Please help me.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Bjorn Gustavsson 2022-10-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1830578-automatically-define-anonymous-functions#answer_1078818

编辑：Bjorn Gustavsson 2022-10-19

在 MATLAB Online 中打开

You can solve this problem using matlabFunction to generate a function-handle:

C=[-1,1,1+1/4*i,-1+1/4*i];
syms x;
f=@(x) x.^2 + 1/4;
df=diff(f,x);
fun = df/f;
fun = matlabFunction(fun);
% Returns a function-handle:
% fun =
% 
%   function_handle with value:
% 
%     @(x)(x.*2.0)./(x.^2+1.0./4.0)
% suitable for integration:
a=abs(integral( fun,1,1,'Waypoints',C))

If you have more parameters in your symbolic expression, then those will be input-arguments to the function-handle, which makes the integration-step a fair bit more klunky. Then you'll have to figure out which order the arguments come in and how to convert that function-handle to a function-handle with the correct fixed and variable input-arguments. But this should solve this problem.

HTH