5/8/02

THIS PACKAGE ENCLOSES SEVERAL MATLAB CODES ASSOCIATED WITH TWO CA MODELS
AND THE DATA ANALYSIS, AS WELL AS THE RAW DATA:

I. SIMULATION MODELS

1) CAoldTCP.M
CA MODEL WITH THE FEEDBACK TRANSMISSION CONTROL PROTOCOL (TCP)

2) CAoldCAC.M
CA MODEL WITH THE CONNECTION-ADMISSION CONTROL (CAC)


II. RAW DATA

1) TCP

* DEFAULT PARAMETERS:
EXPONENTIAL DITRIBUYION PARAMETERS: 
  lamd = 100  ("ON" STATE)
  lamd1 = 500 ("OFF" STATE)
THE THRESHOLD FOR COMPARING THE NORMALIZED ROUND-TRIP TIME (RTT)
  DRTT = 50

* DATA:
CAoldTCP-8-V5.MAT
CAoldTCP-8-V100.MAT
CAoldTCP-16-V20.MAT
CAoldTCP-16-V100.MAT
CAoldTCP-32-V10.MAT
CAoldTCP-32-V100.MAT

* NOTE:
EACH MAT-FILE, FOR EXAMPLE, CAoldTCP-8-V5.MAT, WHERE vmax = 5, 
AND L = 8, includes six matrices:

  Name        Size         Bytes  Class
  XX         64x40000   20480000  double array
  xout        1x40000     320000  double array
  xout1       1x40000     320000  double array
  xrtt        1x17516     140128  double array
  xtld        1x40000     320000  double array
  xtot        1x40000     320000  double array

XX IS A MATRIX (64x40000) OF THE QUEUE LENGTHS OF ALL NODES. xout AND xout1
REPRESENT THE NUMBER OF DATA PACKETS RECEIVED BY TWO NODES DURING EACH vmax. 
xrtt IS THE TIME SERIES OF RTT. xtot IS THE TOTAL NUMBER OF PACKETS TOUCHED
BY A NODE EACH vmax. xtld IS THE TOTAL LOAD OF THE NETWORK EACH vmax.


2) CAC

* DEFAULT PARAMETERS:
EXPONENTIAL DITRIBUYION PARAMETERS: 
  lamd = 100  ("ON" STATE)
  lamd1 = 500 ("OFF" STATE)
THE THRESHOLD FOR COMPARING THE NORMALIZED ROUND-TRIP TIME (RTT)
  DRTT = 50

* DATA:
CAoldCAC-8-V10.MAT
CAoldCAC-16-V10.MAT
CAoldCAC-32-V50.MAT

* NOTE:
EACH MAT-FILE, FOR EXAMPLE, CAoldCAC-16-V10, WHERE vmax = 10, AND L = 16, 
includes six matrices:

  Name        Size         Bytes  Class
  XX        256x20000   40960000  double array
  xout        1x20000     160000  double array
  xout1       1x20000     160000  double array
  xrtt        1x14059     112472  double array
  xtld        1x20000     160000  double array
  xtot        1x20000     160000  double array

XX IS A MATRIX (256x40000) OF THE QUEUE LENGTHS OF ALL NODES. xout AND xout1
REPRESENT THE NUMBER OF DATA PACKETS RECEIVED BY TWO NODES DURING EACH vmax. 
xrtt IS THE TIME SERIES OF RTT. xtot IS THE TOTAL NUMBER OF PACKETS TOUCHED
BY A NODE EACH vmax. xtld IS THE TOTAL LOAD OF THE NETWORK EACH vmax.


III. DATA ANALYSIS

IN OUR PAPER, WE ANALYZE THREE TYPES OF TIME SERIES: THE QUEUE LENGTH (Nr),
xout, AND THE NUMBER OF CONGESTED NODES (y). TO OBSERVE HOW Nr CHANGES AT 
DIFFERENT TIME GRANULARITIES, WE CAN EASILY USE THE COMMAND "plot", E.G., 
plot(XX(1,1:10:10000) WITH THE GRAINSIZE 10.

WE HAVE TWO M-FILES TO EXTRACT TIME SERIES WE WANT FROM RAW DATA:

1) x = taggregat(grainsize,T1,xout)
"taggregat" OUTPUTS A TIME SERIES OF "x" FROM "xout", WITH THE TIME GRANULARITY, 
"grainsize", AND THE LENGTH "T1".  

2) y = congestednode(L,Y,intval,T2,XX)
"congestednode" PRODUCES A TIME SERIES OF "y" FROM "XX", WITH THE TIME GRANULARITY
"intval", AND THE THRESHOLD "Y", WHERE "L" IS THE SYSTEM SIZE, AND "T2" IS THE 
LENGTH OF "XX".

FOR THE POWER SPECTRAL ANALYSIS, WE CAN EASILY USE THE MATLAB COMMANDS, "psd" OR
"fft".


IV. REFERENCE

J. Yuan, K. Mills, "Exploring Collective Dynamics in Communication Networks", 
Journal of Research of the National Institute of Standards and Technology, 
107 (2), pp. 179-191, 2002.