Question

试图求解某些通过循环产生的初始条件不同的方程。

由于“ while”循环，我想获取多个文件，每个文件代表不同的初始条件，但是每个文件都是相同的，数据不会改变。无法看到我哪里出了问题。

from scipy.integrate import odeint

P = 0      # birth rate
d = 0.0001  # natural death percent (per day)
B = 0.0095  # transmission percent  (per day)
G = 0.0001  # resurect percent (per day)
A = 0.0001  # destroy percent  (per day)

# solve the system dy/dt = f(y, t)
def f(y, t):
     Si = y[0]
     Zi = y[1]
     Ri = y[2]
     # the model equations (see Munz et al. 2009)
     f0 = P - B*Si*Zi - d*Si
     f1 = B*Si*Zi + G*Ri - A*Si*Zi
     f2 = d*Si + A*Si*Zi - G*Ri
     return [f0, f1, f2]

Si0 = 50
Zi0 = 20
Ri0= 20

for i in range(20):   
    while Zi0 < 105:
        Zi0 += 5
        y0=[Si0,Zi0,Ri0]
        print(Zi0)


        y = odeint(f, y0, [k*500 for k in range(40)], mxstep=1000)


    print('\n'.join(','.join('%+13.5f'%x for x in row) for row in \
    y),file=open("test_"+str(i)+".txt", "a"))

Answer 1

问题可能是您没有重置Zi0。因此，第一次在for循环中执行while循环，直到Zi0为105。然后第二次for循环开始且Zi0已经为105，因此它将跳过while循环并打印，并继续跳过和打印其余的for循环。

不完全确定要更改每个for循环的内容，因此无法帮助您更多...抱歉。

循环写入多个文件会得到相同的文件

1 个答案: