标签:bm 模型 mi ampl1 约束 range sumnum pow 重力
c = [[0, 0, 0, 0], [0, 2, 4, 4], [0, 4, 1, 2], [0, 4, 2, 2]] o = [0, 8, 14, 10] d = [0, 8, 14, 10] am = [0, 1, 1, 1] ampl1 = [0, 0, 0, 0] bm = [0, 0, 0, 0] bmpl1 = [0, 0, 0, 0] mi = 0.01 e = 0.05 m = 0 for j in range(1, 4): sumnum = 0 for i in range(1, 4): sumnum += am[i] * o[i] * pow(c[i][j], -mi) bm[j] = 1 / sumnum for i in range(1, 4): sumnum = 0 for j in range(1, 4): sumnum += bm[j] * d[j] * pow(c[i][j], -mi) ampl1[i] = 1 / sumnum for j in range(1, 4): sumnum = 0 for i in range(1, 4): sumnum += ampl1[i] * o[i] * pow(c[i][j], -mi) bmpl1[j] = 1 / sumnum def ifif(a, b, c, d): f = [0, 0, 0, 0, 0] for i in range(1, 4): if (a[i] / b[i] > 1 - e and a[i] / b[i] < 1 - e) and (c[i] / d[i] > 1 - e and c[i] / d[i] < 1 + e): f[i] = True if (f[1] and f[2]) and f[3]: return True else: return False for k in range(50): if not ifif(ampl1, am, bmpl1, bm): m += 1 mi += 0.05 for j in range(1, 4): sumnum = 0 for i in range(1, 4): sumnum += am[i] * o[i] * pow(c[j][i], -mi) bm[j] = 1 / sumnum print(bm[j]) for i in range(1, 4): sumnum = 0 for j in range(1, 4): sumnum += bm[j] * d[j] * pow(c[i][j], -mi) ampl1[i] = 1 / sumnum for j in range(1, 4): sumnum = 0 for i in range(1, 4): sumnum += ampl1[i] * o[i] * pow(c[j][i], -mi) bmpl1[j] = 1 / sumnum print(am, ampl1, bm, bmpl1, mi)
2020-05-26
标签:bm,模型,mi,ampl1,约束,range,sumnum,pow,重力 来源: https://www.cnblogs.com/hany-postq473111315/p/12968621.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。