p2.m

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

clf;
x=-1.5:0.1:1.5;
y=x;
[X,Y]=meshgrid(x,y);
Z=0.5.*(1./sqrt((X-0.5).^2+Y.^2)-0.5).^2.*(X.^2+Y.^2)+0.5.*(X.^2+Y.^2);
surf(X,Y,Z)
view(-25,25)
hold on
th=[0:pi/50:2*pi]';
p=1+0.5*cos(th);
q=0.5*sin(th);
plot(p,q)
hold off