How to parameterise a function?

6 visualizaciones (últimos 30 días)
Jennifer Wail
Jennifer Wail el 4 de Jun. de 2015
Editada: Stephen23 el 4 de Jun. de 2015
I was trying to use a multi-variable function in commands that take functions as inputs but kept getting errors and was told to parameterise the function but not sure how. The function is f(x,y,z) = x^2 + 3y^2 + 4z^2 - 2xy + 5x - 3y + 2z. I tried doing:
f = @(x,y,z) = x.^2 + 3.*y.^2 + 4.*z.^2 - 2.*x.*y + 5.*x - 3.*y + 2.*z;
ezplot(f)
fminsearch(f,[0,0,0])
and was told to parameterise my functions so that x(1)=x, x(2)=y... Not sure how to continue though. Please help

Respuesta aceptada

Stephen23
Stephen23 el 4 de Jun. de 2015
Editada: Stephen23 el 4 de Jun. de 2015
The standard answer is to replace each x, y and z with the corresponding element of the input vector:
g = @(V) = V(1).^2 + 3.*V(2).^2 + 4.*V(3).^2 - 2.*V(1).*V(2) + 5.*V(1) - 3.*V(2) + 2.*V(3);
A more readable alternative is to use two functions, which allows you to keep the original function unchanged. This could be useful if it is important to keep exactly the same equation as some reference, or to use exactly those variables names, or just to maintain readability:
f = @(x,y,z) = x.^2 + 3.*y.^2 + 4.*z.^2 - 2.*x.*y + 5.*x - 3.*y + 2.*z;
g = @(V) f(V(1),V(2),V(3));
and then call the function g in all of your following work. It has the disadvantage that it would likely be a smidgen slower, so do not use this if your need to perform millions of iterations.
  1 comentario
Jennifer Wail
Jennifer Wail el 4 de Jun. de 2015
Thank you that did get me further but now having trouble with fminsearch
f = @(x) x(1)^2 + 3*x(2)^2 + 4*x(3)^2 - 2*x(1)*x(2) + 5*x(1) - 3*x(2) + 2*x(3)
[x fval] = fminsearch(f,[0,0,0])
Maybe I'll post a second question on that but thank you again you helped me get further on other problems with that

Iniciar sesión para comentar.

Más respuestas (1)

James Tursa
James Tursa el 4 de Jun. de 2015
Editada: James Tursa el 4 de Jun. de 2015
Do pretty much exactly what the advice says:
f = @(x) x(1)^2 + 3*x(2)^2 + 4*x(3)^2 - 2*x(1)*x(2) + 5*x(1) - 3*x(2) + 2*x(3);
  1 comentario
Jennifer Wail
Jennifer Wail el 4 de Jun. de 2015
I tried that but now fminsearch gives me index out of bounds error

Iniciar sesión para comentar.

Categorías

Más información sobre Programming en Help Center y File Exchange.

Community Treasure Hunt

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

Start Hunting!

Translated by