p1.m

来自「matlab仿真基于相对位置的势场函数参数的选取」· M 代码 · 共 12 行

x=-1.5:0.1:0.5;
u=0.5*x.^2+0.5*(1./(0.5-x)-0.5).^2;
[u_min,i_min]=min(u);
x_text=['x=',num2str(x(i_min))];
u_text=['u=',num2str(u_min)];
min_text=char('minimum',x_text,u_text);
plot(x,u,'b')
hold on
plot(x(i_min),u_min,'r.','MarkerSize',20)
text(x(i_min),u_min+0.4,min_text)
axis([-1.5,0.5,0,5])
hold off