fzero: Not Enough Input Arguments
显示 更早的评论
EDIT: I have edited this section of my question as advised by Ammer Hamza.
Problem Statement:
To minimise an underconstrained linear symbolic system equations in the form of:
a1*X1i + a2*X2i + a3*X3i + a4*X4i + a5*X5i = Yi
(where i=1 to 100)
I have 5 coefficients (X1 to X5), which I have as inputs from a spreadsheet, 5 unknown variables (a1 to a5), and 100 equations.
======================================================================
Hi, I have referred to the other threads on the same error message but I am still unable to solve my problem. I have an underconstrained linear symbolic equation generated from other variables that I would like to minimise. Below is a simplified example of the script I have:
EDITED after comments from Alan Weiss:
a=sym('a',[2,1]);
a1=a(1);a2=a(2);
eqn_a(1)=a(1)+a(2)-3;
eqn_a(2)=2*a(1)+a(2)-4;
eqn_a(3)=a(1)-a(2)+1;
fun=matlabFunction(eqn_a);
x0=[1 10];
fzero(fun,x0)
This is the error message I get:
Error using fzero (line 241)
FZERO cannot continue because user-supplied function_handle ==> @(a1,a2)[a1+a2-3.0,a1.*2.0+a2-4.0,a1-a2+1.0] failed with the error
below.
Not enough input arguments.
Additionally, is this a suitable function or the best function to solve this problem? I have 5 coefficients, 5 variables and 100 equations.
Thank you very much.
3 个评论
Ameer Hamza
2018-6-27
编辑：Ameer Hamza
2018-6-27
Although there are several errors with that use of fminunc() in your code, but it appears that in trying to explain those errors, we will again run into a case of XY problem. fminunc () might be a correct function for this task, but the objective function you are defining is surely not correct. It appears that you are trying to solve some system of equations and are interested in the least square solution. Please explain the actual problem you are trying to solve so that we can suggest the relevant solution.
Ameer Hamza
2018-6-27
Now the problem is clear. This problem has a much easier solution in MATLAB as compared to the using syms, writing all equations and using fminunc(). Please refer to my answer.
采纳的回答
Ameer Hamza
2018-6-27
Since your number of equations (100) is greater than the number of variables (5), the best you can expect is a least square solution, i.e. formulate this as a least square optimization problem as you rightly mentioned in your question. Fortunately solving such system for unconstrained cases require just one line in MATLAB. You don't need to use any function such as fminunc().
X = [x1 x2 x3 x4 x5]; % values from spreadsheet, it will be a [100 x 5] matrix
Y = y; % value of |y| from the spreadsheet.
coefficients = X\Y; % that's it
coefficients will contain the value of a1 to a5 in order.
2 个评论
Ameer Hamza
2018-6-28
fzero can only be used when an exact solution to your system of equation exist. Since you are solving it as a least square problem, fzero is not the right tool. There are several other tools which you can look into, few are listed below with examples for your matrices X and Y,
- lsqlin()
coefficients = lsqlin(X,Y,[],[])
- lsqcurvefit()
f = @(a, X) X*a;
coefficients = lsqcurvefit(f, [0;0;0;0], X, Y)
- fitlm()
coff = fitlm(X,Y,'y~x1+x2+x3+x4+x5-1')
fitlm() will provide much more details than other methods but the all of them will give the same estimate for a since all of them are solving the same problem.
更多回答（1 个）
Alan Weiss
2018-6-27
1 个投票
As formulated, your problem does not make sense.
- You give a scalar value of x0, but you have a 2-D problem.
- You give a 3-output function eqn_a. when objective functions must be scalar-valued.
- You are attempting to solve a set of equations using fminunc, when fsolve is the appropriate solver for a system of equations.
- You do not give the equations in the correct syntax for fsolve.
- For systems of linear equations, you should probably use backslash ( mldivide )
Alan Weiss
MATLAB mathematical toolbox documentation
类别
在 帮助中心 和 File Exchange 中查找有关 Common Operations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
