What procedure I should use to minimize this function?
1 次查看(过去 30 天)
显示 更早的评论
amine&&
2016-7-7
How i should proceed to minimize The objective function?
Thanks!
20 个评论
jgg
2016-7-7
You'll probably want to try a couple of things, but since you have a constraint optimization problem, I suggest fmincon first.
It's also a very, very good idea to move this:
%load the dataset
filename = 'BS_1.xlsx';
y = xlsread(filename);
yS=y(1:288);
Out of your function. You want all the data loading and stuff to be done a single time, not every function evaluation. Pass it into your function:
function X=mFunction(alpha,beta,yS)
instead, then optimize an anonymous version instead:
func = @(a,b)mFunction(a,b,yS)
amine&&
2016-7-7
When I use :
>> func = @(a,b)mFunction(a,b,yS)
func =
@(a,b)mFunction(a,b,yS)
>> x = fmincon(func,xnew,A,b,Aeq,beq,lb,ub)
MATLAB launches the following error :
Error using @(a,b)mFunction(a,b,yS)
Not enough input arguments.
Error in fmincon (line 534)
initVals.f = feval(funfcn{3},X,varargin{:});
Caused by:
Failure in initial user-supplied objective function evaluation. FMINCON cannot continue.
Thanks!
jgg
2016-7-7
Think about it this way; when you minimize your function, you're going to be calling the function thousands, if not millions of times.
Every time you call the function, you have to:
1) Create a filename object
2) Read data from the hard drive into memory
3) Turn the loaded data into a matlab matrix
4) Process the matrix in yS
This is load of un-necessary work which is going to make your minimization much, much slower. If it takes 1 second to read it, and you're going to run this a hundred thousand times, you're wasting 27 hours for no reason.
jgg
2016-7-7
Your func takes two input; you're only passing one.
If you need to optimize over both alpha and beta you probably need to write a function like this:
func = @(x)mFunction(x(1),x(2),yS);
instead. But, this requires these to be scalars; I don't know you function well enough. Just be aware that when you optimize you pass a single vector, so your function needs to take a single vector.
amine&&
2016-7-7
编辑:amine&&
2016-7-7
I totally greement with you. We do not load the data for each use of the function. but I do not know is there a way to recover without being required to load every time. especially when I try to start the minimization procedure I obtien the following error message :
>> x = fmincon(func,xnew,A,b,Aeq,beq,lb,ub)
Undefined function or variable 'y'.
Error in mFunction (line 51)
X = mean((y(26:264)-PREV2(26:264)).^2);
Error in @(x)mFunction(x(1),x(2),yS)
Error in fmincon (line 534)
initVals.f = feval(funfcn{3},X,varargin{:});
Thanks!
jgg
2016-7-7
编辑:jgg
2016-7-7
You're passing yS not y. Your function cannot use something you do not call. You could just pass y instead; change everything above to use y instead of yS, and add the line which creates yS to your function.
This is very straightforward; you need to think carefully about what you're trying to do, and then make sure everything works. The goal is to remove your data loading from your function; you will have to tweak what is passed into the function and how the function is written to get this to work, but it's not difficult to do.
Your workflow, when you do optimization is like this:
data = loadData();
func = @(x)myFunction(x,data)
output = doOptimize(func);
amine&&
2016-7-7
In my case. The following code works perfectly
>> lb = [0,0];
>> ub = [1,1];
>> A = [];
>> b = [];
>> Aeq = [];
>> beq = [];
>> x0 = [0.5,0.5];
>> x = fmincon(@(x)mFunction(x(1),x(2)),x0,A,b,Aeq,beq,lb,ub)
MATLAB gives me as solution :
Local minimum found that satisfies the constraints.
Optimization completed because the objective function is non-decreasing in
feasible directions, to within the default value of the function tolerance,
and constraints are satisfied to within the default value of the constraint tolerance.
<stopping criteria details>
x =
0.0026 0.9964
Whereas when I draw the graph in 3D solution is a little meadows equal to [0.99,0.99]. Thanks!
amine&&
2016-7-7
Even if i change the initial point "e.g (0.99,0.99)", its given't the correct result. For tolerance I do not know how to change it
jgg
2016-7-7
Check the function value near the points; are you sure that it is an optimum? You might have a subtle bug in your function you will need to track down.
amine&&
2016-7-7
I checked and this is what I got :
>> mFunction(0.99,0.99)
ans =
6.6666e-04
>> mFunction(0.8,0.8)
ans =
8.7804e-04
>> mFunction(0.6,0.6)
ans =
0.0014
>> mFunction(0.1,0.1)
ans =
0.0037
jgg
2016-7-7
You said that you ran the optimization and it gave you a point (0.0026 0.9964); what is the value of the function at that point?
jgg
2016-7-7
This value is actually smaller than the 0.99, 0.99 point you found, so the optimization is working properly; this IS the minimum of that function. You probably need to look more closely at the function you are optimizing if this isn't what you expect.
amine&&
2016-7-7
I totally agree with you. I am trying to verify this hypothesis. I rependre you later if I find the solution to this problem.
回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!发生错误
由于页面发生更改,无法完成操作。请重新加载页面以查看其更新后的状态。
您也可以从以下列表中选择网站:
如何获得最佳网站性能
选择中国网站(中文或英文)以获得最佳网站性能。其他 MathWorks 国家/地区网站并未针对您所在位置的访问进行优化。
美洲
- América Latina (Español)
- Canada (English)
- United States (English)
欧洲
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
亚太
- Australia (English)
- India (English)
- New Zealand (English)
- 中国
- 日本Japanese (日本語)
- 한국Korean (한국어)