=0:print(count,'个不满足模型分布的约束条件,将使用福莱特法进行修正。')returncountifcount==0:print("迭代结束")returncount#求qij 求Q转换为oddeffind_Q(i):q=OD[i].iloc[:-1,:-1].values.flatten()foi=Foi[i].repeat(3)fdj=Fdj[i].tolist()*3li=Loi.repeat(3)lj=Ldj.tolist()*3Q=q*...