MATLAB编写Powell算法寻优

这是一个用MATLAB编写的Powell算法的实现。Powell算法是一种非线性优化算法，用于求解无约束的非线性优化问题。这个算法的基本思想是通过组合不同方向的搜索来寻找最小值。

以下是该代码的详细解释：

初始化：函数首先创建一个零矩阵X，用于存储所有的迭代点。然后，创建一个n阶单位矩阵S，用于存储搜索方向。创建一个零向量delta，用于存储每个迭代点的函数下降量。
迭代：在while循环中，函数会进行20次迭代。在每次迭代中，它会遍历每一个搜索方向，计算在当前方向上的最短步长，然后更新迭代点，并计算函数值的下降量。
判断：在每次迭代后，函数会检查新生成的点是否满足精度要求。如果满足，就退出循环。
更新搜索方向：如果新生成的点不满足精度要求，那么函数会生成一个新的搜索方向，并计算反射点。
判别：然后，函数会根据Powell判别法来决定是否采用新的搜索方向。如果采用新的搜索方向，那么旧的搜索方向就会被新的搜索方向替代。否则，函数会选择初始点作为下一次迭代的起点。
这个实现看起来很详细，但是在实际应用中，还需要注意一些细节。例如，该代码没有处理函数不可微或者无解的情况，也没有处理收敛到局部最小值的情况。因此，在实际应用中，可能需要对这个代码进行一些改进。

function [x,fval]=powell(f,x0,e,n)
% Powell 函数f 初始点x0 精度e 未知数个数n

X=zeros(n,n+2); % X矩阵第一列为初始点，第n+1列为Sn方向搜索后的点，最后一列为反射点
X(:,1)=x0; % 将初始点赋值给矩阵X
S=eye(n); % n阶单位矩阵的n列向量分别代表n元方程的n个搜索方向
delta=zeros(1,n); % 预分配内存，用于存放各点的函数下降量

x=x0';
fval=f(x0);
k=0;
while(k<20)  % 最大迭代次数暂定为20
    k=k+1;
    for i=1:n
        syms a;
        a=solve(diff(