Question

我有重复动画的问题。我想使用matplotlib以便可视化矩阵上的每个模拟。这是模拟代码：

import numpy as np
from matplotlib import pyplot as plt

s = np.array([[1,1,1], [1,10,1], [1,1,1]], dtype=np.int8) #Matrix of wages
e = np.zeros((19,), dtype=np.int8) # Vector of rules
e[3]=1
e[12]=1
e[13]=1

mama = np.array([[1, 0, 1], [0, 1, 0], [0, 0, 1]], dtype=np.int8) #The matrix to be tested

def simulation(ma): #Simulation on a given matrix

    n, m = ma.shape
    p = np.zeros((n+2, m+2), dtype=np.int8) #Creates an extended matrix, avoiding conflicts at the edges of the initial matrix. Here I construct a torus
    p[1:-1, 1:-1] = ma #middle
    p[0, 1:-1] = ma[n-1] #the first row of p, the last of ma
    p[-1, 1:-1] = ma[0] #the last row of p, the first of ma
    p[1:-1, 0] = ma[0:, -1] #left col p, right of ma
    p[1:-1, -1] = ma[0:, 0] #right col of p, left of ma
    p[-1, 0] = ma[0, -1] #left bottom corner
    p[-1, -1] = ma[0, 0] #right bottom corner
    p[0, 0] = ma[-1, -1] #left upper corner
    p[0, -1] = ma[-1, 0] #right upper corner

    new = np.zeros(ma.shape, dtype=np.int8)

    v, c = p.shape #verses and columns

    for i in range(1, v):
        for j in range(1, c):
            if p[i-1:i+2, j-1:j+2].shape == (3, 3):
                new[i-1, j-1] = e[np.sum(p[i-1:i+2,j-1:j+2]*s)]
    return new

但是，我想对指定的重复次数运行模拟并可视化模拟的每个步骤，因此我尝试了以下代码：

def rep(fun, mac, ti): #function, matrix, repetitions (time)
    if ti == 1:
        plt.imshow(fun(mac))
        plt.title("Celllar automaton")
        plt.show()

    else:
        f = fun(rep(fun, mac, ti-1))
        plt.imshow(f)
        plt.title("Cellular automaton")
        plt.show()

我得到一个错误：

n, m = ma.shape
AttributeError: 'NoneType' object has no attribute 'shape'

请，你能帮我吗？我真的厌倦了无法可视化我的工作。

可添加：我已将rep替换为：

def shoow(fig):
    plt.imshow(fig)
    plt.title("Cellular automaton")
    plt.show()

def repet(fun, mac, ti):
    c1 = mac
    for i in range(ti):
        f = fun(c1)
        shoow(f)
        c1 = f

但是，它每次都会创建一个新图形。如何获得连续的模拟？

元胞自动机-重复模拟imshow（）

0 个答案: