回顾：系统的能控性、能观性和稳定性及李雅普诺夫方法

能控性反映了输入信号对系统状态的影响和控制能力，能观性反映了输出（量测）对系统状态的识别能力，它们反映了系统（被控对象）本身的内在特性。

这两个概念是卡尔曼在20世纪60年代提出的，是现代控制理论中的两个基本概念。

能控性定义

由于能控性只涉及用外部输入来改变系统状态的问题，故只考虑系统的状态方程：
$\dot x = Ax + Bu$

其中， $x$ 是 $n$ 维的状态向量（点）， $u$ 是 $m$ 维控制输入向量（点）， $A$ 和 $B$ 分别是已知的 $n\times n$ 和 $n\times m$ 维实常数矩阵。

定义：对上述系统的一个状态 $x_0$ ，如果存在一个有限时刻 $T$ 和时间段 $[0,T]$ 上的控制信号 $u(t)$ ，使得在这样一个控制信号（control）作用下，系统状态从 $t=0$ 时刻的初始状态 $x_0$ ，转移到 $t=T$ 时刻的零状态，即 $x(T)=0$ ，则状态 $x_0$ 称为是能控的。若系统的所有状态都是能控的，则称系统是状态完全能控的，简称系统是能控的，有时也记矩阵对 $(A,B)$ 是能控的。

系统的能控性表明：若状态 $x_0$ 是能控的，则一定可以通过设计一个适当的控制律（control law），将系统在有限时间内从 $x_0$ 转移到零状态。

输出能控

在实际的控制系统设计中，我们需要控制的往往是输出量，而不是系统的状态。这种情况下，系统状态能控性对实现输出量的控制既不充分也不必要。

$\left\{ \begin{gathered} \dot x = Ax + Bu \hfill \\ y = Cx \hfill \\ \end{gathered} \right.$

输出能控性定义：对上述系统，若对任意的初始输出 $y_0$ （状态的一个线性组合），存在有限时刻 $T$ 和在时间段 $[0,T]$ 上定义的控制信号 $u(t)$ ，使得在该控制信号作用下，系统的输出从初始输出 $y_0$ 转移到任意给定的最终输出 $y(T)$ ，则系统称为是输出完全能控的，简称输出能控。

系统的能观性

在实际控制系统设计中，我们总希望利用描述系统全部动态行为的状态信息来构造反馈控制器，以使得闭环系统具有尽可能满意的性能。但在一个实际的系统中，并不是所有的状态信息都是直接可测量得到的，能够测量的只是系统的输出，因此如何解决这个矛盾？

状态空间模型中的输出方程建立起了系统的状态变量和输出量之间的关系，从而系统的输出信号中或多或少总包含有系统的状态的信息。那么，是否可以通过观测一段时间内的测量输出信号，或者再结合外加的输入信号（因为输出方程中输出有时也依赖输入信号）来确定出之前某个时刻系统的状态呢？这就是系统状态能否从外部观测或估计的问题，简称系统状态能观性问题。

在讨论能观性条件时，只需要考虑零输入系统：
$\left\{ \begin{gathered} \dot x = Ax \hfill \\ y = Cx \hfill \\ \end{gathered} \right.$

其中， $x$ 是 $n$ 维的系统状态向量， $y$ 是 $p$ 维的测量输出， $A$ 和 $C$ 分别是已知的 $n\times n$ 维和 $p\times n$ 维常数矩阵， $x_0$ 是 $t=0$ 时刻的初始状态向量。

之所以只考虑零输入系统，是因为能观性问题考虑的是用外部的已知信号（如输出信号，控制信号）来估计内部的未知状态。由系统运动分析结论可知，从系统状态空间模型可得系统的状态关于时间的响应在给定初始条件 $x_0$ 下：
$x(t) = {e^{At}}x(0) + \int_0^t {{e^{A(t - \tau )}}Bu(\tau )d\tau }$

从而，可得系统输出响应：
$y(t) = C{e^{At}}x(0) + C\int_0^t {{e^{A(t - \tau )}}Bu(\tau )d\tau }$

由于矩阵 $A,B,C$ 均已知， $u(t)$ 也已知（控制信号由我们设计），所以上式右端的积分项为已知，将它们移到等式的左边：
$y(t) - C\int_0^t {{e^{A(t - \tau )}}Bu(\tau )d\tau } = C{e^{At}}x(0)$

上式左边都是已知信号，而右边是带估计的状态 $x(0)$ 的线性组合。这和上述的零输入系统状态空间模型得到的输出 $y(t) = C{e^{At}}x(0)$ 没有本质区别，即通过左边已知信号来估计右边的未知状态 $x(0)$ 。因此研究系统的状态估计或观测问题只需考虑零输入状态空间模型即可。

定义：对上述系统，若以非零初始状态 $x_0$ 产生的输出响应恒为零，即对所有的时间 $t$ ：
$y(t) = Cx(t) = 0$

则称状态 $x_0$ 是不能观的。若系统中没有不能观的状态（换句话说所有状态都能观），则称系统是完全能观的，简称是能观的，有时也称矩阵对 $(C,A)$ 是能观的。

系统的输出恒为零表明自治系统在非零初始状态 $x_0$ 的激励下仍然是静止的，初始状态对系统输出响应没有任何影响，即在系统输出中不能反映状态 $x_0$ 的任何信息，根据定义这样的状态 $x_0$ 是不能观的。

系统的稳定性

在控制工程中，所设计的系统在受到扰动后，尽管系统会偏离原平衡工作点（稳定点），但在扰动消失后，设计者往往希望系统有能力自动回到并保持在原工作点附近，这就是系统稳定的基本含义。

稳定是一个控制系统能正常工作的基本要求，系统只有在稳定的前提下才能进一步探讨其他性能。因此，稳定性问题一直是控制理论中的一个最基本和最重要的问题，控制系统的稳定性分析是系统分析的首要任务。

1892年，俄国数学力学家李雅普诺夫（A.M. Lyapunov）在他的博士论文《运动稳定性的一般问题》中，提出了著名的李雅普诺夫稳定性理论。该理论作为稳定性判别的一般方法，适用于各类动态系统。李雅普诺夫稳定性理论的核心是提出了判别系统稳定性的两种方法，分别被称为李雅普诺夫第一方法和第二方法。

李雅普诺夫第一方法是通过求解系统的动态方程，然后根据解的性质来判断系统的稳定性，其基本思路和分析方法与古典控制理论是一致的。由于需要求出系统动态方程的解后才能判别系统的稳定性，故也称为判别稳定性的李雅普诺夫间接法。

李雅普诺夫第二方法则是一种定性方法，它无需求解复杂的系统微分方程，而是通过构造一个类似于能量函数的标量李雅普诺夫函数，然后再根据李雅普诺夫函数随时间变化的情况来直接判定系统的稳定性。因此，它特别适合于那些难以求解的非线性系统和时变系统。李雅普诺夫第二方法不仅可以用来分析系统的稳定性，而且还可用于对系统过渡过程特性的评价以及求解参数最优化等问题。李雅普诺夫第二方法最大的优点是它可用于控制系统的设计，从而使得该方法在自动控制的各个分支中都有广泛的应用，是控制理论中最重要的理论和方法之一。

李雅普诺夫意义下的稳定性

图中小球在没有任何外力作用下，它将保持在B点静止不动（稳态）。若给小球一个外力，使之移动到A点，然后让它做自由运动，则小球做震荡运动由于摩擦力的存在，最后在B点稳定下来并静止。

稳定性指的是系统在平衡状态下收到干扰后，系统自由运动的性质。上述描述中，小球最终又稳定在了原平衡点，则这样的系统是稳定的。若小球初始静止状态在D点，则当小球受到干扰后，小球不能再回到D点，这样的系统是不稳定的。

图1 小球运动示意图

在以上小球运动分析中有几个关键的概念。第一个就是平衡状态，图中B点和D点处的状态就是平衡状态，即小球处于静止状态。其次是扰动，小球在受到外部干扰后偏离平衡状态，然后在没有任何外力和扰动作用下做自由运动（自治系统的自由运动）。因此，小球受到的干扰只是初始干扰，而非持续干扰，这就是李雅普诺夫稳定性所处理干扰的特点，从而诸如持续风力干扰等（持续的输入干扰）就不在李雅普诺夫稳定性分析范围之内。最后，系统的稳定与否依赖于小球受干扰前所处的平衡位置，如小球在B点是稳定的，在D点是不稳定的。因此，系统的稳定与否和平衡状态相关，系统稳定性仅仅指的是在某个平衡状态处的稳定性（稳定性都是相对于某个平衡状态而言的）。但若系统只有唯一的平衡状态，则在该平衡状态处的稳定性就可视为整个系统的稳定性；若具有多个平衡状态的系统，其稳定性必须逐个讨论。

平衡状态

由于稳定性是系统在自由运动下的特性，故只需要考虑自治系统：
$\dot x=f(x,t)$

对上述系统而言，若存在状态向量 $x_e$ ，使得对所有时间 $t$ ，都有
$f(x_e,t) \equiv 0$

则称 $x_e$ 为系统的平衡状态或平衡点。事实上，平衡状态指的就是系统的静止状态（稳态）。并不是所有的系统都一定存在平衡状态，有时即使存在也未必是唯一的。

李雅普诺夫意义下的稳定

以下总是假定原点 $x=0$ 是系统 $\dot x=f(x,t)$ 的平衡状态，即 $f(0,t) \equiv 0$ 对所有时间 $t$ 成立。为了分析系统在原点处的稳定性，需要确定系统状态偏离原点的距离。在一般的 $n$ 维实数空间中，点 $x$ 到原点的距离定义为：
$\left\| x \right\| = \sqrt {x_1^2 + x_2^2 + \cdots + x_n^2}$

其中， $\left\| \cdot \right\|$ 称为向量的2范数。

$\left\| x \right\|\leqslant r$ 表示以原点为中心，半径为 $r$ 的球域 $S(r)$ 。当 $r$ 很小时，球域 $S(r)$ 也称为原点的一个邻域。

考虑系统的状态轨迹 $x\left( t \right) = \varphi \left( {t;{x_0},{t_0}} \right)$ ， $\left\| {\varphi \left( {t;{x_0},{t_0}} \right)} \right\| \leqslant r$ 对所有的时间 $t$ 成立表明系统的这一状态轨迹在原点的一个小邻域中。对应于图1，相当于小球始终在B点附近。

定义1：考虑系统 $\dot x=f(x,t)$ 的平衡状态 $x_e=0$ ，如果对任意给定的 $\varepsilon>0$ ，存在一个 $\delta>0$ （与 $\varepsilon$ 和初始时刻 $t_0$ 有关）,使得从球域 $S(\delta)$ 内任一初始状态出发的状态轨迹始终保持在球域 $S(\varepsilon)$ 内，则平衡状态 $x_e=0$ 称为在李雅普诺夫意义下是稳定的。

在图1中，若要使小球运动轨迹不超过A点的高度，则只要初值位置的高度不超过A点高度，就可以保证在以后所有时间内，小球运动时的高度都不会超过A点的高度。对应于定义1，给定A点高度就相当于任意给定 $\varepsilon>0$ ，存在 $\delta$ 的和A点高度相等.

从几何上来看，定义1所定义的系统稳定性意味着：对任意选择的一个球域 $S(\varepsilon)$ ，必存在另一个球域 $S(\delta)$ ，使得对所有的时间 $t$ ，始于球域 $S(\delta)$ 中（初始状态在球域 $S(\delta)$ 内，该球域表示了偏离平衡状态 $x_e=0$ 的界）的状态轨迹总不脱离球域 $S(\varepsilon)$ （注意：球域 $S(\varepsilon)$ 是任意的）。

图2 李雅普诺夫意义下的稳定性

定义2：考虑系统 $\dot x=f(x,t)$ 的平衡状态 $x_e=0$ ，如果平衡状态 $x_e=0$ 是李雅普诺夫意义下稳定的，并且当时 $t \to \infty$ ，始于原点邻域中的轨迹 $x(t) \to 0$ ，则平衡状态 $x_e=0$ 称为在李雅普诺夫意义下是渐进稳定的。

图3 李雅普诺夫意义下的渐进稳定性

在图1中，随着时间 $t$ 趋向于无穷，在B点附近出发运动的小球在摩擦力作用下慢慢回到平衡状态B点，因此B点处的平衡状态是在李雅普诺夫意义下渐进稳定的。

图3和图4表明了所考虑的二阶系统在原点处的渐进稳定性。从图中可以清楚地看出，当时间 $t$ 无限增加时，从球域 $S(\delta)$ 出发的状态轨迹不仅不会超出球域 $S(\varepsilon)$ ，而且最终收敛到原点。图3反映了状态轨迹 $x(t)$ 的有界性和渐进性；图4对状态轨迹 $x(t)$ 随时间变化的状况表示得更为清晰，它反映了初始状态在 $S(\delta)$ 内的状态轨迹随时间的推移，从球域 $S(\varepsilon)$ 范围内被压缩到球域 $S(\mu)$ 范围内。

图4 李雅普诺夫意义下的渐进稳定性

本文讨论的稳定性都是李雅普诺夫意义下的稳定性。在实际应用中，渐进稳定性比稳定性更重要，渐进稳定性表明系统能完全消除扰动的影响。同时需要注意的是，渐进稳定性只是一个局部的概念，它依赖系统的平衡状态。所以简单地确定了系统的渐进稳定性并不意味着系统能正常工作，通常有必要确定系统渐进稳定性的最大范围，即确定在多大范围内出发的状态轨迹将渐进趋向于所考虑的平衡状态。

定义3：若系统是稳定的，且对系统的任意状态，以该状态为初始状态的状态轨迹随着时间的推移都收敛到平衡状态 $x_e=0$ ，则系统称为是大范围渐进稳定的。

由于从状态空间中任意点出发的状态轨迹都要收敛于原点，因此，大范围渐进稳定的系统在整个状态空间中只能有一个平衡状态，这也是系统大范围渐进稳定的必要条件。

定义4：如果存在某个实数 $\varepsilon>0$ ，对不管多么小的 $\delta>0$ ，在球域 $S(\delta)$ 内总存在一个状态 $x_0$ ，使得始于这一状态的状态轨迹最终会离开球域 $S(\varepsilon)$ ，则平衡状态 $x_e=0$ 称为不稳定的。

图5 李雅普诺夫意义下的不稳定性

图1中的平衡状态D点就复合定义4的条件，因此是不稳定的。在图5中，状态轨迹离开了球域，这说明平衡状态是不稳定的。然而，这种情况未必意味着状态轨迹一定将趋于无穷远处。比如图1中的D点虽然不稳定，但随着时间的推移最后可能趋向于另一个平衡点B。

在稳定、渐进稳定和大范围渐进稳定这些定义中的 $\delta$ 一般总是与 $\varepsilon$ 和 $t_0$ 有关。但很多时候 $\delta$ 却与初始时间 $t_0$ 是无关的，此时可进一步称系统为一致稳定、一致渐进稳定和一致大范围渐进稳定。

能量函数

首先分析图1所示小球运动系统B点的稳定性。在一个初始外力的作用下，小球偏离原先的平衡状态到达A点（外力的作用给了小球能量），然后小球做自由运动（不受任何外力）。根据高中物理知识，小球不断做往复运动，能量（动、势能）不断转换。在这个过程中，系统没有从外部吸收能量（无外界输入），故系统总能量从不会增加（单调递减）。其次，在摩擦力的作用下，将消耗系统一定的能量，意味着系统最大的势能在减小，从而小球运动的最高点的高度不断下降。随着时间推移，系统总能量不断减少，最后系统的动能势能都将为零，小球静止在B点。这就是小球在B平衡点处的稳定性。

上面的例子说明系统的能量与系统稳定性之间的密切关系。那么这种能量与系统稳定性之间的关系能不能推广到更一般的系统呢？经过深入分析，李雅普诺夫给出了肯定的答案。然而，一般的系统未必具有那样物理意义清晰的能量函数。为此，李引入了虚拟能量函数的概念，并根据该虚拟能量函数沿系统状态轨迹随时间的变化情况，提出了一般系统基于能量函数的李雅普诺夫稳定性分析方法。

定理1：考虑非线性系统：
$\dot x=f(x(t),t)$

原点是该系统的平衡状态，即 $f(0,t) \equiv 0$ 。如果存在一个具有连续一阶偏导数的标量函数 $V(x,t)$ ，且满足一下条件：
（1） $V(x,t)$ 是正定的（标量函数恒大于0）；
（2）沿系统任意轨迹， $V(x,t)$ 关于时间 $t$ 的导数 $dV(x,t)/dt$ 是负定的；
则系统在原点处的平衡状态是渐进稳定的。满足以上条件（1）和（2）的标量函数 $V(x,t)$ 称为是系统的一个李雅普诺夫函数。进而，若 $\left\| x \right\| \to \infty$ 时，有 $V(x,t) \to \infty$ （径向无穷大），则在原点处的平衡状态 $x_e=0$ 是大范围渐进稳定的。

定理说明：

（1）定理1给出的系统稳定性条件仅仅是充分的，即如果找到一个李雅普诺夫函数 $V(x,t)$ ，则系统一定是渐进稳定的。但若找不到这样的李雅普诺夫函数，并不能说明系统是不稳定的。
（2）对于非线性系统，通过构造具体的李雅普诺夫函数，可以证明系统在某个稳定域内是渐进稳定的，但这并不意味着稳定域外的运动就是不稳定的。然后，可以证明：对于线性系统，如果存在渐进稳定的平衡点，则它必定是大范围渐进稳定的。
（3）若定理1条件（2）中的 $dV(x,t)/dt$ 是半负定的，则系统在平衡状态是稳定的。
（4）定理1既适用于线性系统、非线性系统，也适用于定常系统、时变系统。

注意：定理1中关于 $dV(x,t)/dt$ 必须是负定的条件还是比较苛刻的，因为它要求对所有非零的 $x$ ， $dV(x,t)/dt$ 都小于零。而事实上，只需要在系统的状态轨迹上 $V(x,t)$ 是减少的，即在系统的状态轨迹上 $V(x,t)$ 小于零就可以了。

定理2：考虑非线性系统：
$\dot x=f(x(t),t)$

原点是系统的平衡状态。若存在具有连续一阶偏导数的标量函数 $V(x,t)$ ，满足以下条件：
（1） $V(x,t)$ 是正定的；
（2）沿系统的任意轨迹， $V(x,t)$ 关于时间的导数 $dV(x,t)/dt$ 是半负定的；
（3）在系统的任意轨迹上， $dV(x,t)/dt$ 不恒等于零；
（4）当 $\left\| x \right\| \to \infty$ 时，有 $V(x,t) \to \infty$ ，则在原点处的平衡状态 $x_e=0$ 是大范围渐进稳定的。

关于系统的不稳定性

定理3：考虑非线性系统：
$\dot x=f(x(t),t)$

原点是系统的平衡状态。若存在具有连续一阶偏导数的标量函数 $V(x,t)$ ，满足以下条件：
（1） $V(x,t)$ 在原点附件的某一领域内是正定的；
（2） $dV(x,t)/dt$ 在同样的领域内也是正定的。
则系统在原点处的平衡状态是不稳定的。

显然，当 $dV(x,t)/dt$ 正定时，表示系统的能量在不断增大（不可能保持在原点的小邻域内），故系统的状态必将发散，远离原点。所以，系统是不稳定的。

李雅普诺夫稳定性方法在控制系统分析中的应用

李雅普诺夫稳定性方法在控制系统分析和设计中有着广泛的应用。它不仅可以用来判别一个系统（可以是非线性、时变）的稳定性，或者确定系统中某些参数的取值范围，使得系统保持稳定，还可以用于设计使得闭环系统稳定的控制器，即稳定化控制器的设计；线性系统时间常数的估计；确定系统的最优参数等。

参数优化问题

考虑线性系统：
$\dot x = A(\alpha)x，x(0) = x_0$

其中， $x$ 是系统的 $n$ 维状态向量，系统矩阵中含有可调参数 $\alpha$ 。一般的，参数 $\alpha$ 不仅可以影响系统的稳定性，而且还可以影响系统的动态特性。因此，希望选取最优参数 $\alpha$ ，使得系统不仅是渐近稳定的，同时还使得系统的性能指标：
$J{\text{ = }}\int_0^\infty {{x^{\rm T}}Qxdt}$

最小化，其中 $Q$ 是对称正定加权矩阵，可以根据任务设计。这样一个问题称为参数优化问题，其目的在于保证系统稳定的前提下，使得系统具有较好的过渡特性（动态过程）。性能指标 $J$ 越小，系统状态衰减到零的速度就越快，调节时间越短，震荡幅度也越小，故动态性能就越好。

我们可以用李雅普诺夫稳定性分析方法有效地解决这个问题，这种方法不仅能保证所求得的参数使系统渐近稳定，而且可以避免求解系统状态的微分方程和性能指标积分。

由于选择的参数 $\alpha$ 要保证系统是渐近稳定的，则必须对任意给定的对称正定矩阵 $R$ ，李雅普诺夫方程：
${A^{\rm T}}\left( \alpha \right)P + PA\left( \alpha \right) = - R$

存在唯一对称正定矩阵 $P$ （该矩阵需要我们去找）。此时， $V\left( x \right) = {x^{\rm T}}Px$ 是系统的一个李雅普诺夫函数，且沿着系统任意轨迹，李雅普诺夫函数的导数 $\dot V\left( x \right)$ ：
$dV(x)/dt = - x^{\rm T}Rx$

两边分别对时间 $t$ 积分，并利用系统的渐近稳定性：
$- \int_0^\infty {{x^{\rm T}}Rxdt} = \int_0^\infty {\frac{{dV\left( x \right)}}{{dt}}dt} = - V\left( {x\left( 0 \right)} \right) = - x_0^{\rm T}P{x_0}$

因此， $\int_0^\infty {{x^{\rm T}}Rxdt} = x_0^{\rm T}P{x_0}$

由于矩阵 $R$ 为任意对称正定阵，我们选 $R = Q$ ，则可得： $\int_0^\infty {{x^{\rm T}}Qxdt} = x_0^{\rm T}P{x_0}$

其中矩阵 $P$ 满足李雅普诺夫方程：
${A^{\rm T}}\left( \alpha \right)P + PA\left( \alpha \right) = - Q$

因此，我们得到了系统的性能指标可以通过求解以上一个静态的李雅普诺夫矩阵方程来计算，显然这比求解一个微分方程和积分要简单得多。从李雅普诺夫方程可以看出，李雅普诺夫矩阵 $P$ 依赖参数 $\alpha$ 。因此，
$J = {x_0}^{\rm T}P\left( \alpha \right){x_0}$

从而，原来的参数优化问题，转为求解参数 $\alpha$ 使得 $J$ 最小。由于 $J$ 是凸的，我们可以通过求一个无约束的极值问题，得到参数 $\alpha$ 的解析解，一般情况下参数 $\alpha$ 的取值与初始状态 $x_0$ 有关。
$\frac{{dJ}}{{d\alpha }} = 0$

特别地，从LQR的视角我们还可以知道，求得的最优李雅普诺夫矩阵 $P$ 就是在给定加权矩阵 $Q$ 和初始状态 $x_0$ 下系统的最优值函数。

最后编辑于：2018.12.13 22:14:01

人面猴
序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
沈念sama阅读 204,293评论 6赞 478
死咒
序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
沈念sama阅读 85,604评论 2赞 381
救了他两次的神仙让他今天三更去死
文/潘晓璐我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
开封第一讲书人阅读 150,958评论 0赞 337
道士缉凶录：失踪的卖姜人
文/不坏的土叔我叫张陵，是天一观的道长。经常有香客问我，道长，这世上最难降的妖魔是什么？我笑而不...
开封第一讲书人阅读 54,729评论 1赞 277
港岛之恋（遗憾婚礼）
正文为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
茶点故事阅读 63,719评论 5赞 366
恶毒庶女顶嫁案：这布局不是一般人想出来的
文/花漫我一把揭开白布。她就那样静静地躺着，像睡着了一般。火红的嫁衣衬着肌肤如雪。梳的纹丝不乱的头发上，一...
开封第一讲书人阅读 48,630评论 1赞 281
城市分裂传说
那天，我揣着相机与录音，去河边找鬼。笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。我是一名探鬼主播，决...
沈念sama阅读 38,000评论 3赞 397
双鸳鸯连环套：你想象不到人心有多黑
文/苍兰香墨我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
开封第一讲书人阅读 36,665评论 0赞 258
万荣杀人案实录
序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
沈念sama阅读 40,909评论 1赞 299
护林员之死
正文独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋以下内容为张勋视角年9月15日...
茶点故事阅读 35,646评论 2赞 321
白月光启示录
正文我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
茶点故事阅读 37,726评论 1赞 330
活死人
序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
沈念sama阅读 33,400评论 4赞 321
日本核电站爆炸内幕
正文年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
茶点故事阅读 38,986评论 3赞 307
男人毒药：我在死后第九天来索命
文/蒙蒙一、第九天我趴在偏房一处隐蔽的房顶上张望。院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
开封第一讲书人阅读 29,959评论 0赞 19
一桩弑父案，背后竟有这般阴谋
文/苍兰香墨我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。一阵脚步声响...
开封第一讲书人阅读 31,197评论 1赞 260
情欲美人皮
我被黑心中介骗来泰国打工，没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。一个月前我还...
沈念sama阅读 44,996评论 2赞 349
代替公主和亲
正文我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
茶点故事阅读 42,481评论 2赞 342