CN107563801A - 一种基于隐马尔可夫模型的消费者线下消费行为预测方法 - Google Patents
一种基于隐马尔可夫模型的消费者线下消费行为预测方法 Download PDFInfo
- Publication number
- CN107563801A CN107563801A CN201710728972.6A CN201710728972A CN107563801A CN 107563801 A CN107563801 A CN 107563801A CN 201710728972 A CN201710728972 A CN 201710728972A CN 107563801 A CN107563801 A CN 107563801A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msubsup
- consumer
- hidden markov
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于隐马尔可夫模型的消费者线下消费行为预测方法，包括如下步骤：步骤一、数据预处理得到观测序列；步骤二、初始化隐马尔可夫模型；步骤三、使用Baum‑Welch算法和观测序列(o1o2...ot)对模型进行迭代；步骤四、使用Viterbi算法和训练好的模型，计算观测序列(o1o2...ot)对应的最优的状态序列，然后预测t+1时刻的最有可能的状态yt+1，然后计算t+1时刻概率最大的观测值；步骤五、通过计算的t+1时刻概率最大的观测值ot+1并对照商家编号即可得到消费者在t+1时刻最有可能消费的线下商店shopt+1。本发明的有益效果是：提供了一种基于隐马尔可夫模型的消费者线下消费行为的预测方法，使用概率图模型更加准确高效地对不同用户的线下消费行为进行预测。
Description
技术领域
本发明涉及一种基于隐马尔可夫模型的消费者线下消费行为的预测方法，主要是一种对消费者历史线下消费数据进行处理，通过隐马尔可夫模型预测该用户未来最有可能消费的线下商店的方法。
背景技术
随着电子商务的不断发展，O2O模式和手机移动支付也在迅速的发展。用户在线下消费时会产生大量的消费数据，这些数据中包含了很多有用的信息。通过分析这些信息，可以分析各个消费者的消费习惯从而预测消费者未来可能的消费行为，这对给用户提供更好的个性化推荐服务和商家优化经营和管理策略都是非常重要的。
如何更好的利用历史消费数据从而有效快速地预测消费者未来的消费行为，是本领域技术人员急需解决的问题。
发明内容
本发明的目的是克服现有技术中的不足，提供一种基于隐马尔可夫模型的消费者线下消费行为的预测方法。
基于隐马尔可夫模型的消费者线下消费行为预测方法，包括如下步骤：
步骤一、数据预处理得到观测序列：
抽取消费者的长度为t的历史商家消费数据{shop1，shop2，…，shopt}并对消费者的历史消费线下商家进行编号；然后根据时间顺序把消费数据按照商家编号转换成观测值序列(o1o2…ot)；
步骤二、初始化隐马尔可夫模型：
通过一个三元组参数确定一个隐马尔可夫模型，λ＝(A，B，π)指代，分别包含转移矩阵A、混淆矩阵B和初始状态矩阵π；给参数A，B，π置随机初始值λ＝（A1，B1，π1)；
步骤三、使用Baum-Welch算法和观测序列(o1o2…ot)对模型进行迭代：
已知观测序列O＝(o1o2…oT)，估计参数模型λ＝（A，B，π)，使得该模型下观测概率P＝(O|λ)最大；
Baum-Welch算法具体步骤：
输入：观测数据O＝(o1,o2…oT)
输出：隐马尔可夫模型参数
(1)初始化：对于n＝0,选取得到模型λ(0)＝(A(0)，B(0)，π(0))
(2)递推，对于n＝1，2，…，
其中，αij＝P(yt+1＝Sj|yt＝Si) bij＝P(xt＝Oj|yt＝Si) πi＝P(y1＝Si)，
αt(i)为前向变量，即输出为(o1，o2，…ot)，在t时刻位于状态Si的概率；βt(i)为后向变量，即时刻t状态为Si后输出为(ot+1，ot+2…oT)的概率；
(3)终止，得到模型参数λ(n+1)＝(A(n+1)，B(n+1)，π(n+1))
步骤四、使用Viterbi算法和训练好的模型，计算观测序列（o1o2…ot)对应的最优的状态序列；然后预测t+1时刻的最有可能的状态yt+1，然后计算t+1时刻概率最大的观测值；
Viterbi算法具体步骤：
输入：观测数据o＝(o1，o2…oT)
输出：隐马尔可夫模型参数
输出：最优隐状态概率I＝(i1，i2…iT)
(1)初始化：δ1(i)＝πibi(o1),1≤i≤N ψ1(i)＝0
(2)递推，对于t＝2，3，…，T
其中，δt(i)是Viterbi变量，ψt(i)是路径记忆变量；
(3)终止，
(4)回溯最优路径，t＝T-1，T-2，…，1
it＝ψt+1(it+1)；
步骤五、通过计算的t+1时刻概率最大的观测值ot+1并对照商家编号即可得到消费者在t+1时刻最有可能消费的线下商店shopt+1。
作为优选:步骤二中，隐马尔可夫模型参数包括三元组λ＝(A，B，π)，具体含义：
转移矩阵
其中aij＝P(yt+1＝Sj|yt＝Si)(1≤i,j≤N)表示t时刻状态是Si且t+1时刻状态是Sj的概率；
混淆矩阵
其中bij＝P(xt＝Oj|yt＝Si)(1≤i≤N,1≤j≤M)表示t时刻状态是Si且观测值是时刻状态是Oj的概率；
初始状态矩阵πi＝P(y1＝Si)，1≤i≤N
其中πi表示在初始t＝1时刻状态是Si的概率。
作为优选：步骤三中，隐马尔可夫模型中的变量分为2组；第一组是状态变量Y＝{y1,y2,…,yn}，其中yi表示第i时刻的系统状态；假定状态是隐藏的、不可被观测的；第二组是观测值变量X＝{x1,x2,…,xn}，其中xi表示第i时刻的观测值；在隐马尔可夫模型中，系统在多个状态S＝{s1,s2,…,sN}之间转换，因此状态变量yi的取值范围为N个取值的离散空间。
本发明的有益效果是:提供了一种基于隐马尔可夫模型的消费者线下消费行为的预测方法，使用概率图模型更加准确高效地对不同用户的线下消费行为进行预测。
附图说明
图1是本方法流程图；
图2是隐马尔可夫模型示意图；
图3是实验性能比较图。
具体实施方式
下面结合实施例对本发明做进一步描述。下述实施例的说明只是用于帮助理解本发明。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。
如图1所示，本实施方式中的基于隐马尔可夫模型的消费者线下消费行为预测方法，其步骤如下：
1.数据预处理得到观测序列：
抽取消费者的长度为t的历史商家消费数据{shop1，shop2，…，shopt}并对消费者的历史消费线下商家进行编号。然后根据时间顺序把消费数据按照商家编号转换成观测值序列(o1o2…ot)
2.初始化隐马尔可夫模型：
通过一个三元组参数就能确定一个隐马尔可夫模型，通常λ＝(A，B，π)指代，分别包含转移矩阵A、混淆矩阵B和初始状态矩阵π。初始化一个隐马尔可夫模型即给参数A，B，π置随机初始值λ＝(A1，B1，π1)。
3.使用Baum-Welch算法和观测序列(o1o2…ot)对模型进行迭代：
如图2所示，隐马尔可夫模型中的变量可分为2组。第一组是状态变量Y＝{y1,y2,…,yn}，其中yi表示第i时刻的系统状态。通常假定状态是隐藏的、不可被观测的，因此状态变量也被称为隐变量。第二组是观测值变量X＝{x1,x2,…,xn}，其中xi表示第i时刻的观测值。在隐马尔可夫模型中，系统通常在多个状态S＝{s1,s2,…,sN}之间转换，因此状态变量yi的取值范围通常是有N个取值的离散空间。
Baum-Welch算法是为了解决隐马尔可夫模型的参数估计问题而提出的。具体来说就是已知观测序列O＝(o1o2…oT)，估计参数模型λ＝(A，B，π)，使得该模型下观测概率P＝(O|λ)最大。
Baum-Welch算法具体步骤：
输入：观测数据O＝(o1，o2…oT)
输出：隐马尔可夫模型参数
(1)初始化：对于n＝0,选取得到模型λ(0)＝(A(0)，B(0)，π(0))
(2)递推，对于n＝1，2，…，
其中，aij＝P(yt+1＝Sj|yt＝Si) bij＝P(xt＝Oj|yt＝Si) πi＝P(y1＝Si)，
αt(i)为前向变量，即输出为(o1，o2…ot)，在t时刻位于状态Si的概率；βt(i)为后向变量，即时刻t状态为Si后输出为(ot+1，ot+2…oT)的概率。
(3)终止，得到模型参数λ(n+1)＝(A(n+1)，B(n+1)，π(n+1))
4.使用Viterbi算法和训练好的模型，计算观测序列(o1o2…ot)对应的最优的状态序列。然后预测t+1时刻的最有可能的状态yt+1，然后计算t+1时刻概率最大的观测值。
Viterbi算法是一种动态规划算法，可用于为给定的观测序列和模型找到最有可能的隐藏状态序列。
Viterbi算法具体步骤：
输入：观测数据O＝(o1，o2…oT)
输出：隐马尔可夫模型参数
输出：最优隐状态概率I＝(i1，i2…iT)
(1)初始化：δ1(i)＝πibi(o1),1≤i≤N ψ1(i)＝0
(2)递推，对于t＝2，3，…，T
其中，δt(i)是Viterbi变量，ψt(i)是路径记忆变量。
(3)终止，
(4)回溯最优路径，t＝T-1，T-2，…，1
it＝ψt+1(it+1)；
5.通过计算的t+1时刻概率最大的观测值ot+1并对照商家编号即可得到消费者在t+1时刻最有可能消费的线下商店shopt+1。
实验及结果：
为了验证该方法的预测效果，我们在一个真实的数据集上进行了实验。该数据集包含1057个不同的用户和2000个商家，每个用户拥有120次以上的消费记录。并且为了验证不同训练长度对预测结果的影响，我们选取了(20，40，60，80，100)五种不同的训练长度对这1057个用户进行实验。最终实验结果表明，我们的预测方法在训练长度为60时准确度最高，达到了82.02％，在训练长度为20时准确度最低，也有80.86％。同时我们也与最频繁项预测方法(most frequent,简称MF)和仅使用转移矩阵的隐马尔可夫预测方法(简称HMM-PTP)进行对比(如图3所示)。MF方法最高的预测准确率为81.23％，最低为80.24％；HMM-PTP方法的最高预测准确率为72.89％，最低为71.91％。由此可见，相对于其他两种预测方法，我们的方法能取得更好的预测效果。
Claims (3)
1.一种基于隐马尔可夫模型的消费者线下消费行为预测方法，其特征在于：包括如下步骤：
步骤一、数据预处理得到观测序列：
抽取消费者的长度为t的历史商家消费数据{shop1，shop2，...，shopt}并对消费者的历史消费线下商家进行编号；然后根据时间顺序把消费数据按照商家编号转换成观测值序列(o1o2...ot)；
步骤二、初始化隐马尔可夫模型：
通过一个三元组参数确定一个隐马尔可夫模型，λ＝(A，B，π)指代，分别包含转移矩阵A、混淆矩阵B和初始状态矩阵π；给参数A，B，π置随机初始值λ＝(A1，B1，π1)；
步骤三、使用Baum-Welch算法和观测序列(o1o2...ot)对模型进行迭代：
已知观测序列O＝(o1o2...oT)，估计参数模型λ＝(A，B，π)，使得该模型下观测概率P＝(O|λ)最大；
Baum-Welch算法具体步骤：
输入：观测数据O＝(o1，o2...oT)
输出：隐马尔可夫模型参数
1)初始化：对于n＝0,选取得到模型λ(0)＝(A(0)，B(0)，π(0))
2)递推，对于n＝1，2，...，
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>a</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>&xi;</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>&gamma;</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>b</mi>
<mi>j</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<msub>
<mi>o</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
</mrow>
<mrow>
<mi>T</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>&gamma;</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</msubsup>
<msub>
<mi>&gamma;</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>&pi;</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中，aij＝P(yt+1＝Sj|yt＝Si) bij＝P(xt＝Oj|yt＝Si) πi＝P(y1＝Si)，
<mrow>
<msub>
<mi>&xi;</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>i</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mi>i</mi>
<mo>,</mo>
<msub>
<mi>i</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mi>j</mi>
<mo>,</mo>
<mi>O</mi>
<mo>|</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>O</mi>
<mo>|</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msub>
<mi>&alpha;</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>b</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>o</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>P</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>i</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mi>i</mi>
<mo>,</mo>
<mi>O</mi>
<mo>|</mo>
<mi>&lambda;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>&beta;</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msub>
<mi>&alpha;</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>b</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>o</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
αt(i)为前向变量，即输出为(o1，o2...ot)，在t时刻位于状态Si的概率；βt(i)为后向变量，即时刻t状态为Si后输出为(ot+1，ot+2...oT)的概率；
3)终止，得到模型参数λ(n+1)＝(A(n+1)，B(n+1)，π(n+1))；
步骤四、使用Viterbi算法和训练好的模型，计算观测序列(o1o2...ot)对应的最优的状态序列；然后预测t+1时刻的最有可能的状态yt+1，然后计算t+1时刻概率最大的观测值；
Viterbi算法具体步骤：
输入：观测数据O＝(o1，o2...oT)
输出：隐马尔可夫模型参数
输出：最优隐状态概率I＝(i1，i2...iT)
1)初始化：δ1(i)＝πibi(o1),1≤i≤N ψ1(i)＝0
2)递推，对于t＝2，3，...，T
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&delta;</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mi>max</mi>
<mrow>
<mn>1</mn>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<mi>N</mi>
</mrow>
</munder>
<mo>&lsqb;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
<msub>
<mi>b</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>o</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&psi;</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mrow>
<mi>arg</mi>
<mi>max</mi>
</mrow>
<mrow>
<mn>1</mn>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<mi>N</mi>
</mrow>
</munder>
<mo>&lsqb;</mo>
<msub>
<mi>&delta;</mi>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
<msub>
<mi>b</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>o</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中，δt(i)是Viterbi变量，ψt(i)是路径记忆变量；
3)终止，
4)回溯最优路径，t＝T-1，T-2，...，1
it＝ψt+1(it+1)；
步骤五、通过计算的t+1时刻概率最大的观测值ot+1并对照商家编号即可得到消费者在t+1时刻最有可能消费的线下商店shopt+1。
2.根据权利要求1所述的基于隐马尔可夫模型的消费者线下消费行为预测方法，其特征在于：步骤二中，隐马尔可夫模型参数包括三元组λ＝(A，B，π)，具体含义：
转移矩阵
其中aij＝P(yt+1＝Sj|yt＝Si)(1≤i,j≤N)表示t时刻状态是Si且t+1时刻状态是Sj的概率；混淆矩阵
其中bij＝P(xt＝Oj|yt＝Si)(1≤i≤N,1≤j≤M)表示t时刻状态是Si且观测值是时刻状态是Oj的概率；
初始状态矩阵πi＝P(y1＝Si)，1≤i≤N
其中πi表示在初始t＝1时刻状态是Si的概率。
3.根据权利要求1所述的基于隐马尔可夫模型的消费者线下消费行为预测方法，其特征在于：步骤三中，隐马尔可夫模型中的变量分为2组；第一组是状态变量Y＝{y1,y2,…,yn}，其中yi表示第i时刻的系统状态；假定状态是隐藏的、不可被观测的；第二组是观测值变量X＝{x1,x2,…,xn}，其中xi表示第i时刻的观测值；在隐马尔可夫模型中，系统在多个状态S＝{s1,s2,…,sN}之间转换，因此状态变量yi的取值范围为N个取值的离散空间。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710728972.6A CN107563801A (zh) | 2017-08-23 | 2017-08-23 | 一种基于隐马尔可夫模型的消费者线下消费行为预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710728972.6A CN107563801A (zh) | 2017-08-23 | 2017-08-23 | 一种基于隐马尔可夫模型的消费者线下消费行为预测方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107563801A true CN107563801A (zh) | 2018-01-09 |
Family
ID=60975627
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710728972.6A Pending CN107563801A (zh) | 2017-08-23 | 2017-08-23 | 一种基于隐马尔可夫模型的消费者线下消费行为预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107563801A (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108573399A (zh) * | 2018-02-28 | 2018-09-25 | 中国银联股份有限公司 | 基于转移概率网络的商户推荐方法及其系统 |
| CN108681923A (zh) * | 2018-05-16 | 2018-10-19 | 浙江大学城市学院 | 一种基于改进型隐马尔可夫模型的消费者消费行为预测方法 |
| CN109086306A (zh) * | 2018-06-22 | 2018-12-25 | 国网浙江省电力有限公司 | 基于混合隐马尔可夫模型的原子事件标签的提取方法 |
| CN109606288A (zh) * | 2018-12-17 | 2019-04-12 | 长安大学 | 一种公交车车厢智能监测灭火系统及控制方法 |
| CN109685564A (zh) * | 2018-12-21 | 2019-04-26 | 苏州易泰勒电子科技有限公司 | 一种基于项目反应理论的线下消费者决策行为预测方法 |
| CN110727593A (zh) * | 2019-10-12 | 2020-01-24 | 四川新网银行股份有限公司 | 基于历史需求数据挖掘的测试用例智能生成方法及介质 |
| CN112069045A (zh) * | 2020-08-14 | 2020-12-11 | 西安理工大学 | 基于隐马尔可夫的云平台软件性能预测方法 |
| CN113065431A (zh) * | 2021-03-22 | 2021-07-02 | 浙江理工大学 | 一种基于隐马尔可夫模型和循环神经网络的人体违规行为预测方法 |
| CN115481812A (zh) * | 2022-09-30 | 2022-12-16 | 广州天长信息技术有限公司 | 一种基于隐马尔科夫模型的绿通车混装预测方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070004385A1 (en) * | 2005-06-29 | 2007-01-04 | Microsoft Corporation | Principals and methods for balancing the timeliness of communications and information delivery with the expected cost of interruption via deferral policies |
| CN105930921A (zh) * | 2016-04-13 | 2016-09-07 | 石河子大学 | 一种棉叶螨的预测方法 |
-
2017
- 2017-08-23 CN CN201710728972.6A patent/CN107563801A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070004385A1 (en) * | 2005-06-29 | 2007-01-04 | Microsoft Corporation | Principals and methods for balancing the timeliness of communications and information delivery with the expected cost of interruption via deferral policies |
| CN105930921A (zh) * | 2016-04-13 | 2016-09-07 | 石河子大学 | 一种棉叶螨的预测方法 |
Non-Patent Citations (1)
| Title |
|---|
| 宋涛等: "基于二次聚类和隐马尔可夫链的持卡消费行为预测", 《计算机应用》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108573399A (zh) * | 2018-02-28 | 2018-09-25 | 中国银联股份有限公司 | 基于转移概率网络的商户推荐方法及其系统 |
| CN108681923A (zh) * | 2018-05-16 | 2018-10-19 | 浙江大学城市学院 | 一种基于改进型隐马尔可夫模型的消费者消费行为预测方法 |
| CN109086306A (zh) * | 2018-06-22 | 2018-12-25 | 国网浙江省电力有限公司 | 基于混合隐马尔可夫模型的原子事件标签的提取方法 |
| CN109606288A (zh) * | 2018-12-17 | 2019-04-12 | 长安大学 | 一种公交车车厢智能监测灭火系统及控制方法 |
| CN109685564A (zh) * | 2018-12-21 | 2019-04-26 | 苏州易泰勒电子科技有限公司 | 一种基于项目反应理论的线下消费者决策行为预测方法 |
| CN110727593A (zh) * | 2019-10-12 | 2020-01-24 | 四川新网银行股份有限公司 | 基于历史需求数据挖掘的测试用例智能生成方法及介质 |
| CN112069045A (zh) * | 2020-08-14 | 2020-12-11 | 西安理工大学 | 基于隐马尔可夫的云平台软件性能预测方法 |
| CN113065431A (zh) * | 2021-03-22 | 2021-07-02 | 浙江理工大学 | 一种基于隐马尔可夫模型和循环神经网络的人体违规行为预测方法 |
| CN115481812A (zh) * | 2022-09-30 | 2022-12-16 | 广州天长信息技术有限公司 | 一种基于隐马尔科夫模型的绿通车混装预测方法 |
| CN115481812B (zh) * | 2022-09-30 | 2023-09-08 | 广州天长信息技术有限公司 | 一种基于隐马尔科夫模型的绿通车混装预测方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107563801A (zh) | 一种基于隐马尔可夫模型的消费者线下消费行为预测方法 | |
| Cai et al. | Solving nonlinear water management models using a combined genetic algorithm and linear programming approach | |
| Mongillo et al. | Online learning with hidden Markov models | |
| CN101305509B (zh) | 用于确定估计电池状态向量的系统和方法 | |
| Nannen et al. | A method for parameter calibration and relevance estimation in evolutionary algorithms | |
| CN108681923A (zh) | 一种基于改进型隐马尔可夫模型的消费者消费行为预测方法 | |
| US20150254554A1 (en) | Information processing device and learning method | |
| CN108090788B (zh) | 基于时间信息集成模型的广告转化率预估方法 | |
| CN105471631A (zh) | 基于流量趋势的网络流量预测方法 | |
| Wang et al. | Cell selection with deep reinforcement learning in sparse mobile crowdsensing | |
| CN105760952A (zh) | 一种基于卡尔曼滤波与自适应模糊神经网络的负荷预测方法 | |
| Ma et al. | A novel kernel regularized nonlinear GMC (1, n) model and its application | |
| Tang et al. | Bayesian augmented Lagrangian algorithm for system identification | |
| Branisavljevic et al. | Fuzzy approach in the uncertainty analysis of the water distribution network of Becej | |
| Al‐Shehri | Artificial neural network for forecasting residential electrical energy | |
| Guorong et al. | Traffic flow forecasting based on PCA and wavelet neural network | |
| Baur et al. | Predicting high-dimensional heterogeneous time series employing generalized local states | |
| CN113452473B (zh) | 在反向散射中预测最佳信道的方法、装置及存储介质 | |
| Ghomi et al. | A new analytical algorithm and generation of Gaussian quadrature formula for stochastic network | |
| CN108053266A (zh) | 一种专利价值预估方法以及装置 | |
| CN107330726A (zh) | 基于自动拟合均值函数的泊松过程用户到店行为预测方法 | |
| Cui | On asymptotics of t-type regression estimation in multiple linear model | |
| Yu et al. | A combination forecasting model of extreme learning machine based on genetic algorithm optimization | |
| Nazemi et al. | Fuzzy-stochastic linear programming in water resources engineering | |
| Liu et al. | Echo state networks with double-reservoir for time-series prediction |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180109 |
|
| RJ01 | Rejection of invention patent application after publication |