Not enough input arguments in embedded function when using fsolve!

Question

Mohmmad Teymouri 2015-7-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/228813-not-enough-input-arguments-in-embedded-function-when-using-fsolve

编辑： Walter Roberson 2015-7-15

I have a code which itself is a function but I have also many embedded functions in it. The main function defines a system of non-linear equations. Input arguments to this code are [w, del, v, Pg, Qg]. When I define input arguments and pass them to the function there is no problem and I have output. But when I want to solve it using fsolve (with trust-region-dogleg, because my non-linear system of equations is square) I receive "not enough input arguments" error for one of my embedded functions. What should I do? Why I don't get any errors in case of not using fsolve, but I get the mentioned error while using fsolve?

4 个评论
显示 2更早的评论隐藏 2更早的评论

Mohmmad Teymouri 2015-7-8

编辑：Walter Roberson 2015-7-15

在 MATLAB Online 中打开

The code is as follows:

function output = problem( w, del, v, Pg, Qg )
%UNTITLED2 forms the PQ, and Droop equations which are then used in Newton-Trust-Region method t be solved
%   It should be mentioned that since still there are no PV buses in the
%   problem, this method lacks the equations for PV-buses, But the adding
%   procedure is straightforward and can be made.
inputdata;
l=length(input(:,1));  %----number of lines
for m=1:l
q=input(m,2);
for ind=1:13
mat(q,ind)=input(m,ind);         %----Rearranging inputdata
end
end
Ppq(1,:)=zeros(1,3);
Qpq(1,:)=zeros(1,3);
cntrpq=2;
cntrdrp=1;
for m=1:mxndnm
    if mat(m,13)==1
        for n=1:3
            Pl=Pload(w,v);
            P=Pbus(del,v);
            Ppq(cntrpq,n)=Pl(m,n)+P(m,n);
%             -Pgspec(m,n);
              Ql=Qload(w,v);
              Q=Qbus(del,v);
              Qpq(cntrpq,n)=Ql(m,n)+Q(m,n);
  %             -Qgspec(m,n);
          end
          cntrpq=cntrpq+1;
      else 
          if mat(m,13)==3
          for n=1:3
              Pl=Pload(w,v);
              P=Pbus(del,v);
              mp=data(m,4);
              pgen=pg(w,mp);
              Pdrp(cntrdrp,n)=Pl(m,n)+P(m,n)-pgen;
              Ql=Qload(w,v);
              Q=Qbus(del,v);
              nq=data(m,5);
              qgen=qg(v(m,n),nq);
              Qdrp(cntrdrp,n)=Ql(m,n)+Q(m,n)-qgen;
              end
              vab(cntrdrp)=v(m,1)-v(m,2);
              vac(cntrdrp)=v(m,1)-v(m,3);
              delab(cntrdrp)=del(m,1)-del(m,2)-(2*pi/3);
              delac(cntrdrp)=del(m,1)-del(m,3)+(2*pi/3);
              Pgdrp(cntrdrp)=Pg(cntrdrp,1)+Pg(cntrdrp,2)+Pg(cntrdrp,3)-pgen;
              Qgdrp(cntrdrp)=Qg(cntrdrp,1)+Qg(cntrdrp,2)+Qg(cntrdrp,3)-qgen;
          cntrdrp=cntrdrp+1;
          end
      end
  end
  Pdrp
      Qdrp
      Ppq
      Qpq
      vab
      vac
      delab
      delac
      Pgdrp
      Qgdrp
  output=[Pdrp;
      Qdrp;
      Ppq;
      Qpq;
      vab;
      vac;
      delab;
      delac;
      Pgdrp;
      Qgdrp];
  end

Also I provide

w=377
del=zeros(25,3)
v=ones(25,3)
Pg=ones(3,3)
Qg=ones(3,3)

as my initial inputs. The error I get is:

Error using problem (line 26)
Not enough input arguments.
Error in fsolve (line 217)
            fuser = feval(funfcn{3},x,varargin{:});
Caused by:
    Failure in initial user-supplied objective function evaluation. FSOLVE cannot continue.

Torsten 2015-7-8

What are the variables you want fsolve to solve for ?

What are the equations you want fsolve to solve ?

Where is the call to fsolve ?

Best wishes

Torsten.

Mohmmad Teymouri 2015-7-8

fsolve should solve for [w,del,v,Pg,Qg]

It should be noted that w is a scalar. size(del)=size(v)=[25 3], size(Pg)=size(Qg)=[3 3].

The equations are mentioned in output in end of the mentioned codes.

The point is that for this case which is a 25-bus power grid I have 56*3 equations. And also regarding the variables I have 56*3+1 variables. w is unknown and also del is unknown. But for making a reference regarding which the other elements of del will be measured, I set del(1,1)=0.

I don't know also how should I feed this to fsolve, to know that del(1,1)=0 is a known variable and not one of the unknown variables, therefore I have a square problem in which I have 56*3 equations and 56*3 variables.

I call fsolve in command line in MATLAB as follows:

OPTIONS=optimoptions('fsolve','algorithm','trust-region-dogleg');

X=fsolve(@problem, [w,del,v,Pg,Qg],OPTIONS)

请先登录，再进行评论。

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

Answer 1

Brendan Hamm 2015-7-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/228813-not-enough-input-arguments-in-embedded-function-when-using-fsolve#answer_185435

编辑：Brendan Hamm 2015-7-8

在 MATLAB Online 中打开

First off your function to be evaluated should accept only a vector as input. I notice you are providing an initial point:

x0 = [w,del,v,Pg,Qg];

If size(del) = [25 3], and w is a scalar, then you cannot concatenate them into a vector as it seems you are doing:

del = rand(25,3);
w = 4;
x0 = [w, del]
Error using horzcat
Dimensions of matrices being concatenated are not consistent.

What you need to do is rearrange all of the design variables into a vector to be passed into your objective function. Inside your objective function you can then go ahead and extract them into their respective matrices etc..

x0 = [w;del(:);v(:);Pg(:),Qg(:)]; % A big long vector
X = fsolve(@problem,x0,OPTIONS);
function xnew = problem(x)
w = x(1);
del = x(2:76);
del = reshape(del,25,3);
v = x(77:151);
v = reshape(v,25,3);
Pg = x(152:160);
Pg = reshape(Pg,3,3);
Qg = x(161:169);
Qg = reshape(Qg,3,3);
% Rest of your function here
end

4 个评论
显示 2更早的评论隐藏 2更早的评论

Mohmmad Teymouri 2015-7-9

Dear sir, your comments are really useful and helpful for me.

I have done your last comment and till now I have no problem running my code.

But there is a matter which makes my work still incomplete.

When I feed the variables as a (168*1) matrice, and after that sort it in my function (for example to reshape del, v, and so on) this codes are the first lines of my function (after declaring the function name of course).

Because of this in every iteration (when using fsolve) the assigned values are updated to their initial value, and therefore I have no progress in my convergence.

The question is that how could I make fsolve update my special variables (w, del, v, Pg, Qg) in every iteration to their new value? Also I mention that I need to keep del(1,1)=0 (always unchanged).

With best regards and thanks in advance

Brendan Hamm 2015-7-14

在 MATLAB Online 中打开

I should mention that the output of the function problem should be the function evaluation F(x) at the value of x passed in. fsolve will then attempt to find the x which satisfies F(x) = 0. So the definition of problem I give is not complete and it should actually compute the function evaluated at x, I just provide the method by which the objective function needs to be defined (1 input of design variables).

Once the function definition is complete and the anonymous function objFun created, you just need to pass objFun to the fsolve function.

x0 = rand(168,1);
x = fsolve(objFun,x0);
disp(objFun(x)) % Should be close to zero.

So fsolve first passes the values in x0 to your objective function and evaluates it. It then shifts x0 by a small amount to approximate the gradient and choose a direction to move each element of x in. It passes this new value of x back to the objective function and repeats until a termination criteria is met when it returns to you the vector which caused termination.

请先登录，再进行评论。