matlab计算速度很慢,有没有大神帮忙看看该怎么修改？

Question

本人matlab新手，现在这个matlab运行不出来结果，运行速度太慢，有没有大神帮忙看看，非常感谢！
function [best_x,best_fx,count]=DFP(x0,ess)
colormap Jet
% ###########################
syms x1 x2 t;
f=0.5*2000*sqrt((3-x1)^2+(8-x2)^2)+0.5*3000*sqrt((8-x1)^2+(2-x2)^2)+0.75*2500*sqrt((2-x1)^2+(5-x2)^2)+0.75*1000*sqrt((6-x1)^2+(4-x2)^2)+0.75*1500*sqrt((8-x1)^2+(8-x2)^2);
fx=diff(f,x1);%求表达式f对x1的一阶求导
fy=diff(f,x2);%求表达式f对x2的一阶求导
fi=[fx fy];%构造函数f的梯度函数
%初始点的梯度和函数值
g0=subs(fi,[x1 x2],x0);
f0=subs(f,[x1 x2],x0);
H0=eye(2); %输出x0,f0,g0
x0
f0
g0
xk=x0;
fk=f0;
gk=g0;
Hk=H0;
k=1;
while(norm(gk)>ess)%迭代终止条件||gk||<=ess
disp('************************************************************')
disp(['第' num2str(k) '次寻优'])
%确定搜索方向
pk=-Hk*gk';
%由步长找到下一点x(k+1)
xk=xk+t*pk';
f_t=subs(f,[x1 x2],xk); %构造一元搜索的一元函数φ(t) %由一维搜索找到最优步长
df_t=diff(f_t,t);
tk=solve(df_t);
if tk~=0
tk=double(tk);
else
break;
end
%计算下一点的函数值和梯度
xk = subs(xk,t,tk)
fk=subs(f,[x1 x2],xk)
gk0=gk;
gk=subs(fi,[x1 x2],xk)
%DPF校正公式，找到修正矩阵
yk=gk-gk0;
sk=tk*pk';
Hk=Hk;
Hk=Hk-(Hk*yk'*yk*Hk)/(yk*Hk*yk')+sk'*sk/(yk*sk')%修正公式
k=k+1;
end

disp('结果如下：')
best_x=xk;%最优点
best_fx=fk;%最优值
count=k-1;
end

Answer 1

首先，我运行了你的程序，你的程序subs函数处无法运行

大概看了下你的代码，计算不出结果，有一个可能，迭代没收敛，一直在循环中
最后，一点建议，对于计算函数在某点的导数不一定只用符号求导，可以用有限差分法里的高精度离散求导逼近，速度会快很多