[发明专利]Ising图模型的区间传播推理方法

权利要求书

1.一种Ising图模型的区间传播推理方法，通过定义Ising均值场计算树，并在计算树上实现期望区间传播推理过程，其特征在于，包括下列步骤：

(1)定义Ising均值场计算树：在Ising均值场推理下，变量簇c_γ的计算树模型为一四元组T(D_γ，R，M，Q)，其中：

1)D_γ：以c_γ为根节点的变量簇节点集D_γ＝{c_γ}∪Ch(c_γ)∪Ch(Ch(c_γ))∪…，其中，Ch(c_γ)表示c_γ的子节点集，Ch(c_γ)＝{c_β|x_i∈c_γ，x_j∈c_β，(i，j)∈E，γ≠β}，Ch(Ch(c_γ))表示变量集Ch(c_γ)的子节点集合：Ch(Ch(cγ))=∪cβ∈Ch(cγ)Ch(cβ),]]>

2)R：关系集R＝{<c_α，c_β>|c_α∈D_γ，c_β∈Ch(c_α)}，其中，关系<c_α，c_β>表示c_β是c_α的子节点，

3)M：计算树上的消息自底向上单向传播，记Mα-out={μi|i∈Vcα}]]>表示c_α输出消息集，

Mα-in=∪cβ∈Ch(cα)Mβ-out]]>表示c_α输入消息集，则M＝{M_α-out|c_α∈D_γ}，

4)Q：概率分布集Q＝{q_α(c_α；θ′_α)|c_α∈D_γ}，且θijxixj∝exp{Σi∈Vcαθ′ixi+Σ(i,j)∈Ecαθijxixj},]]>

其中，θ′i=θi+Σμt∈Mα-inθitμt,]]>

(2)定义Ising均值场截枝计算树：Ising均值场推理下，变量簇c_γ的剪枝计算树模型是一四元组T_c(D_γ，R，M，Q)，其中：

1)D_γ：以c_γ为根节点的变量簇节点集D_γ＝{c_γ}∪CCh(c_γ)∪CCh(CCh(c_γ))∪…，其中，CCh(c_γ)表示c_γ的子节点集，CCh(c_γ)＝Ch(c_γ)/An(c_γ)，CCh(CCh(c_γ))表示变量集CCh(c_γ)的截枝子节点集：CCh(CCh(cγ))=∪cβ∈CCh(cγ)CCh(cβ),]]>

2)R：关系集R＝{<c_α，c_β>|c_α∈D_γ，c_β∈CCh(c_α)}，其中，关系<c_α，c_β>表示c_β是c_α的子节点，

3)M：计算树上的消息自底向上单向传播，记Mα-out={μi|i∈Vcα}]]>表示c_α输出消息集，

Mα-in=∪cβ∈Ch(cα)Mβ-out]]>表示c_α输入消息集，则M＝{M_α-out|c_α∈D_γ}，

4)Q：概率分布集Q＝{q_α(c_α；θ′_α)|c_α∈D_γ}，且θijxixj∝exp{Σi∈Vcαθ′ixi+Σ(i,j)∈Ecαθijxixj},]]>其中，θ′i=θi+Σμt∈Mα-inθitμt;]]>

(3)在截枝计算树上，根据c_α的输入消息区间计算该变量簇的概率分布区间，令M_α-in^int.表示变量簇c_α输入消息区间的集合，即Mα-inint.={[μtl,μtu]|t∈Vcβ,]]>c_β∈Ch(c_α)}，变量簇c_α概率分布的参数区间为

θi′l=min{θi+Σμt∈Mα-inθitμt|μtl≤μt≤μtu},]]>

θi′u=max{θi+Σμt∈Mα-inθitμt|μtl≤μt≤μtu},]]>

θα′l={θi′l,θij|i,j∈Vcα,(i,j)∈Ecα},]]>

θα′u={θi′u,θij|i,j∈Vcα,(i,j)∈Ecα}.]]>

令A＝{a₁，a₂，…}，B＝{b₁，b₂，…}表示元素一一对应的等势集合，且A≤B＝{a_i≤b_i|i＝1，2，…}，则变量簇节点c_α的概率分布区间为

{qα(cα;θ′α)|θα′l≤θ′α≤θα′u};]]>

(4)基于变量簇的概率分布区间，利用和积算法计算变量簇c_α的输出消息区间：在概率分布q_α(c_α；θ′_α)上运行和积算法计算变量x_i∈c_α的期望μ_i，即μ_i＝Sum-Prod(q_α(c_α；θ′_α))，当给定参数θ′_i取值区间时，μ_i的极值取在参数区间的端点处，即：

μil=min{Sum-Prod(qα(cα;θ′α))|(θ′1,θ′2,···)∈{θ1′l,θ1′u}×{θ2′l,θ2′u}×···},]]>

μiu=max{Sum-Prod(qα(cα;θ′α))|(θ′1,θ′2,···)∈{θ1′l,θ1′u}×{θ2′l,θ2′u}×···},]]>

则变量簇节点c_α输出消息区间集合为Mα-outint.={[μil,μiu]i∈Vcα},]]>计算时，根据底层变量簇节点输入消息的取值区间，自底向上逐层进行消息区间传播计算出根变量簇变量期望区间。