matlab最速下降法代码

最速下降法（Steepest Descent Method）是一种求解无约束优化问题的基本方法之一，是一种迭代算法。其基本思想是：在每一步中取使函数值下降最快的方向作为搜索方向。下面给出matlab的最速下降法代码：

function [f, x] = steepestDescent(f,x0,epsilon)
% f: 目标函数，x0: 初始点，epsilon: 精度
% f: 目标函数，x0: 初始点，epsilon: 精度
% fprime: 目标函数的梯度
% x，f：迭代过程中产生的最优点和函数值

MAX_ITER = 10000; % 迭代次数上限
alpha = 0.2; % 步长

for i = 1 : MAX_ITER
fprime = gradient(f);
d = -fprime; % 梯度方向
falpha = @(alpha) f(x0+alpha*d); % 定义目标函数在alpha处的取值
alpha = fminunc(falpha,0); % 用Fminunc函数求解alpha
x = x0 + alpha*d;
if abs(f(x)-f(x0)) < epsilon
break; % 达到精度要求
end
x0 = x; % 上一个迭代的x
end
end

该代码利用matlab函数gradient计算目标函数的梯度，并用fminunc函数计算在当前迭代点处的步长alpha。在达到精度要求后，返回当前产生的最优点和最优解。