WebRTC拥塞控制算法

　　背景
　　网络拥塞是基于IP协议的数据报交换网络中常见的一种网络传输问题，它对网络传输的质量有严重的影响，网络拥塞是导致网络吞吐降低，网络丢包等的主要原因之一，这些问题使得上层应用无法有效的利用网络带宽获得高质量的网络传输效果。特别是在通信领域，网络拥塞导致的丢包，延迟，抖动等问题，严重的影响了通信质量，如果不能很好的解决这些问题，一个通信产品就无法在现实环境中正常使用。丢包带宽估计
　　基于丢包的拥塞控制比较简单，其基本思想是根据丢包的多少来判断网络的拥塞程度。丢包越多则认为网络越拥塞，那么我们就要降低发送速率来缓解网络拥塞；如果没有丢包，这说明网络状况很好，这时候就可以提高发送码率，向上探测是否有更多的带宽可用。实现该算法有两点：一是获得接收端的丢包率，一是确定降低码率和提升码率的阈值。
　　Ashat（i）Ashat（i1），2丢包率10，保持不变Ashat（i）Ashat（i1）（10。5p），丢包率10，下降，p是丢包率Ashat（i）1。05Ashat（i1），丢包率2，上升
　　webrtc中发送端收到RTCPRR报文并解析得到丢包率后，根据下图公式计算发送端码率：当丢包率大于0。1时，说明网络发生拥塞，此时降低发送端码率；当丢包率小于0。02时，说明网络状况良好，此时增大发送端码率；其他情况下，发送端码率保持不变。延迟带宽估计GoogREMB
　　基于接收端的延迟算法，利用延迟值，通过卡尔曼滤波器估计出下一时刻的网络带宽趋势，效果的准确性和公平性不如TransportCC，目前已经被webrtc的新版本所淘汰TransportCC
　　基于发送端的延迟算法，也是利用区间延迟值，通过TrendLine滤波器（最小二乘法滤波器），通过斜率的增加或者减小来判断当前网络的拥塞状况时延梯度计算
　　首先我们通过一张图看下时延梯度的计算：
　　对于两个包组：i以及i1，它们的时延梯度：
　　在WebRTC中，延迟梯度不是一个个包来计算的，而是通过将包分组，然后计算这些包组之间的延迟，这样做可以减少计算次数，同时减少误差。包组的划分原则为在一个bursttime间隔内的一系列的包构成一个组。建议bursttime为5ms。
　　C音视频开发WebRTC学习资料：点击领取音视频开发（资料文档视频教程面试题）（FFmpegWebRTCRTMPRTSPHLSRTP）
　　TrendLine滤波器（最小二乘法滤波器）
　　使用到了线性回归的方法进行时延趋势预测，通过最小二乘法求得拟合的直线斜率
　　对于一堆样本点（x，y），拟合直线方程yaxb的斜率a按如下公式计算
　　具体方式就是对距离平方和的最小值，根据求导推出a的值
　　代码实现路径，googcctrendlineestimator。cc，实现和上述推导公式一致absl：：optionaldoubleLinearFitSlope（conststd：：dequeTrendlineEstimator：：PacketTimingpackets）｛RTCDCHECK（packets。size（）2）；Computethecenterofmass。doublesumx0；doublesumy0；for（constautopacket：packets）｛计算x，y的总值sumxpacket。arrivaltimems；sumypacket。smootheddelayms；｝计算x，y的平均值doublexavgsumxpackets。size（）；doubleyavgsumypackets。size（）；Computetheslopeksum（xixavg）（yiyavg）sum（xixavg）2doublenumerator0；doubledenominator0；for（constautopacket：packets）｛doublexpacket。arrivaltimems；doubleypacket。smootheddelayms；numerator（xxavg）（yyavg）；denominator（xxavg）（xxavg）；｝if（denominator0）returnabsl：：nullopt；returnnumeratordenominator；｝
　　整体调用流程为LinearFitSlopeUpdateTrendlineLinearFitSlopeLinearFitSlopeTrendlineEstimator：：DetectTrendlineEstimator：：UpdateThresholdUpdateTrendline
　　这个方法主要是做一些数据的统计计算之后，调用LinearFitSlope计算斜率voidTrendlineEstimator：：UpdateTrendline（doublerecvdeltams，doublesenddeltams，int64tsendtimems，int64tarrivaltimems，sizetpacketsize）｛时延变化：接收时间差发送时间差constdoubledeltamsrecvdeltamssenddeltams；numofdeltas；numofdeltasstd：：min（numofdeltas，kDeltaCounterMax）；if（firstarrivaltimems1）firstarrivaltimemsarrivaltimems；Exponentialbackofffilter。accumulateddelaydeltams；BWETESTLOGGINGPLOT（1，accumulateddelayms，arrivaltimems，accumulateddelay）；根据1）中的时延累积值计算得到平滑后的时延：smoothingcoef为0。9smootheddelaysmoothingcoefsmootheddelay（1smoothingcoef）accumulateddelay；BWETESTLOGGINGPLOT（1，smootheddelayms，arrivaltimems，smootheddelay）；if（delayhist。size（）settings。windowsize）delayhist。popfront（）；Simplelinearregression。doubletrendprevtrend；当队列delayhist大小等于设定的窗口大小时，开始进行时延变化趋势计算，得到直线斜率if（delayhist。size（）settings。windowsize）｛trendLinearFitSlope（delayhist）。valueor（trend）；if（settings。enablecap）｛absl：：optionaldoublecapComputeSlopeCap（delayhist，settings）；Weonlyusethecaptofilteroutoverusedetections，nottodetectadditionalunderuses。if（trend0cap。hasvalue（）trendcap。value（））｛trendcap。value（）；｝｝｝BWETESTLOGGINGPLOT（1，trendlineslope，arrivaltimems，trend）；传入trend来计算当前网络状态Detect（trend，senddeltams，arrivaltimems）；｝Detect
　　该函数主要根据时延变化增长趋势计算当前网络状态。在Detect函数内部，会根据前面计算得到的斜率得到一个调整后的斜率值：modifiedtrend：
　　【小编强力推荐】音视频开发教程视频：【免费】FFmpegWebRTCRTMPNDKAndroid音视频流媒体高级开发学习视频教程腾讯课堂
　　C音视频开发WebRTC学习资料：点击领取音视频开发（资料文档视频教程面试题）（FFmpegWebRTCRTMPRTSPHLSRTP）
　　voidTrendlineEstimator：：Detect（doubletrend，doubletsdelta，int64tnowms）｛if（numofdeltas2）｛hypothesisBandwidthUsage：：kBwNormal；return；｝constdoublemodifiedtrendstd：：min（numofdeltas，kMinNumDeltas）trendthresholdgain；prevmodifiedtrendmodifiedtrend；BWETESTLOGGINGPLOT（1，T，nowms，modifiedtrend）；BWETESTLOGGINGPLOT（1，threshold，nowms，threshold）；1。if（modifiedtrendthreshold）｛if（timeoverusing1）｛Initializethetimer。Assumethatwevebeenoverusinghalfofthetimesincetheprevioussample。timeoverusingtsdelta2；｝else｛Incrementtimertimeoverusingtsdelta；｝overusecounter；23if（timeoverusingoverusingtimethresholdoverusecounter1）｛4if（trendprevtrend）｛timeoverusing0；overusecounter0；hypothesisBandwidthUsage：：kBwOverusing；｝｝｝elseif（modifiedtrendthreshold）｛timeoverusing1；overusecounter0；hypothesisBandwidthUsage：：kBwUnderusing；｝else｛timeoverusing1；overusecounter0；hypothesisBandwidthUsage：：kBwNormal；｝prevtrendtrend；UpdateThreshold（modifiedtrend，nowms）；｝
　　constdoublemodifiedtrendstd：：min（numofdeltas，kMinNumDeltas）trendthresholdgain；
　　numofdeltas表示包组间延迟梯度计算的次数，取值范围是〔2，60〕，thresholdgain是Trendline滤波器的增益参数，其默认值为4，这两个变量都会对延迟梯度趋势值进行放大。这里的关键就在于为何要放大的意义trend在于是一个斜率值，取值范围（1，1），这个值很小，而阈值要跟随的变化而变化，这可能导致很容易的超过从而触发连续的过载信号。放大可以使得算法不会因为的波动而过于敏感。
　　判断为过载的条件延迟梯度斜率值当前阈值过载总时长kOverUsingTimeThreshold（10ms）过载次数1当前延迟梯度斜率值上一次的延迟梯度斜率值，即延迟在不断恶化。UpdateThreshold
　　阈值threshold动态调整为了改变算法对时延梯度的敏感度。根据主要是，时延梯度是变化的，有时很大，有时很小，如果阈值是固定的，对于时延梯度来说可能太大或者太小，这样就会出现检测不够敏感，无法检测到网络拥塞，或者过于敏感，导致一直检测为网络拥塞；voidTrendlineEstimator：：UpdateThreshold（doublemodifiedtrend，int64tnowms）｛根据斜率上升还是下降定义K系数kdown（0。039）kup（0。0087），constdoublekfabs（modifiedtrend）threshold？kdown：kup；constint64tkMaxTimeDeltaMs100；计算上一次更新time的时间，和100ms取最小值int64ttimedeltamsstd：：min（nowmslastupdatems，kMaxTimeDeltaMs）；将k和timedeltams，进行计算，更新阈值thresholdk（fabs（modifiedtrend）threshold）timedeltams；thresholdrtc：：SafeClamp（threshold，6。f，600。f）；lastupdatemsnowms；｝