savefig.asv

一个计算CHUA混沌系统Lyapunov指数MATLAB程序。

        % Plot U Region Devision Planes
        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        vertex_matrix = [-xlimit -ylimit sys_variables(7);-xlimit ylimit sys_variables(7);xlimit -ylimit sys_variables(7);xlimit ylimit sys_variables(7)];
        G_ABC_pp(5) = patch('Vertices',vertex_matrix,'Faces',[1 2 4 3],'EdgeColor','b','FaceColor','none','FaceAlpha',0.2,'EdgeAlpha',0.2);
        text(xlimit,ylimit,sys_variables(7),'  U_{1}','Color',[0 0 1],'FontWeight','bold');
        
        %Plot P+ Outer EigenPlane
        %%%%%%%%%%%%%%%%%%%%%% here sys_variables(7) = E
        inner_x1 = pointplane(Eigen_Planes(2,1),Eigen_Planes(2,2),Eigen_Planes(2,3),-Eigen_Planes(2,4),0,ylimit,sys_variables(7),'x');
        inner_x2 = pointplane(Eigen_Planes(2,1),Eigen_Planes(2,2),Eigen_Planes(2,3),-Eigen_Planes(2,4),0,-ylimit,sys_variables(7),'x');
        inner_y1 = pointplane(Eigen_Planes(2,1),Eigen_Planes(2,2),Eigen_Planes(2,3),-Eigen_Planes(2,4),xlimit,0,sys_variables(7),'y');
        inner_y2 = pointplane(Eigen_Planes(2,1),Eigen_Planes(2,2),Eigen_Planes(2,3),-Eigen_Planes(2,4),-xlimit,0,sys_variables(7),'y');
        inner_z = sys_variables(7);
        
        point_number = 1;
        if (inner_x1 >= -xlimit)&&(inner_x1 <= xlimit)
            point(point_number,:) = [inner_x1 ylimit inner_z];
            P_ABC_PP(plot_number) = plot3(point(point_number,1),point(point_number,2),point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            mirror_point(point_number,:)=[(2*Points(1,3)-point(point_number,1)) (2*Points(2,3)-point(point_number,2)) (2*Points(3,3)-point(point_number,3))];
            P_ABC_PP(plot_number) = plot3(mirror_point(point_number,1),mirror_point(point_number,2),mirror_point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            point_number = point_number+1;
        end
        if (inner_x2 >= -xlimit)&&(inner_x2 <= xlimit)
            point(point_number,:) = [inner_x2 -ylimit inner_z];
            P_ABC_PP(plot_number) = plot3(point(point_number,1),point(point_number,2),point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            mirror_point(point_number,:)=[(2*Points(1,3)-point(point_number,1)) (2*Points(2,3)-point(point_number,2)) (2*Points(3,3)-point(point_number,3))];
            P_ABC_PP(plot_number) = plot3(mirror_point(point_number,1),mirror_point(point_number,2),mirror_point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            point_number = point_number+1;
        end
        if (inner_y1 >= -ylimit)&&(inner_y1 <= ylimit)
            point(point_number,:) = [xlimit inner_y1 inner_z];
            P_ABC_PP(plot_number) = plot3(point(point_number,1),point(point_number,2),point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            mirror_point(point_number,:)=[(2*Points(1,3)-point(point_number,1)) (2*Points(2,3)-point(point_number,2)) (2*Points(3,3)-point(point_number,3))];
            P_ABC_PP(plot_number) = plot3(mirror_point(point_number,1),mirror_point(point_number,2),mirror_point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            point_number = point_number+1;
        end
        if (inner_y2 >= -ylimit)&&(inner_y2 <= ylimit)
            point(point_number,:) = [-xlimit inner_y2 inner_z];
            P_ABC_PP(plot_number) = plot3(point(point_number,1),point(point_number,2),point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            mirror_point(point_number,:)=[(2*Points(1,3)-point(point_number,1)) (2*Points(2,3)-point(point_number,2)) (2*Points(3,3)-point(point_number,3))];
            P_ABC_PP(plot_number) = plot3(mirror_point(point_number,1),mirror_point(point_number,2),mirror_point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            point_number = point_number+1;
        end
        % Plot P+ outer Plane
        if point_number ~= 1
            vertex_matrix = [point(1,:); point(2,:); mirror_point(1,:); mirror_point(2,:)];
            G_ABC_pp(6) = patch('Vertices',vertex_matrix,'Faces',[1 2 3 4],'EdgeColor','c','FaceColor','none','FaceAlpha',0.2,'EdgeAlpha',0.2);
            text(vertex_matrix(4,1),vertex_matrix(4,2),vertex_matrix(4,3),'  E^{c}(P^{+})','Color','c','FontWeight','bold');
           
            %text((vertex_matrix(2,1)+vertex_matrix(4,1))/1,(vertex_matrix(2,2)+vertex_matrix(4,2))/1,(vertex_matrix(2,3)+vertex_matrix(4,3))/1,'  E^{c}P_{+}(0)','Color','c','FontWeight','bold');
            %text((xlimit/5),(ylimit/5),pointplane(Eigen_Planes(1,1),Eigen_Planes(1,2),Eigen_Planes(1,3),Eigen_Planes(1,4),(xlimit/5),(ylimit/5),0,'z'),'  E^{c}P_{+}(0)','Color','m','FontWeight','bold');
        end
        
        %Plot P- Outer EigenPlane
        %%%%%%%%%%%%%%%%%%%%%% here sys_variables(7) = E
        inner_x1 = pointplane(Eigen_Planes(2,1),Eigen_Planes(2,2),Eigen_Planes(2,3),Eigen_Planes(2,4),0,ylimit,-sys_variables(7),'x');
        inner_x2 = pointplane(Eigen_Planes(2,1),Eigen_Planes(2,2),Eigen_Planes(2,3),Eigen_Planes(2,4),0,-ylimit,-sys_variables(7),'x');
        inner_y1 = pointplane(Eigen_Planes(2,1),Eigen_Planes(2,2),Eigen_Planes(2,3),Eigen_Planes(2,4),xlimit,0,-sys_variables(7),'y');
        inner_y2 = pointplane(Eigen_Planes(2,1),Eigen_Planes(2,2),Eigen_Planes(2,3),Eigen_Planes(2,4),-xlimit,0,-sys_variables(7),'y');
        inner_z = -sys_variables(7);
        
        point_number = 1;
        if (inner_x1 >= -xlimit)&&(inner_x1 <= xlimit)
            point(point_number,:) = [inner_x1 ylimit inner_z];
            P_ABC_PP(plot_number) = plot3(point(point_number,1),point(point_number,2),point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            mirror_point(point_number,:)=[(2*Points(1,1)-point(point_number,1)) (2*Points(2,1)-point(point_number,2)) (2*Points(3,1)-point(point_number,3))];
            P_ABC_PP(plot_number) = plot3(mirror_point(point_number,1),mirror_point(point_number,2),mirror_point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            point_number = point_number+1;
        end
        if (inner_x2 >= -xlimit)&&(inner_x2 <= xlimit)
            point(point_number,:) = [inner_x2 -ylimit inner_z];
            P_ABC_PP(plot_number) = plot3(point(point_number,1),point(point_number,2),point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            mirror_point(point_number,:)=[(2*Points(1,1)-point(point_number,1)) (2*Points(2,1)-point(point_number,2)) (2*Points(3,1)-point(point_number,3))];
            P_ABC_PP(plot_number) = plot3(mirror_point(point_number,1),mirror_point(point_number,2),mirror_point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            point_number = point_number+1;
        end
        if (inner_y1 >= -ylimit)&&(inner_y1 <= ylimit)
            point(point_number,:) = [xlimit inner_y1 inner_z];
            P_ABC_PP(plot_number) = plot3(point(point_number,1),point(point_number,2),point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            mirror_point(point_number,:)=[(2*Points(1,1)-point(point_number,1)) (2*Points(2,1)-point(point_number,2)) (2*Points(3,1)-point(point_number,3))];
            P_ABC_PP(plot_number) = plot3(mirror_point(point_number,1),mirror_point(point_number,2),mirror_point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            point_number = point_number+1;
        end
        if (inner_y2 >= -ylimit)&&(inner_y2 <= ylimit)
            point(point_number,:) = [-xlimit inner_y2 inner_z];
            P_ABC_PP(plot_number) = plot3(point(point_number,1),point(point_number,2),point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            mirror_point(point_number,:)=[(2*Points(1,1)-point(point_number,1)) (2*Points(2,1)-point(point_number,2)) (2*Points(3,1)-point(point_number,3))];
            P_ABC_PP(plot_number) = plot3(mirror_point(point_number,1),mirror_point(point_number,2),mirror_point(point_number,3),'.', ...
                'MarkerSize',10,'Color',[0 0 0]);hold on;
            plot_number = plot_number + 1;
            point_number = point_number+1;
        end

        % Plot P- outer Plane
        if point_number ~= 1
            vertex_matrix = [point(1,:); point(2,:); mirror_point(1,:); mirror_point(2,:)];
            G_ABC_pp(7) = patch('Vertices',vertex_matrix,'Faces',[1 2 3 4],'EdgeColor','c','FaceColor','none','FaceAlpha',.2,'EdgeAlpha',0.2);
            
            text(vertex_matrix(4,1),vertex_matrix(4,2),vertex_matrix(4,3),'  E^{c}(P^{-})','Color','c','FontWeight','bold');
            %text((vertex_matrix(2,1)+vertex_matrix(4,1))/1.5,(vertex_matrix(2,2)+vertex_matrix(4,2))/1.5,(vertex_matrix(2,3)+vertex_matrix(4,3))/1.5,'  E^{c}P_{-}(0)','Color','c','FontWeight','bold');
            %text((xlimit/5),(ylimit/5),pointplane(Eigen_Planes(1,1),Eigen_Planes(1,2),Eigen_Planes(1,3),Eigen_Planes(1,4),(xlimit/5),(ylimit/5),0,'z'),'  E^{c}(0)','Color','m','FontWeight','bold');
        end
        
        axis([-xlimit xlimit -ylimit ylimit -zlimit zlimit]);
        

        
        % Color Axes - Necessary after reset it would seem
        %set(G_ABC_pp(1),'XColor',[1 0 0],'YColor',[0 1 0],'ZColor',[0 0 1]);