旋转的太极图，使用matlab打开绘图新世界

r=2;k=0.25;
delta_theta = pi/60;
theta0=0:pi/60:2*pi;%外面的大圆
theta1=0:pi/30:pi;%中间的分割线需要的半圆
x1=r*cos(theta1);
y1=r*sin(theta1);%x1,y1画最外圈的半个圆
x2=-r*cos(theta1);
y2=-r*sin(theta1);%x2,y2画最外圈的半个圆

x3=cos(0)+k*cos(theta0);
y3=-sin(0)+k*sin(theta0);
x4=-cos(0)+k*cos(theta0);
y4=sin(0)-k*sin(theta0);%x3,y3,x4,y4画中间的两个小圆

x5=cos(theta1)-cos(0);
y5=sin(theta1)+sin(0);

x6=cos(theta1)+cos(0);
y6=-sin(theta1)-sin(0);%x5,y5,x6,y6画中间的分割线

h01=plot(x3,y3,'w',x4,y4,'k');%小圆
hold on;
h02=plot(x5,y5,'k',x6,y6,'k');%分割线   %生成图像
axis equal

f01=fill(x2,y2,'k'); %一个半圆涂成黑色
f02=fill(x5,y5,'k');
f03=fill(x6,y6,'w');%分割线部分那里涂成黑白
line([0,x6],[0,y6],'color','w');%画一条白线遮住黑线
f04=fill(x3,y3,'k');
f05=fill(x4,y4,'w');%中间小圆涂成黑白

h03=plot(x1,y1,'k');%大圆
h04=plot(x2,y2,'k');
%涂色的时候不一定就是想要的地方涂成想要的颜色，多调试就可以啦

F(1)=getframe;%F=getframe 捕获显示在屏幕上的当前坐标区作为影片帧

delete(h01); 
delete(h02);

delete(f01);
delete(f02);
delete(f03);
delete(f04);
delete(f05);

for i=1:120
  x1=r*cos(theta1+delta_theta*i);
  y1=r*sin(theta1+delta_theta*i);
  x2=r*cos(theta1+delta_theta*i);
  y2=r*sin(theta1+delta_theta*i);
  x3=cos(0-delta_theta*i)+k*cos(theta0-delta_theta*i);
  y3=-sin(0-delta_theta*i)+k*sin(theta0-delta_theta*i);
  x4=-cos(0-delta_theta*i)+k*cos(theta0-delta_theta*i);
  y4=sin(0-delta_theta*i)+k*sin(theta0-delta_theta*i);

  x5=-cos(theta1-delta_theta*i)-cos(0-delta_theta*i);
  y5=sin(theta1-delta_theta*i)+sin(0-delta_theta*i);

  x6=cos(theta1-delta_theta*i)+cos(0-delta_theta*i);
  y6=-sin(theta1-delta_theta*i)-sin(0-delta_theta*i);  %和上面一样，设置参数，为画图和涂色做准备
   
  h11=plot(x3,y3,'w',x4,y4,'k');%小圆
  hold on;
  h12=plot(x5,y5,'k',x6,y6,'k');%分割线   %生成图像
  axis equal

  f11=fill(x2,y2,'k'); %一个半圆涂成黑色
  f12=fill(x5,y5,'w');
  f13=fill(x6,y6,'k');%分割线部分那里涂成黑白
  line([0,x5],[0,y5],'color','w');
  f14=fill(x3,y3,'w');
  f15=fill(x4,y4,'k');%中间小圆涂成黑白

  h13=plot(x1,y1,'k');%大圆
  h14=plot(x2,y2,'k');
%涂色的时候不一定就是想要的地方涂成想要的颜色，多调试就可以啦
  F(i+1) = getframe;

  delete(h11); 
  delete(h12); 
  
  delete(f11);
  delete(f12);
  delete(f13);
  delete(f14);
  delete(f15);
end
movie(F);  %movie函数用来播放动画帧