让我们指出，之前的工作（例如，[23]和[18]）表明，这种依赖性实际上是依赖于由i（t）=vb（t）定义的队列大小的不平衡性- va（t）vb（t）+va（t）=qb（t）- qa（t）qb（t）+qa（t）（21），其中va/b（t）表示在时间t的最佳询价/投标价格下的可用量，我们假设订单具有与之前定义的AES相对应的恒定量。这种不平衡是解释即时购买压力的最简单代表。在这方面，在图9和图10中，我们以对数标度绘制了每个订单类型的参数f`（qia，qjb），作为平衡I的函数，计算为与分位数（qia，qjb）相关的状态区间的中值不平衡。让我们注意到，在这些图中，点大小表示相应qib的大小。通过查看图9和图10，我们可以得出以下观察结果：首先，订单上的不平衡变化捕获了密集参数f`（qa，qb）的大部分变化，这在大型蜱虫资产（Bund）上更为明显，因为f可以跨越几乎三个数量级（对于P和M事件），因为I的范围为-1至+1。小勾号资产（DAX）的f`（qa，qb）变化较小，因此不平衡的影响不太明显，尽管P和L事件仍然可见，但程度较低。这与观察结果一致，即不平衡是Largettick资产中价变化的一个很好的预测因子，而它的预测能力对小刻度资产的预测能力则不那么显著（参见例[12]）。

在分析这些数据时需要牢记的一个方面是，特别是对于一个小的tick资产，在本书的更深层次上也有大量的信息，这里没有考虑到这些信息。在与外滩结果相对应的图9中，我们观察到，对于中间价格变化（P事件），似乎存在一种阈值效应，即对于低于I=-1/2，向上的价格增量被显著抑制，而较大的正失衡只会略微增加向上价格变化的可能性。这是很自然的，因为不平衡小于-1/2对应于这样的情况，即任务大小至少比出价大小大三倍，这使得中间价不太可能向上移动。更令人惊讶的是，代理似乎强烈地根据订单簿的状态来决定使用攻击性订单（M事件）。也许他们急于在向中价移动之前获得最后可用的流动性（事实上，当不平衡较大时，f对于M事件变得非常大，即当ask大小几乎没有任何剩余时）。限制和取消订单似乎对该书的状态不那么敏感。正如我们已经指出的那样，对于DAX，对不平衡的依赖性不太明显（见图10），除中间价格变化外，所有因子f`（qa，qb）都随I而变弱。但是，请注意，对于市场（M）和取消（C）订单，我们观察到一个制度转换，即I=1/2。这些机制分别对应于非常大和非常小的ask队列大小。

在后一种情况下（当价差等于一个滴答声时），不太可能出现不改变中间价的市场和取消订单。最后但并非最不重要的一点是，让我们注意到，不平衡未捕获到的强度变化主要位于I=0附近，其中各个队列大小（几乎相等的qa或qb）似乎具有重要影响。因此，例如，对于外滩，fP+（q，q）似乎随着q的增加而增加。实际上，这可以通过更深入的分析来解释，该分析表明，平均价差随着q的增加而增加，当q较大时，它非常接近2。对于市场订单也可以这样说。或者，根据等式（19），在[14]1.0 0.5 0.0 0.5 1.0不平衡21012log（f）=P+1.0 0.5 0.0 0.5 1.0不平衡21012log（f）=La1.0 0.5 0.0 0.5 1.0不平衡21012log（f）=Ca1.0 0.5 0.0 0 0.5 1.0不平衡21012log（f）=图9：从左到右，从上到下，l=P+，La，Ca的log（fl Maof QRHmodel作为不平衡21的函数，外滩未来。bid侧和ask侧的分位数相同，分别对应于qa=qb=]0，80]，qa=qb=]80，165]，qa=qb=]165，258]，qa=qb=]258，386]和qa=qb=]386+∞[.由于该示例说明了其他微观结构变量，尤其是扩散，因此应在更完整的模型中加以考虑。然而，这超出了当前工作的范围。估计强度和经验强度的比较。为了进一步验证QRH-II模型，我们测试了重现经验强度的能力。平均强度的MLE由ate q=（qa，qb）15次写入（在非参数框架中）^∧`（q）=^∧（q）Pt'kq（tl-k） =qPtkq（t-k） =qwith^∧（q）=平均值（tk- tk公司-1 | q（t-k） =q）-1，（22）其中运算符表示（…）对应于经验平均值。这个估计仅仅基于对过程{Nlt}l的观察。

可以使用∧'''QRH-II（t`）的QRH-II参数形式计算相应的估计量'∧''''QRH-II（q），-k | q），这导致∧` QRH（q）=平均值（λ` QRH-II（t`），-k） | q（t`，-k） =q）（23）为了综合所获得的结果，我们选择不对所有类型的阶和所有状态q的∧` QRH-II（q）与∧`（q）进行比较。相反，对于每种类型的阶，我们报告加权相对误差`定义为：`=Pq^∧`（q）-^∧` QRH-II（q）N`（q）Pq∧`（q）N`（q）（24）Bund和DAX的加权相对误差如表6所示。我们观察到`约为10%，这与经验观察到的强度相符。1.0 0.5 0.0 0.5 1.0不平衡1.00.50.00.51.0log（f）=P+1.0 0.5 0.0 0.5 1.0不平衡1.00.50.00.51.0log（f）=La1.0 0.5 0.0 0 0.5 1.0不平衡1.00.50.00.51.0log（f）=Ca1.0 0 0 0.5 0 0 0.5 1.0不平衡1.00.50.00.51.0log（f）=图10：从左到右，从上到下l=P+、La、Ca、Maof QRH模型的对数（fl（qia、qjb）），作为不平衡21、DAX指数期货的函数。bid侧和ask侧的分位数相同，分别对应于qa=qb=]0,2]，qa=qb=]2,3]，qa=qb=]3,5]，qa=qb=]5,8]和qa=qb=[8+∞[.P+P-LaLbCaCbMaMbBund 14.2%10.5%6.7%6.0%7.4%8.1%4.0%12.0%DAX 8.2%5.9%0.5%4.7%7.1%1.1%1.6%5.9%表6：核范数矩阵顺序类型分析的平均强度误差。最后，为了完成我们结果的分析，在图11中，我们显示了Bund和DAX的估计范数{Rφ\'m（t）dt}\'mf的矩阵。当忽略队列依赖性时，这些矩阵提供了不同事件类型之间的平均交互信息。因此，它们是[4]中图4所示的核范数矩阵的对应项。

我们恢复了[4]中强调的许多功能，例如限制、市场和取消订单的强大对角线组件（订单拆分的签名），以及市场订单和价格变动似乎对限制和取消的影响要远远大于相反的事实。值得注意的是，对于外滩而言，当中间价上升（或下降）时，它会降低（或上升）ask（或最佳）侧的限价订单率，因为有效价格接近最佳ask，所以相对于市场而言，发送限价没有任何收益。这也解释了为什么它会增加（或减少）最佳询问（或出价）端的取消订单发送率。在DAX上或多或少可以看到同样的影响（但减弱了很多）。由于DAX上的滴答声很小，所以它会减弱，因此有效的价格论证并不强烈。让我们指出，完全相同的论点可以用来解释市场订单对限制和取消订单的影响。市场秩序对价格变化的影响主要是因为在我们的设置下，市场秩序以最佳价格消耗流动性，最终可能产生新的价格。由于DAX是一种小型的交易资产，而最高价的排队规模比外滩小，因此市场订单更有可能产生新的价格。因此，市场秩序对价格变化的影响更为明显。我们可以看到，对于DAX未来，出于同样的原因，限价订单和取消订单也会对价格变化产生强烈影响。我们还可以看到，除了订单分割带来的自我激励效应外，对于外滩而言，一方的限价订单与另一方的取消订单相结合。这可以看作是一种简单的做市商策略（重新平衡头寸）。

ItP+PMAMBLALBCABP+PMAMBLALBCAB0.90.60.30.00.30.60.9P+PMAMBLALBCABP+PMAMBLACAB0.90.60.30.00.30.60.9图11：使用最小二乘估计QRH-II模型的估计矩阵normsRφlm（t）dt。Bundfuture位于左侧，DAX future位于右侧。1.0 0.5 0.0 0.5 1.0I0.00.10.20.30.40.5密度1.0 0.5 0.0 0.5 1.0I0.000.050.100.150.20密度图12：外滩（左侧面板）和DAX（右侧面板）P+事件后观察到的不平衡值的经验频率。请注意两种情况下负不平衡值的大分量。当然不会出现在DAX上，因为刻度大小要小得多，相同的再平衡策略不会对双方都产生必要的影响。关于这些特征的解释，我们请读者参考第4节【4】。平均反转率bund 0.65DAX 0.56表7：平均价格反转的度量：两个连续的中间价格变化事件具有相反方向的经验概率。虽然对于大多数特征，QRH-II模型（图11）和纯霍克斯模型（图4中的图4）是相似的，但我们可以观察到在P→ P和P→ T子矩阵。我们无法在QRH-II中发现P+到P+之间没有强激发-和P-P+（图11中左上角的2×2子矩阵），这应该是高频价格均值回归的特征（如[4]所述）。对于纯霍克斯过程模型，如工作[4]，价格的平均回归反映在P→ P核范数子矩阵（即强P+→ P-和P-→ P+条款）。这是因为（平均而言）P+事件会产生更多的P-事件多于P+事件，反之亦然。如表7所示，实际的中间价序列是强线性回归，因此，在QRH-II模型中，这一特征可能应该用队列反应函数f（qa，qb）来解释。

事实上，如图12所示，在中间价格变化之后，不平衡发展为有利于价格朝相反方向移动。例如，价格上涨P+在大多数情况下会导致负不平衡，要么是因为百思买队列QA满了，要么仅仅是因为在价差内发送了单笔限购订单。根据图的结果。9和10，在这种情况下，我们有fP+（qa，qb） 1和fP-（qa，qb）>1。相反，价格下跌后，P-, 我们会有fP-（qa、qb） 1和fP+（qa，qb）>1。2×2子矩阵P→ 在纯霍克斯模型中，P估计值很可能对应于QRH-II霍克斯子矩阵乘以反对角线上的一个大因子和对角线上的一个小因子。基于不平衡影响的同一种论证可以用来解释P的高对角线值→ T子矩阵在反对角线上观察到了[4]中的最高值。4总结与展望在本文中，我们介绍了两种“队列反应性霍克斯模型”（QRH-I和QRH-II），旨在分别改进Huang等人[16]关于劳动力动力学的队列反应性和Bacry等人[4]的模型的方法。我们表明，shuch模型可以通过参数化方法轻松校准。我们对欧洲期货交易所（Eurex）的两种不同未来资产（namelyBund和DAX）订单簿数据的经验发现表明，队列反应性和过去订单流量依赖性都与未来订单簿事件发生的可能性有关。事实上，这两种模型都在优度方面表现为纯霍克斯模型和纯队列反应模型。就QRH-I模型而言，我们的框架允许人们保持在马尔可夫过程的框架内，马尔可夫过程具有遍历性，因此，与Huang et al。

方法中，我们可以定义并估计与模型不变度量相关的队列大小分布。QRH-IImodel引导我们通过考虑最佳出价和最佳ASK队列的状态来确定Bacry等人的发现【4】。事实证明，容量不平衡可以解释大多数队列依赖性。QRH-II模型可以通过考虑队列大小（以及不平衡）和订单流量之间的相互作用来改进，就像在QRH-I模型中一样，还可以包括对小tick资产利差的明确依赖关系。我们所做的一些实质性简化也可以删除，以获得更现实的模型，例如，通过采用与[19]类似的方法放弃单一订单大小的假设。除了考虑这些模型的各种应用来设计和优化高频交易策略外，从数值的角度来看，仍需确定的方法允许处理等式（2）中定义的完整QR模型，其中不仅外生核和霍克斯核的队列依赖性是任意的，而且可以同时考虑到给定深度的所有订单队列。也许最近在统计学习中发展的一些方法和技术将有助于处理此类高维问题。从数学的角度来看，对于依赖队列的Hawkes模型的稳定性和平稳性条件还有待进一步的理解。

更根本的是，对于观察到的外生强度的形状以及订单流量到达率在不同市场参与者（理性）行为方面的不平衡依赖性的明确理解仍然是一个悬而未决的问题。附录A。1不变分布存在性的证明我们有动机证明第一部分中提出的QRH-I模型的不变分布的存在性，因为如果存在这种不变分布，我们可以通过长期模拟QRH-I模型来近似队列大小的经验分布。利用lyapunov函数进行了证明。首先，让我们定义：o\'mu（t）=Ztα\'mue-βu（t-s） dNms（25）通过采用定义（23）和（6），强度函数现在的形式为：λ`（t）=u`（q（t-)) +DXm=1UXu=1o\'mu（t），（26）到总系数q（t）6=0，对于使用队列大小的订单，必须考虑到这一点。我们注意到J增加队列大小的订单类型集，I减少队列大小的订单类型集。我们进一步假设P`∈Iu`（q）≥ c`q和um（q）≤ c*, m级∈ J、 队列大小由订单流量之和决定：q（t）=Xm∈JNm公司-X个`∈INm（27）设~o（t）为所有（`，m，u）的o\'mu（t）分量垂直叠加得到的向量∈{1，…，D}×{1，…，U}。要验证向量过程（q，~ o）是否是马尔可夫过程并不困难。然后，我们的目标是为（q，~ o）T构造一个Lyapunov函数。

首先，我们分别为q和~o构造Lyapunov函数，然后将它们组合在一起。A、 1.1 oLet us的Lyapunov函数首先写出o（t）的微分形式：do\'mu（t）=-βuo ` mu（t）dt+α` mudNmt（28）然后对于任何映射R2D+Uto R的任意合适函数F，微型生成器的形式为：LF（~ o）=Xmλm（F（~ o+m（~ o））- F（~ o））-X`，m，uβuo`muFo` mu（29），其中λmis事件到达点过程N的维数m的概率，以及m（~o）是由该事件引起的~o跳跃。然后，我们将矩阵A定义为：A ` m=Xuα` muβu（30）。我们假设QRH-I模型的Hawkes部分描述的序流动力学对应于稳定的Hawkes过程。根据Perron-Frobenius定理，Asatis fies的最大特征值κ为0<κ<1，以及~ κsaties的相关特征向量l、，l> 0。然后我们注意到δ\'mu：=δ\'u=`βu（31）我们选择函数Vof~o为：V（~o）=X\'，m，uδ\'muo\'mu（32），使用上面定义的符号，我们可以验证lv（~o）=Xm（um+Xp，qompq）`=0∨q（t）>0X\'，uδ\'muα\'mu-X`，m，uβuo`muδ`mu≤Xm（um+Xp，qompq）X\'，uδ\'muα\'mu-X`，m，uβuo` muδ\'mu=X`，m，uumδ\'muα\'mu+Xm（Xp，qompq）（X`，u`βuα\'mu）-X`，m，uβuo` muδ\'mu=C+Xm（Xp，qompq）（X``Xuα\'muβu）-X`，m，uβuo` muδ\'mu=C+Xm（Xp，qompq）mκ-X`，m，u`o\'mu=C- (1 - κ） X`，m，u`o\'mu=C- (1 - κ） X`，m，uβuδ\'muo\'mu≤ -ρV+C（33），其中常数ρ选择为ρ=（1- κ） infβu（34）A.1.2 Lyapunov函数对于q类似，让我们首先写出q（t）的微分形式：dq（t）=Xm∈JdNmt公司-X个`∈IdN\'t（35）我们一直使用相同的符号 如等式29所示。

对于将R+映射到R的任意合适函数F，最小生成元的形式为：LF（q）=Xmλm（F（q+m（q））- F（q））（36）接下来，我们选择q:V（q）：=q（37）的恒等式函数Vas，那么可以很容易地验证V（q）是q:LV（q）=Xm的Lyapunov函数∈Jλm（q）-X个`∈Iλ`（q）≤Xm公司∈Jum（q）-X个`∈Iu`（q）+X`，m，uo`mu≤Xm公司∈Jc公司*-X个`∈Ic\'q+X\'，m，uo\'mu≤ -X个`∈集成电路q+C≤ -ρV+C（38）提醒常数C*, 之前定义了c。常数ρ简单地选择为ρ=max\'c\'。此外，必须注意的是，此处c：=Xm∈Jc公司*+X`，m，uo`mu（39），它依赖于~ o。另一个注释是因为m和c*如果固定，则必须存在ρ*谁满意c<ρ*X`，m，uo` mu（40）A.1.3（q，~ o）的Lyapunov函数最后一步是为q和~ o建立一个函数V：V（q，~ o）=V（q）+ηV（~ o）（41），然后我们应用先前对Vand V得出的结论。直接计算表明：LV（q，~ o）=LV（q）+ηLV o（~ o）=-ρV+C-ρηV+Cη≤ -ρV+ρ*X`，m，uo`mu-ρηV+C（42）注意，由于系数δ\'muin和o\'mu均为正，因此必须存在一个满足ρ2ηV（~o）>ρ的η*X`，m，uo`mu（43）通过将该不等式代入等式42，我们最终证明V是（q，~o）：LV（q，~o）的李雅普诺夫函数≤ -ρV-ρ2ηV+C≤ - min（ρ，ρ2η）V+C（44）在假设A的谱半径小于1的情况下，假设Lyapunov函数的存在和上述几何漂移条件，文献[20]中的定理7.1保证过程qis是遍历的。此外，它以指数速度收敛到其唯一的平稳分布。A、 2 QRH模型对数似然函数及其梯度的计算对于QRH模型，对数似然是u和α的函数。让我们注意t“k类型的第k个事件的时间戳”，N“类型的事件总数”。

用这种表示法，L（~α，~u）=DX`=1NlXk=1logu’（q（t’，-k） +DXm=1UXu=1α\'muβuZtke-βu（t-s） dNms公司-DX`=1ZTu`（q（t））+DXm=1UXu=1α\'muβuZse-βu（s-v） dNmvds（45）我们首先定义g和g，其值不取决于参数u和α。有趣的是，G和G只需要计算一次，然后就可以重用它们来加速G似然函数及其梯度的计算。gmu（t）=Xtmk<tβue-βu（t-tmk）（46）和gmu（t）=Ztgmu（s）ds=ZtZsβue-βu（s-v） dNmvds（47）可以使用以下递推关系计算g和g：计算ggmu（tk）=Xtmk<tkβue-βu（tk）-tmk）（48）=Xtmk<tk-1βue-βu（tk-1.-tmk）e-βu（tk-tk公司-1） +Xtk-1.≤tmk<tkβue-βu（tk-tmk）（49）=e-βu（tk-tk公司-1） 总经理（tk-1） +βue-βu（tk-tk公司-1） 类型（t+k-1） =m（50）计算GGmu（tk）- Gmu（tk-1） =Ztktk-1gmu（s）ds（51）=1- e-βu（tk-tk公司-1） βugmu（tk-1) +1.- e-βu（tk-tk公司-1)类型（t+k-1） =m（52）我们稍后将看到，对于每个tk，我们只需要g和g at的值。A、 2.1对数似然函数通过利用g和g的定义，对数似然函数可以用以下方式重写：L（α，u）=DX`=1N` Xk=1logu’（q（t’，-k） ）+DXm=1UXu=1α\'mugmu（tk）-DX`=1NXk=1u`（q（t-k） ）（塔卡- tk公司-1) -DX`=1u\'q（t+N）（t- tN）-DX`=1NXk=1DXm=1UXu=1α`μ（Gmu（tk））- Gmu（tk-1)) -DX`=1DXm=1UXu=1α`μ（Gmu（T））- Gmu（tN））（53）A.2.2梯度稍微滥用符号q和q（t），直接计算表明：Lu`（q）=N`Xk=1q（t`，-k） =qu`（q（t`），-k） ）+PDm=1PUu=1α\'mugmu（t\'k）-NXk=1q（t-k） =q（tk- tk公司-1) - 1q（t+N）=q（t- tN）（54）Lα\'mu=N\'Xk=1gmu（t\'k）u\'（q（t`），-k） ）+PDm=1PUu=1α\'mugmu（t\'k）-NXk=1（Gmu（tk）-Gmu（tk-1))-（Gmu（T）-Gmu（tN））（55）A.3计算QRH模型的最小二乘函数及其梯度A。3.1第二部分所述QRH模型的最小二乘函数。

让我们通过组合ask端和bide端的状态来记录LOB的状态q。q=qa×qbf（q（t））：=f（qa（t-), qb（t-)) （56）那么在时间t，维数的强度函数为：λ\'t=f\'（q（t））u\'+XmXuα\'muZtβue-βu（t-s） dNms公司（57）与对数似然函数及其梯度的计算类似，我们首先确定一些中间值，其值不依赖于u、α或f。它们只需在预处理阶段计算一次。然后，在计算最小二乘函数及其梯度时可以重用它们。g\'u（t）=Xt\'k<tβue-βu（t-t`k）（58）Gmu（q）=ZTgmu（s）1type（q（t））=qdt（59）Hmmuu（q）=ZTYPE（q（t））=qgmu（t）gm0u（t）dt（60）D（q）=ZTYPE（q（t））=qdt（61）Cm（q）=NmXk=1type（tm，-k） =q（62）计算ggmu（tk）=Xtmk<tkβue-βu（tk-tmk）（63）=Xtmk<tk-1βue-βu（tk-1.-tmk）e-βu（tk-tk公司-1） +Xtk-1.≤tmk<tkβue-βu（tk-tmk）（64）=e-βu（tk-tk公司-1） 总经理（tk-1） +βue-βu（tk-tk公司-1） 类型（t-k） =m（65）计算GGmu（tk）- Gmu（tk-1） =Ztktk-1gmu（s）ds（66）=1- e-βu（tk-tk公司-1） βugmu（tk-1) +1.- e-βu（tk-tk公司-1)类型（t+k-1） =m（67）计算HHmmuu（q）=XkZtktk-1类型（q（t））=qgmu（t）gm0u（t）dt（68）=Xktype（q（t-k） ）=qgmu（tk-1） gm0u（tk-1） e类-（βu+βu）（tk-tk公司-1） （69）A.3.2最小二乘函数对于前一部分中定义的中间变量，我们可以将20中定义的最小二乘函数改写为更有效的计算形式。在尺寸`，R`（θ）=ZTλ`（t；θ| Ft）dt中-N`Xk=1λ`（t`k；θ| Ft`k）=ZTf`（q（t））u\'+XmXuα\'muZtβue-βu（t-s） dNms公司dt公司-N\'Xk=1f\'（q（t\'k））u\'+XmXuα\'muZt\'kβue-βu（tmk-s） dNms公司（70）让我们分别命名上述第I项和第II项公式中的第一项和第二项。

我们可以进一步将项I分解为新项项项I.1、项I.2等：项I=ZTf`（q（t））u`+Xmα\'muXuZtβue-βu（t-s） dNms公司dt=ZTf`（q（t））u`dt+ZT2f`（q（t））u`XmXuα\'muZtβue-βu（t-s） dNms公司dt+ZT2f`（q（t））XmXuα\'muuZtβue-βu（t-s） dNms公司dt（71）对于第I.1项，可根据上述中间变量计算：对于第I.2项，ZTf`（q（t））u\'dt=u\'XqD（q）f`（q）（72），ZT2f`（q（t））u`XmXuα\'muZtβue-βu（t-s） dNms公司dt=ZT2f`（q（t））u`XmXuα\'mugmu（t）dt=2u\'Xqf\'（q）XmXuα\'muXuGmu（q）（73）对于术语I.3，ZT2f`（q（t））XmXuα\'muZtβue-βu（t-s） dNms公司dt=ZT2f`（q（t））XmXmXuXuα\'muα\'mugmu（t）gmu（t）dt=XmXmXuXuα\'muα\'muXqf`（q）Hmmuu（q）（74）为了便于记法，我们使用术语II的符号。然后我们把它分解成第二项。1和第二学期。2：II=N\'Xk=1f\'（q（t\'k））u\'+XmXuα\'muZt\'kβue-βu（tmk-s） dNms公司=N\'Xk=1f\'（q（t\'k））u\'+N\'Xk=1f\'（q（t\'k））XmXuα\'muZtmkβue-βu（tmk-s） dNms公司（75）第二期。1，对于第II项，N\'Xk=1f\'（q（t\'k））u`=u\'Xqf（q）C\'（q）（76）。2，N\'Xk=1f\'（q（t\'k））XmXuα\'muZtmkβue-βu（tmk-s） dNms公司=N\'Xk=1f\'（q（t\'k））XmXuα\'mugmu（t\'k）（77）20中定义的最小二乘函数可以通过将所有这些项相加来计算。A、 3.3梯度利用上述中间变量和结果，直接计算表明：Ru`=2u\'XqD（q）f`（q）+2Xqf`（q）Xmα\'muXuGmu（q）- 2Xqf（q）Cm（q）（78）Rα\'mu=2u\'Xqf\'（q）Gmu（q）+2XmXuα\'muXqf`（q）Hmmuu（q）- 2NmXk=1f`（q（tmk））gmu（tmk-)(79)Rf`（q）=Xq2f`（q）u\'D（q）+4u\'f`（q）Xmα\'muXuGmu（q）+ 2XmXmXuXuα\'muα\'muXqf`（q）Hmmuu（q）- u\'C\'（q）- 2N\'Xk=1类型（t\'k-)=qXmXuα\'mugmu（t\'k）（80）参考文献[1]F.Abergel和Jedidi A.基于Hawkes过程的极限指令簿的长期行为。SSRNElectronic杂志，2015年。[2] F.Abergel、M.Anane、A.Chakraborti、A.Jedidi和I.Muni Toke。限制订单簿。剑桥大学出版社，2016年。[3] 弗雷德里克·阿伯格尔和艾门·杰迪。订单建模的数学方法。16(5):1–40,2010.[4] E.Bacry、T.Jaisson和Muzy J。

缓慢递减Hawkes核的估计：高频订单动态的应用。《定量金融》，16（8）：1179–120120016。[5] E.Bacry、I.Mastromatteo和J.F.Muzy。Hawkes流程融资。《市场微观结构与流动性》，01（01）：15500052015。[6] Emmanuel Bacry和Jean-Fran,cois Muzy。霍克斯过程的一阶和二阶统计特征和非参数估计。IEEE信息论学报，62（4）：2184–2202202202016。[7] 理查德·C·布拉德利。强混合条件的基本性质。调查和一些开放性问题。概率。《调查》，2:107–144，2005年。[8] Anirban Chakraborti、Ioane Muni Toke、Marco Patriarca和Fr'ed'eric Abergel。经济物理学评论：I.经验事实。数量。财务。，11(7):991–1012, 2011.[9] R.Cont、S.Stoikov和R.Talreja。订单动态的随机模型。运营研究，58:549–56320010。[10] D.J.Daley。点过程理论导论：基本理论和方法，第1卷。斯普林格，2008年。[11] 安德鲁·道和杰莫·彭德。队列霍克斯过程：短暂的自我兴奋。arXive印刷品，第arXiv页：1811.042821918年11月。[12] 马丁·古尔德和朱利叶斯·博纳特。在limitorder book中，队列不平衡是一个提前一点的价格预测因素。《市场微观结构与流动性》，2（02）：16500062016。[13] 马丁·古尔德、梅森·阿波特、斯泰西·威廉姆斯、马克·麦克唐纳、丹尼尔·杰芬和萨姆·道维森。限制订单簿。数量。财务。，13(11):1709–1742, 2013.[14] N.R.Hansen、P.Reynaud Bouret和V.Rivoirard。”多元点过程的lasso和概率不等式”。伯努利，21（1）：83–1432015。[15] A.G.霍克斯。一些自激和互激点过程的谱。Biometrika，58:83-901971年。[16] W.Huang、Lehalle C.A.和Rosenbaum M.模拟和分析订单数据：队列反应模型。

《美国统计协会杂志》，110（509）：107–1222015。[17] 艾门·杰迪和弗里德·埃里克·阿伯格尔。基于derbook模型的hawkes过程的稳定性和价格标度极限。SSRN电子期刊，第1-22页，2013年5月5日。[18] Charles Albert Lehalle、Othmane Mounjid和Mathieu Rosenbaum。基于流动性的最佳交易策略。2018年第1-39页。[19] 陆晓飞和埃里克·阿伯格尔。订单建模和做市策略。第1-232018页。[20] S.Meyn和R.L.Tweedie。马尔可夫过程的稳定性iii：连续时间过程的Foster-lyapunov准则。S、 应用概率的进展，25:518–548。[21]Maxime Morariu Patrichi和Mikko S.Pakkanen。状态相关的Hawkes过程及其在限制订单建模中的应用。arXiv电子打印，第arXiv页：1809.080602018年9月。[22]Christine A Palour和Duane J Seppi。限价订单市场：一项调查。《金融中介与银行手册》，2008年5:63–95。马塞洛·兰巴迪、埃曼纽尔·巴克里和法布里齐奥·利洛。订单量在订单动态中的作用：多变量霍克斯过程分析。《定量金融》，17（7）：999–1020，2017年。[24]马塞洛·兰巴迪、埃曼纽尔·巴克里和让·弗兰·科伊斯·穆齐。分离和量化市场参与者的波动性贡献。arXiv预印本arXiv:1807.070362018。[25]E.Smith、J.D.Farmer、L.Gillemot和S.Krishnamurthy。连续双重拍卖的统计理论。《定量金融》，6:481–5142003。