function aa

clear;clc;close all
load matlab3
for i=1:3;
    hold on;
    if i==1   %必须是用两个等号。
          figure1=figure('color','w')  
        eta12=eta1(1:5:end,:);

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

        theta=pi/2-eta12(:,3);
        
        myr=linspace(0,0,length(x));
        myg=linspace(0,1,length(x));
        myb=linspace(0,0,length(x));
        for jj=1:length(x)   
           plot(x,y);hold on;
            plotSquare2(x(jj),y(jj),theta(jj),1/30,[myr(jj),myg(jj),myb(jj)]);
            axis equal
            pause(0.1)
    hold off;
        end
        colormap([myr', myg', myb'])
        colorbar
        


    elseif i==2
         figure1=figure('color','w')
        eta22=eta2(1:5:end,:);

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

        theta=pi/2-eta22(:,3);
        
        myr=linspace(0,0,length(x));
        myg=linspace(0,1,length(x));
        myb=linspace(0,0,length(x));
        for jj=1:length(x)
            plot(x,y);hold on;
            plotSquare2(x(jj),y(jj),theta(jj),1/30,[myr(jj),myg(jj),myb(jj)]);
            axis equal
            pause(0.1)
    hold off;
        end
        colormap([myr', myg', myb'])
        colorbar
        
        
    elseif i==3        
          figure1=figure('color','w')
        eta32=eta1(1:5:end,:);

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

        theta=pi/2-eta32(:,3);
        
        myr=linspace(0,0,length(x));
        myg=linspace(0,1,length(x));
        myb=linspace(0,0,length(x));
        for jj=1:length(x)   
            plot(x,y);hold on;
            plotSquare2(x(jj),y(jj),theta(jj),1/30,[myr(jj),myg(jj),myb(jj)]);
            axis equal
            pause(0.1)
    hold off;
        end
        colormap([myr', myg', myb'])
        colorbar
        
    end
   
    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)