What am I doing wrong!?

2012 2 7
4 个回答
回答已采纳
2 次查看（30 天）

0 个投票

Hi,
I have the following code, it's part of a larger file that is trying (trying being the crucial word here) to perform optimisation using powell's method:
% testing
clear
a = zeros(1,1);
x0 = [0.2, 0.4, 0.6];
s1 = [0.4, 0.4, 1.6];
dx1 = a*x0;
x1 = x0 + dx1;
func = (x1(1)-x1(2))^2+2*(x1(2)-x1(3))^2+3*(x1(3)-1)^2;
[a] = feval(func, x1(1), x1(2), x1(3));
What I am trying to do is find the value of 'a' that minimises 'func'?
Thanks

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

 采纳的回答

Kevin Holst
Kevin Holst 2012-2-7

0 个投票

Ah I see the problem now. In your original post you had
s1 = [0.4, 0.4, 1.6];
but nothing using s1 in it. Now that I see what s1 is used for here's the solution:
x0 = [0.2, 0.4, 0.6];
s1 = [0.4, 0.4, 1.6];
x1 = @(n,a) x0(n) + s1(n)*a;
func = @(a)(x1(1,a)-x1(2,a))^2+2*(x1(2,a)-x1(3,a))^2+3*(x1(3,a)-1)^2;
a = fminsearch(func, 0);
a = 0.1364

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Will
Will 2012-2-7
Thank you so much! That's exactly what I wanted. Are you able to explain what having the (n) does? I would understand if it was in a 'for' loop.

请先登录，再进行评论。

更多回答（3 个）

Sean de Wolski
Sean de Wolski 2012-2-7

0 个投票

Syntax issues:
a = zeros(1,1);
x0 = [0.2, 0.4, 0.6];
s1 = [0.4, 0.4, 1.6];
dx1 = a*x0;
x1 = x0 + dx1;
func = @(x)(x(1)-x(2))^2+2*(x(2)-x(3))^2+3*(x(3)-1)^2;
a = fminsearch(func, x1(1:3));

2 个评论
显示 无 隐藏 无

Will
Will 2012-2-7
Thank you very much!
Is 'feval' still used?
What does @(x) do to the 'func' line?
Will
Will 2012-2-7
Also, that code works for finding the values of x that minimise func. I am trying to find the value of 'a' which is a scaler multiplier.

请先登录，再进行评论。

Matt Kindig
Matt Kindig 2012-2-7

0 个投票

This is not the approach you want to follow. First, you are calling feval, which will not minimize anything, only evaluate a specified. Instead, you will need one of the solvers in the Optimization Toolbox that is designed to minimize a function, such as fminunc. You will pass a function handle into fminunc that defines your function with the variable you wish to change (a) as the input parameter. Something like this should work:
x0 = [0.2, 0.4, 0.6];
s1 = [0.4, 0.4, 1.6]; %note that you don't use this anywhere
func = @(a) ((a+1)*(x0(1)-x0(2))).^2 + 2*((a+1)*(x0(2)-x0(3))).^2 + 3*( (a+1)*(x0(3)-1)).^2;
a = fminunc(func, 0);

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Will
Will 2012-2-7
Thanks for the replies but neither of those suggestions give me the answer I am trying to obtain, probably because I haven't explained it. This is the maths of what I want to achieve.
Original function f1=(x1-x2)^2+2*(x2-x3)^2+3*(x3-1)^2 X1=X0+a1*S1 X1=[0.2 0.4 0.6]+a1*[0.4 0.4 0.6]=[0.2+0.4a1 0.4+0.4a1 0.6+1.6a1]
Substituting the new values of x1,x2,x3 into the original function yields:
fa1=-0.2^2+2*(-0.2-1.2a1)^2+3*(-0.4+1.6a1)^2
Minimisation of this function should return a value of a1=0.1364
Thanks

请先登录，再进行评论。

Kevin Holst
Kevin Holst 2012-2-7

0 个投票

I'd suggest using the following code, it's not pretty, but it gets the job done without any toolboxes.
x0 = [0.2, 0.4, 0.6];
x1 = @(n,a) (1+a)*x0(n);
func = @(a)(x1(1,a)-x1(2,a))^2+2*(x1(2,a)-x1(3,a))^2+3*(x1(3,a)-1)^2;
a = fminsearch(func, 0);

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Will
Will 2012-2-7
Thanks for the replies but neither of those suggestions give me the answer I am trying to obtain, probably because I haven't explained it. This is the maths of what I want to achieve.
Original function f1=(x1-x2)^2+2*(x2-x3)^2+3*(x3-1)^2 X1=X0+a1*S1 X1=[0.2 0.4 0.6]+a1*[0.4 0.4 0.6]=[0.2+0.4a1 0.4+0.4a1 0.6+1.6a1]
Substituting the new values of x1,x2,x3 into the original function yields:
fa1=-0.2^2+2*(-0.2-1.2a1)^2+3*(-0.4+1.6a1)^2
Minimisation of this function should return a value of a1=0.1364
Thanks

请先登录，再进行评论。

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by