NLP之句子相似度之入门篇
Posted AI蜗牛之家
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了NLP之句子相似度之入门篇相关的知识,希望对你有一定的参考价值。
文章目录
如下在师兄的博文基础上修改: 静觅 » 自然语言处理中句子相似度计算的几种方法
1.基于统计的方法
1.1.编辑距离计算
编辑距离,英文叫做 Edit Distance,又称 Levenshtein 距离,是指两个字串之间,由一个转成另一个所需的最少编辑操作次数,如果它们的距离越大,说明它们越是不同。许可的编辑操作包括将一个字符替换成另一个字符,插入一个字符,删除一个字符。
例如我们有两个字符串:string 和 setting,如果我们想要把 string 转化为 setting,需要这么两步:
- 第一步,在 s 和 t 之间加入字符 e。
- 第二步,把 r 替换成 t。
所以它们的编辑距离差就是 2,这就对应着二者要进行转化所要改变(添加、替换、删除)的最小步数。
那么用 Python 怎样来实现呢,我们可以直接使用 distance 库:
import distance
def edit_distance(s1, s2):
return distance.levenshtein(s1, s2)
s1 = 'string'
s2 = 'setting'
print(edit_distance(s1, s2))
这里我们直接使用 distance 库的 levenshtein() 方法,传入两个字符串,即可获取两个字符串的编辑距离了。
运行结果如下:
2
这里的 distance 库我们可以直接使用 pip3 来安装:pip3 install distance
这样如果我们想要获取相似的文本的话可以直接设定一个编辑距离的阈值来实现,如设置编辑距离为 2,下面是一个样例:
import distance
def edit_distance(s1, s2):
return distance.levenshtein(s1, s2)
strings = [
'你在干什么',
'你在干啥子',
'你在做什么',
'你好啊',
'我喜欢吃香蕉'
]
target = '你在干啥'
results = list(filter(lambda x: edit_distance(x, target) <= 2, strings))
print(results)
这里我们定义了一些字符串,然后定义了一个目标字符串,然后用编辑距离 2 的阈值进行设定,最后得到的结果就是编辑距离在 2 及以内的结果,运行结果如下:
['你在干什么', '你在干啥子']
通过这种方式我们可以大致筛选出类似的句子,但是发现一些句子例如“你在做什么” 就没有被识别出来,但他们的意义确实是相差不大的,因此,编辑距离并不是一个好的方式,但是简单易用。
1.2.杰卡德系数计算
杰卡德系数,英文叫做 Jaccard index, 又称为 Jaccard 相似系数,用于比较有限样本集之间的相似性与差异性。Jaccard 系数值越大,样本相似度越高。
实际上它的计算方式非常简单,就是两个样本的交集除以并集得到的数值,当两个样本完全一致时,结果为 1,当两个样本完全不同时,结果为 0。
算法非常简单,就是交集除以并集,下面我们用 Python 代码来实现一下:
from sklearn.feature_extraction.text import CountVectorizer
import numpy as np
def jaccard_similarity(s1, s2):
def add_space(s):
return ' '.join(list(s))
# 将字中间加入空格
s1, s2 = add_space(s1), add_space(s2)
# 转化为TF矩阵
cv = CountVectorizer(tokenizer=lambda s: s.split())
corpus = [s1, s2]
vectors = cv.fit_transform(corpus).toarray()
# 求交集
numerator = np.sum(np.min(vectors, axis=0))
# 求并集
denominator = np.sum(np.max(vectors, axis=0))
# 计算杰卡德系数
return 1.0 * numerator / denominator
s1 = '你在干嘛呢'
s2 = '你在干什么呢'
print(jaccard_similarity(s1, s2))
这里我们使用了 Sklearn 库中的 CountVectorizer 来计算句子的 TF 矩阵,然后利用 Numpy 来计算二者的交集和并集,随后计算杰卡德系数。
这里值得学习的有 CountVectorizer 的用法,通过它的 fit_transform() 方法我们可以将字符串转化为词频矩阵,例如这里有两句话“你在干嘛呢”和“你在干什么呢”,首先 CountVectorizer 会计算出不重复的有哪些字,会得到一个字的列表,结果为:
['么', '什', '你', '呢', '嘛', '在', '干']
这个其实可以通过如下代码来获取,就是获取词表内容:cv.get_feature_names()
接下来通过转化之后,vectors 变量就变成了:
[[0 0 1 1 1 1 1]
[1 1 1 1 0 1 1]]
它对应的是两个句子对应词表的词频统计,这里是两个句子,所以结果是一个长度为 2 的二维数组,比如第一句话“你在干嘛呢”中不包含“么”字,那么第一个“么”字对应的结果就是0,即数量为 0,依次类推。
后面我们使用了 np.min() 方法并传入了 axis 为 0,实际上就是获取了每一列的最小值,这样实际上就是取了交集,np.max() 方法是获取了每一列的最大值,实际上就是取了并集。
二者分别取和即是交集大小和并集大小,然后作商即可,结果如下:
0.5714285714285714
这个数值越大,代表两个字符串越接近,否则反之,因此我们也可以使用这个方法,并通过设置一个相似度阈值来进行筛选。
1.3.TF 计算
第三种方案就是直接计算 TF 矩阵中两个向量的余弦相似度了,实际上就是求解两个向量夹角的余弦值,就是点乘积除以二者的模长,公式如下:
c
o
s
θ
=
a
⋅
b
∣
a
∣
∣
b
∣
cosθ=\\frac a \\cdot b|a||b|
cosθ=∣a∣∣b∣a⋅b
更多关于余弦相似度:
上面我们已经获得了 TF 矩阵,下面我们只需要求解两个向量夹角的余弦值就好了,代码如下:
from sklearn.feature_extraction.text import CountVectorizer
import numpy as np
from scipy.linalg import norm
def tf_similarity(s1, s2):
def add_space(s):
return ' '.join(list(s))
# 将字中间加入空格
s1, s2 = add_space(s1), add_space(s2)
# 转化为TF矩阵
cv = CountVectorizer(tokenizer=lambda s: s.split())
corpus = [s1, s2]
vectors = cv.fit_transform(corpus).toarray()
# 计算TF系数
return np.dot(vectors[0], vectors[1]) / (norm(vectors[0]) * norm(vectors[1]))
s1 = '你在干嘛呢'
s2 = '你在干什么呢'
print(tf_similarity(s1, s2))
在在这里我们使用了 np.dot() 方法获取了向量的点乘积,然后通过 norm() 方法获取了向量的模长,经过计算得到二者的 TF 系数,结果如下:
0.7302967433402214
1.4.TFIDF 计算
另外除了计算 TF 系数我们还可以计算 TFIDF 系数,TFIDF 实际上就是在词频 TF 的基础上再加入 IDF 的信息,IDF 称为逆文档频率,不了解的可以看下阮一峰老师的讲解:http://www.ruanyifeng.com/blog/2013/03/tf-idf.html,里面对 TFIDF 的讲解也是十分透彻的。
关于idf(t)应该这样认识:一个词语在文档集合中出现了n次,文档集合总数为N。idf(t)来源于信息论。那么每篇文档出现这个词语的概率为:
n
/
N
n/N
n/N,所以这篇文档出现这个词语的信息量为:
−
l
o
g
(
n
/
N
)
-log(n/N)
−log(n/N)。这个和信息熵有些类似
(
−
P
(
x
)
l
o
g
P
(
x
)
)
(-P(x)logP(x))
(−P(x)logP(x)),在数据挖掘的过滤法进行特征选择时,需要用到互信息,其实是计算信息增益,还有决策树。把
−
l
o
g
(
n
/
N
)
-log(n/N)
−log(n/N)变换一下,
l
o
g
(
N
/
n
)
log(N/n)
log(N/n),为了避免0的出现,进行平滑处理,就是上面的公式(就像朴素贝叶斯需要拉普拉斯平滑处理一样)。
下面我们还是借助于 Sklearn 中的模块 TfidfVectorizer 来实现,代码如下:
from sklearn.feature_extraction.text import TfidfVectorizer
import numpy as np
from scipy.linalg import norm
def tfidf_similarity(s1, s2):
def add_space(s):
return ' '.join(list(s))
# 将字中间加入空格
s1, s2 = add_space(s1), add_space(s2)
# 转化为TF矩阵
cv = TfidfVectorizer(tokenizer=lambda s: s.split())
corpus = [s1, s2]
vectors = cv.fit_transform(corpus).toarray()
# 计算TF系数
return np.dot(vectors[0], vectors[1]) / (norm(vectors[0]) * norm(vectors[1]))
s1 = '你在干嘛呢'
s2 = '你在干什么呢'
print(tfidf_similarity(s1, s2))
这里的 vectors 变量实际上就对应着 TFIDF 值,内容如下:
[[0. 0. 0.4090901 0.4090901 0.57496187 0.4090901 0.4090901 ]
[0.49844628 0.49844628 0.35464863 0.35464863 0. 0.35464863 0.35464863]]
运行结果如下:
0.5803329846765686
所以通过 TFIDF 系数我们也可以进行相似度的计算。
1.5.BM25
BM25 is a bag-of-words retrieval function that ranks a set of documents based on the query terms appearing in each document, regardless of the inter-relationship between the query terms within a document (e.g., their relative proximity). It is not a single function, but actually a whole family of scoring functions, with slightly different components and parameters. One of the most prominent instantiations of the function is as follows.
BM25算法,通常用来作搜索相关性平分。一句话概况其主要思想:对Query进行语素解析,生成语素qi;然后,对于每个搜索结果D,计算每个语素qi与D的相关性得分,最后,将qi相对于D的相关性得分进行加权求和,从而得到Query与D的相关性得分。
BM25算法的一般性公式如下:
S
c
o
r
e
(
Q
,
d
)
=
∑
i
n
W
i
⋅
R
(
q
i
,
d
)
Score(Q,d) = \\sum _i ^n W_i \\cdot R(q_i,d)
Score(Q,d)=i∑nWi⋅R(qi,d)
其中,
Q
Q
Q表示Query,
q
i
q_i
qi表示
Q
Q
Q解析之后的一个语素(对中文而言,我们可以把对Query的分词作为语素分析,每个词看成语素qi。);
d
d
d表示一个搜索结果文档;
W
i
W_i
Wi表示语素
q
i
q_i
qi的权重;
R
(
q
i
,
d
)
R(q_i,d)
R(qi,d)表示语素
q
i
q_i
qi与文档
d
d
d的相关性得分。
下面我们来看如何定义
W
i
W_i
Wi。判断一个词与一个文档的相关性的权重,方法有多种,较常用的是IDF。这里以IDF为例,公式如下:
I
D
F
(
q
i
)
=
l
o
g
N
−
n
(
q
i
)
+
0.5
n
(
q
i
)
+
0.5
IDF(q_i) = log \\frac N-n(q_i)+0.5 n(q_i)+0.5
IDF(qi)=logn(qi)+0.5N−n(qi)+0.5
其中, N N N为索引中的全部文档数, n ( q i ) n(q_i) n(qi)为包含了qi的文档数。
根据IDF的定义可以看出,对于给定的文档集合,包含了 q i q_i qi的文档数越多, q i q_i qi的权重则越低。也就是说,当很多文档都包含了 q i q_i qi时, q i q_i qi的区分度就不高,因此使用 q i q_i qi来判断相关性时的重要度就较低。
我们再来看语素
q
i
q_i
qi与文档d的相关性得分
R
(
q
i
,
d
)
R(q_i,d)
R(qi,d)。首先来看BM25中相关性得分的一般形式:
R
(
q
i
,
d
)
=
f
i
(
k
1
+
1
)
f
i
+
K
⋅
d
f
i
(
k
2
+
1
)
q
f
i
+
k
2
R(q_i,d) = \\frac f_i(k_1+1) f_i+K \\cdot \\fracdf_i(k_2+1)qf_i+k_2
R(qi,d)=fi+Kfi(k1+1)以上是关于NLP之句子相似度之入门篇的主要内容,如果未能解决你的问题,请参考以下文章
NLP之gensim库python实现文本相似度/匹配/查重