以下是MATLAB程序,很简单.我就不写注释了
clc
clear
axis([-2 2 -2 2])
hold on
fai0=pi/12
fai2=pi*(2/3)
fai3=pi/3
fai4=pi
fai5=1.5*pi
fai6=pi*(4/3)
fai7=pi*(-1/6)
fai8=0.01*pi*(0.00002)
u=0.5
r=0.001
t=0.01:0.01:20
d=sqrt(r*r+(abs(u)-r*r)*exp(-2*u*t))
p=(sqrt(abs(u))*r)./d
x=p.*cos(t)
y=p.*sin(t)
plot(x,y,'k')
hold on
x=p.*cos(t+fai4)
y=p.*sin(t+fai4)
plot(x,y,'c')
u=0.5
r=2
d=sqrt(r*r+(abs(u)-r*r)*exp(-2*u*t))
p=(sqrt(abs(u))*r)./d
x=p.*cos(t)
y=p.*sin(t)
plot(x,y,'r')
hold on
%
x=p.*cos(t+fai5)
y=p.*sin(t+fai5)
plot(x,y,'m')
x=p.*cos(t+fai2)
y=p.*sin(t+fai2)
plot(x,y,'k')
x=p.*cos(t+fai6)
y=p.*sin(t+fai6)
plot(x,y,'g')
x=p.*cos(t+fai7)
y=p.*sin(t+fai7)
plot(x,y,'y')
x=p.*cos(t+fai4)
y=p.*sin(t+fai4)
plot(x,y)
x=p.*cos(t+fai3)
y=p.*sin(t+fai3)
plot(x,y,'m')
