13-六月
-2011
17:55 456029056 456029056 620781440 620781440 456029056 620781440
核心出口H3C-GigabitEthernet0/3 _TO_学生
Min/Max/Average bps of Recv 1.0 Gbps Xmit 1.0 Gbps
LAST 7 DAYS
AVERAGE
MIN/MAX TRANSMIT RECEIVE
DATE / TIME MIN/MAX RECEIVE BPS BPS BPS 13-六
月-2011
AVERAGE TRANSMIT
BPS
0:05 99511848 99511848 64362632 64362632 99511848 64362632
13-六
月-2011
0:15 85477992 85477992 80723752 80723752 85477992 80723752
13-六
月-2011
0:25 21495884 21495884 10103505 10103505 21495884 10103505
13-六
月-2011
0:35 62327008 62327008 50474480 50474480 62327008 50474480
13-六
月-2011
0:45 3301606.75 3301606.75 49534832 49534832 3301607 49534832
13-六
月-2011
0:55 59641512 59641512 34217256 34217256 59641512 34217256
附录3 模型3的程序：
for i=1:length(n)
% 计算这个子序列 a=floor(N/n(i));
15
% 仓健这个子序列的矩阵
X=reshape(x(1:a*n(i)),n(i),a);
% 估算这个子序列的平均值 ave=mean(X);
% 给这个序列的每一个值除以平均值 cumdev=X-ones(n(i),1)*ave;
% 估算累计离差
cumdev=cumsum(cumdev);
% 估算这个标准偏差 switch method case 'Hurst'
% Hurst-Mandelbrot 参数 stdev=std(X); case 'Lo' % Lo 参数 for j=1:a sq=0; for k=0:q
v(k+1)=sum(X(k+1:n(i),j)'*X(1:n(i)-k,j))/(n(i)-1); if k>0
sq=sq+(1-k/(q+1))*v(k+1); end end
stdev(j)=sqrt(v(1)+2*sq); end case 'MW'
% Moody-Wu 参数 for j=1:a sq1=0; sq2=0; for k=0:q
v(k+1)=sum(X(k+1:n(i),j)'*X(1:n(i)-k,j))/(n(i)-1); if k>0
sq1=sq1+(1-k/(q+1))*(n(i)-k)/n(i)/n(i); sq2=sq2+(1-k/(q+1))*v(k+1); end end
stdev(j)=sqrt((1+2*sq1)*v(1)+2*sq2); end
16
case 'Parzen' % Parzen 参数 if mod(q,2)~=0
error('在\参数中q 必须是2.'); end for j=1:a sq1=0; sq2=0; for k=0:q
v(k+1)=sum(X(k+1:n(i),j)'*X(1:n(i)-k,j))/(n(i)-1); if k>0 & k<=q/2
sq1=sq1+(1-6*(k/q)^2+6*(k/q)^3)*v(k+1); elseif k>0 & k>q/2
sq2=sq2+(1-(k/q)^3)*v(k+1); end end
stdev(j)=sqrt(v(1)+2*sq1+2*sq2); end otherwise
error('你应该付给 \另一个值.'); end
% 估算(R(t)/s(t))
rs=(max(cumdev)-min(cumdev))./stdev;
clear stdev
% 取这个平均值 R/S的对数
logRS(i,1)=log10(mean(rs));
if nargout>1
% 开始计算log(E(R/S)) j=1:n(i)-1;
s=sqrt((n(i)-j)./j); s=sum(s);
% 估算log(E(R/S))
logERS(i,1)=log10(s/sqrt(n(i)*pi/2));
%其它估算log(E(R/S))
%logERS(i,1)=log10((n(i)-0.5)/n(i)*s/sqrt(n(i)*pi/2)); %logERS(i,1)=log10(sqrt(n(i)*pi/2));
17
end
if nargout>2
% 估算 V
V(i,1)=mean(rs)/sqrt(n(i)); end end
附录4分析周期性时和突发性的图
此图是核心学生服务器上每天学生宿舍的平均接收量同时间的关系图
第
一
天
第二天:
18
第
三
天
第
四
天
19