【手把手制作三阶魔方模拟器】用MATLAB让你的魔方动起来

已于 2023-06-26 09:02:47 修改
阅读量1.6w 收藏 25
点赞数 6
分类专栏： MATLAB 其他 文章标签： 魔方模拟 MATLAB
于 2021-02-09 15:09:33 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xsz591541060/article/details/113728493
版权
MATLAB 同时被 2 个专栏收录
80 篇文章 116 订阅
订阅专栏
其他
21 篇文章 3 订阅
订阅专栏

【手把手制作三阶魔方模拟器】用MATLAB让你的魔方动起来

by 今天不飞了

有一个酷爱魔方的朋友，托我给他定制一个专门用于“训练魔方观察和预判能力”的程序。听完需求之后觉得很有趣，就答应了，并决定把整个制作过程公开。

今天不飞了，一起写代码吧！

上次已经画好三阶魔方，是时候让它动起来了。话不多说，直接上代码。

1 定义魔方初始状态

function [axisBlock,cornerBlock, edgeBlock] = init(axisPosList,cornerPosList,edgePosList)

% 轴块属性初始化
axisBlock = cell(6,1);
for n = 1:6
    axisBlock{n}.RotateMatrix = [1,0,0;0,1,0;0,0,1];
    axisBlock{n}.position = axisPosList(n,:);
    axisBlock{n}.color = ones(1,6)*7;
    axisBlock{n}.color(n) = n;
end
% 角块属性初始化
cornerBlock = cell(8,1);
for n = 1:8
    cornerBlock{n}.RotateMatrix = [1,0,0;0,1,0;0,0,1];
    cornerBlock{n}.position = cornerPosList(n,:);
    cornerBlock{n}.color = ones(1,6)*7;
end
cornerBlock{1}.color([1,2,5]) = [1,2,5];
cornerBlock{2}.color([1,2,3]) = [1,2,3];
cornerBlock{3}.color([1,3,4]) = [1,3,4];
cornerBlock{4}.color([1,4,5]) = [1,4,5];
cornerBlock{5}.color([6,2,5]) = [6,2,5];
cornerBlock{6}.color([6,2,3]) = [6,2,3];
cornerBlock{7}.color([6,3,4]) = [6,3,4];
cornerBlock{8}.color([6,4,5]) = [6,4,5];
% 棱块属性初始化
edgeBlock = cell(12,1);
for n = 1:12
    edgeBlock{n}.RotateMatrix = [1,0,0;0,1,0;0,0,1];
    edgeBlock{n}.position = edgePosList(n,:);
    edgeBlock{n}.color = ones(1,6)*7;
end
edgeBlock{1}.color([1,2]) = [1,2];
edgeBlock{2}.color([1,3]) = [1,3];
edgeBlock{3}.color([1,4]) = [1,4];
edgeBlock{4}.color([1,5]) = [1,5];
edgeBlock{5}.color([2,3]) = [2,3];
edgeBlock{6}.color([3,4]) = [3,4];
edgeBlock{7}.color([4,5]) = [4,5];
edgeBlock{8}.color([5,2]) = [5,2];
edgeBlock{9}.color([6,2]) = [6,2];
edgeBlock{10}.color([6,3]) = [6,3];
edgeBlock{11}.color([6,4]) = [6,4];
edgeBlock{12}.color([6,5]) = [6,5];

end

2 操作模块

2.1 全操作定义

就是所谓的“U U’ R R’ f r M x” 等等，直接上成品，对推导过程感兴趣的，可以上B站看录播视频（见文末）

function [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(opt, ...
    blockIdx, axisBlock, cornerBlock, edgeBlock)

switch opt
    
    case 1 % F
        R = getRotateMatrix(1,0,0);
        
        for n = blockIdx.cornerIdx([1,2,6,5])
            cornerBlock{n}.RotateMatrix = cornerBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([1,5,9,8])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.cornerIdx([1,2,6,5]) = blockIdx.cornerIdx([2,6,5,1]) ;
        blockIdx.edgeIdx([1,5,9,8]) = blockIdx.edgeIdx([5,9,8,1]);
        
    case -1 % F'
        R = getRotateMatrix(-1,0,0);
        
        for n = blockIdx.cornerIdx([1,2,6,5])
            cornerBlock{n}.RotateMatrix = cornerBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([1,5,9,8])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.cornerIdx([1,2,6,5]) = blockIdx.cornerIdx([5,1,2,6]) ;
        blockIdx.edgeIdx([1,5,9,8]) = blockIdx.edgeIdx([8,1,5,9]);
        
    case 2 % R
        R = getRotateMatrix(0,1,0);
        
        for n = blockIdx.cornerIdx([2,3,7,6])
            cornerBlock{n}.RotateMatrix = cornerBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([2,6,10,5])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.cornerIdx([2,3,7,6]) = blockIdx.cornerIdx([3,7,6,2]) ;
        blockIdx.edgeIdx([2,6,10,5]) = blockIdx.edgeIdx([6,10,5,2]);
        
    case -2 % R'
        R = getRotateMatrix(0,-1,0);
        
        for n = blockIdx.cornerIdx([2,3,7,6])
            cornerBlock{n}.RotateMatrix = cornerBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([2,6,10,5])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.cornerIdx([2,3,7,6]) = blockIdx.cornerIdx([6,2,3,7]) ;
        blockIdx.edgeIdx([2,6,10,5]) = blockIdx.edgeIdx([5,2,6,10]);
        
    case 3 % U
        R = getRotateMatrix(0,0,1);
        
        for n = blockIdx.cornerIdx([5,6,7,8])
            cornerBlock{n}.RotateMatrix = cornerBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([9,10,11,12])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.cornerIdx([5,6,7,8]) = blockIdx.cornerIdx([6,7,8,5]) ;
        blockIdx.edgeIdx([9,10,11,12]) = blockIdx.edgeIdx([10,11,12,9]);
        
    case -3 % U'
        R = getRotateMatrix(0,0,-1);
        
        for n = blockIdx.cornerIdx([5,6,7,8])
            cornerBlock{n}.RotateMatrix = cornerBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([9,10,11,12])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.cornerIdx([5,6,7,8]) = blockIdx.cornerIdx([8,5,6,7]) ;
        blockIdx.edgeIdx([9,10,11,12]) = blockIdx.edgeIdx([12,9,10,11]);
    case 4 % B
        R = getRotateMatrix(-1,0,0);
        
        for n = blockIdx.cornerIdx([3,4,8,7])
            cornerBlock{n}.RotateMatrix = cornerBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([3,7,11,6])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.cornerIdx([3,4,8,7]) = blockIdx.cornerIdx([4,8,7,3]) ;
        blockIdx.edgeIdx([3,7,11,6]) = blockIdx.edgeIdx([7,11,6,3]);
        
    case -4 % B'
        R = getRotateMatrix(1,0,0);
        
        for n = blockIdx.cornerIdx([3,4,8,7])
            cornerBlock{n}.RotateMatrix = cornerBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([3,7,11,6])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.cornerIdx([3,4,8,7]) = blockIdx.cornerIdx([7,3,4,8]) ;
        blockIdx.edgeIdx([3,7,11,6]) = blockIdx.edgeIdx([6,3,7,11]);
        
    case 5 % L
        R = getRotateMatrix(0,-1,0);
        
        for n = blockIdx.cornerIdx([1,5,8,4])
            cornerBlock{n}.RotateMatrix = cornerBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([4,8,12,7])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.cornerIdx([1,5,8,4]) = blockIdx.cornerIdx([5,8,4,1]) ;
        blockIdx.edgeIdx([4,8,12,7]) = blockIdx.edgeIdx([8,12,7,4]);
        
    case -5 % L'
        R = getRotateMatrix(0,1,0);
        
        for n = blockIdx.cornerIdx([1,5,8,4])
            cornerBlock{n}.RotateMatrix = cornerBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([4,8,12,7])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.cornerIdx([1,5,8,4]) = blockIdx.cornerIdx([4,1,5,8]) ;
        blockIdx.edgeIdx([4,8,12,7]) = blockIdx.edgeIdx([7,4,8,12]);
        
        
        
    case 6 % D
        R = getRotateMatrix(0,0,-1);
        
        for n = blockIdx.cornerIdx([1,2,3,4])
            cornerBlock{n}.RotateMatrix = cornerBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([1,2,3,4])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.cornerIdx([1,2,3,4]) = blockIdx.cornerIdx([4,1,2,3]) ;
        blockIdx.edgeIdx([1,2,3,4]) = blockIdx.edgeIdx([4,1,2,3]);
        
    case -6 % D'
        R = getRotateMatrix(0,0,1);
        
        for n = blockIdx.cornerIdx([1,2,3,4])
            cornerBlock{n}.RotateMatrix = cornerBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([1,2,3,4])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.cornerIdx([1,2,3,4]) = blockIdx.cornerIdx([2,3,4,1]) ;
        blockIdx.edgeIdx([1,2,3,4]) = blockIdx.edgeIdx([2,3,4,1]);
        
    case 13 % S
        R = getRotateMatrix(1,0,0);
        
        for n = blockIdx.axisIdx([1,3,6,5])
            axisBlock{n}.RotateMatrix = axisBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([2,10,12,4])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.axisIdx([1,3,6,5]) = blockIdx.axisIdx([3,6,5,1]) ;
        blockIdx.edgeIdx([2,10,12,4]) = blockIdx.edgeIdx([10,12,4,2]);
    case -13 % S'
        R = getRotateMatrix(-1,0,0);
        
        for n = blockIdx.axisIdx([1,3,6,5])
            axisBlock{n}.RotateMatrix = axisBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([2,10,12,4])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.axisIdx([1,3,6,5]) = blockIdx.axisIdx([5,1,3,6]) ;
        blockIdx.edgeIdx([2,10,12,4]) = blockIdx.edgeIdx([4,2,10,12]);
        
    case 14 % M
        R = getRotateMatrix(0,-1,0);
        
        for n = blockIdx.axisIdx([1,4,6,2])
            axisBlock{n}.RotateMatrix = axisBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([1,3,11,9])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.axisIdx([1,4,6,2]) = blockIdx.axisIdx([2,1,4,6]) ;
        blockIdx.edgeIdx([1,3,11,9]) = blockIdx.edgeIdx([9,1,3,11]);
    case -14 % M'
        R = getRotateMatrix(0,1,0);
        
        for n = blockIdx.axisIdx([1,4,6,2])
            axisBlock{n}.RotateMatrix = axisBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([1,3,11,9])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.axisIdx([1,4,6,2]) = blockIdx.axisIdx([4,6,2,1]) ;
        blockIdx.edgeIdx([1,3,11,9]) = blockIdx.edgeIdx([3,11,9,1]);
        
    case 15 % E
        R = getRotateMatrix(0,0,-1);
        
        for n = blockIdx.axisIdx([2,3,4,5])
            axisBlock{n}.RotateMatrix = axisBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([5,6,7,8])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.axisIdx([2,3,4,5]) = blockIdx.axisIdx([5,2,3,4]) ;
        blockIdx.edgeIdx([5,6,7,8]) = blockIdx.edgeIdx([8,5,6,7]);
    case -15 % E'
        R = getRotateMatrix(0,0,1);
        
        for n = blockIdx.axisIdx([2,3,4,5])
            axisBlock{n}.RotateMatrix = axisBlock{n}.RotateMatrix*R;
        end
        for n = blockIdx.edgeIdx([5,6,7,8])
            edgeBlock{n}.RotateMatrix = edgeBlock{n}.RotateMatrix*R;
        end
        blockIdx.axisIdx([2,3,4,5]) = blockIdx.axisIdx([3,4,5,2]) ;
        blockIdx.edgeIdx([5,6,7,8]) = blockIdx.edgeIdx([6,7,8,5]);
        
    case 7 % f
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(1, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % F
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(13, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % S
        
    case -7 % f'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-1, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % F'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-13, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % S'
        
    case 8 % r
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(2, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % R
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-14, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % M'
        
    case -8 % r'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-2, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % R'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(14, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % M
    case 9 % u
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(3, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % U
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-15, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % E'
        
    case -9 % u'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-3, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % U'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(15, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % E
        
    case 10 % b
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-4, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % B
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-13, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % S'
        
    case -10 % b'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-4, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % B'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(13, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % S
        
    case 11 % l
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(5, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % L
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(14, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % M
        
    case -11 % l'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-5, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % L'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-14, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % M'
    case 12 % d
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(6, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % D
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(15, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % E
        
    case -12 % d'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-6, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % D'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-15, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % E'
        
    case 16 % x (y)
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(8, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % r
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-5, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % L'
        
    case -16 % x'(y')
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-8, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % r'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(5, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % L
        
    case 17 % y (z)
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(9, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % u
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-6, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % D'
        
    case -17 % y'(z')
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-9, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % u'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(6, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % D
        
    case 18 % z (x)
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(7, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % f
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-4, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % B'
        
    case -18 % z'(x')
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(-7, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % f'
        [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(4, ...
            blockIdx, axisBlock, cornerBlock, edgeBlock); % B
        
end

end

2.2 旋转公式小函数


function R = getRotateMatrix(alpha,beta,gamma)
% alpha,beta,gamma 为1或-1，分别表示顺时针转和逆时针转

alpha = -pi/2*alpha;
beta = -pi/2*beta;
gamma = -pi/2*gamma;

Rx = [1,0,0;0,cos(alpha),-sin(alpha);0,sin(alpha),cos(alpha)];
Ry = [cos(beta),0,sin(beta);0,1,0;-sin(beta),0,cos(beta)];
Rz = [cos(gamma),-sin(gamma),0;sin(gamma),cos(gamma),0;0,0,1];
R = (Rx*Ry*Rz)';

end

2.3操作识别

咱们直接用魔方手法代号，让它自动转为上面的编号把
注：由于字符串里面不能再写 ’这个符号，所以我们用，代替。比如RUR’U’ ，代码里用RUR,U,。

function [optIdx,optNum] = getOptIdx(optList)

optListNum = length(optList);
optNum = 0;
optIdx = zeros(optListNum,1);
for k = 1:optListNum
    tmp = strfind('FRUBLDfrubldSMExyz',optList(k));
    if ~isempty(tmp)
        optNum = optNum+1;
        optIdx(optNum) = tmp;
    else
        if optNum >0
            optIdx(optNum) = -optIdx(optNum);
        end
    end
end
optIdx(optNum+1:end) = [];

end

3 显示模块

3.1 初始化

画一个未打乱的魔方

function H = initshowcube(axisBlock, cornerBlock, edgeBlock, ...
    blockVertices, blockFace, color, colorAlpha)

% 轴块
for n = 1:6
    for f = 1:6
        V = (axisBlock{n}.position+blockVertices(blockFace{f},:))*axisBlock{n}.RotateMatrix;
        C = color{axisBlock{n}.color(f)};
        A = colorAlpha(axisBlock{n}.color(f));
        H.axis{n}{f} = patch('Faces',[1 2 3 4],'Vertices',V,'FaceColor',C,...
            'FaceAlpha',A);
    end
end


% 棱
for n = 1:12
    for f = 1:6
        V = (edgeBlock{n}.position+blockVertices(blockFace{f},:))*edgeBlock{n}.RotateMatrix;
        C = color{edgeBlock{n}.color(f)};
        A = colorAlpha(edgeBlock{n}.color(f));
        H.edge{n}{f} = patch('Faces',[1 2 3 4],'Vertices',V,'FaceColor',C,...
            'FaceAlpha',A);
    end
end

% 角块
for n = 1:8
    for f = 1:6
        
        V = (cornerBlock{n}.position+blockVertices(blockFace{f},:))*cornerBlock{n}.RotateMatrix;
        C = color{cornerBlock{n}.color(f)};
        A = colorAlpha(cornerBlock{n}.color(f));
        H.corner{n}{f} = patch('Faces',[1 2 3 4],'Vertices',V,'FaceColor',C,...
            'FaceAlpha',A);
    end
end

end

3.2 显示操作后的魔方

function showcube(H, axisBlock, cornerBlock, edgeBlock, ...
    blockVertices, blockFace, color, colorAlpha)

% 轴块
for n = 1:6
    for f = 1:6
        V = (axisBlock{n}.position+blockVertices(blockFace{f},:))*axisBlock{n}.RotateMatrix;
        C = color{axisBlock{n}.color(f)};
        A = colorAlpha(axisBlock{n}.color(f));
        H.axis{n}{f}.Vertices = V;
        H.axis{n}{f}.FaceColor = C;
        H.axis{n}{f}.FaceAlpha = A;
    end
end

% 棱
for n = 1:12
    for f = 1:6
        V = (edgeBlock{n}.position+blockVertices(blockFace{f},:))*edgeBlock{n}.RotateMatrix;
        C = color{edgeBlock{n}.color(f)};
        A = colorAlpha(edgeBlock{n}.color(f));
        H.edge{n}{f}.Vertices = V;
        H.edge{n}{f}.FaceColor = C;
        H.edge{n}{f}.FaceAlpha = A;
    end
end

% 角块
for n = 1:8
    for f = 1:6
        V = (cornerBlock{n}.position+blockVertices(blockFace{f},:))*cornerBlock{n}.RotateMatrix;
        C = color{cornerBlock{n}.color(f)};
        A = colorAlpha(cornerBlock{n}.color(f));
        H.corner{n}{f}.Vertices = V;
        H.corner{n}{f}.FaceColor = C;
        H.corner{n}{f}.FaceAlpha = A;
    end
end

end

4 测试

简单测试

clear; close all; clc

%% 基本定义（魔方几何信息） 尺寸颜色透明度之类的可以根据需要自己修改
blockVertices = [ 1,-1,-1; 1, 1,-1;-1, 1,-1;-1,-1,-1;...
    1,-1, 1; 1, 1, 1;-1, 1, 1;-1,-1, 1]*0.95;
blockFace{1} = [1,2,3,4];
blockFace{2} = [1,2,6,5];
blockFace{3} = [2,3,7,6];
blockFace{4} = [3,4,8,7];
blockFace{5} = [4,1,5,8];
blockFace{6} = [5,6,7,8];

% 6轴
axisPosList = [0, 0,-1; 1, 0, 0; 0, 1, 0; -1, 0, 0; 0,-1, 0; 0, 0, 1]*2;
% 8角
cornerPosList = [1,-1,-1; 1, 1,-1; -1, 1,-1; -1,-1,-1;...
    1,-1, 1; 1, 1, 1; -1, 1, 1; -1,-1, 1]*2;
% 12棱
edgePosList = [1, 0,-1; 0, 1,-1;-1, 0,-1; 0,-1,-1;...
    1, 1, 0;-1, 1, 0;-1,-1, 0; 1,-1, 0;
    1, 0, 1; 0, 1, 1;-1, 0, 1; 0,-1, 1]*2;

% 颜色定义 白 蓝 红 绿 橙 黄
color = {[1,1,1],[0,0.2,0.8],[0.8,0,0],[0,0.8,0.2],[1,0.4,0],[1,1,0],[0.1,0.1,0.1]};
colorAlpha = [1,1,1,1,1,1,0.9];

% 初始槽位
blockIdx.axisIdx = [1,2,3,4,5,6];
blockIdx.cornerIdx = [1,2,3,4,5,6,7,8];
blockIdx.edgeIdx = [1,2,3,4,5,6,7,8,9,10,11,12];

%% 初始化
[axisBlock,cornerBlock,edgeBlock] = init(axisPosList,cornerPosList,edgePosList);

% 显示
figure('Position',[100,100,600,700])
axis off
axis equal
axis([-1,1,-1,1,-1,1]*4)
xlabel('x')
ylabel('y')
zlabel('z')
view([5,2,2])

H = initshowcube(axisBlock, cornerBlock, edgeBlock, ...
    blockVertices, blockFace, color, colorAlpha);

%% 操作
% 设置操作
optList = 'RUUR,U,RUR,U,RU,R,'; 
% 转为编号
[optIdx,optNum] = getOptIdx(optList);
% 循环转动
for k = 1:optNum
	% 执行
    [blockIdx, axisBlock, cornerBlock, edgeBlock] = operation(optIdx(k), blockIdx, axisBlock, cornerBlock, edgeBlock);
    % 显示
    showcube(H, axisBlock, cornerBlock, edgeBlock, ...
        blockVertices, blockFace, color, colorAlpha)
    pause(0.1)
end

效果图

在这里插入图片描述
大家肯定发现了，这样实现的话缺少“动感”，给人感觉就是颜色在变，没有转动的感觉。
没错，下一期就一起来实现“动态旋转”。

其他

  1. 程序有bug，或有更好的提议，欢迎留言告诉我
  2. 代码录播视频见B站【三阶魔方】手把手从零实现，魔方观察训练程序（二）让魔方动起来
今天不飞了
关注
专栏目录
标题：使用MATLAB生成Excel表格
TechInk的博客
09-03 1610
我们设置了表格的表头，并使用MATLAB的数据填充了表格的内容。在MATLAB中，我们可以使用内置的函数和工具来生成和编辑Excel表格。接下来，我们可以使用MATLAB中的数据来填充表格的内容。假设我们有一个包含学生信息的结构体数组，我们可以遍历该数组并将数据逐行写入Excel表格中。我们可以使用Excel对象的方法来创建新的工作簿和工作表。在Excel表格中，通常会有一个表头，用于描述数据的各个列的含义。当我们完成数据填充后，可以使用以下代码保存Excel文件并关闭Excel对象。
读取电子表格数据并创建表格（Matlab）
CyberBladeX的博客
09-19 951
在上面的代码中，我们假设导入的数据存储在名为"ImportedData"的变量中。然后，我们使用Matlab的"table"函数创建一个表格，指定每一列的变量名。如果没有安装，可以通过在Matlab命令窗口中输入"ver"来检查"Import Tool"工具箱是否已安装。在"Import Tool"窗口中，单击"Open"按钮，并浏览到包含要读取的电子表格文件的位置。当您完成预览并进行必要的调整后，单击"Import Selection"按钮将数据导入到Matlab工作环境中。如有任何问题，请随时提问。
3D魔方_3D3阶魔方_
10-03
模拟真实三阶魔方的功能，实现在电脑上玩魔方。
用Matlab生成表格
03-12
这个一个m文件，利用Matlab，直接生成latex代码，可直接插入latex中，生成表格，方便快捷
Robotarium MATLAB 模拟器使用教程
最新发布
gitblog_00067的博客
09-13 781
Robotarium MATLAB 模拟器使用教程 robotarium-matlab-simulator MATLAB simulator for the Robotarium! 项目地址: https://gitcode.co...
MATLAB生成excel表格
热门推荐
HPC_ZY
11-13 5万+
利用excel对数据进行复杂处理有一定局限性，于是选择导入MATLAB进行处理。此贴针对MATLAB将数据导出至excel，列举一些简单的例子，帮助大家学习。
matlab生成word和excel
09-16
里面有详细的程序，而且带有注释方遍学习，易于掌握,简单
matlab 表格生成
m_zhangJingDong的博客
05-17 7363
figure(3); set(3,'Name','Calibration results','NumberTitle','off'); subplot(2,1,1); ssmy=[-2.224834e+02 ;0.000000e+00 ;1.990932e-03; -4.004428e-06; 1.218167e-08]; ssit=[-2.169126e+02 ;0.000000e+00; 1....
【手把手制作三阶魔方模拟器】用MATLAB绘制一个一阶魔方
HPC_ZY
02-04 1553
有一个酷爱魔方的朋友，托我给他定制一个专门用于“训练魔方观察和预判能力”的程序。听完需求之后觉得很有趣，就答应了，并决定把整个制作过程公开。
Matlab之Simulink菜鸟级经典教程（手把手滴让你成为仿真高手）
05-16
本教程“Matlab之Simulink菜鸟级经典教程（手把手滴让你成为仿真高手）”针对初学者，旨在帮助他们快速掌握Simulink的基本操作和核心概念，从而成为一名熟练的仿真工程师。 Simulink是一种基于图形用户界面的建模...
【手把手制作三阶魔方模拟器】用MATLAB绘制一个三阶魔方
HPC_ZY
02-05 1万+
【手把手制作三阶魔方模拟器】用MATLAB绘制一个一阶魔方绘制三阶魔方每个块的中心为每个块六面上色# 其他 by 今天不飞了 有一个酷爱魔方的朋友，托我给他定制一个专门用于“训练魔方观察和预判能力”的程序。听完需求之后觉得很有趣，就答应了，并决定把整个制作过程公开。 今天不飞了，一起写代码吧！ 上次画了一个一阶魔方，这次就把它怼成三阶。话不多说，直接上代码。 绘制三阶魔方每个块的中心 以魔方中心轴为原点，魔方每个块尺寸为222单位， 则整个魔方尺寸为666，那么每个块的中心坐标易推导，我直接上代码，你们自
手把手教你用EXCEL制作日历.pdf
11-16
手把手教你用EXCEL制作日历.pdf
手把手教你天线设计——用MATLAB仿真天线方向图
07-12
手把手教你天线设计——用MATLAB仿真天线方向图 本篇文章主要介绍了天线设计的基本概念和MATLAB仿真天线方向图的方法。文章首先介绍了天线的基本概念，包括天线的定义、分类、基本单元和互易定理等。然后，对于半波...
matlab生成excel文档源码
猪猪侠
02-16 7271
1. 创建一个本地的Excel服务端,其实就是打开Excel try Excel = actxGetRunningServer('Excel.Application'); catch Excel = actxserver('Excel.Application'); end 表达意思如下 try 如果Excel 服务器已经打开，返回其句柄 catch ...
利用MATLAB生成Word和Excel文档
06-11
组件对象模型（COM） MATLAB中的ActiveX控件接口技术 利用MATLAB生成Word文档 利用MATLAB生成Excel文档
matlab 建立excel表格----例子
yyywww666的专栏
09-08 1万+
function ceshi_Excel %利用MATLAB生成Excel文档 %   ceshi_Excel % %   Copyright 2009 - 2010 xiezhh.  %   $Revision: 1.0.0.0 $  $Date: 2009/12/14 20:40:00 $ filespec_user
Matlab Excel开发（VBA）
ChanMon的博客
02-29 1128
Range 对象存储了对指定工作表内某个数据范围的引用。
3阶魔方模拟
11-30
很好用，Scramble是打乱。界面很好，而且还显示已转步数
手把手教程：零基础使用MATLAB完成基于深度学习U-Net模型的遥感影像分类
07-05
所以，本教程手把手教你使用MATLAB中的深度学习框架，完成遥感影像分类的具体操作步骤。本教程给出完整的代码、操作手擦、原始训练及测试样本，旨在最大限度的简化操作步骤，让深度学习零基础的学生也可以轻松跑通...
【matlab】生成excel表格，并读写excel内容
weixin_30292745的博客
12-05 1524
本段程序的主要功能是：按行或者按列依次循环往表格里写数据。 1 function main() 2 clc; 3 close all; 4 clear; 5 6 year=2017;month=12; 7 daynum=eomday(year,month); 8 % date=201712; 9 % daynum=cal_day(date)...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值