[發(fā)明專利]一種基于聚類算法的移動(dòng)用戶一人多卡識(shí)別方法

在審

• 申請(qǐng)?zhí)枺?/td>: 201811593501.X
• 申請(qǐng)日: 2018-12-25
• 公開(kāi)（公告）號(hào): CN111436017A
• 公開(kāi)（公告）日: 2020-07-21
• 發(fā)明（設(shè)計(jì)）人: 楊占軍;張廣志;秦星星;于海薇
• 申請(qǐng)（專利權(quán)）人: 北京融信數(shù)聯(lián)科技有限公司
• 主分類號(hào): H04W4/02
• 分類號(hào): H04W4/02;H04W4/029;H04W8/18
• 代理公司: 北京兆君聯(lián)合知識(shí)產(chǎn)權(quán)代理事務(wù)所(普通合伙) 11333
• 代理人: 初向慶

1.一種基于聚類算法的移動(dòng)用戶一人多卡識(shí)別方法，其特征在于：

(1)獲取電信運(yùn)營(yíng)商信令數(shù)據(jù)，對(duì)每個(gè)IMSI識(shí)別號(hào)所在的基站位置、進(jìn)出小區(qū)時(shí)間信息進(jìn)行采集并存儲(chǔ)；

(2)將城市區(qū)域劃分為若干個(gè)網(wǎng)格，對(duì)同一網(wǎng)格區(qū)間，計(jì)算區(qū)間內(nèi)包含的所有基站的平均經(jīng)緯度，作為該網(wǎng)格的經(jīng)緯度；

(3)對(duì)每個(gè)基站，計(jì)算其到各個(gè)網(wǎng)格之間的經(jīng)緯度，將其歸為某個(gè)網(wǎng)格區(qū)間，實(shí)現(xiàn)區(qū)間的重新劃分；

(4)重復(fù)步驟(3)，直到算法收斂；

(5)每次讀取一對(duì)IMSI識(shí)別號(hào)，設(shè)為a IMSI和b IMSI，分別找到其所對(duì)應(yīng)的網(wǎng)格經(jīng)緯度，計(jì)算兩個(gè)IMSI在各時(shí)間點(diǎn)下的距離，計(jì)算方法如下：

(51)先找到a IMSI的一個(gè)時(shí)間點(diǎn)t，再找到b IMSI在t時(shí)間點(diǎn)前的最后一個(gè)時(shí)間點(diǎn)t1以及在t時(shí)間點(diǎn)后的第一個(gè)時(shí)間點(diǎn)t2；

(52)對(duì)應(yīng)出a IMSI在時(shí)間點(diǎn)t所在位置的經(jīng)緯度，以及b IMSI在時(shí)間點(diǎn)t1和t2時(shí)的位置經(jīng)緯度；

(53)分別計(jì)算出a IMSI在t時(shí)間點(diǎn)和b IMSI在t1時(shí)間點(diǎn)時(shí)的距離，以及a IMSI在t時(shí)間點(diǎn)和b IMSI在t2時(shí)間點(diǎn)時(shí)的距離，計(jì)算公式如下：

Dis1[t]＝R*acos(sinpi(y1/180)*sinpi(y2/180)+cospi(y1/180)*cospi(y2/180)*cospi((x1-x2)/180))；

公式中R表示地球半徑；Dis1[t]表示a IMSI在t時(shí)刻與b IMSI在該時(shí)刻前的最后一個(gè)時(shí)間點(diǎn)t1的距離；x1、x2分別表示兩個(gè)loc-ci的經(jīng)度，y1、y2分別表示兩個(gè)loc-ci的緯度；

同理可以計(jì)算出a IMSI在t時(shí)刻與b IMSI在該時(shí)刻后的第一個(gè)時(shí)間點(diǎn)t2的距離Dis2[t]；

同時(shí)，計(jì)算出兩個(gè)時(shí)間差：

Time1[t]＝t-t1；Time2[t]＝t2-t；

(54)計(jì)算a IMSI在t時(shí)間點(diǎn)下和b IMSI在t1時(shí)間點(diǎn)的距離閾值Disthr1[t]和最大距離Dismax1[t]：

Disthr1[t]＝Speedthr*Time1[t]；

Dismax1[t]＝Speedmax*Time1[t]；

這里Speedthr表示速度閾值，設(shè)為10公里/小時(shí)；Speedmax表示最大速度，設(shè)為20公里/小時(shí)；

同理，計(jì)算出a IMSI在t時(shí)間點(diǎn)下和b IMSI在t2時(shí)間點(diǎn)的距離閾值Disthr2[t]和最大距離Dismax2[t]；

(6)按照上述方法進(jìn)行循環(huán)計(jì)算，分別計(jì)算出a IMSI在各個(gè)時(shí)間點(diǎn)i下與b IMSI最接近的兩個(gè)時(shí)間點(diǎn)的距離Dis1[i]和Dis2[i]，以及兩個(gè)距離閾值Disthr1[i]和Disthr2[i],兩個(gè)最大距離值Dismax1[i]和Dismax2[i]；

(7)根據(jù)上述計(jì)算得到的距離值、距離閾值、和最大距離值，進(jìn)行一人多號(hào)判斷：

(71)如果在任意時(shí)間點(diǎn)i，距離Dis1[i]和Dis2[i]分別大于Dismax1[i]和Dismax2[i]，則直接判斷該對(duì)IMSI為非同一軌跡，即非一人多號(hào)；否則，進(jìn)一步按以下步驟(72)進(jìn)行置信度判斷；

(72)如果距離Dis1[i]小于或等于距離閾值Disthr1[i]時(shí)，按以下公式計(jì)算似然值L1[i]：

L1[i]＝1-(Dis1[i]/Disthr1[i])*α；

如果距離Dis1[i]大于距離閾值Disthr1[i]，而小于或等于最大距離Dismax1[i]時(shí)，則似然值L1[i]按以下公式計(jì)算：

L1[i]＝(1-α)-(Dis1[i]/Dismax1[i])*0.1；

上述公式中，α為經(jīng)驗(yàn)值，默認(rèn)為α＝0.6；

同理，計(jì)算出似然值L2[i]；

計(jì)算出一對(duì)IMSI在各時(shí)間點(diǎn)上似然值的均值；

進(jìn)一步，設(shè)定一個(gè)閾值β，β＝0.8；

如果一對(duì)IMSI在各時(shí)間點(diǎn)上的似然值的均值小于閾值β，則認(rèn)為該對(duì)IMSI為非同一軌跡，認(rèn)定為非一人多號(hào)；

如果一對(duì)IMSI在各時(shí)間點(diǎn)上的似然值的均值大于等于閾值β，則認(rèn)為該對(duì)IMSI在所述的時(shí)間段內(nèi)為同一軌跡；設(shè)定一個(gè)包含有多個(gè)連續(xù)的所述時(shí)間段的時(shí)間周期，如果該對(duì)IMSI在一個(gè)時(shí)間周期內(nèi)有大于設(shè)定數(shù)量的時(shí)間段被判定為同一軌跡，則認(rèn)定該對(duì)IMSI為一人多號(hào)。