使用while循环重复循环的步骤,以达到正确的答案(Using while loop to repe

2019-10-28 14:24发布

我有以下的代码,我只写它的最后一部分。 出于某些原因,我们使用random.normal值。 因为很明显的func功能有三个变量。

我的问题是:计算的步骤是101和正好我们有100个输出的每个参数值。 我该如何告诉Python的,如果你看到Wo[i]>0Omo[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])

Answer 1:

您可以使用一般的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


Answer 2:

插入while循环重复你所需要的计算。

L > R WoOmo与更新OmnWn值。 因此,如果这个值是Omn<0Wn>0 ,你需要重新calulate他们。 当L <= RWoOmo从先前迭代计算。 由于前一次迭代已经是Wo[i]<=0Omo[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])


文章来源: Using while loop to repeat a step of loop to reach the proper answer