Minimizing a multivariable anonymous function that contains vectors using fminsearch

2017 7 22
1 个回答
更新时间：2017 7 23
3 次查看（30 天）

0 个投票

I am trying to minimize I0 and a in the function handle below. rBes and rCutFit are vectors with the same dimensions. The variables, wave, a, C, u, and r are constants that have already been defined in the script.
R2=@(I0,a) (rBes(u)-I0*((2*rCutFit((((2*pi)/wave)*a*sin(atan((u*C)/r))))/(((2*pi)/wave)*a*sin(atan((u*C)/r))))^.2)).^2; MIN=fminsearch(R2,[1,10^-6]);
I keep getting the error below everytime I run the code. Please let me know if you have any advice for fixing this.
"Not enough input arguments.
Error in @(I0,a)(rBes(u)-I0*((2*rCutFit((((2*pi)/wave)*a*sin(atan((u*C)/r))))/(((2*pi)/wave)*a*sin(atan((u*C)/r))))^.2)).^2
Error in fminsearch (line 189) fv(:,1) = funfcn(x,varargin{:});"

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

回答（1 个）

James Tursa
James Tursa 2017-7-22

0 个投票

Make R2 accept a vector input. E.g., use x as the argument and replace I0 with x(1) and a with x(2)
R2=@(x) (rBes(u)-x(1)*((2*rCutFit((((2*pi)/wave)*x(2)*sin(atan((u*C)/r))))/(((2*pi)/wave)*x(2)*sin(atan((u*C)/r))))^.2)).^2;
MIN=fminsearch(R2,[1,10^-6]);

9 个评论
显示 7更早的评论 隐藏 7更早的评论

Tana Si
Tana Si 2017-7-22
编辑：Tana Si 2017-7-22
I have tried your suggestion. This is my code:
wave=473*10^-9; r=0.143; C=4.022*10^-5; u=30;
R2=@(x) (rBes(u)-x(2)*((2*rCutFit((((2*pi)/wave)*x(2)*sin(atan((u*C)/r))))/(((2*pi)/wave)*x(2)*sin(atan((u*C)/r))))^.2)).^2; MIN=fminsearch(R2,[1,10^-6]);
This is the error I received:
"Subscript indices must either be real positive integers or logicals.
Error in @(x)(rBes(u)-x(2)*((2*rCutFit((((2*pi)/wave)*x(2)*sin(atan((u*C)/r))))/(((2*pi)/wave)*x(2)*sin(atan((u*C)/r))))^.2)).^2
Error in fminsearch (line 189) fv(:,1) = funfcn(x,varargin{:});"
Please let me know if you have any further suggestions. Thanks.
James Tursa
James Tursa 2017-7-22
You didn't replace that 2nd "a" with x(2)
Also, the error you are getting typically indicates you are shadowing a function with a variable. E.g., maybe you have a variable named "sin" or "atan" or "rCutFit" or "rBes" in your workspace that is shadowing those functions.
Tana Si
Tana Si 2017-7-23
So sorry. I replaced the 2nd a and got the same error. rCutFit is a column vector with the same dimensions as rBes.
Tana Si
Tana Si 2017-7-23
I do not have a variable named sin and atan, but I do have rCutFit and rBes as defined column vectors before. rBes is from experimental data and rCutFit is generated through MATLAB's besselj.
Walter Roberson
Walter Roberson 2017-7-23
rCutFit((((2*pi)/wave)*x(2)*sin(atan((u*C)/r))))
mean that you want to calculate (((2*pi)/wave)*x(2)*sin(atan((u*C)/r))) and use that as an index into rCurFit . That is going to fail except under unlikely circumstances that the computation produces a positive integer.
Did you need a multiplication between rCurFit and the following expression?
Tana Si
Tana Si 2017-7-23
I need to express the vector rCutFit, which was produced by the besselj function in MATLAB and reshaped by removing a section of its first values, as a function of ((2*pi)/wave)*sin(atan((u*C)/r)). rCutFit is J1 in the mathematical function attached with its input, kasin(theta)), expressed in a different form.
rBes is an experimentally derived vector and was not generated by MATLAB.
Tana Si
Tana Si 2017-7-23
Please note that in the numerator, kasin(theta) is not a multiplicative factor, but an input into the J1 function (the bessel function of the first kind, first order).
Walter Roberson
Walter Roberson 2017-7-23
Arrays (and vectors) have particular sizes at any one time, and are indexed with () notation and what follows in the () must be a list of positive integers, possibly together with colons.
Functions also use () notation to invoke them, and can have non-integral arguments that are function parameters.
You cannot use non-integral values (other than colon) in the () after a vector name (except in intermediate calculations provided the overall result of the calculation is positive integer).
If you need to perform a calculation on a non-integral argument and get a vector or array output, then you should make rCutFit into a function that does the appropriate bessel calculation and removes whatever and reshapes as is suitable.

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by