How to Call a function inside for loop?

Question

Vaswati Biswas 2021-11-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1590654-how-to-call-a-function-inside-for-loop

编辑： Stephen23 2021-11-19

I am writing a function like this

function Eqn = eff(neff,hf)
nf=vpa(2.1511);
ns=vpa(1.5264);
nc=vpa(1.3354);
rho=1;
lambda=532.3;
Eqn = (((2*pi./lambda)*(sqrt(nf^2-neff.^2)*hf))-atan((nf/nc)^(2*rho)*sqrt ((neff.^2-nc^2)/(nf^2-neff.^2)))-atan((nf/ns)^(2*rho)*sqrt ((neff^2-ns^2)/(nf^2-neff.^2))) -(m*pi))==0;% equation 8 from paper
end

then I am calling it inside a for loop using the following code

hf = 75:20:350;
m= 0;
NEFF = zeros(size(hf)); 
for ii=1:numel(hf)
options = optimoptions('fsolve','Display','none','PlotFcn',@optimplotfirstorderopt);
neff = [0,5];
NEFF(ii)=fsolve(@(neff,hf(ii))eff(neff),neff,options
end
plot(hf, NEFF)

But I am getting an error

Invalid expression. When calling a function or indexing a variable, use parentheses. Otherwise, check for
mismatched delimiters.

what mistake I am making here? Kindly help

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Stephen23 2021-11-19

编辑：Stephen23 2021-11-19

在 MATLAB Online 中打开

Your code is missing severay parentheses and that anonymous function is not defined correctly.

Replace this

NEFF(ii)=fsolve(@(neff,hf(ii))eff(neff),neff,options

with something like

fnh  = @(n)eff(n,hf(ii));
NEFF(ii) = fsolve(fnh,neff,options);

But I cannot run your code as the objective function needs to be debugged first.

请先登录，再进行评论。

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

Answer 1

Sulaymon Eshkabilov 2021-11-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1590654-how-to-call-a-function-inside-for-loop#answer_835549

在 MATLAB Online 中打开

Here is the corrected code, but still you should check the equation. The found roots are complex valued:

hf = 75:20:350;
NEFF = zeros(size(hf)); 
for ii=1:numel(hf)
options = optimoptions('fsolve','Display','none','PlotFcn',@optimplotfirstorderopt);
neff = [0; 5];
NEFF(ii)=fsolve(@(hf) eff(hf(ii), neff), neff,options);
end
%plot(hf, NEFF)
function F = eff(hf, neff)
nf=(2.1511);
ns=(1.5264);
nc=(1.3354);
rho=(1);
lambda=(532.3);
m=(0);
F = (((2*pi./lambda)*(sqrt(nf^2-neff.^2)*hf))-atan((nf/nc)^(2*rho)*sqrt ((neff.^2-nc^2)/(nf^2-neff.^2)))-atan((nf/ns)^(2*rho)*sqrt ((neff.^2-ns^2)/(nf^2-neff.^2))) -(m*pi));% equation 8 from paper
end