Not enough input arguments in embedded function when using fsolve!

2015 7 7
2 个回答
回答已采纳
更新时间：2015 7 15
3 次查看（30 天）

0 个投票

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
Mohmmad Teymouri 2015-7-8
编辑：Walter Roberson 2015-7-15
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
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
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)

请先登录，再进行评论。

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

 采纳的回答

Brendan Hamm
Brendan Hamm 2015-7-8
编辑：Brendan Hamm 2015-7-8

0 个投票

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
Mohmmad Teymouri 2015-7-8
Thanks a lot for your very useful recommendation.
I did it, and I am just one step until I run my code perfectly true.
The problem is that I have 169 variables, and 168 equations. Regarding trust-region-dogleg method which is one of three ways fsolve uses for solving non-linear system of equations, I should have a square problem.
Also I should mention that the first element of my del(25*3) is a known variable. But I don't know how to declare it as a known variable so fsolve can recognize it and just consider 168 variables out of 169 variables as unknown.
By known I mean it has fixed value of zero ==> del(1,1)=0 (always)
Brendan Hamm
Brendan Hamm 2015-7-8
Any parameters which are used in the objective function must be passed as separate arguments.
function xnew = problem(x,del1)
% Shift the previous extraction
w = x(1);
del = [del1;x(2:75)]; % Concatenate the fixed value
del = reshape(del,25,3);
v = x(76:150);
v = reshape(v,25,3);
Pg = x(151:159);
Pg = reshape(Pg,3,3);
Qg = x(160:168);
Qg = reshape(Qg,3,3);
% Rest of your function here
end
Now when we call fsolve it expects the only input to the objective function to be the vector of design variables. To get around this we use what we refer to as anonymous function handles (denoted by @).
objFun = @(x) problem(x,0); % Create anonymous function handle
Now objfun behaves just like the function problem with the second inputs value (0) hard coded into it.
x = rand(168,1);
xNew = objFun(x);
xNew2 = problem(x,0);
all(xNew == xnew2);
ans =
1
Mohmmad Teymouri
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
Brendan Hamm 2015-7-14
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.

请先登录，再进行评论。

更多回答（1 个）

Walter Roberson
Walter Roberson 2015-7-8

0 个投票

fsolve passes in a vector of values. Your routine expects 4 values. Try
fsolve(@(x) problem(x(1),x(2),x(3),x(4)),......

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by