TCP Fluid Model

 没有数学理论指导的研究，没有意义。本篇主要介绍Tcp拥塞控制中的流模型[1]，现在假设，一条链路上只有一个tcp流。

符号	含义
$W(t)$	拥塞窗口
$R(t)$	RTT
$Q(t)$	网络队列长度
$X(t)$	发送速率
$C$	瓶颈链路带宽

 需要点背景，经典的TCP的拥塞控制算法遵从AIMD原则.TCP在拥塞避免阶段，每收到一个ACK，窗口大小增加$\frac{1}{W}$,经过一个RTT，窗口大小增加1.这句话用数学公式表示：
$$dW=\frac{1}{R(t)}dt\text{\hspace{1em}(1)}$$
 (1)当$X(t)<C$时，就是发送速率小于网络中的带宽处理能力，不存在排队现象，有如下结论：
$$\dot{W}=\frac{1}{R_{min}}\text{\hspace{1em}(2)}$$
$$R(t)=R_{min}\text{\hspace{1em}(3)}$$
$$Q(t)=0\text{\hspace{1em}(4)}$$
$$X(t)=\frac{W(t)}{R_{min}}\text{\hspace{1em}(5)}$$
 (2)当$X(t)>C$时，发送速率超出了网络的处理能力，网络中开始出现排队现象，有如下结论：
$$\dot{W}(t)=\frac{1}{R(t)}\text{\hspace{1em}(6)}$$
$$R(t)=R_{min}+\frac{Q(t)}{C}\text{\hspace{1em}(7)}$$
$$Q=\int (X(t)-C)dt\Rightarrow \dot{Q}(t)=X(t)-C\text{\hspace{1em}(8)}$$
$$X(t)=(1+\frac{1}{W(t)})C\text{\hspace{1em}(9)}$$
 公式(9)不是很明显，按照[1]的说法，发送速率已经超出了瓶颈链路的处理能力，ACK返回的速率$ack\text{\_}rate$为瓶颈链路的处理能力$C$,所以$X(t)=(1+\frac{1}{W(t)})\ast ack\text{\_}rate$.
 以上几个公式可以说是Tcp研究中的麦克斯韦方程组。
 求解，对(7)求导：
KaTeX parse error: \tag works only in display equations
 根据(6)(10),得到：
$$1=\dot{W}(t)R(t)=W(t)\dot{R}(t)\Rightarrow \frac{d(R(t)W(t))}{dt}=2\text{\hspace{1em}(11)}$$
 并且,对(12)求导可以得到(11)左边等式的结论，说明在达到瓶颈链路的处理能力后，窗口的增加，带来时延的增加，发送速率并没有增益。
$$\frac{W(t)}{R(t)}=C\text{\hspace{1em}(12)}$$
 将(12)带入(11):
$$CR\dot{R}=1\Rightarrow {\int }_{R_{min}}^{R}RdR={\int }_0^t\frac{1}{C}dt\Rightarrow R(t)=\sqrt{\frac{2}{C}t+R_{min}^2}\text{\hspace{1em}(13)}$$
 (13)代入(12):
$$W(t)=CR(t)=\sqrt{2C\ast t+(C{R}_{min}{)}^2}=\sqrt{2C\ast t+(W(0){)}^2}\text{\hspace{1em}(14)}$$
$$\dot{Q}(t)=\frac{C}{W(t)}=\dot{W}(t)\Rightarrow Q(t)=W(t)-W(0)\text{\hspace{1em}(15)}$$
$$X(t)=C+\frac{1}{R(t)}\text{\hspace{1em}(16)}$$
 而且，还可以得到：
$$\dot{Q}(t)=\frac{C}{W(t)}=\frac{1}{R(t)}=\dot{W}(t)\text{\hspace{1em}(17)}$$

 推广到多条流：

符号	含义
$W_i(t)$	第i个tcp流拥塞窗口
$R_i(t)$	第i个tcp流RTT
$Q(t)$	网络队列长度
$X(t)$	所有数据流发送速率和
$X_i(t)$	第i个tcp流发送速率
$C$	瓶颈链路带宽
$Y_i(t)$	第i个tcp流返回ack的速率

$$\dot{W_i}(t)=\frac{1}{R_i(t)}\text{\hspace{1em}(18)}$$
$$R_i(t)=R_{i,min}+\frac{Q(t)}{C}\text{\hspace{1em}(19)}$$
$$\dot{Q}(t)=X(t)-C\Rightarrow \dot{Q}(t)=-C+\frac{N}{R_i(T)}{W}_i(t)\text{\hspace{1em}(20)}$$
$$X_i(t)=(1+\frac{1}{W_i(t)})Y_i(t)\text{\hspace{1em}(21)}$$
 根据(17)可以得到：
$$\dot{Q}(t)=\sum _{i=1}^N\frac{1}{R_i(t)}=\sum _{i=1}^N\frac{1}{R_{i,min}+\frac{Q(t)}{C}}\text{\hspace{1em}(22)}$$
 公式(22)还有另外的表达形式，参考[2].
[1]A Note on Tcp Fluid Model
[2]Hollot C V, Misra V, Towsley D, et al. Analysis and design of controllers for AQM routers supporting TCP flows[J]. IEEE Transactions on automatic control, 2002, 47(6): 945-959.