谓词逻辑的归结反演(人工智能python)

一、 前言

        计算机专业人工智能课程的实验,笔者耗时接近15个小时完成,特写一篇博客记录艰辛历程,如有错误、建议和优化欢迎指出

二、 算法原理

(1) 归结反演

        将目标子句的反~a加入前提子句集KB得到初始子句集KB`,选择两个子句作为亲本子句,将它们的互补文字进行合一操作,然后将亲本子句中的参数用合一进行替换,最后合并两个亲本子句并删除掉互补文字,得到一个新的归结子句。反复进行上述操作,直到得到一个空的归结子句,说明KB`不可满足,即可得KB|=a

(2) 支持集策略

        亲本子句至少有一个来自目标子句及其归结式,即通过仅选择与目标相关的子句进行归结来减少搜索空间,从而最大程度地减少无用的归结

(3) 回溯化简

        从空子句出发,一步一步往回推,找到与归结出空子句相关的有用归结式,忽略掉无用的归结式

 

(4) 算法描述

        先用支持集策略和归结原理得到全部归结式,再从空子句开始回溯化简把用不到的归结式删除并更新有用归结式的编号

(5) 伪代码:

将初始子句集加入到归结式列表

令支持集只有目标子句

循环:

        令归结集为空

        遍历子句集每个子句1和每个子句2

                if子句2在支持集

                        遍历子句1里的文字1和子句2里的文字2

                                if文字1和文字2互补

                                        得到索引和最一般合一

                                                if最一般合一存在

                                                        置换后归结得到新子句

                                                        if新子句不重复

                                                                得到归结式

                                                                更新归结集、更新归结式列表

                                                        if新子句是空子句

                                                                退出外循环

        更新子句集、更新支持集

将得到空子句的所有有关归结式汇总得到有用归结式列表

新归结式列表 = 初始子句列表 + 有用归结式列表

遍历新归结列表的每个归结式

        重新编号

        重新写归结式

输出最终的归结列表

 

三、算法实现

(1)数据结构

  • 谓词:首字母大写,不含∃,∀量词符号
  • 参数:用逗号“,”间隔,不加空格,在()内
  • 文字:用字符串存储,符号‘¬’用‘~’代替
  • 子句:用元组存储
  • 子句集:用列表存储
  • 合一置换:用字典存储
  • 常量:单词a, b, c等指定小写字母表示
  • 变量:u、v、w、x、y、z等指定小写字母表示
  • 归结式:用字符串存储,格式为:步骤编号 R[子句编号1,子句编号2]{变量=常量} = 子句 (如果一个字句包含多个公式就用编号 a,b,c 区分,否则不用加编号区分)
  • 归结过程:用列表存储

(2)变量含义

  • clauseset:子句集(列表)
  • clause:子句(元组)
  • literal:文字(字符串)
  • para:参数(列表)
  • unification:合一(字典)
  • index:编号(字符串)
  • sequence:归结式(字符串)
  • result:归结过程(列表)
  • number:子句的索引(整型)
  • useful_process:有用的归结式(列表)
  • KB:初始子句集(列表)

(3)自定义函数

1. Complementary():判断两个文字是否互补

        输入:两个文字

        输出:若存在则返回两个文字索引,若不存在则返回None

#判断两个文字是否互补
def Complementary(literal1,literal2):
    #谓词:从文字开头到'('的子字符串
    end1 = literal1.find('(')
    end2 = literal2.find('(')
    #互补:文字间只有~的区别
    if literal2[0] == '~' and literal2[1:end2] == literal1[:end1]:
        return True
    if literal1[0] == '~' and literal1[1:end1] == literal2[:end2]:
        return True 
    return False

2. Index():得到归结所用的编号

        输入:子句索引、文字索引、子句长度

        输出:归结式的子句编号

#得到题目要求的归结编号
def Index(literal_index,clause_index,length):
    if length == 1: #如果子句只有一个元素,则文字索引不再需要
        index = str(clause_index+1)
    else:           #否则将文字索引变为字母
        index = str(clause_index+1) + chr(ord('a')+literal_index)
    return index

3. Parameter():得到文字的参数列表

        输入:文字

        输出:参数列表

#得到文字的参数列表
def Parameter(literal):
    #参数是在()内且以,分隔的多个子字符串
    start = literal.find('(')
    end = literal.find(')')
    para = literal[start+1:end].split(',')
    return para

4. Unify():得到最一般合一(不完全正确,注意甄别!!)

        输入:两个文字的参数列表

        输出:合一字典

#得到合一
def Unify(para1,para2):
    variable = {'u','v','w','x','y','z','uu','vv','ww','xx','yy','zz'}  #变量集
    unification = {}
    if para1 == para2:  #若一开始就相等,直接返回空合一
        return unification 
    while True:
        for i,j in zip(para1,para2):
            if i in variable and j in variable: #若都是变量,不可合一
                return None
            elif i not in variable and j not in variable and i != j:    #若都是常量且不相等,不可合一
                return None
            elif i in variable and j not in variable:
                unification[i] = j
                break
            elif j in variable and i not in variable:
                unification[j] = i
                break
        #合一复合
        para1 = [unification[item] if item in unification else item for item in para1]
        para2 = [unification[item] if item in unification else item for item in para2]
        if para1 == para2:
            break
    return unification

5. Substitute():对子句置换

        输入:子句、合一

        输出:置换后的新子句

#置换旧子句得到新子句
def Substitute(unification,clause):
    if unification == {}:   #字典为空,说明无需置换,直接返回原子句
        return clause
    newclause = []
    for literal in clause:
        start = literal.find('(')
        end = literal.find(')')
        para = literal[start + 1:end].split(',')   #得到文字的参数列表
        para = [unification[item] if item in unification else item for item in para]
        newliteral = literal[:start + 1] + ','.join(para) + literal[end:]  #合并成新子句
        newclause.append(newliteral)
    return tuple(newclause) #返回新子句,注意是元组格式

6. Resolve():归结

        输入:两个子句、两个互补文字索引

        输出:归结后的新子句

#归结:消除互补项
def Resolve(clause1,clause2,literal_index1,literal_index2):
    newclause = list(clause1) + list(clause2)
    newclause.remove(clause1[literal_index1])
    newclause.remove(clause2[literal_index2])
    newclause = list(OrderedDict.fromkeys(newclause))   #消除相同的子句
    return tuple(newclause)

7. Sequence():得到归结式

        输入:归结子句、子句编号、合一

        输出:归结式

#得到归结式
def Sequence(newclause,unification,index1,index2):
    string = ''
    if unification == {}:    #如果字典为空,说明不需要输出合一
        string += 'R[' + index1 + ',' + index2 + '] = '
    else:
        string += 'R[' + index1 + ',' + index2 + ']{'
        for key,value in unification.items():
            string += key + '=' + value + ','
        string = string[:-1]
        string += '} = '
    string += str(newclause)
    return string

8. Refutation():得到归结过程

        输入:初始子句集

        输出:归结过程

#反演:支持集策略
def Refutation(KB):
    clauseset = list(KB)            #拷贝一份,防止更改原初始子句集
    supportset = [KB[-1]]           #支持集,默认KB最后一个元素是目标子句的否定
    result = ['归结顺序:'] + KB      #将0位置补充元素,确保编号和列表索引对应
    while True:
        new_clauseset = []  #此级得到的全部的新的子句
        clause_index1 = 0
        for clause1 in clauseset:
            if clause1 in KB:
                clause_index2 = 0
                for clause2 in clauseset:
                    if clause1 != clause2 and clause2 in supportset: #其中一个亲本子句必须来自支持集
                        literal_index1 = 0
                        for literal1 in clause1:
                            literal_index2 = 0
                            for literal2 in clause2:
                                #判断是否互补
                                if Complementary(literal1,literal2):
                                    #得到索引
                                    index1 = Index(literal_index1,clause_index1,len(clause1))
                                    index2 = Index(literal_index2,clause_index2,len(clause2))
                                    #得到互补文字
                                    literal1 = clause1[literal_index1]
                                    literal2 = clause2[literal_index2]
                                    #得到参数列表
                                    para1 = Parameter(literal1)
                                    para2 = Parameter(literal2)
                                    #得到合一置换
                                    unification = Unify(para1,para2)
                                    #合一置换不存在则退出
                                    if unification == None:
                                        break
                                    #得到子句
                                    newclause1 = Substitute(unification,clause1)    
                                    newclause2 = Substitute(unification,clause2)  
                                    #得到归结子句
                                    newclause = Resolve(newclause1,newclause2,literal_index1,literal_index2)
                                    #归结子句存在于原子句集则退出
                                    if any([set(newclause)==(set(item)) for item in clauseset]):
                                        break
                                    #归结子句存在于新子句集则退出
                                    if any([set(newclause)==(set(item)) for item in new_clauseset]):
                                        break
                                    #得到新归结式
                                    sequence = Sequence(newclause,unification,index1,index2)
                                    #加入结果列表和新子句列表
                                    result.append(sequence)
                                    new_clauseset.append(newclause)
                                    #当且仅当新子句为空,则可退出函数,返回结果列表
                                    if newclause == ():
                                        return result
                                literal_index2 += 1
                            literal_index1 += 1
                    clause_index2 += 1
            clause_index1 += 1
        clauseset += new_clauseset      #更新子句列表
        supportset += new_clauseset     #更新支持集

9. Number():得到归结式的子句索引

        输入:归结式

        输出:两个子句索引

#得到归结式的子句索引
def Number(clause):
    start = clause.find('[')
    end = clause.find(']')
    number = clause[start+1:end].split(',')
    #将文字索引去掉
    num1 = int(''.join(item for item in number[0] if not item.isalpha()))
    num2 = int(''.join(item for item in number[1] if not item.isalpha()))
    return num1,num2

10. Renumber():得到新归结式的子句索引

        输入:原归结过程、归结式、待更新索引、初始子句集大小

        输出:更新后的索引

#得到新归结式的子句索引
def Renumber(num,result,useful_process,size):
    if num <= size: #如果是初始子句集的,直接返回
        return num
    #找到亲本子句
    sequence = result[num]
    begin = sequence.find('(')
    aim_clause = sequence[begin:]
    #找到亲本子句在化简子句集的编号
    for i in range(size+1,len(useful_process)):
        begin = useful_process[i].find('(')
        if useful_process[i][begin:] == aim_clause:
            return i 

11. Resequence():更新归结式

        输入:原归结式、原索引、新索引

        输出:新归结式

#更新归结式
def Resequence(sequence,num1,num2,newnum1,newnum2):
    # 第一次替换:替换第一个编号
    start = sequence.find(num1)
    end = start + len(num1)
    sequence = sequence[:start] + newnum1 + sequence[end:]
    # 第二次替换:替换第二个编号
    end = start + len(newnum1)
    start = sequence.find(num2, end)
    end = start + len(num2)
    sequence = sequence[:start] + newnum2 + sequence[end:]
    return sequence

12. Simplify():化简归结过程

        输入:原归结过程、初始子句集大小

        输出:新归结过程

#化简归结过程
def Simplify(result,size):
    base_process = result[0:size+1] #初始子句集
    useful_process = []             #有用子句集
    number = [len(result)-1]        #用作队列,先将空子句的索引入列
    while number != []:
        number0 = number.pop(0)                 #提取队列首元素,即有用子句的索引
        useful_process.append(result[number0])  #将有用子句加入到有用子句集            
        num1,num2 = Number(result[number0])     #得有用子句用到的亲本子句索引
        #如果是初始子句集就无需加入
        if num1 > size:
            number.append(num1)
        if num2 > size:
            number.append(num2)
    #得到新的归结过程
    useful_process.reverse()
    useful_process = base_process + useful_process
    #将归结过程重新编号
    for i in range(size+1,len(useful_process)): 
        num1,num2 = Number(useful_process[i])
        newnum1 = str(Renumber(num1,result,useful_process,size))
        newnum2 = str(Renumber(num2,result,useful_process,size))
        useful_process[i] = Resequence(useful_process[i],str(num1),str(num2),newnum1,newnum2)
    return useful_process

13. Print():输出归结过程

        输入:归结过程

        输出:编号后的归结过程

#打印结果
def Print(result):
    print(result[0])
    for i in range(1,len(result)):
        print(i,result[i])

14. ResolutionFOL():归结反演(汇总)

        输入:初始子句集

        输出:归结反演

#归结反演
def ResolutionFOL(KB):
    result = Refutation(KB)
    new_result = Simplify(result,len(KB))
    Print(new_result)

 

四、完整代码

from collections import OrderedDict

#判断两个文字是否互补
def Complementary(literal1,literal2):
    #谓词:从文字开头到'('的子字符串
    end1 = literal1.find('(')
    end2 = literal2.find('(')
    #互补:文字间只有~的区别
    if literal2[0] == '~' and literal2[1:end2] == literal1[:end1]:
        return True
    if literal1[0] == '~' and literal1[1:end1] == literal2[:end2]:
        return True 
    return False

#得到题目要求的归结编号
def Index(literal_index,clause_index,length):
    if length == 1: #如果子句只有一个元素,则文字索引不再需要
        index = str(clause_index+1)
    else:           #否则将文字索引变为字母
        index = str(clause_index+1) + chr(ord('a')+literal_index)
    return index

#得到文字的参数列表
def Parameter(literal):
    #参数是在()内且以,分隔的多个子字符串
    start = literal.find('(')
    end = literal.find(')')
    para = literal[start+1:end].split(',')
    return para

#得到合一
def Unify(para1,para2):
    variable = {'u','v','w','x','y','z','uu','vv','ww','xx','yy','zz'}  #变量集
    unification = {}
    if para1 == para2:  #若一开始就相等,直接返回空合一
        return unification 
    while True:
        for i,j in zip(para1,para2):
            if i in variable and j in variable: #若都是变量,不可合一
                return None
            elif i not in variable and j not in variable and i != j:    #若都是常量且不相等,不可合一
                return None
            elif i in variable and j not in variable:
                unification[i] = j
                break
            elif j in variable and i not in variable:
                unification[j] = i
                break
        #合一复合
        para1 = [unification[item] if item in unification else item for item in para1]
        para2 = [unification[item] if item in unification else item for item in para2]
        if para1 == para2:
            break
    return unification


#置换旧子句得到新子句
def Substitute(unification,clause):
    if unification == {}:   #字典为空,说明无需置换,直接返回原子句
        return clause
    newclause = []
    for literal in clause:
        start = literal.find('(')
        end = literal.find(')')
        para = literal[start + 1:end].split(',')   #得到文字的参数列表
        para = [unification[item] if item in unification else item for item in para]
        newliteral = literal[:start + 1] + ','.join(para) + literal[end:]  #合并成新子句
        newclause.append(newliteral)
    return tuple(newclause) #返回新子句,注意是元组格式

#归结:消除互补项
def Resolve(clause1,clause2,literal_index1,literal_index2):
    newclause = list(clause1) + list(clause2)
    newclause.remove(clause1[literal_index1])
    newclause.remove(clause2[literal_index2])
    newclause = list(OrderedDict.fromkeys(newclause))   #消除相同的子句
    return tuple(newclause)

#得到归结式
def Sequence(newclause,unification,index1,index2):
    string = ''
    if unification == {}:    #如果字典为空,说明不需要输出合一
        string += 'R[' + index1 + ',' + index2 + '] = '
    else:
        string += 'R[' + index1 + ',' + index2 + ']{'
        for key,value in unification.items():
            string += key + '=' + value + ','
        string = string[:-1]
        string += '} = '
    string += str(newclause)
    return string

#反演:支持集策略
def Refutation(KB):
    clauseset = list(KB)            #拷贝一份,防止更改原初始子句集
    supportset = [KB[-1]]           #支持集,默认KB最后一个元素是目标子句的否定
    result = ['归结顺序:'] + KB      #将0位置补充元素,确保编号和列表索引对应
    while True:
        new_clauseset = []  #此级得到的全部的新的子句
        clause_index1 = 0
        for clause1 in clauseset:
            if clause1 in KB:
                clause_index2 = 0
                for clause2 in clauseset:
                    if clause1 != clause2 and clause2 in supportset: #其中一个亲本子句必须来自支持集
                        literal_index1 = 0
                        for literal1 in clause1:
                            literal_index2 = 0
                            for literal2 in clause2:
                                #判断是否互补
                                if Complementary(literal1,literal2):
                                    #得到索引
                                    index1 = Index(literal_index1,clause_index1,len(clause1))
                                    index2 = Index(literal_index2,clause_index2,len(clause2))
                                    #得到互补文字
                                    literal1 = clause1[literal_index1]
                                    literal2 = clause2[literal_index2]
                                    #得到参数列表
                                    para1 = Parameter(literal1)
                                    para2 = Parameter(literal2)
                                    #得到合一置换
                                    unification = Unify(para1,para2)
                                    #合一置换不存在则退出
                                    if unification == None:
                                        break
                                    #得到子句
                                    newclause1 = Substitute(unification,clause1)    
                                    newclause2 = Substitute(unification,clause2)  
                                    #得到归结子句
                                    newclause = Resolve(newclause1,newclause2,literal_index1,literal_index2)
                                    #归结子句存在于原子句集则退出
                                    if any([set(newclause)==(set(item)) for item in clauseset]):
                                        break
                                    #归结子句存在于新子句集则退出
                                    if any([set(newclause)==(set(item)) for item in new_clauseset]):
                                        break
                                    #得到新归结式
                                    sequence = Sequence(newclause,unification,index1,index2)
                                    #加入结果列表和新子句列表
                                    result.append(sequence)
                                    new_clauseset.append(newclause)
                                    #当且仅当新子句为空,则可退出函数,返回结果列表
                                    if newclause == ():
                                        return result
                                literal_index2 += 1
                            literal_index1 += 1
                    clause_index2 += 1
            clause_index1 += 1
        clauseset += new_clauseset      #更新子句列表
        supportset += new_clauseset     #更新支持集

#得到归结式的子句索引
def Number(clause):
    start = clause.find('[')
    end = clause.find(']')
    number = clause[start+1:end].split(',')
    #将文字索引去掉
    num1 = int(''.join(item for item in number[0] if not item.isalpha()))
    num2 = int(''.join(item for item in number[1] if not item.isalpha()))
    return num1,num2

#得到新归结式的子句索引
def Renumber(num,result,useful_process,size):
    if num <= size: #如果是初始子句集的,直接返回
        return num
    #找到亲本子句
    sequence = result[num]
    begin = sequence.find('(')
    aim_clause = sequence[begin:]
    #找到亲本子句在化简子句集的编号
    for i in range(size+1,len(useful_process)):
        begin = useful_process[i].find('(')
        if useful_process[i][begin:] == aim_clause:
            return i 

#更新归结式
def Resequence(sequence,num1,num2,newnum1,newnum2):
    # 第一次替换:替换第一个编号
    start = sequence.find(num1)
    end = start + len(num1)
    sequence = sequence[:start] + newnum1 + sequence[end:]
    # 第二次替换:替换第二个编号
    end = start + len(newnum1)
    start = sequence.find(num2, end)
    end = start + len(num2)
    sequence = sequence[:start] + newnum2 + sequence[end:]
    return sequence

#化简归结过程
def Simplify(result,size):
    base_process = result[0:size+1] #初始子句集
    useful_process = []             #有用子句集
    number = [len(result)-1]        #用作队列,先将空子句的索引入列
    while number != []:
        number0 = number.pop(0)                 #提取队列首元素,即有用子句的索引
        useful_process.append(result[number0])  #将有用子句加入到有用子句集            
        num1,num2 = Number(result[number0])     #得有用子句用到的亲本子句索引
        #如果是初始子句集就无需加入
        if num1 > size:
            number.append(num1)
        if num2 > size:
            number.append(num2)
    #得到新的归结过程
    useful_process.reverse()
    useful_process = base_process + useful_process
    #将归结过程重新编号
    for i in range(size+1,len(useful_process)): 
        num1,num2 = Number(useful_process[i])
        newnum1 = str(Renumber(num1,result,useful_process,size))
        newnum2 = str(Renumber(num2,result,useful_process,size))
        useful_process[i] = Resequence(useful_process[i],str(num1),str(num2),newnum1,newnum2)
    return useful_process

#打印结果
def Print(result):
    print(result[0])
    for i in range(1,len(result)):
        print(i,result[i])

#归结反演
def ResolutionFOL(KB):
    result = Refutation(KB)
    new_result = Simplify(result,len(KB))
    Print(new_result)

 

五、测试

测试代码

KB1 = [('GradStudent(sue)',),('~GradStudent(x)','Student(x)'),('~Student(x)','HardWorker(x)'),
     ('~HardWorker(sue)',)]

KB2 = [('On(tony,mike)',),('On(mike,john)',),('Green(tony)',),('~Green(john)',),
      ('~On(xx,yy)','~Green(xx)','Green(yy)')]

KB3 = [('A(tony)',), ('A(mike)',), ('A(john)',), ('L(tony,rain)',), ('L(tony,snow)',),
      ('~A(x)', 'S(x)', 'C(x)'), ('~C(y)', '~L(y,rain)'), ('L(z,snow)', '~S(z)'),
      ('~L(tony,u)', '~L(mike,u)'), ('L(tony,v)', 'L(mike,v)'), ('~A(w)', '~C(w)', 'S(w)')]

ResolutionFOL(KB1)
ResolutionFOL(KB2)
ResolutionFOL(KB3)

结果:在ResolutionFOL函数中加了一步,输出还没化简的归结过程,用于比较分析-->Print(result)

 

 

 

 

六、总结

(1) 评价

  • 时间复杂度主要取决于得到未化简归结过程
  • 完备的

(2) 不足之处

  • 归结过程中每次迭代,选择clause1和clause2都是从clauseset的头到尾,这会产生许多重复判定,造成不必要的时间开销
  • 支持集策略只能指导其中一个子句的选择,无法同时指导两个子句的选择
  • 空间复杂度较高

(3) 改进方法

  • 结合多个搜索策略:可以将多个策略联合使用,判断什么情况下用什么策略更好
  • 添加删除策略:应该删除纯文字、重言式、类含子句等
  • 引入启发式方法:可以根据子句的复杂度、覆盖的变量范围等因素进行评估,并优先选择可能导致有效归结的子句进行处理

七、参考资料

  • 47
    点赞
  • 36
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值