我们可以在python中实现4个变量的K-Map算法,如下图。该函数接受 SOP(乘积之和)形式的布尔函数和变量名称,并返回简化的简化表示。基本上,您需要创建包含 8、4、2 等 2 次幂的总项的矩形组,并尝试在一组中覆盖尽可能多的元素(我们需要覆盖所有元素)。
例如,函数可以表示为 F(w,x,y,z) = xy' + x'z' + wxz + wx'y SOP 形式为 f(w,x,y,z)= ∑(0,2,4,5,8,10,11,12,13,15),由下表可知:
从下一段代码sn-p的输出可以看出,程序输出简化形式x¬y + ¬x¬z + wyz,其中布尔变量x的取反在代码中表示为¬x。
from collections import defaultdict
from itertools import permutations, product
def kv_map(sop, vars):
sop = set(sop)
not_covered = sop.copy()
sop_covered = set([])
mts = [] # minterms
# check for minterms with 1 variable
all_3 = [''.join(x) for x in product('01', repeat=3)]
for i in range(4):
for v_i in [0,1]:
if len(not_covered) == 0: continue
mt = ('' if v_i else '¬') + vars[i]
s = [x[:i]+str(v_i)+x[i:] for x in all_3]
sop1 = set(map(lambda x: int(x,2), s))
if len(sop1 & sop) == 8 and len(sop_covered & sop1) < 8: # if not already covered
mts.append(mt)
sop_covered |= sop1
not_covered = not_covered - sop1
if len(not_covered) == 0:
return mts
# check for minterms with 2 variables
all_2 = [''.join(x) for x in product('01', repeat=2)]
for i in range(4):
for j in range(i+1, 4):
for v_i in [0,1]:
for v_j in [0,1]:
if len(not_covered) == 0: continue
mt = ('' if v_i else '¬') + vars[i] + ('' if v_j else '¬') + vars[j]
s = [x[:i]+str(v_i)+x[i:] for x in all_2]
s = [x[:j]+str(v_j)+x[j:] for x in s]
sop1 = set(map(lambda x: int(x,2), s))
if len(sop1 & sop) == 4 and len(sop_covered & sop1) < 4: # if not already covered
mts.append(mt)
sop_covered |= sop1
not_covered = not_covered - sop1
if len(not_covered) == 0:
return mts
# check for minterms with 3 variables similarly (code omitted)
# ... ... ...
return mts
mts = kv_map([0,2,4,5,8,10,11,12,13,15], ['w', 'x', 'y', 'z'])
mts
# ['x¬y', '¬x¬z', 'wyz']
下面的动画展示了上面的代码如何(贪婪地)简化了以 SOP 形式给出的布尔函数(基本目标是用最少的 2 次方块覆盖所有 1)。由于算法是贪心的,它可能会卡在某个局部最小值,我们需要小心。