三阶魔方自动求解及动态可视化matlab代码

最新推荐文章于 2022-06-01 09:51:48 发布
最新推荐文章于 2022-06-01 09:51:48 发布
阅读量8.1k 收藏 127
点赞数 41
文章标签： 程序人生 经验分享 机器学习 算法 matlab
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xiaoquanquanxiao/article/details/113056028
版权
这篇博客介绍了使用MATLAB实现三阶魔方的自动求解和动态可视化过程，包括魔方状态表示、旋转行为定义、随机打乱和逆序恢复、公式恢复算法。作者通过编写函数对魔方的每一步旋转进行建模，并实现了动态显示魔方旋转的动画。虽然目前的算法仍需改进以减少旋转步数，但已能成功完成魔方的还原。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

三阶魔方自动求解及动态可视化matlab代码


第一次写博客，想总结分享下以前做过的一些有趣的东西，目的是为了回望过去与展望未来，同时为了提高自己的写作表达能力。

思路与步骤

  1. 三阶魔方有6个面，每个面有 3 × 3 3\times3 3×3小块，用一个 6 × 3 × 3 6\times3\times3 6×3×3的矩阵来保存魔方的状态；
  2. 定义4+12种魔方旋转行为：整体旋转（左右上下共4种方式AaCc），侧面（6个侧面顺逆时针旋转共12种方式LlRrUuDdFfBb），旋转可视化界面方便调试（三维小面作图fill3创建 6 × 3 × 3 = 54 6\times3\times3=54 6×3×3=54个小面对象，对象成员属性包括颜色和4个三维坐标点，定义了全局函数对小面进行绕轴旋转操作，如此可以可视化调试，确保程序是按人想法运行），旋转行为改变矩阵的状态矩阵；
  3. 定义随机打乱和逆序恢复的函数，测试确保人为给出旋转公式，魔方能正确旋转；
  4. 编写按公式恢复魔方的函数，记录恢复过程的旋转过程公式（取一个魔方一边拧一边编程，考虑所有情况）
    1)：顶层拼十字：简化为4次恢复顶面棱中小块（不破坏已恢复的结果）
    2)：顶层拼四角：简化为4次恢复顶面的角小块（不破坏已恢复的结果）
    第一层恢复完成，剩下可按公式恢复魔方，把第一层置底
    3)：二层拼棱角：简化为4次恢复二层棱角小块（不破坏已恢复的结果）
    第二层恢复完成
    4)：按公式顶层拼十字
    5)：按公式顶层四角对应恢复
    6)：按公式顶面同色恢复
    7)：按公式顶面棱边恢复
    第三层恢复完成
  5. 编写一个公式简化的函数，如消除相邻正反操作、消除4次相同操作、3次同向旋转替换为一次反向旋转操作、不相关跳跃旋转整合等，未能考虑所有可化简的情况，还未能化简为最少旋转步数公式。时间精力有限，待以后有新思路的时候再玩。

总结

很多指令都是现查现用的，编程模块化，过程分解，编程周期较长，编完以后，整体思路清晰，但单独模块缺少注释，不好看懂，特别是手动旋转对应编程部分，时间久了容易忘记。这是2019年1月编的代码（忘记了具体参考过哪些代码了），现在整理下思路，这里三阶魔方恢复算法是普通的按公式还原算法，主要难点在于判断魔方的状态，我编程的过程中是实际拿着个魔方一边转一边进行状态观察遍历，方法比较笨，可以实现自动魔方还原，但是还原步骤较多，程序中通过增加一些函数进行了步骤的简化
function TurnManu=simplifyTurnManu(TurnManu)
但是简化结果并不彻底，需要更高级的简化函数，或高级的魔方恢复算法。

思考

如果能够用快速相机记录并识别魔方高手旋转魔方的方法，通过机器学习、聚类、迁移等思想，改善算法，数据即算法？如果没有如果，给定初始的算法，能不能让算法自己通过运行自学习，更新算法，或者产生模拟高级的数据用于机器学习？感觉这会不会是一个很好的案例？数据为王？算法为王？Alpha Zero从0开始？有人实现了吗？

参考链接

没作记录，忘了，对参考过代码的博主表示歉意，如有发现雷同的地方，欢迎提醒，以作补充。

魔方

源代码

%% 主调函数
% ---定义魔方的初始状态
face0=ones(3,3);
MFstate={
   face0;face0+1;face0+2;face0+3;face0+4;face0+5};
% % MFstate=faceTurn(MFstate,'r');
% % MFstate=sideTurn(MFstate,'r'); 
[MFstate,TurnManu]=disorganize(MFstate,10);
% DispTurnManu(TurnManu,1);
figure;pltMoFang(MFstate);
set(gcf,'color','k');axis square;axis off;view(160,30)

TurnManu=simplifyTurnManu(TurnManu,0);
TurnManu0=DispTurnManu(TurnManu);
disp({
   '魔方打乱操作：';TurnManu0.'})

% MFstate=disorganize(MFstate,TurnManu);

[MFstate1,TurnManu1]=InvRecover(MFstate,TurnManu);
TurnManu0=DispTurnManu(TurnManu1);
disp({
   '魔方逆序恢复操作：';TurnManu0.'})
figure;pltMoFang(MFstate1);
set(gcf,'color','k');axis square;axis off;view(160,30)

 [MFstate,TurnManu2]=formulaRecover(MFstate);
 TurnManu2=simplifyTurnManu(TurnManu2,0);
TurnManu0=DispTurnManu(TurnManu2);
disp({
   '魔方公式恢复操作：';TurnManu0.'})
figure;pltMoFang(MFstate);
set(gcf,'color','k');axis square;axis off;view(160,30)
-----------------------------------------------------------------------------------------------------------------------------------
%% 动画显示
face0=ones(3,3);
MFstate={
   face0;face0+1;face0+2;face0+3;face0+4;face0+5};
limL=-4;limU=1;
close all;fig=figure; fig.Color=[1,1,1];  fig.Position=[400 200 550 600];
fig.ToolBar='none';  fig.MenuBar='none';  fig.NumberTitle='off';  fig.Name='魔方'; 
hds=pltMoFang(MFstate); 
set(gcf,'color','k');axis([limL,limU,limL,limU,limL,limU]);axis off;view(160,30);
pause(4);
TurnManu_all=[TurnManu;TurnManu2];
for n=1:length(TurnManu_all)
    pause(0.00001);
    turnAroundP(hds,TurnManu_all(n),10,0.01,limL,limU);
    MFstate=disorganize(MFstate,TurnManu_all(n));
    hold off;hds=pltMoFang(MFstate);axis([limL,limU,limL,limU,limL,limU]);axis off;view(160,30);
end
%% 部分子函数编写
%% 公式恢复魔方
function [MFstate,TurnManu]=formulaRecover(MFstate)
TurnManu=[];
        %---第一层恢复
        %将顶层拼十字
        num=zeros(1,4);
        for n=1:4
            [MFstate,TnMn,num]=recoverF10(MFstate,num);TurnManu=[TurnManu;TnMn];
            if sum(num)==4;break;end
            MFstate=faceTurn(MFstate,'r');num=circshift(num,-1);TurnManu=[TurnManu;'a'];
        end
        TurnManu=[TurnManu;('- ').'];
        %将顶层4个角恢复，完成第一层恢复
        num=zeros(1,4);
        for n=1:4
            [MFstate,TnMn,num]=recoverF4(MFstate,num);TurnManu=[TurnManu;TnMn];
            if sum(num)==4;break;end
            MFstate=faceTurn(MFstate,'r');num=circshift(num,-1);TurnManu=[TurnManu;'a'];
        end
        

        TurnManu=[TurnManu;('- ').'];
        %第二层恢复 4个角恢复,完成第二层恢复
        [MFstate,TnMn]=disorganize(MFstate,'CC');TurnManu=[TurnManu;TnMn];
        num=zeros(1,4);
        for n=1:4
            [MFstate,TnMn,num]=recover2F4(MFstate,num);TurnManu=[TurnManu;TnMn];
            if sum(num)==4;break;end
            MFstate=faceTurn(MFstate,'l');num=circshift(num,-1);TurnManu=[TurnManu;'A'];
        end
        
        TurnManu=[TurnManu;('- ').'];
        %第三层公式恢复
        %---步骤1：顶上拼十字
         [MFstate,TnMn]=recover3F10(MFstate);TurnManu=[TurnManu;TnMn];
                 
        TurnManu=[TurnManu;('- ').'];
        %---步骤2：四个角对应上
         [MFstate,TnMn]=recover3F4(MFstate);TurnManu=[TurnManu;TnMn];
         
         TurnManu=[TurnManu;('- ').'];
         %---步骤3：顶面恢复
         [MFstate,TnMn]=recover3Face(MFstate);TurnManu=[TurnManu;TnMn];
         while isempty(TnMn)==0
             [MFstate,TnMn]=recover3Face(MFstate);TurnManu=[TurnManu;TnMn];
         end
         
         TurnManu=[TurnManu;('- ').'];
         %---步骤4：顶棱中恢复
         [MFstate,TnMn]=recover3FArris(MFstate);TurnManu=[TurnManu;TnMn];
         while isempty(TnMn)==0
             [MFstate,TnMn]=recover3FArris(MFstate);TurnManu=[TurnManu;TnMn];
         end
         
end
% ---定义魔方整体旋转行为---改变朝前面
% r,l,u,d 共4种行为
function MFstate=faceTurn(MFstate,manu)
    if manu=='r'
            MFstate={
   MFstate{
   5};rot90(MFstate{
   3},2);MFstate{
   1};fliplr(MFstate{
   4}.');rot90(MFstate{2},2);fliplr(MFstate{6}).'};
    elseif manu=='l'
            MFstate={
   MFstate{
   3};rot90(MFstate{
   5},2);rot90(MFstate{
   2},2);fliplr(MFstate{
   4}).';MFstate{1};fliplr(MFstate{6}.')};
    elseif manu=='u'
            MFstate={
   MFstate{
   4};MFstate{
   6};fliplr(MFstate{
   3}.');MFstate{2};fliplr(MFstate{5}).';MFstate{
   1}};
    elseif manu=='d'
            MFstate={
   MFstate{
   6};MFstate{
   4};fliplr(MFstate{
   3}).';MFstate{1};fliplr(MFstate{5}.');MFstate{
   2}};
    else 
        return;
    end
end
----------------------------------------------------------------------------------------------------------------------------------------------------------
% ---定义一个旋转过程的函数
function turnAroundP(hds,TurnManu,stepNum,tSec,limL,limU)
if nargin==4
    limL=-3;limU=0;
end
if TurnManu=='a'
    [az,el]=view();
    azStep=linspace(0,-90,stepNum);
    for n=1:stepNum
        view(az+azStep(n),el);pause(tSec);
    end

elseif TurnManu=='A'
    [az,el]=view();
    azStep=linspace(0,90,stepNum);
    for n=1:stepNum
        view(az+azStep(n),el);pause(tSec);
    end
  
elseif TurnManu=='c' || TurnManu=='C'
    azStep=90/(stepNum+1);
    if TurnManu=='C'
         azStep=-azStep;
    end
    direct=[1,0,0];
    origin=-[1.5,1.5,1.5];
    for n=1:stepNum
        for k=1:6
            for kk=1:3
                for kkk=1:3
                    h=hds{
   k}{
   kk,kkk};
                    rotate(h,direct,azStep,origin);
                end
            end
        end
        axis([limL,limU,limL,limU,limL,limU]);axis off;view(160,30);pause(tSec);
    end
       
elseif TurnManu=='r' || TurnManu=='R'
    azStep=90/(stepNum+1);
    if TurnManu=='R'
         azStep=-azStep;
    end
    direct=[1,0,0];
    origin=-[1.5,1.5,1.5];
    for n=1:stepNum
        for k=3
            for kk=1:3
                for kkk=1:3
                    h=hds{
   k}{
   kk,kkk};
                    rotate(h,direct,azStep,origin);
                end
            end
        end
        
        for k=[1,2,4,6]
            for kk=1:3
                for kkk=3
                    h=hds{
   k}{
   kk,kkk};
                    rotate(h,direct,azStep,origin);
                end
            end
        end
        axis([limL,limU,limL,limU,limL,limU]);axis off;view(160,30);pause(tSec);
    end
       
elseif TurnManu=='l' || TurnManu=='L'
    azStep=-90/(stepNum+1);
    if TurnManu=='L'
         azStep=-azStep;
    end
    direct=[1,0,0];
    origin=-[1.5,1.5,1.5];
    for n=1:stepNum
        for k=5
            for kk=1:3
                for kkk=1:3
                    h=hds{
   k}{
   kk,kkk};
                    rotate(h,direct,azStep,origin);
                end
            end
        end
        
        for k=[1,2,4,6]
            for kk=1:3
                for kkk=1
                    h=hds{
   k}{
   kk,kkk};
                    rotate(h,direct,azStep,origin);
                end
            end
        end
        axis([limL,limU,limL,limU,limL,limU]);axis off;view(160,30);pause(tSec);
    end
       
elseif TurnManu=='u' || TurnManu=='U'
    azStep=-90/(stepNum+1);
    if TurnManu=='U'
         azStep=-azStep;
    end
    direct=[0,0,1];
    origin=-[1.5,1.5,1.5];
    for n=1:stepNum
        for k=6
            for kk=1:3
                for kkk=1:3
                    h=hds{
   k}{
   kk,kkk};
                    rotate(h,direct,azStep,origin);
                end
            end
        end
        
        for k=[1,3,5]
            for kk=1
                for kkk=1:3
                    h=hds{
   k}{
   kk,kkk};
                    rotate(h,direct,azStep,origin);
                end
            end
        end
 
        for k=2
            for kk=3
                for kkk=1:3
                    h=hds{
   k}{
   kk,kkk};
                    rotate(h,direct,azStep,origin);
                end
            end
        end
        axis([limL,limU,limL,limU,limL,limU]);axis off;view(160,30);pause(tSec);
    end    
       
elseif TurnManu=='d' || TurnManu=='D'
    azStep=90/(stepNum+1);
    if TurnManu=='D'
         azStep=-azStep;
    end
    direct=[0,0,1];
    origin=-[1.5,1.5,1.5];
    for n=1:stepNum
        for k=4
            for kk=1:3
                for kkk=1:3
                    h=hds{
   k}{
   kk,kkk};
                    rotate(h,direct,azStep,origin);
                end
            end
        end
        
        for k=[1,3,5]
            for kk=3
                for kkk=1:3
                    h=hds{
   k}{
   kk,kkk};
                    rotate(h,direct,azStep,origin);
                end
            end
        end
 
        for k=2
            for kk=1
                for kkk=1:3
                    h=hds{
   k}{
   kk,kkk};
                    rotate(h,direct,azStep,origin);
                end
            end
        end
        axis([limL,limU,limL,limU,limL,limU]);axis off;view(160,30);pause(tSec);
    end    
       
elseif TurnManu=='f' || TurnManu=='F'
    azStep=-90/(stepNum+1);
    if TurnManu=='F'
         azStep=-azStep;
    end
    direct=[0,1,0];
    origin=-[1.5,1.5,1.5];
    for n=1:stepNum
        for k=1
            for kk=1:3
                for kkk=1:3
                    h=hds{
   k}{
   kk,kkk};
                    rotate(h,direct,azStep,origin);
                end
            end
        end
            for kk=1:3
                    h=hds{
   3}{
   kk,1};
                    rotate(h,direct,azStep,origin);
            end
            for kk=1:3
                    h=hds{
   5}{
   kk,3};
                    rotate(h,direct,azStep,origin);
            end
            
            for kkk=1:3
                    h=hds{
   4}{
   1,kkk};
                    rotate(h,direct,azStep,origin);
            end
            
            for kkk=1:3
                    h=hds{
   6}{
   3
最低0.47元/天 解锁文章
matlab3阶魔方阵,3*3三阶魔方阵解决方法
weixin_30074929的博客
03-20 2969
C/C++ codeint[][] MagicMatrix ( int n ){//Odd order.if ( mod ( n, 2 ) )[ mj, mi ] = MakeMeshGridSquareMatrix ( n ); // mj : mesh along row,// while mi : mesh along column// both mj and mi are square ...
10月22日matlab学习记录
zaofan93的博客
10-22 310
语法知识 3.13奇数阶魔方矩阵的编程 matlab中魔法矩阵的程序是隐藏的，用type命令显示magic.m。魔方矩阵的设计思路之一，是构造一个扩展的零矩阵A ，随后用顺序数矩阵v填入A的相应对角线，再用行、列叠加使矩阵A的中心n阶矩阵零元素填满。最后取出中心n阶方阵，即得魔方矩阵 魔方阵，古代又称“纵横图”，是指组成元素为自然数1、2、…、n2的平方的n×n的方阵，其中每个元素值都不相等，且每...
三阶,三阶魔方七步还原法,matlab
09-10
简单滑模控制实验
三阶魔方的的智能求解源代码及魔方的数学问题
07-16
魔方求解涉及到数学上群的问题，任意组合的三阶魔方可以在20步内还原，源代码采用AI算法复原魔方。相关的数学论述请看附件。
matlab魔方复原源代码-Rubik-s-Cube:3D魔方解算器和GUI演示
06-13
matlab魔方光源代码魔方 3D Rubik's Cube Solver and GUI Presenting 我的第一个MATLAB程序项目和清晰的操作步骤在'solving-rubiks-cube.pdf'文件中给出。 视频演示显示在“RubiksCube.zip”文件中 Scrambled Rubik's Cube: 魔方动作: Finished Rubik's Cube:
matlab自己还原魔方,复原魔方，原来这么简单（内附详细图文解说）
weixin_32312081的博客
03-22 955
转自：有点儿闲这是一篇魔方复原的教程，看了两个小时视频学会的，然后又自己练习了一遍，果然大神们说得对，三阶魔方很简单，知道规律剩下的就交给练习好了。基础介绍(学习视频教授的叫法，挺形象的，这里就延续原视频的叫法，按照自己的理解复述一遍教程)基础介绍第一步：复原底部十字底部十字第二步：复原第一层复原第一层1①还原底部四个角，以右上角为基准，白色冲右。寻找与白色角块其他面相同颜色的中心块，如图，则白色...
python实现三阶魔方还原
weixin_44946961的博客
04-27 6034
python实现三阶魔方还原 思路 复原魔方困难问题的分解： ​ 1、用合适的数据结构表示出三阶魔方的六个面以及每一面的颜色 ​ 2、每一次不同旋转操作的实现 ​ 3、复原时如何判断当前魔方应该使用哪一种公式 本次实验实现了前两个操作，具体思路是： 用numpy库中的矩阵将六个面分别表示出来，按上下左右前后的顺序放入列表中。再依据流行公式里的方法编写对每一个面进行旋转操作的函数，调用函数实现魔方的旋转。最后输入指令可得到旋转之后的魔方，以及执行逆序指令后验证魔方还原。 预备知识 矩阵：使用nu
MATLAB对三阶魔方建模并进行旋转操作
葬骸缤纷的专栏
09-27 8952
对于魔方的表示：先看一串天书般的字母：UF UR UB UL DF DR DB DL FR FL BR BL UFR URB UBL ULF DRF DFL DLB DBR这种表示法是由一个叫Mike Reid的兄弟首先使用的，它表示一个已经被解好的魔方。先不要被这串字母吓倒，看算法就像追mm一样，要迎难而上。仔细观察，你会发现其中只有六种字母：U: UpF: FrontR: RightL: Le...
MATLAB画三维动态魔方/旋转魔方/旋转立方体
qq_42053235的博客
04-24 9117
以下是我的思路 先了解几个重要的函数 patch()函数 原理：点按顺序连成封闭多边形。 使用： point_sequence=[1,2,3,4]; %点连接的顺序 square_xyz=[3,1,3;1,1,3;1,-1,3;3,-1,3]; %三维坐标4x3矩阵 patch('Faces',point_sequence,'Vertices',square_xyz,'FaceColor','blue'); axis([-3,3,-3,3,-3,3]);%坐标系范围 view(3); rotx
使用matlab实现的3D魔方功能
01-06
里面是使用matlab实现的三维魔方，有按钮可以对魔方进行旋转，魔方的六个面颜色都不同，类似现实中玩的三阶魔方，源代码可以直接运行（.m文件）。
3D自动解魔方源代码
05-01
标题 "3D自动解魔方源代码" 涉及的核心技术是计算机图形学和算法设计，特别是针对3D魔方解决的算法。这个程序能够模拟并解决魔方的复杂问题，展示了C++编程的强大功能。 描述中提到的“可以自动解魔方的程序”，...
MATLAB与C++结合实现三阶魔方三维还原仿真的技术解析
最新发布
04-04
使用场景及目标：适用于希望深入理解魔方求解算法及其可视化实现的学习者；也可用于教学演示，帮助学生更好地掌握相关概念和技术。 其他说明：作者分享了许多实用技巧，包括如何避免连续转动同一面的问题、提高解算...
MATLAB工具箱大全- 魔方模拟器与规划求解
02-09
1. **三维可视化**：真实地展示魔方的外观，允许用户通过点击或键盘命令来旋转各个面。 2. **状态记录**：保存和恢复魔方的任何状态，便于分析和比较。 3. **算法实施**：允许用户输入自定义的解魔方步骤，或者自动...
【手把手制作三阶魔方模拟器】用MATLAB绘制一个三阶魔方
HPC_ZY
02-05 1万+
【手把手制作三阶魔方模拟器】用MATLAB绘制一个一阶魔方绘制三阶魔方每个块的中心为每个块六面上色# 其他 by 今天不飞了 有一个酷爱魔方的朋友，托我给他定制一个专门用于“训练魔方观察和预判能力”的程序。听完需求之后觉得很有趣，就答应了，并决定把整个制作过程公开。 今天不飞了，一起写代码吧！ 上次画了一个一阶魔方，这次就把它怼成三阶。话不多说，直接上代码。 绘制三阶魔方每个块的中心 以魔方中心轴为原点，魔方每个块尺寸为222单位， 则整个魔方尺寸为666，那么每个块的中心坐标易推导，我直接上代码，你们自
MATLAB生成excel表格
热门推荐
HPC_ZY
11-13 5万+
利用excel对数据进行复杂处理有一定局限性，于是选择导入MATLAB进行处理。此贴针对MATLAB将数据导出至excel，列举一些简单的例子，帮助大家学习。
MATALB虚拟魔方构建
yjh_2019的博客
06-01 3556
利用matlab搭建了一个虚拟魔方，并对魔方进行简单的操作。
自己动手做一个虚拟魔方——建模、还原、可视化
qq_45362415的博客
04-02 2832
今天为了教妹妹玩儿魔方，又拿起了初中学会的几个公式。其实我拼魔方很菜的，只是知道几个公式(准确的说是三个…)突发奇想，能不能用python建立一个存在于电脑中的数字化普通三阶魔方。这篇文章记录了我从初步构思到最终实现的过程(如果真能实现…) 由于代码量比较大，所以特意写下这篇博文作为我创建虚拟魔方过程的目录。 虚拟魔方创建过程 建模——对虚拟魔方完成最基本的建模工作，实现各种操作 可视化...
评论 19
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值