<p>好的,所以我基本上修复了你的代码,因为,没有冒犯,但它之前有点混乱(见注释)。你知道吗</p>
<pre><code>import random
import matplotlib.pyplot as plt
fig, axs = plt.subplots(nrows=2, sharex=True)
L = 100 # Number of lattice sites
alpha = 0.2 # Rate of entry
beta = 0.4 # Rate of exit
n_time = 200000 # Simulation steps
record_each = 2000
state = [0]*(L + 1)
record = [] # store a record of the total number of states occupied
for itr in range(n_time):
rand_int = random.randint(0, L)
if rand_int == 0:
if random.random() < alpha:
state[1] = 1
elif rand_int == L:
if random.random() < beta:
state[L] = 0
elif state[rand_int] == 1 and state[rand_int + 1] == 0:
state[rand_int + 1] = 1
state[rand_int] = 0
if itr % record_each == 0:
yaxis = [i for i in range(L) if state[i + 1] == 1]
axs[1].plot([itr]*len(yaxis), yaxis, 'r.')
record.append(sum(state)) # add the total number of states to the record
axs[0].plot(range(0, n_time, record_each), record) # plot the record
axs[1].set_xlabel('Number of steps')
axs[1].set_ylabel('System configuration')
axs[0].set_ylabel('Number of states occupied')
plt.show()
</code></pre>
<p>这个输出</p>
<p><a href="https://i.stack.imgur.com/G40rt.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/G40rt.png" alt="enter image description here"/></a></p>