我有以下的代码,我只写它的最后一部分。 出于某些原因,我们使用random.normal
值。 因为很明显的func
功能有三个变量。
我的问题是:计算的步骤是101
和正好我们有100
个输出的每个参数值。 我该如何告诉Python的,如果你看到Wo[i]>0
和Omo[i]<0
重新计算步骤(这将使用另一种随机数),直到我们有正确的答案(Wo[i]<0 and Omo[i]>0
。
我们必须这样做,因为 ,如果我们打印输出,然后删除那些不符合我们的条件值,我们有例如,60个值,而不是100.我知道我们可以使用while loop
,但我不知道怎么办。
N=101
for i in range (1,N):
R=np.random.uniform(0,1)
Omn[i] = Omo[i-1] + 0.05 * np.random.normal()
Wn[i] = Wo[i-1] + 0.05 * np.random.normal()
Mn[i] = Mo[i-1] + 0.1 * np.random.normal()
L = exp(-0.5 * ( func(Omn[i], Wn[i], Mn[i] ) - func( Omo[i-1], Wo[i-1], Mo[i-1] )))
if L>R:
Wo[i]=Wn[i]
Mo[i]=Mn[i]
Omo[i]=Omn[i]
else:
Wo[i]=Wo[i-1]
Mo[i]=Mo[i-1]
Omo[i]=Omo[i-1]
print(Wo[i],Mo[i],Omo[i])
您可以使用一般的while循环和计算之后,例如检查条件:
for i in range (1,N):
R=np.random.uniform(0,1)
while True:
Omn[i] = Omo[i-1] + 0.05 * np.random.normal()
Wn[i] = Wo[i-1] + 0.05 * np.random.normal()
Mn[i] = Mo[i-1] + 0.1 * np.random.normal()
if Wo[i] <= 0 and Omo[i] >= 0:
break
# rest of code
插入while
循环重复你所需要的计算。
当L > R
Wo
和Omo
与更新Omn
和Wn
值。 因此,如果这个值是Omn<0
或Wn>0
,你需要重新calulate他们。 当L <= R
, Wo
和Omo
从先前迭代计算。 由于前一次迭代已经是Wo[i]<=0
和Omo[i]>=0
则不需要重复这些计算。
因此,如下面的代码,你只需要重新计算OMN和Wn的变量:
Omo[0] = 0.24
Wo[0] = -0.2
Mo[0] = 1.0
N = 101
for i in range (1,N):
is_valid = False
if __debug__:
print "Calculating position " + str(i)
while (not is_valid):
Omn[i] = Omo[i-1] + 0.05 * np.random.normal()
Wn[i] = Wo[i-1] + 0.05 * np.random.normal()
Mn[i] = Mo[i-1] + 0.1 * np.random.normal()
if __debug__:
print "- is_valid iteration values: " + str(Wn[i]) + " " + str(Omn[i])
print "- is_valid previous values: " + str(Wo[i-1]) + " " + str(Omo[i-1])
is_valid = Omn[i] >= 0 and Wn[i] <= 0
R = np.random.uniform(0,1)
L = exp(-0.5 * ( func(Omn[i], Wn[i], Mn[i] ) - func( Omo[i-1], Wo[i-1], Mo[i-1] )))
if L > R:
Omo[i] = Omn[i]
Wo[i] = Wn[i]
Mo[i] = Mn[i]
else:
Omo[i] = Omo[i-1]
Wo[i] = Wo[i-1]
Mo[i] = Mo[i-1]
print(Wo[i], Mo[i], Omo[i])