p3.m

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

x=-1.5:0.1:1;
u=0.5*x.^2+0.5*((1./(0.5-x)-0.5).^2).*((-x).^0.5);
plot(x,u);
hold on
for n=1:3;
u=0.5*x.^2+0.5*((1./(0.5-x)-0.5).^2).*((-x).^n);
plot(x,u);
end
hold off