旋转的太极

   今天突发奇想，想用matlab实现一个旋转的太极。说干就干。写程序，调程序，费劲周折，终于实现了这个2B的想法。
 

%绘制一个旋转的太极
function taiji
close all
clear
clc
r=1;%太极的大圆半径
figure;
title('太极','FontSize',25)
axis([-r r -r r])  %限制绘图范围
axis vis3d off         %三维坐标修正,关闭坐标轴箭头
hold on
 
[x,y]=circle_draw(0,0,r,0,pi);%上半个大圆
fill(x,y,'w','edgealpha',0)%填充颜色，去掉边框
[x,y]=circle_draw(0,0,r,pi,2*pi);%下半个大圆
fill(x,y,'k','edgealpha',0)
[x,y]=circle_draw(r/2,0,r/2,pi,2*pi);%右面下半个圆
fill(x,y,'w','edgealpha',0)
[x,y]=circle_draw(-r/2,0,r/2,0,pi);%左面上半个圆
fill(x,y,'k','edgealpha',0)
[x,y]=circle_draw(r/2,0,0.2*r,0,2*pi);%右边小圆
fill(x,y,'k','edgealpha',0)
[x,y]=circle_draw(-r/2,0,0.2*r,0,2*pi);%左边小圆
fill(x,y,'w','edgealpha',0)
 
for j=1:5%旋转5次
    for i=0:360
        view(i,90)
        pause(0.01)%调节旋转速度
    end
end
 


function [x,y]=circle_draw(x0,y0,r,alpha,beta)%返回一个弧心为(x0,y0)，半径为r，起始角度为alpha,终止角度为beta的一段弧的坐标
angle=linspace(alpha,beta,100);%将弧的范围等分100份
x=r*cos(angle)+x0;
y=r*sin(angle)+y0;

这是最终的结果，还算可以