小人表白MATLAB代码上课讲义.docxVIP

小人表

白MATLAB代

码

functiongaobai

sc=get(0,ScreenSize);

figure(toolbar,none,...position,[sc(3)/450sc(3)*2/4sc(4)-100],...

name,告白CodebyQIAOXINGXING);

set(gca,position,[0011]);

set(gca,xtick,[],ytick,[]);

holdon;

axisequal;

axis([-08-.33.3]);

degree=pi/180;

deg1=[120:-2:60].*degree;

deg2=-deg1;

t1=linspace(-120,-60,length(deg1))*degree;

t2=-pi-t1;

u2=[0,0.5,.5*exp(pi/6*i)+.5];

h1=plot(0);

h2=plot(0);

h3=plot(0);

h4=plot(1);

A=0;

forc=0:5

forb=1:length(deg1)delete(h1);delete(h2);delete(h3);delete(h4);B=exp(i*deg1(b))+A;

C=exp(i*deg2(b))+B;u=[ABC];

E=[.3*exp(i.*(-90:270).*degree),-1.3*i];E=E+1.3*i+B;

u1=u2*exp(i*t1(b))+.8*i+B;u11=u2*exp(i*t2(b))+.8*i+B;

h1=plot(u,linewidth,3);h2=plot(E,linewidth,3);

h3=plot(u1,linewidth,3);h4=plot(u11,linewidth,3);pause(.001)

end

CC=A;A=C;

end

pause(.2)

t11=[-120:2:-15]*degree;

t22=linspace(120,195,length(t11))*degree;

ford=1:length(t11)delete(h3);delete(h4);

u1=u2*exp(i*t11(d))+.8*i+B;

u11=[.5*exp(i*t22(d)),0,.5*exp(-15*i*degree)]-.5*exp(i*t22(d))+.8*i+B;h3=plot(u1,linewidth,3);h4=plot(u11,linewidth,3);

pause(.001)end

hua_t=[-180:9:120]*degree;

hua_u=exp(i*hua_t)+2/sqrt(3)*exp(-i*pi/6);

hua_u1=hua_u*exp(i*120*degree);

hua_u2=hua_u1*exp(i*120*degree);

hua_uu=.09*([hua_u,hua_u1,hua_u2,-i*5]+i*5);

hua_uu1=u1(3)+hua_uu;

t33=[-25:15,15:-1:-25]*degree;

foree=1:4

ifee==4t33=[-25:3:15,15:-3:-20]*degree;end;fore=1:length(t33)

delete(h3);

u1=u2*exp(i*t33(e))+.8*i+B;h3=plot(u1,linewidth,3);pause(.001)

endend

h5=fill(real(hua_uu1),imag(hua_uu1),r);

h6=plot(hua_uu1,r,linewidth,3);

pause(1)

C=CC;

t1=[120:-.5:100]*degree;

t2=linspace(120,180,length(t1))*degree;

t3=linspace(60,0,length(t1))*degree;

forf=1:length(t1)delete(h1);delete(h2);delete(h3);delete(h4);delete(h5);delete(h6);A1=.5*exp(i*t1(f))+A;

B=.5*exp(i*t2(f))+A