function aa

clear;clc;close all
load matlab3
        eta12=eta1(1:5:end,:);

        x=eta12(:,2);
        y=eta12(:,1);
        % theta=eta1(:,3); % ?

        theta1=pi/2-eta12(:,3);
         
        eta22=eta2(1:5:end,:);

        x=eta22(:,2);
        y=eta22(:,1);
        % theta=eta1(:,3); % ?

        theta2=pi/2-eta22(:,3);
         
        eta32=eta1(1:5:end,:);

        x=eta32(:,2);
        y=eta32(:,1);
        % theta=eta1(:,3); % ?

        theta3=pi/2-eta32(:,3);
        
        myr=linspace(0,0,length(x));
        myg=linspace(0.5,1,length(x));
        myb=linspace(0,0,length(x));
        for jj=1:length(x)   
            plot(x,y);hold on;colormap([myr', myg', myb'])
            plotSquare2(x(jj),y(jj),theta1(jj),1/30,[myr(jj),myg(jj),myb(jj)]);
            plotSquare2(x(jj),y(jj),theta2(jj),1/30,[myr(jj),myg(jj),myb(jj)]);
            plotSquare2(x(jj),y(jj),theta3(jj),1/30,[myr(jj),myg(jj),myb(jj)]);
             colorbar
            axis([min(x)-1 max(x)+1 min(y)-1 max(y)+1])
            pause(0.05)
            hold off;
        end

function plotSquare2(x0,y0,theta,a,mycolor)
xs=[-2*a,-2*a,2*a,3*a,2*a,-2*a];
ys=[a,-a,-a,0,a,a];

xy=[cos(theta),-sin(theta);
    sin(theta),cos(theta)]*[xs;ys] ;
fill(xy(1,:)+x0,xy(2,:)+y0,mycolor)

3楼: Originally posted by sogooda at 2012-04-09 23:09:03:
或者自动一点儿用
pause(1)