數(shù)值計算day3-求解線性方程組(上)

上節(jié)課主要介紹了非線性方程的幾種數(shù)值解法解恰，其中包括交叉法（二分法、線性插值法）和開放法（牛頓法浙于、割線法护盈、固定點法）。本節(jié)課主要介紹線性方程組的數(shù)值求解方法羞酗，主要分為直接法和迭代法兩類腐宋。直接法包括高斯消去法(Gauss elimination)、高斯約當法(Gauss-Jordan)以及LU分解法，迭代法包括Jacobi法和高斯塞德爾法(Gauss-Seidel) 胸竞。

1. 高斯消去法（Gauss Elimination Method）

高斯消去法是通過將線性方程組轉(zhuǎn)化為上三角的形式欺嗤，然后一步一步回代(back substitution)求解的方法。
例： $\begin{cases}x_1+2x_2=5 \\ 3x_1+4x_2=6\end{cases}$

系統(tǒng)法： $\begin{cases}x_1+2x_2=5 \\ 3x_1+4x_2=6\end{cases} \rightarrow \begin{cases}x_1+2x_2=5 \\ 0-2x_2=-9\end{cases} \rightarrow \begin{cases}x_1+2x_2=5 \\ x_2=4.5\end{cases} \rightarrow \begin{cases}x_1=-4 \\ x_2=4.5\end{cases}$
系統(tǒng)法的增廣矩陣形式：將方程寫為矩陣的形式 $\begin{bmatrix}1 & 2 \\3 & 4\end{bmatrix} \begin{bmatrix}x_1 \\ x_2 \end{bmatrix}=\begin{bmatrix}5\\6\end{bmatrix}$ 系統(tǒng)法的求解過程可以寫成如下增廣矩陣的變換形式 $\begin{bmatrix}1 & 2 & 5 \\3 & 4 & 6\end{bmatrix} \rightarrow \begin{bmatrix}1 & 2 & 5 \\0 & -2 & -9\end{bmatrix} \rightarrow \begin{bmatrix}1 & 2 & 5 \\0 & 1 & 4.5\end{bmatrix} \rightarrow \begin{bmatrix}1 & 0 & -4 \\0 & 1 & 4.5\end{bmatrix}$

例：
$\begin{cases}2x_1+x_2-x_3=1\\x_1+2x_2+x_3=8\\-x_1+x_2-x_3=-5\end{cases}$

matrix:
$\begin{bmatrix}2 & 1 & -1 & 1\\ 1 &2 & 1 & 8\\ -1 & 1 & -1 & -5 \end{bmatrix} \rightarrow \begin{bmatrix}2 & 1 & -1 & 1\\ 0 &1.5 & 1.5 & 7.5\\ 0 & 1.5 & -1.5 & -4.5 \end{bmatrix}$

$\rightarrow \begin{bmatrix}2 & 1 & -1 & 1\\ 0 &1.5 & 1.5 & 7.5\\ 0 & 0 & -3 & -12\end{bmatrix} \rightarrow \begin{cases}-3x_2 = -12\\x_3=4\end{cases} \rightarrow \begin{cases}x_1 =2 \\x_2=1\\x_3=4\end{cases}$

類似的卫枝， $n \times n$ 形式的方程組: $\begin{cases}a_{11}x_1+a_{12}x_2+\cdots+a_{1n}x_n=b_1 \\ a_{21}x_1+a_{22}x_2+\cdots+a_{2n}x_n=b_2\\ \vdots\\ a_{n1}x_1+a_{n2}x_2+\cdots+a_{nn}x_n=b_n\\ \end{cases}$
可以寫作如下矩陣形式： $A=\begin{bmatrix} a_{11} & a_{12} & \cdots & a_{1n}\\ a_{21} & a_{22} & \cdots & a_{2n}\\ \vdots\\ a_{n1} & a_{n2} & \cdots & a_{nn} \end{bmatrix}, x=\begin{bmatrix}x_1 \cdots x_n\end{bmatrix}, b=\begin{bmatrix}b_1 \cdots b_n\end{bmatrix}, Ax=b$

高斯消去法的核心：

方程組寫為矩陣形式
矩陣行變換：將一行中所有元素乘上一個標量煎饼，減去（加上）另外一行，該操作不會改變方程的解
逐行變換校赤，直至將矩陣轉(zhuǎn)換為便于求解的形式

一般而言吆玖，如下三種矩陣形式方便求解：

對角形式： $\begin{bmatrix}1 & & 0\\ & \ddots & \\0 & & 1 \end{bmatrix}$ ， $x_n=b_n$
上三角形式： $\begin{bmatrix} a_{11} & a_{12} & \cdots & a_{1n}\\ & a_{22} & \cdots & a_{2n}\\ & & \ddots&\vdots\\ & & & a_{nn}\\ \end{bmatrix}$ 痒谴，使用回代法求解(back substitution): $a_{nn}x_n=b_n \rightarrow x_n = \frac{b_n}{a_{nn}}$ $a_{n-1,n-1}x_{n-1}+ a_{n-1,n}x_n=b_{n-1}\rightarrow x_{n-1}=\frac{b_{n-1}-a_{n-1,n}x_n}{a_{n-1,n-1}}$ 重復下去 $x_i=\frac{b_i-\sum^n_{j=i+1}a_{ij}x_j}{a_{ii}}$
下三角形式： $\begin{bmatrix} a_{11} & & &\\ a_{21} & a_{22} & &\\ \vdots & \vdots & \ddots&\\ a_{n1}& a_{n2}& \cdots & a_{nn}\\ \end{bmatrix}$ 使用前向替換法(forward substitution): $a_{11}x_1=b_1 \rightarrow x_1 = \frac{b_1}{a_{11}}$ $a_{2,1}x_{1}+ a_{2,2}x_2=b_{2}\rightarrow x_{2}=\frac{b_{2}-a_{2,1}x_1}{a_{2,2}}$ 重復下去 $x_i=\frac{b_i-\sum^{i-1}_{j=1}a_{ij}x_j}{a_{ii}}$

Matlab實現(xiàn)（求解 $Ax=b$ , 使用左除 $x=A\backslash b$ ）:

>> A = [1 2;3 4]

A =

     1     2
     3     4

>> b = [5 6]'

b =

     5
     6

>> x = A\b

x =

   -4.0000
    4.5000

>> format long
>> x = A\b

x =

  -3.999999999999999
   4.499999999999999

MATLAB實現(xiàn)（高斯消去法）：

function x=Gauss(a,b)
  ab=[a,b];
  [R,C]=size(ab);
  for j=1:R-1
      for i=j+1:R
          ab(i,j:C)=ab(i,j:C)-ab(i,j)/ab(j,j)*ab(j,j:C);
      end
  end
  x=zeros(R,1);
  x(R)=ab(R,C)/ab(R,R);
  for i=R-1:-1:1
      x(i)=(ab(i,C)-ab(i,i+1:R)*x(i+1:R))/ab(i,i);
  end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
>> A = [1 2;3 4]

A =

     1     2
     3     4

>> b = [5 6]'

b =

     5
     6

>> format long
>> Gauss(A,b)

ans =

  -4.000000000000000
   4.500000000000000

MATLAB實現(xiàn)（回代法與前向替換法）

function y = BackwardSub(a,b)
% The function solves a system of linear equations ax=b 
% where a is upper triangular by using backward substitution.
% Input variables:
% a  The matrix of coefficients.
% b  A column vector of constants.
% Output variable:
% y  A colum vector with the solution.

n = length(b);
y(n,1) = b(n)/a(n,n);
for i = n-1:-1:1
    y(i,1)=(b(i)-a(i,i+1:n)*y(i+1:n,1))./a(i,i);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y = ForwardSub(a,b)
% The function solves a system of linear equations ax=b 
% where a is lower triangular by using forward substitution.
% Input variables:
% a  The matrix of coefficients.
% b  A column vector of constants.
% Output variable:
% y  A colum vector with the solution.

n = length(b);
y(1,1) = b(1)/a(1,1);
for i = 2:n
    y(i,1)=(b(i)-a(i,1:i-1)*y(1:i-1,1))./a(i,i);
end

2. 高斯主元消去(Gauss Pivoting)

可以發(fā)現(xiàn)衰伯，當主元為 $0$ 或非常小時铡羡，使用高斯消去法會失去精度积蔚，因此，在必要時可以交換行的順序: $\begin{bmatrix} 0 & 2 & 5\\ 3 & 4& 5 \end{bmatrix},$ 高斯主元法就是通過交換行的順序烦周，來保證主元不為0尽爆，并且盡可能取絕對值較大的值，下面通過一個案例读慎，說明主元很小時漱贱，高斯消去法會造成計算精度上出現(xiàn)較大誤差，而高斯主元法的計算結(jié)果會更為精確夭委。
例: $\begin{cases}0.0003x_1 + 12.34x_2=12.343 \\0.4321 x_1+x_2=5.321\end{cases}$ 使用高斯消去:

$\frac{0.4321}{0.0003}=1440$ (round-off error, since $0.0003*1440 = 0.4320$ )
$\begin{cases}0.0003x_1 + 12.34x_2=12.343 \\ 0.0001 x_1+17770x_2=-17760\end{cases} \rightarrow x_2 = -\frac{17760}{17770}=0.9994 \rightarrow x_1 = \frac{12.343-12.34*0.9994}{0.0003}=33.33$

MATLAB運算結(jié)果：

>> A = [0.0003 12.34;0.4321 1]

A =

   0.000300000000000  12.340000000000000
   0.432100000000000   1.000000000000000

>> b = [12.343 5.321]'

b =

  12.343000000000000
   5.321000000000000

>> A\b

ans =

  10.000000000000000
   1.000000000000000

若交換行序（高斯主元消去）：

$\frac{0.0003}{0.4321}=0.0006943$ (round-off error, since $0.4321*0.0006943 = 0.0003$ )
$\begin{cases}0.4321 x_1+x_2=5.321 \\ 0x_1 + 12.34x_2=12.34\end{cases} \rightarrow x_2 = -\frac{12.34}{12.34}=1\rightarrow x_1 = \frac{5.321-1}{0.4321}=10$

MATLAB實現(xiàn)（高斯主元消去）

function x=GaussPivot(a,b)
  ab=[a,b];
  [R,C]=size(ab);
  for j=1:R-1
      if ab(j,j)==0
          for k=j+1:R
              if ab(k,j)~=0
                  abTemp=ab(j,:);
                  ab(j,:)=ab(k,:);
                  ab(k,:)=abTemp;
                  break
              end
          end          
      end
      for i=j+1:R
          ab(i,j:C)=ab(i,j:C)-ab(i,j)/ab(j,j)*ab(j,j:C);
      end
  end
  x=zeros(R,1);
  x(R)=ab(R,C)/ab(R,R);
  for i=R-1:-1:1
      x(i)=(ab(i,C)-ab(i,i+1:R)*x(i+1:R))/ab(i,i);
  end
end

3. 高斯約當法（Gauss-Jordan Method幅狮，同時適用于求矩陣的逆）

高斯約當法一般采用高斯主元消去法，不同的是株灸，高斯約當法每次會將主元規(guī)范化為 $1$ 崇摄，且最終將原矩陣轉(zhuǎn)化為單位矩陣。
例： $\begin{bmatrix}4 & 1 & 2 & 21\\ 2 & -2 & 2 & 8\\ 1 & -2 & 4 & 16\end{bmatrix} \rightarrow \begin{bmatrix}1 & -2 & 4 & 16\\ 2 & -2 & 2 & 8 \\4 & 1 & 2 & 21 \end{bmatrix} \rightarrow \begin{bmatrix}1 & -2 & 4 & 16\\ 0 & 2 & -6 & -24 \\0 & 9 & -14 & -43 \end{bmatrix}$ $\rightarrow \begin{bmatrix}1 & -2 & 4 & 16\\ 0 & 1 & -3 & -12 \\0 & 9 & -14 & -43 \end{bmatrix} \rightarrow \begin{bmatrix}1 & 0 & -2 & -8\\ 0 & 1 & -3 & -12 \\0 & 0 & 13 & 65 \end{bmatrix} \rightarrow \begin{bmatrix}1 & 0 & -2 & -8\\ 0 & 1 & -3 & -12 \\0 & 0 & 1 & 5 \end{bmatrix}$ $\rightarrow \begin{bmatrix}1 & 0 & 0 & 2\\ 0 & 1 & 0 & 3 \\0 & 0 & 1 & 5 \end{bmatrix}$

可以看到慌烧，高斯約當法可以很方便地用來求解一個矩陣的逆： $AB=I$ , $A\begin{bmatrix}b_1 & b_2 & b_3 \end{bmatrix} \rightarrow \begin{bmatrix}a_{11} & a_{12} & a_{13} & 1 & 0 & 0\\ a_{21} & a_{22} & a_{23} & 0 & 1 & 0\\ a_{31} & a_{32} & a_{33} & 0 & 0 & 1\end{bmatrix} \rightarrow \begin{bmatrix}1 & 0 & 0 &c_{11} & c_{12} & c_{13} & \\ 0 & 1 & 0 & c_{21} & c_{22} & c_{23} &\\ 0 & 0 & 1& c_{31} & c_{32} & c_{33} \end{bmatrix}$

4. LU分解法（LU-decomposition Method）

LU分解是將一個矩陣分解為一個下三角矩陣和上三角矩陣的乘積： $A=LU$ , 其中 $L$ 表示下三角矩陣, $U$ 表示上三角矩陣逐抑。LU分解之后，求解線性方程 $Ax=b$ 就等價于求解如下兩個方程： $Ux=y,Ly=b$

LU分解主要有兩種方法來實現(xiàn)：

高斯消去(Gauss Elimination);
Crout's Method

4.1 Crout's Method

$\begin{bmatrix} a_{11} & a_{12} & \cdots & a_{1n}\\ a_{21} & a_{22} & \cdots & a_{2n}\\ \vdots\\ a_{n1} & a_{n2} & \cdots & a_{nn} \end{bmatrix}=\begin{bmatrix}L_{11} & & & & \\ L_{21} & L_{22} & & & \\ \vdots & \vdots & \ddots & \\ L_{n1} & L_{n2} & \cdots & L_{nn} \end{bmatrix} \begin{bmatrix}1 & U_{12} & U_{13} & \cdots & U_{1n} \\ & 1 & U_{23} & \cdots &U_{2n}\\ & & \ddots & & \vdots \\ & & & & 1\end{bmatrix}$

例:
$\begin{bmatrix} a_{11} & a_{12} & a_{13}\\ a_{21} & a_{22} & a_{23}\\ a_{31} & a_{32} & a_{33} \end{bmatrix}=\begin{bmatrix} L_{11} & 0 & 0\\ L_{21} & L_{22} & 0\\ L_{31} & L_{32} & L_{33} \end{bmatrix}\begin{bmatrix} 1 & U_{12} & U_{13}\\ 0 &1 & U_{23}\\ 0 & 0 & 1 \end{bmatrix}=\begin{bmatrix} L_{11} & L_{11}U_{12} & L_{11}U_{13}\\ L_{21} & L_{21}U_{12}+L_{22} & L_{21}U_{13}+ L_{22}U_{23}\\ L_{31} & L_{31}U_{12}+L_{32} & L_{33}+ L_{31}U_{13}+ L_{32}U_{23} \end{bmatrix}$
采用待定系數(shù)法求解： $L_{11}=a_{11},L_{21}=a_{21},L_{31}=a_{31}$ $U_{12} = \frac{a_{12}}{L_{11}}, U_{13}=\frac{a_{13}}{L_{11}}$ $L_{22}=a_{22}-L_{21}U_{12},L_{32}=a_{32}-L_{31}U_{12}$ $U_{23}=\frac{a_{23}-L_{21}U_{13}}{L_{22}}$ $L_{33}=a_{33}-L_{31}U_{13}-L_{32}U_{23}$
結(jié)合計算過程屹蚊，可以推導出厕氨，對于 $n\times n$ 矩陣，使用Crout's Method進行LU分解的通用形式為：

$L$ 矩陣的第一列 $L_{i1}=a_{i1},i=1,2,...,n$
$U$ 矩陣的對角元素 $U_{ii}=1,i=1,2,...,n$
$U$ 矩陣的第 $1$ 行 $U_{1j} = \frac{a_{1j}}{L_{11}},j=2,...,n$
$L$ 矩陣的其他元素 $L_{ij}=a_{ij}-\sum^{j-1}_{k=1}L_{ik}U_{kj},i=2,3,...,n,j=2,3,...,i$
$U$ 矩陣的其他元素 $U_{ij} = \frac{a_{ij}-\sum^{j-1}_{k=1}L_{ik}U_{kj}}{L_{ii}},i=2,3,...,n,j=i+1,i+2,...,n$

MATLAB實現(xiàn)(LU分解汹粤，Crout's Method)

function [L, U] = LUdecompCrout(A)
% The function decomposes the matrix A into a lower triangular matrix L 
% and an upper triangular matrix U, using Crout's method such that A=LU.
% Input variables:
% A  The matrix of coefficients.
% Output variable:
% L  Lower triangular matrix.
% U  Upper triangular matrix.

[R, C] = size(A);
for i = 1:R
    L(i,1) = A(i,1);
    U(i,i) = 1;
end
for j = 2:R
    U(1,j)= A(1,j)/L(1,1);
end
for i = 2:R
    for j = 2:i
        L(i,j)=A(i,j)-L(i,1:j-1)*U(1:j-1,j);
    end
    for j = i+1:R
        U(i,j)=(A(i,j)-L(i,1:i-1)*U(1:i-1,j))/L(i,i);
    end
end

4.2 高斯消去

$\begin{bmatrix} a_{11} & a_{12} & \cdots & a_{1n}\\ a_{21} & a_{22} & \cdots & a_{2n}\\ \vdots\\ a_{n1} & a_{n2} & \cdots & a_{nn} \end{bmatrix}=\begin{bmatrix}1 & & & & \\ m_{21} & 1 & & & \\ m_{31} & m_{32} & 1 & &\\ \vdots & \vdots & \vdots & \ddots & \\ m_{n1} & m_{n2} & \cdots & m_{n,n-1} & 1\end{bmatrix} \begin{bmatrix}a'_{11} & a'_{12} & \cdots & a'_{1n} \\ & a'_{22} & \cdots & a'_{2n}\\ & & \ddots & \vdots\\ & & & a'_{nn}\end{bmatrix}$

general formulation: $m_{ij}=\frac{a'_{ij}}{a'_{jj}}$

例: $\begin{bmatrix}1 & 2\\3 & 4\end{bmatrix}=\begin{bmatrix}1 & 0\\3 & 1\end{bmatrix}\begin{bmatrix}1 & 2\\0 & -2\end{bmatrix}$

高斯消去法進行LU分解也可以通過待定系數(shù)法求解出矩陣參數(shù)：

$U$ 矩陣第一行： $U_{1j}=a_{1j}命斧，j=1,2,...,n$
$L$ 矩陣對角元素： $L_{ii}=1,i=1,2,...,n$
$L$ 矩陣第一列： $L_{i1}=\frac{a_{i1}}{U_{11}}$
$U$ 矩陣其他元素： $U_{ij}=a_{ij}-\sum^{i-1}_{k=1}L_{ik}U_{kj},j=2,...,n,i=2,...,j$
$L$ 矩陣其他元素： $L_{ij}=\frac{a_{i,j}-\sum^{j-1}_{k=1}L_{ik}U_{kj}}{L_{jj}},j=2,...,n,i=j+1,...,n$

MATLAB實現(xiàn)(LU分解，高斯法嘱兼，待定系數(shù))

function [L,U] = LUdecopGauss(A)
[R,C] = size(A);
for i = 1:C
    L(i,i) = 1;
    U(1,i) = A(1,i);
end
for i = 2:R
    L(i,1) = A(i,1)/U(1,1);
end

for j = 2:C
    for i = 2:j
        U(i,j) = A(i,j)-L(i,1:i-1)*U(1:i-1,j);
    end
    for i= (j+1):R
        L(i,j) = (A(i,j)-L(i,1:j-1)*U(1:j-1,j))/U(j,j);
    end
end
end

高斯消去的待定系數(shù)法與Crout's 方法的待定系數(shù)法相比冯丙，要難以理解，可讀性差，主要在于通解形式與常規(guī)的求解順序不一致胃惜。為便于理解泞莉，還有一種拉格朗日形式來實現(xiàn)高斯消去法：
回顧高斯消去的過程：

第一步運算 $A = \begin{bmatrix} a_{11} & a_{12} & \cdots & a_{1n}\\ a_{21} & a_{22} & \cdots & a_{2n}\\ \vdots\\ a_{n1} & a_{n2} & \cdots & a_{nn} \end{bmatrix}\rightarrow A'=\begin{bmatrix} a_{11} & a_{12} & \cdots & a_{1n}\\ 0 & a'_{22} & \cdots & a'_{2n}\\ \vdots\\ 0 & a'_{n2} & \cdots & a'_{nn} \end{bmatrix}=L_1A,L_1=\begin{bmatrix} 1 & 0 & \cdots & 0\\ -\frac{a_{21}}{a_{11}} & 1 & \cdots & 0\\ \vdots\\ -\frac{a_{n1}}{a_{11}} & 0 & \cdots & 1 \end{bmatrix}$ 這是一步行變換
第二步運算 $A'=\begin{bmatrix} a_{11} & a_{12} & \cdots & a_{1n}\\ 0 & a'_{22} & \cdots & a'_{2n}\\ \vdots\\ 0 & a'_{n2} & \cdots & a'_{nn} \end{bmatrix}\rightarrow A''=\begin{bmatrix} a_{11} & a_{12} & \cdots & a_{1n}\\ 0 & a'_{22} & \cdots & a'_{2n}\\ \vdots\\ 0 & 0 & \cdots & a''_{nn} \end{bmatrix}=L_2A‘,L_2=\begin{bmatrix} 1 & 0 & \cdots & 0\\ 0 & 1 & \cdots & 0\\ \vdots\\ 0 & -\frac{a'_{n2}}{a'_{22}} & \cdots & 1 \end{bmatrix}$
第 $n-1$ 步運算后 $\begin{bmatrix}U_{11} & U_{12} & \cdots & U_{1n} \\ & U_{22} & \cdots & U_{2n}\\ & & \ddots & \vdots\\ & & & U_{nn}\end{bmatrix}=L_{n-1}L_{n-2}...L_1A$ 此時即為LU分解中的U矩陣，而 $L=L^{-1}_{1}...L^{-1}_{n-2}L^{-1}_{n-1}=\begin{bmatrix} 1 & 0 & \cdots & 0\\ \frac{a_{21}}{a_{11}} & 1 & \cdots & 0\\ \vdots\\ \frac{a_{n1}}{a_{11}} & 0 & \cdots & 1 \end{bmatrix}\begin{bmatrix} 1 & 0 & \cdots & 0\\ 0 & 1 & \cdots & 0\\ \vdots\\ 0 & \frac{a'_{n2}}{a'_{22}} & \cdots & 1 \end{bmatrix}...=\begin{bmatrix} 1 & 0 & \cdots & 0\\ \frac{a_{21}}{a_{11}} & 1 & \cdots & 0\\ \vdots\\ \frac{a_{n1}}{a_{11}} & \frac{a'_{n2}}{a'_{22}} & \cdots & 1 \end{bmatrix}$ 可以看到 $L_{ij}=\frac{a_{ij}}{a'_{jj}}$

MATLAB實現(xiàn)(LU分解船殉，高斯消去鲫趁，拉格朗日形式)

function [L,U]=LUGauss2(A)
  [R,C]=size(A);
   L = zeros(size(A));
  for i = 1:C
    L(i,i) = 1;
  end
  for j=1:R-1
      for i=j+1:R
          L(i,j) = A(i,j)/A(j,j);
          A(i,j:C)=A(i,j:C)-A(i,j)/A(j,j)*A(j,j:C);
      end
  end
  U = A;
end

5. 總結(jié)

本節(jié)課主要介紹了求解線性方程組的幾種直接法，包括高斯消去利虫、高斯主元消去挨厚、高斯約當、LU分解法等糠惫，其中高斯約當法可以用來求矩陣的逆疫剃。矩陣的LU分解又可以通過高斯消去和Crout 方法來實現(xiàn)，對參數(shù)矩陣進行LU分解后硼讽，很容易可以求解方程巢价。下節(jié)課介紹求解線性方程組的迭代法。

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者

人面猴
序言：七十年代末固阁，一起剝皮案震驚了整個濱河市壤躲，隨后出現(xiàn)的幾起案子，更是在濱河造成了極大的恐慌备燃，老刑警劉巖碉克，帶你破解...
沈念sama閱讀 212,080評論 6贊 493
死咒
序言：濱河連續(xù)發(fā)生了三起死亡事件，死亡現(xiàn)場離奇詭異并齐，居然都是意外死亡漏麦，警方通過查閱死者的電腦和手機，發(fā)現(xiàn)死者居然都...
沈念sama閱讀 90,422評論 3贊 385
救了他兩次的神仙讓他今天三更去死
文/潘曉璐我一進店門况褪，熙熙樓的掌柜王于貴愁眉苦臉地迎上來撕贞，“玉大人，你說我怎么就攤上這事窝剖÷榈В” “怎么了？”我有些...
開封第一講書人閱讀 157,630評論 0贊 348
道士緝兇錄：失蹤的賣姜人
文/不壞的土叔我叫張陵赐纱，是天一觀的道長脊奋。經(jīng)常有香客問我，道長疙描，這世上最難降的妖魔是什么诚隙？我笑而不...
開封第一講書人閱讀 56,554評論 1贊 284
?港島之戀（遺憾婚禮）
正文為了忘掉前任，我火速辦了婚禮起胰，結(jié)果婚禮上久又，老公的妹妹穿的比我還像新娘巫延。我一直安慰自己，他們只是感情好地消，可當我...
茶點故事閱讀 65,662評論 6贊 386
惡毒庶女頂嫁案：這布局不是一般人想出來的
文/花漫我一把揭開白布炉峰。她就那樣靜靜地躺著，像睡著了一般脉执。火紅的嫁衣襯著肌膚如雪疼阔。梳的紋絲不亂的頭發(fā)上，一...
開封第一講書人閱讀 49,856評論 1贊 290
城市分裂傳說
那天半夷，我揣著相機與錄音婆廊，去河邊找鬼。笑死巫橄，一個胖子當著我的面吹牛淘邻，可吹牛的內(nèi)容都是我干的。我是一名探鬼主播湘换，決...
沈念sama閱讀 39,014評論 3贊 408
雙鴛鴦連環(huán)套：你想象不到人心有多黑
文/蒼蘭香墨我猛地睜開眼宾舅，長吁一口氣：“原來是場噩夢啊……” “哼！你這毒婦竟也來了枚尼？” 一聲冷哼從身側(cè)響起贴浙，我...
開封第一講書人閱讀 37,752評論 0贊 268
萬榮殺人案實錄
序言：老撾萬榮一對情侶失蹤砂吞，失蹤者是張志新（化名）和其女友劉穎署恍，沒想到半個月后，有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體蜻直，經(jīng)...
沈念sama閱讀 44,212評論 1贊 303
?護林員之死
正文獨居荒郊野嶺守林人離奇死亡盯质，尸身上長有42處帶血的膿包…… 初始之章·張勛以下內(nèi)容為張勛視角年9月15日...
茶點故事閱讀 36,541評論 2贊 327
?白月光啟示錄
正文我和宋清朗相戀三年，在試婚紗的時候發(fā)現(xiàn)自己被綠了概而。大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片呼巷。...
茶點故事閱讀 38,687評論 1贊 341
活死人
序言：一個原本活蹦亂跳的男人離奇死亡，死狀恐怖赎瑰，靈堂內(nèi)的尸體忽然破棺而出王悍，到底是詐尸還是另有隱情，我是刑警寧澤餐曼，帶...
沈念sama閱讀 34,347評論 4贊 331
?日本核電站爆炸內(nèi)幕
正文年R本政府宣布压储，位于F島的核電站，受9級特大地震影響源譬，放射性物質(zhì)發(fā)生泄漏集惋。R本人自食惡果不足惜，卻給世界環(huán)境...
茶點故事閱讀 39,973評論 3贊 315
男人毒藥：我在死后第九天來索命
文/蒙蒙一踩娘、第九天我趴在偏房一處隱蔽的房頂上張望刮刑。院中可真熱鬧，春花似錦、人聲如沸雷绢。這莊子的主人今日做“春日...
開封第一講書人閱讀 30,777評論 0贊 21
一樁弒父案，背后竟有這般陰謀
文/蒼蘭香墨我抬頭看了看天上的太陽翘紊。三九已至胶惰，卻和暖如春，著一層夾襖步出監(jiān)牢的瞬間霞溪，已是汗流浹背孵滞。一陣腳步聲響...
開封第一講書人閱讀 32,006評論 1贊 266
情欲美人皮
我被黑心中介騙來泰國打工，沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留鸯匹，地道東北人坊饶。一個月前我還...
沈念sama閱讀 46,406評論 2贊 360
代替公主和親
正文我出身青樓，卻偏偏與公主長得像殴蓬，于是被迫代替她去往敵國和親匿级。傳聞我的和親對象是個殘疾皇子，可洞房花燭夜當晚...
茶點故事閱讀 43,576評論 2贊 349