python 分类_python解多面体分类

最新推荐文章于 2023-02-23 08:00:00 发布
最新推荐文章于 2023-02-23 08:00:00 发布
阅读量120 收藏
点赞数
文章标签： python 分类

[0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][0, 0, 0][0, 1, 3][0, 1, 4][0, 1, 5][0, 1, 6][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][0, 2, 3][0, 2, 4][0, 2, 5][0, 2, 6][0, 2, 7][0, 0, 0][0, 3, 1][0, 3, 2][8, 8, 8][0, 3, 4][0, 3, 5][0, 3, 6][0, 3, 7][0, 0, 0][0, 4, 1][0, 4, 2][0, 4, 3][8, 8, 8][0, 4, 5][0, 4, 6][0, 4, 7][0, 0, 0][0, 5, 1][0, 5, 2][0, 5, 3][0, 5, 4][8, 8, 8][0, 5, 6][0, 5, 7][0, 0, 0][0, 6, 1][0, 6, 2][0, 6, 3][0, 6, 4][0, 6, 5][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][0, 7, 2][0, 7, 3][0, 7, 4][0, 7, 5][0, 0, 0][8, 8, 8][0, 0, 0][8, 8, 8][0, 0, 0][1, 0, 3][1, 0, 4][1, 0, 5][1, 0, 6][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][0, 0, 0][8, 8, 8][1, 3, 4][1, 3, 5][1, 3, 6][0, 0, 0][0, 0, 0][8, 8, 8][0, 0, 0][1, 4, 3][8, 8, 8][1, 4, 5][1, 4, 6][0, 0, 0][0, 0, 0][8, 8, 8][0, 0, 0][1, 5, 3][1, 5, 4][8, 8, 8][1, 5, 6][0, 0, 0][0, 0, 0][8, 8, 8][0, 0, 0][1, 6, 3][1, 6, 4][1, 6, 5][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][2, 0, 3][2, 0, 4][2, 0, 5][2, 0, 6][2, 0, 7][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][8, 8, 8][2, 3, 4][2, 3, 5][2, 3, 6][2, 3, 7][0, 0, 0][0, 0, 0][8, 8, 8][2, 4, 3][8, 8, 8][2, 4, 5][2, 4, 6][2, 4, 7][0, 0, 0][0, 0, 0][8, 8, 8][2, 5, 3][2, 5, 4][8, 8, 8][2, 5, 6][2, 5, 7][0, 0, 0][0, 0, 0][8, 8, 8][2, 6, 3][2, 6, 4][2, 6, 5][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][2, 7, 3][2, 7, 4][2, 7, 5][0, 0, 0][8, 8, 8][0, 0, 0][3, 0, 1][3, 0, 2][8, 8, 8][3, 0, 4][3, 0, 5][3, 0, 6][3, 0, 7][0, 0, 0][8, 8, 8][0, 0, 0][8, 8, 8][3, 1, 4][3, 1, 5][3, 1, 6][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][8, 8, 8][3, 2, 4][3, 2, 5][3, 2, 6][3, 2, 7][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][3, 4, 1][3, 4, 2][8, 8, 8][8, 8, 8][3, 4, 5][3, 4, 6][3, 4, 7][0, 0, 0][3, 5, 1][3, 5, 2][8, 8, 8][3, 5, 4][8, 8, 8][3, 5, 6][3, 5, 7][0, 0, 0][3, 6, 1][3, 6, 2][8, 8, 8][3, 6, 4][3, 6, 5][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][3, 7, 2][8, 8, 8][3, 7, 4][3, 7, 5][0, 0, 0][8, 8, 8][0, 0, 0][4, 0, 1][4, 0, 2][4, 0, 3][8, 8, 8][4, 0, 5][4, 0, 6][4, 0, 7][0, 0, 0][8, 8, 8][0, 0, 0][4, 1, 3][8, 8, 8][4, 1, 5][4, 1, 6][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][4, 2, 3][8, 8, 8][4, 2, 5][4, 2, 6][4, 2, 7][0, 0, 0][4, 3, 1][4, 3, 2][8, 8, 8][8, 8, 8][4, 3, 5][4, 3, 6][4, 3, 7][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][4, 5, 1][4, 5, 2][4, 5, 3][8, 8, 8][8, 8, 8][4, 5, 6][4, 5, 7][0, 0, 0][4, 6, 1][4, 6, 2][4, 6, 3][8, 8, 8][4, 6, 5][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][4, 7, 2][4, 7, 3][8, 8, 8][4, 7, 5][0, 0, 0][8, 8, 8][0, 0, 0][5, 0, 1][5, 0, 2][5, 0, 3][5, 0, 4][8, 8, 8][5, 0, 6][5, 0, 7][0, 0, 0][8, 8, 8][0, 0, 0][5, 1, 3][5, 1, 4][8, 8, 8][5, 1, 6][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][5, 2, 3][5, 2, 4][8, 8, 8][5, 2, 6][5, 2, 7][0, 0, 0][5, 3, 1][5, 3, 2][8, 8, 8][5, 3, 4][8, 8, 8][5, 3, 6][5, 3, 7][0, 0, 0][5, 4, 1][5, 4, 2][5, 4, 3][8, 8, 8][8, 8, 8][5, 4, 6][5, 4, 7][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][5, 6, 1][5, 6, 2][5, 6, 3][5, 6, 4][8, 8, 8][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][5, 7, 2][5, 7, 3][5, 7, 4][8, 8, 8][0, 0, 0][8, 8, 8][0, 0, 0][6, 0, 1][6, 0, 2][6, 0, 3][6, 0, 4][6, 0, 5][8, 8, 8][0, 0, 0][0, 0, 0][8, 8, 8][0, 0, 0][6, 1, 3][6, 1, 4][6, 1, 5][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][6, 2, 3][6, 2, 4][6, 2, 5][8, 8, 8][0, 0, 0][0, 0, 0][6, 3, 1][6, 3, 2][8, 8, 8][6, 3, 4][6, 3, 5][8, 8, 8][0, 0, 0][0, 0, 0][6, 4, 1][6, 4, 2][6, 4, 3][8, 8, 8][6, 4, 5][8, 8, 8][0, 0, 0][0, 0, 0][6, 5, 1][6, 5, 2][6, 5, 3][6, 5, 4][8, 8, 8][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][7, 0, 2][7, 0, 3][7, 0, 4][7, 0, 5][0, 0, 0][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][7, 2, 3][7, 2, 4][7, 2, 5][0, 0, 0][8, 8, 8][0, 0, 0][0, 0, 0][7, 3, 2][8, 8, 8][7, 3, 4][7, 3, 5][0, 0, 0][8, 8, 8][0, 0, 0][0, 0, 0][7, 4, 2][7, 4, 3][8, 8, 8][7, 4, 5][0, 0, 0][8, 8, 8][0, 0, 0][0, 0, 0][7, 5, 2][7, 5, 3][7, 5, 4][8, 8, 8][0, 0, 0][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0]---------------------------

[0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][0, 0, 0][0, 1, 3][0, 1, 4][0, 1, 5][0, 1, 6][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][0, 2, 3][0, 2, 4][0, 2, 5][0, 2, 6][0, 2, 7][0, 0, 0][0, 3, 1][0, 3, 2][8, 8, 8][0, 3, 4][0, 3, 5][0, 3, 6][0, 3, 7][0, 0, 0][0, 4, 1][0, 4, 2][0, 4, 3][8, 8, 8][0, 4, 5][0, 4, 6][0, 4, 7][0, 0, 0][0, 5, 1][0, 5, 2][0, 5, 3][0, 5, 4][8, 8, 8][0, 5, 6][0, 5, 7][0, 0, 0][0, 6, 1][0, 6, 2][0, 6, 3][0, 6, 4][0, 6, 5][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][0, 7, 2][0, 7, 3][0, 7, 4][0, 7, 5][0, 0, 0][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][1, 0, 3][1, 0, 4][1, 0, 5][1, 0, 6][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][1, 3, 4][1, 3, 5][1, 3, 6][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][1, 4, 3][8, 8, 8][1, 4, 5][1, 4, 6][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][1, 5, 3][1, 5, 4][8, 8, 8][1, 5, 6][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][1, 6, 3][1, 6, 4][1, 6, 5][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][2, 0, 3][2, 0, 4][2, 0, 5][2, 0, 6][2, 0, 7][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][2, 3, 4][2, 3, 5][2, 3, 6][2, 3, 7][0, 0, 0][0, 0, 0][0, 0, 0][2, 4, 3][8, 8, 8][2, 4, 5][2, 4, 6][2, 4, 7][0, 0, 0][0, 0, 0][0, 0, 0][2, 5, 3][2, 5, 4][8, 8, 8][2, 5, 6][2, 5, 7][0, 0, 0][0, 0, 0][0, 0, 0][2, 6, 3][2, 6, 4][2, 6, 5][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][2, 7, 3][2, 7, 4][2, 7, 5][0, 0, 0][8, 8, 8][0, 0, 0][3, 0, 1][3, 0, 2][0, 0, 0][3, 0, 4][3, 0, 5][3, 0, 6][3, 0, 7][0, 0, 0][8, 8, 8][0, 0, 0][0, 0, 0][3, 1, 4][3, 1, 5][3, 1, 6][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][0, 0, 0][3, 2, 4][3, 2, 5][3, 2, 6][3, 2, 7][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][3, 4, 1][3, 4, 2][0, 0, 0][8, 8, 8][3, 4, 5][3, 4, 6][3, 4, 7][0, 0, 0][3, 5, 1][3, 5, 2][0, 0, 0][3, 5, 4][8, 8, 8][3, 5, 6][3, 5, 7][0, 0, 0][3, 6, 1][3, 6, 2][0, 0, 0][3, 6, 4][3, 6, 5][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][3, 7, 2][0, 0, 0][3, 7, 4][3, 7, 5][0, 0, 0][8, 8, 8][0, 0, 0][4, 0, 1][4, 0, 2][4, 0, 3][0, 0, 0][4, 0, 5][4, 0, 6][4, 0, 7][0, 0, 0][8, 8, 8][0, 0, 0][4, 1, 3][0, 0, 0][4, 1, 5][4, 1, 6][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][4, 2, 3][0, 0, 0][4, 2, 5][4, 2, 6][4, 2, 7][0, 0, 0][4, 3, 1][4, 3, 2][8, 8, 8][0, 0, 0][4, 3, 5][4, 3, 6][4, 3, 7][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][4, 5, 1][4, 5, 2][4, 5, 3][0, 0, 0][8, 8, 8][4, 5, 6][4, 5, 7][0, 0, 0][4, 6, 1][4, 6, 2][4, 6, 3][0, 0, 0][4, 6, 5][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][4, 7, 2][4, 7, 3][0, 0, 0][4, 7, 5][0, 0, 0][8, 8, 8][0, 0, 0][5, 0, 1][5, 0, 2][5, 0, 3][5, 0, 4][0, 0, 0][5, 0, 6][5, 0, 7][0, 0, 0][8, 8, 8][0, 0, 0][5, 1, 3][5, 1, 4][0, 0, 0][5, 1, 6][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][5, 2, 3][5, 2, 4][0, 0, 0][5, 2, 6][5, 2, 7][0, 0, 0][5, 3, 1][5, 3, 2][8, 8, 8][5, 3, 4][0, 0, 0][5, 3, 6][5, 3, 7][0, 0, 0][5, 4, 1][5, 4, 2][5, 4, 3][8, 8, 8][0, 0, 0][5, 4, 6][5, 4, 7][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][5, 6, 1][5, 6, 2][5, 6, 3][5, 6, 4][0, 0, 0][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][5, 7, 2][5, 7, 3][5, 7, 4][0, 0, 0][0, 0, 0][8, 8, 8][0, 0, 0][6, 0, 1][6, 0, 2][6, 0, 3][6, 0, 4][6, 0, 5][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][0, 0, 0][6, 1, 3][6, 1, 4][6, 1, 5][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][6, 2, 3][6, 2, 4][6, 2, 5][0, 0, 0][0, 0, 0][0, 0, 0][6, 3, 1][6, 3, 2][8, 8, 8][6, 3, 4][6, 3, 5][0, 0, 0][0, 0, 0][0, 0, 0][6, 4, 1][6, 4, 2][6, 4, 3][8, 8, 8][6, 4, 5][0, 0, 0][0, 0, 0][0, 0, 0][6, 5, 1][6, 5, 2][6, 5, 3][6, 5, 4][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][7, 0, 2][7, 0, 3][7, 0, 4][7, 0, 5][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][8, 8, 8][7, 2, 3][7, 2, 4][7, 2, 5][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][7, 3, 2][8, 8, 8][7, 3, 4][7, 3, 5][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][7, 4, 2][7, 4, 3][8, 8, 8][7, 4, 5][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][7, 5, 2][7, 5, 3][7, 5, 4][8, 8, 8][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0][0, 0, 0]

import numpy as np

v=np.mat([[0,0,0,0],[1,0,0,0],[0,1,0,0],[0,0,1,0],[0,0,0,1],[2,-1,-1,-1],[1,1,0,0],[-1,0,0,0]]);

A=np.mat([[1,7],[1,2],[6,7]]);

C=np.mat([[2,3,4,5],[3,4,5,6]]);

z=7;

def f(x,y,z):

if x==y or y==z:

return 0

else:

return x

def g(x,y,z):

if x==y:

return z

else:

return y

def h(x,y,z,k):

if x==y:

return [z,k]

elif x==z:

return [y,k]

else:

return [y,z]

def K(w,i,j):

t=v[g(w,i,j)]

#   print(t)

#   print(v[w])

#   print(np.dot(v[w],np.transpose(t))*t)

return v[w]-np.dot(v[w],np.transpose(t))*t/np.dot(t,np.transpose(t))

def F(i,j):

import numpy as np

A=np.mat([[1,7],[1,2],[6,7]]);

aa=False;

a=[0]*3;

b=[0]*3;

for u in range(3):

a[u]=np.array_equal(np.mat([i,j]),np.array(A[u,:])) ;

b[u]=np.array_equal(np.mat([j,i]),np.array(A[u,:])) ;

aa=aa or a[u];

aa=aa or b[u];

if aa==True or i==j:

return [0,0,0]

else:

return [i,j,0]

#for i in range(8):

#    for j in range(8)    :

#        print(F(i,j));

def V(i,j,w):

a=i==j or j==w or w==i ;

b=F(i,j)==[0,0,0] or F(j,w)==[0,0,0] or F(i,w)==[0,0,0];

if a or b==True:

return [0,0,0]

else:

return [i,j,w]

#for i in range(8):

#    for j in range(8)    :

#        for w in range(8):

#          print(V(i,j,w),end=""),

#

def Q(w,i,j):

if w==0:

return [0,0,0]

elif w!=i and w!=j or F(i,j)==[0,0,0]:

return V(i,j,w)

else:

return [8,8,8]

#for i in range(8):

#    for j in range(8)    :

#        for w in range(8):

#          print(Q(w,i,j),end="");

print("---------------------------")

for i in range(8):

for j in range(8)    :

for k in range(8):

#          print(Q(f(k,i,j),i,j),end="");

if Q(f(k,i,j),i,j)==[8,8,8]:

print([i,j,k])

---------------------------

[0, 1, 1]

[0, 2, 2]

[0, 3, 3]

[0, 4, 4]

[0, 5, 5]

[0, 6, 6]

[0, 7, 7]

[1, 3, 3]

[1, 4, 4]

[1, 5, 5]

[1, 6, 6]

[2, 3, 3]

[2, 4, 4]

[2, 5, 5]

[2, 6, 6]

[2, 7, 7]

[3, 1, 1]

[3, 2, 2]

[3, 4, 4]

[3, 5, 5]

[3, 6, 6]

[3, 7, 7]

[4, 1, 1]

[4, 2, 2]

[4, 3, 3]

[4, 5, 5]

[4, 6, 6]

[4, 7, 7]

[5, 1, 1]

[5, 2, 2]

[5, 3, 3]

[5, 4, 4]

[5, 6, 6]

[5, 7, 7]

[6, 1, 1]

[6, 2, 2]

[6, 3, 3]

[6, 4, 4]

[6, 5, 5]

[7, 2, 2]

[7, 3, 3]

[7, 4, 4]

[7, 5, 5]

def q(w,i,j):

if Q(w,i,j)!=[8,8,8]:

return Q(w,i,j)

else:

return [Q(f(1,i,j),i,j),Q(f(2,i,j),i,j),Q(f(3,i,j),i,j),Q(f(4,i,j),i,j),Q(f(5,i,j),i,j),Q(f(6,i,j),i,j),Q(f(7,i,j),i,j)]

import numpy as np

v=np.mat([[0,0,0,0],[1,0,0,0],[0,1,0,0],[0,0,1,0],[0,0,0,1],[2,-1,-1,-1],[1,1,0,0],[-1,0,0,0]]);

A=np.mat([[1,7],[1,2],[6,7]]);

C=np.mat([[2,3,4,5],[3,4,5,6]]);

z=7;

def f(x,y,z):

if x==y or y==z:

return 0

else:

return x

def g(x,y,z):

if x==y:

return z

else:

return y

def h(x,y,z,k):

if x==y:

return [z,k]

elif x==z:

return [y,k]

else:

return [y,z]

def K(w,i,j):

t=v[g(w,i,j)]

#   print(t)

#   print(v[w])

#   print(np.dot(v[w],np.transpose(t))*t)

return v[w]-np.dot(v[w],np.transpose(t))*t/np.dot(t,np.transpose(t))

def F(i,j):

import numpy as np

A=np.mat([[1,7],[1,2],[6,7]]);

aa=False;

a=[0]*3;

b=[0]*3;

for u in range(3):

a[u]=np.array_equal(np.mat([i,j]),np.array(A[u,:])) ;

b[u]=np.array_equal(np.mat([j,i]),np.array(A[u,:])) ;

aa=aa or a[u];

aa=aa or b[u];

if aa==True or i==j:

return [0,0,0]

else:

return [i,j,0]

#for i in range(8):

#    for j in range(8)    :

#        print(F(i,j));

def V(i,j,w):

a=i==j or j==w or w==i ;

b=F(i,j)==[0,0,0] or F(j,w)==[0,0,0] or F(i,w)==[0,0,0];

if a or b==True:

return [0,0,0]

else:

return [i,j,w]

#for i in range(8):

#    for j in range(8)    :

#        for w in range(8):

#          print(V(i,j,w),end=""),

#

def Q(w,i,j):

if w==0:

return [0,0,0]

elif w!=i and w!=j or F(i,j)==[0,0,0]:

return V(i,j,w)

else:

return [8,8,8]

#for i in range(8):

#    for j in range(8)    :

#        for w in range(8):

#          print(Q(w,i,j),end="");

#print("---------------------------")

#for i in range(8):

#    for j in range(8)    :

#        for k in range(8):

##          print(Q(f(k,i,j),i,j),end="");

#          if Q(f(k,i,j),i,j)==[8,8,8]:

#              print([i,j,k])

def q(w,i,j):

if Q(w,i,j)!=[8,8,8]:

return Q(w,i,j)

else:

return [Q(f(1,i,j),i,j),Q(f(2,i,j),i,j),Q(f(3,i,j),i,j),Q(f(4,i,j),i,j),Q(f(5,i,j),i,j),Q(f(6,i,j),i,j),Q(f(7,i,j),i,j)]

q(5,7,5)

Out[16]: [[0, 0, 0], [7, 5, 2], [7, 5, 3], [7, 5, 4], [8, 8, 8], [0, 0, 0], [0, 0, 0]]

import numpy as np

v=np.mat([[0,0,0,0],[1,0,0,0],[0,1,0,0],[0,0,1,0],[0,0,0,1],[2,-1,-1,-1],[1,1,0,0],[-1,0,0,0]]);

A=np.mat([[1,7],[1,2],[6,7]]);

C=np.mat([[2,3,4,5],[3,4,5,6]]);

z=7;

def f(x,y,z):

if x==y or y==z:

return 0

else:

return x

def g(x,y,z):

if x==y:

return z

else:

return y

def h(x,y,z,k):

if x==y:

return [z,k]

elif x==z:

return [y,k]

else:

return [y,z]

def K(w,i,j):

t=v[g(w,i,j)]

#   print(t)

#   print(v[w])

#   print(np.dot(v[w],np.transpose(t))*t)

return v[w]-np.dot(v[w],np.transpose(t))*t/np.dot(t,np.transpose(t))

def F(i,j):

import numpy as np

A=np.mat([[1,7],[1,2],[6,7]]);

aa=False;

a=[0]*3;

b=[0]*3;

for u in range(3):

a[u]=np.array_equal(np.mat([i,j]),np.array(A[u,:])) ;

b[u]=np.array_equal(np.mat([j,i]),np.array(A[u,:])) ;

aa=aa or a[u];

aa=aa or b[u];

if aa==True or i==j:

return [0,0,0]

else:

return [i,j,0]

#for i in range(8):

#    for j in range(8)    :

#        print(F(i,j));

def V(i,j,w):

a=i==j or j==w or w==i ;

b=F(i,j)==[0,0,0] or F(j,w)==[0,0,0] or F(i,w)==[0,0,0];

if a or b==True:

return [0,0,0]

else:

return [i,j,w]

#for i in range(8):

#    for j in range(8)    :

#        for w in range(8):

#          print(V(i,j,w),end=""),

#

def Q(w,i,j):

if w==0:

return [0,0,0]

elif w!=i and w!=j or F(i,j)==[0,0,0]:

return V(i,j,w)

else:

return [8,8,8]

#for i in range(8):

#    for j in range(8)    :

#        for w in range(8):

#          print(Q(w,i,j),end="");

#print("---------------------------")

#for i in range(8):

#    for j in range(8)    :

#        for k in range(8):

##          print(Q(f(k,i,j),i,j),end="");

#          if Q(f(k,i,j),i,j)==[8,8,8]:

#              print([i,j,k])

def q(w,i,j):

if Q(w,i,j)!=[8,8,8]:

return Q(w,i,j)

else:

return [Q(f(1,i,j),i,j),Q(f(2,i,j),i,j),Q(f(3,i,j),i,j),Q(f(4,i,j),i,j),Q(f(5,i,j),i,j),Q(f(6,i,j),i,j),Q(f(7,i,j),i,j)]

q(5,7,5)

q(7,5,5)

q(5,5,7)

Out[17]: [[0, 0, 0], [5, 7, 2], [5, 7, 3], [5, 7, 4], [0, 0, 0], [0, 0, 0], [8, 8, 8]]

q(5,5,7)

Out[18]: [[0, 0, 0], [5, 7, 2], [5, 7, 3], [5, 7, 4], [0, 0, 0], [0, 0, 0], [8, 8, 8]]

q(7,5,5)

Out[19]: [0, 0, 0]

q(5,7,5)

Out[20]: [[0, 0, 0], [7, 5, 2], [7, 5, 3], [7, 5, 4], [8, 8, 8], [0, 0, 0], [0, 0, 0]]

weixin_39779928
关注
十二面体搜索算法python的实现
pannhao的博客
08-28 385
一种三维追踪算法
利用python在abaqus中画Voronoi多面体简单示例
最新发布
Niklas4088的博客
12-19 431
为Voronoi多面体各面上的顶点序号。或 CAE 新一个part 后删除几何。为Voronoi多面体顶点坐标。从rpy文件读取响应指令。从rpy文件读取响应指令。从rpy文件读取响应指令。
pypoman:多面体几何的Python模块
05-23
Python中的多面体操作 该库在实现了常见的操作，例如和。 有关详细信息，请参见完整的。 安装 通过以下方法确保您具有所有系统范围的依赖项： sudo apt-get install cython libglpk-dev python python-dev python-pip python-scipy 然后，安装模块本身： pip install pypoman 例子 顶点枚举 我们可以通过A * x <= b来计算以半空间表示形式描述的多边形的顶点列表： from numpy import array from pypoman import compute_polytope_vertices A = array ([ [ - 1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ], [ 0 ,
Python绘制正二十面体
微小冷的学习笔记
12-10 3438
正20面体的12个顶点刚好可以分为三组，每一组都是一个符合黄金分割比例的长方形，而且这三个长方形是互相正交的。
python三维曲面合并_绘制多面体的三维曲面
weixin_35993037的博客
01-29 839
我尝试使用Python和Matplotlib来渲染一个多面体的3D曲面，由但是我的代码(如下所示)似乎没有正确地绘制它。应该怎么做呢？在失败的尝试：%matplotlib inlineimport numpy as npimport matplotlib.pyplot as pltfrom mpl_toolkits.mplot3d import Axes3Dfrom matplotlib impo...
pplpy:Parma多面体库（PPL）的Python包装器
05-15
PPL Python包装器（移至gitlab） 该项目已移至gitlab。 请考虑参观
newmesh_3D_混凝土_python_python_骨料_混凝土骨料_
10-02
本主题聚焦于“三维voronoi多面体混凝土土骨料网格划分”，这是一个利用Python编程语言进行的复杂几何计算过程。Voronoi图是一种在几何学中广泛使用的工具，用于构建一个区域内最近点的集合。在混凝土骨料建模中，...
3D-voronoi_Voronoi3D_dirichlet图_python_泰森多边形_voronoi_
09-29
在3D-Voronoi图中，这些多边形被扩展为多面体，每个多面体的边界由站点之间的中垂线决定，这些中垂线在三维空间中形成面的边界。 在Python中实现3D-Voronoi图，我们可以利用如`scipy.spatial.Voronoi`这样的库。`...
vor3d_混凝土_随机骨料_骨料_3Dvoronoi多面体生成_
10-03
本教程将详细讲解如何利用Python编程语言中的vor3d库来生成混凝土中的随机骨料，即3D Voronoi多面体。 首先，让我们了解一下Voronoi图的基本概念。Voronoi图，又称Dirichlet分区，是一种几何分割方法，它将空间分割...
matlab开发-多面体和多面体
08-24
matlab开发-多面体和多面体。在多面体和多面体上执行操作的一套程序
Python：利用python编程实现三维图像绘制展示(六面体旋转、三维球柱状体、下雪场景等)
近期请国内外头部出版社可尽快私信博主！——心比天高，仗剑走天涯，保持热爱，奔赴向梦想！低调，谦虚，自律，反思，成长，还算是比较正能量的博主，公益免费传播……内心特别想在AI界做出一些可以推进历史进程影响力的东西(兴趣使然，有点小情怀，也有点使命感呀)…
11-08 5683
Python：利用python编程实现三维图像绘制展示(六面体旋转、三维球柱状体、下雪场景等) 目录 利用python编程实现三维图像绘制展示(六面体旋转、三维球柱状体、下雪场景等) 1、实现六面体旋转 2、三维球柱状体 3、下雪场景 利用python编程实现三维图像绘制展示(六面体旋转、三维球柱状体、下雪场景等) 1、实现六面体旋转 2、三维球柱状体 3、下雪场景 ...
【python】用plotly绘制正二十面体
微小冷的学习笔记
02-23 1934
通过plotly绘制正二十面体及其顶点和棱，通过这个案例掌握三维散点图、三维曲线图以及三角剖分图。
在Python里应用Openscad实现3D建模(修改简化版)之3D螺纹建模初探1
伪民科小白随意折腾随心记
06-11 1376
在Python里应用Openscad实现3D建模(修改简化版)-2 –SolidPython学习笔记2 –form https://github.com/SolidCode/SolidPython 参考Parkinbotshortcuts.scad做了部分修改使得代码更为简洁易读 以下是对solidpython的readme的学习笔记(其中的部分code进行了简化)： –3D螺纹建模初...
python matplatlib 画边长为1的正四面体和正四棱锥
小米mm修仙路
07-27 3309
import matplotlib.pyplot as plt from mpl_toolkits.mplot3d.art3d import Poly3DCollection, Line3DCollection fig = plt.figure() ax = fig.gca(projection='3d') ax.set_xlabel('X') ax.set_xlim3d(0, 2)...
python 3D绘制立体几何
majunfu
05-13 9427
直接复制就能用，写的简单，请勿吐槽import numpy as np import mpl_toolkits.mplot3d import matplotlib.pyplot as pltx=[0,3,0,3,1.5] y=[0,0,3,3,1.5] z=[0,0.8,1.5,2.3,3]ax=plt.subplot(111,projection='3d') for i in range(len(
OpenGL（二）图形绘制之平面多面体的绘制
爱拼才会赢
07-07 5451
通过绘制一个三菱锥初步了解绘制平面多面体。 #include #include #include #include #include //窗口的大小 GLsizei windowWidth; GLsizei windowHeight; //旋转角度参数 static GLfloat xRot = 0.0f; static GLfloat yRot = 0.0f; //确定多边形的绕法的方
python三角斜边_多面体的Delaunay三角化(Python)
weixin_39996134的博客
12-01 868
我试图在python中获得Delaunay Triangulation的多面体,以便我可以计算质心.我看到scipy.spatial中有一个Delaunay函数,它在n维中起作用.问题是文档显示了二维使用,并没有告诉我如何处理更高的尺寸.能够将这个对象分解成数组可能会为我解决这个问题,但我不知道该怎么做.我遇到的问题是我不知道如何在输出对象时验证这是否正常工作.我在Google上找不到关于如何绘制...
ABAQUS中Python脚本实现Voronoi多面体制作教程
ABAQUS作为一个强大的有限元分析软件，虽然主要以几何建模和求解力学问题为主，但通过Python脚本功能，可以实现Voronoi多面体的自动生成。 通常，人们会采用外部工具如MATLAB或其他专门软件来创建Voronoi结构，然后...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值