Multicast-based inference of temporal delay characteristics in packet data networks

Abstract

Disclosed are method and apparatus for characterizing the temporal delay characteristics of a packet data network by multicast-based inference. Multicast probes are transmitted from a source node to a plurality of receiver nodes, which record the delays of the multicast probes. From the aggregate data comprising recorded delays of the end-to-end paths from the source node to each receiver node, temporal delay characteristics of individual links within the network may be calculated. In a network with a tree topology, the complexity of calculations may be reduced through a process of subtree partitioning.

Description

BACKGROUND OF THE INVENTION

The present invention relates generally to network characterization of packet delay, and more particularly to network characterization of packet delay by multicast-based inference.

Packet data networks, such as Internet Protocol (IP) networks, were originally designed to transport basic data in a packetized format. Increasingly, however, other services, such as voice over IP (VoIP) and video on demand (VOD), are utilizing packet data networks. These services, in general, have more stringent requirements for network quality of service (QoS) than basic data transport. Depending on the application, QoS is characterized by different parameters. In addition to packet loss, an important parameter is packet delay. Services such as VoIP, for example, operate in real time (or, at least, near-real time). Excessive delay will result in poor voice quality. Even if only data is being transported, competing services using the same transport network may have different QoS requirements. For example, near-real time system control will have more stringent delay requirements than download of music files. In some instances, QoS requirements are set by service level agreements between a network provider and a customer.

Measurement of various network parameters is essential for network planning, architecture, administration, and diagnostics. Some parameters may be measured directly by network equipment, such as routers and switches. Since different network providers typically do not share this information with other network providers and with end users, however, system-wide information is generally not available to a single entity. Additionally, the measurement capabilities of a piece of network equipment are typically dependent on proprietary network operation systems of equipment manufacturers. The limitations of internal network measurements are especially pronounced in the public Internet, which comprises a multitude of public and private networks, often stitched together in a haphazard fashion.

A more general approach to network characterization, therefore, needs to be independent of measurements captured by equipment internal to the transport network. That is, the measurements need to be performed by user-controlled hosts attached to the network. One approach is for one host to send a test message to another host to characterize the network link between them. A standard message widely utilized in IP networks is a “ping”. Host A sends a ping to Host B. Assuming that Host B is operational, if the network connection between Host A and Host B is operational, Host A will receive a reply message from Host B. A field in the reply message records the round-trip time (RTT). If Host A does not receive a reply within a user-defined timeout interval, it declares the message to have been lost. Pings are examples of point-to-point messages between two hosts. As the number of hosts connected to the network increases, the number of point-to-point test messages increases to the level at which they are difficult to administer. They may also produce a significant load on both the hosts and the transport network. A key requirement of any test tool is that it must not corrupt the system under test. In addition to the above limitations, in some instances, pings may not provide the level of network characterization required for adequate network planning, architecture, administration, and diagnostics.

What is needed is a network characterization tool which provides detailed parameters on the network, runs on hosts controlled by end users, and has minimal disturbance on the operations of the hosts and transport network.

BRIEF SUMMARY OF THE INVENTION

Temporal delay characteristics in packet data networks are characterized by multicast-based inference. A packet data network comprises a set of nodes connected by a set of paths. Each path may comprise a set of individual links. In multicast-based inference, multiple test messages (probes) are multicast from a source node to a set of receiver nodes. Each receiver node records the delays of the probes transmitted along an end-to-end path from the source node to the receiver node. From the aggregate delay data collected by the set of receiver nodes, temporal delay characteristics of individual links may be calculated. In addition to average delay per unit time, temporal delay characteristics comprise parameters such as the number of probes with delays less than a specified value and the number of probes with delays greater than a specified value. Probes with delays greater than a threshold value may be declared to be lost probes. In embodiments in which the topology of the packet data networks are trees, calculations may be simplified by a process of subtree partitioning.

These and other advantages of the invention will be apparent to those of ordinary skill in the art by reference to the following detailed description and the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows a schematic of a packet data communications system;

FIG. 2 shows a schematic of a tree model of a network;

FIG. 3 shows a schematic of a network test architecture;

FIG. 4 shows a flow chart of a multicast-based temporal network characterization process;

FIG. 5 is a schematic for subtree partitioning of a binary tree;

FIG. 6 is a schematic for subtree partitioning of an arbitrary tree;

FIG. 7 shows a flowchart of a multicast-based temporal delay network characterization process; and,

FIG. 8 is a schematic of a computer for performing a multicast-based network characterization process.

DETAILED DESCRIPTION

FIG. 1 shows a network architecture schematic of an example of a communications system comprising packet data network 102 and end-user nodes 104-110. Within packet data network 102 are edge nodes 112-118 and intermediate nodes 120 and 122. An edge node connects an end-user node to a packet data network. An intermediate node connects nodes within a network. In some instances, a node may serve as both an edge node and an intermediate node. Herein, an edge node and an intermediate node are considered to be logically equivalent, and “intermediate nodes” comprise both intermediate nodes and edge nodes. Herein, “nodes” comprise both physical nodes and logical nodes. An example of a physical end-user node is a host computer. An example of a logical end-user node is a local area network. An example of a physical intermediate node is a router. An example of a logical intermediate node is a subnetwork of routers, switches, and servers. In all instances, an end user may access and control an end-user node. Access and control policies for an intermediate node, however, are set by a network provider, which, in general, is a different entity from an end user. In general, an end user may not have permission to access and control an intermediate node.

Nodes are connected via network links, which comprise physical links and logical links. In FIG. 1, links 124-140 represent physical links. Examples of physical links include copper cables and optical fiber. Links 142 and 144 represent logical links. For example, logical link 142 represents an end-to-end network link along which data is transmitted between end-user node 106 and end-user node 104. Logical link 142 comprises physical links 126, 134, 132, and 124. A logical link may further comprise segments which are also logical links. For example, if intermediate node 114 is a router, there is both a physical link for signal transport across the router and a logical link for data transport across the router. Logical links may span multiple combinations of end-user nodes and intermediate nodes. Logical links may span multiple networks. Since a physical link may also be considered a logical link, a network link may also be referred to herein simply as a “link”. Herein, an end-to-end network link connecting one node to another node may also be referred to as a “path”. A path may comprise multiple links.

In an embodiment, characterization of packet data network 102 is performed by multicasting test messages from a source node (for example, end-user node 106) to receiver nodes (for example, end-user nodes 104, 108, and 110). Analysis of the test messages transmitted from source node 106 and received by a specific receiver node (for example, node 104) yields characteristics of the path from the source node 106 to the specific receiver node 104. In addition, test messages received at all the receiver nodes may be aggregated to infer characteristics of internal network links. For example, in FIG. 1, test messages transmitted from the source node 106 to receiver nodes 104, 108, and 110 all must pass through the common network link defined by (link 126-node 114-link 134-node 120). Thus, if test messages are received at any one of the receiver nodes 104, 108, and 110, then source 106, link 126, node 114, link 134, and node 120 are all operational. (The assumption here is that if a node is operational for one link passing through it, it is operational for all links passing through it. In some instances, this assumption may not hold.) If receiver node 108 receives the test messages, but receiver node 104 does not, then it can be inferred that transmission failed along the network link defined by (link 132-node 112-link 124).

The process of characterizing a packet data network by multicasting test messages from a source node and analyzing the aggregate of test messages received by multiple receiver nodes is referred to herein as “multicast-based inference of network characteristics (MINC)”. Previous applications of MINC have characterized average packet loss. (Herein, “packet loss” will be referred to simply as “loss”.) See, for example, R. Caceres et al., “Multicast-Based Inference of Network-Internal Loss Characteristics,” IEEE Transactions in Information Theory, vol. 45, pp. 26-45, 2002. Average loss, however, provides only coarse characterization of network loss characteristics. It is well known, for example, that packet data networks are susceptible to noise (for example, electromagnetic interference), which may cause packets to be lost. Losses may be much greater during a noise burst than during quasi-quiet periods. It is also well known, for example, that traffic in packet data networks is bursty. Traffic congestion may cause packets to be lost. Losses may be much greater during heavy traffic load than during light traffic load. Simple average values of loss, therefore, do not adequately capture network characteristics. Advantageous procedures for MINC described herein expand the range of network characterization to include temporal loss characteristics and temporal delay characteristics of packet data networks. Herein, “temporal loss characteristics” refers to values of network loss as a function of time. Examples of temporal loss characteristics are discussed below.

Advantageous procedures for MINC are illustrated herein for packet data networks with a tree topology. FIG. 2 shows a graphical representation of a network viewed as a logical multicast tree T 200 comprising a set V of nodes 202-216, V={0, k, b-g), and a set L of links 218-230, L={link k, link b-link g}. In the tree model, node 0202 is the root node; node k 204, node b 206, and node c 208 are branch nodes; and node d 210-node g 216 are leaf nodes. Herein, a branch node in a tree model, as illustrated in FIG. 2, is equivalent to an intermediate node in a network architecture model, as illustrated in FIG. 1. Herein, the following genealogical terminology is also used: Node k 204 is a “child” of root node 0202; node b 206 and node c 208 are “children” of node k 204; node d 210 and node e 212 are children of node b 206; and node f 214 and node g 216 are children of node c 208. Other examples of genealogical terminology used herein include: node b 206 is the “father” of node d 210; and node b 206, node k 204, and root node 0202 are all “ancestors” of node d 210.

In the tree model illustrated in FIG. 2, the network topology is known to the end user. One skilled in the art may develop other embodiments which apply to networks in which the network topology is not a priori known to the end user. See, for example, N. G. Duffield et al., “Multicast Topology Inference from Measured End-to-End Loss,” IEEE Transactions in Information Theory, vol. 48, pp. 26-45, 2002. In a tree model, a non-root node node l receives messages from one and only one node, referred to as the unique father node f(l) of node l. One skilled in the art may develop other embodiments to characterize networks in which a specific node may have more than one father node. See, for example, T. Bu et al., “Network Tomography on General Topologies,” Proceedings ACM Sigmetrics 2002, Marina Del Rey, Calif., Jun. 15-19, 2002.

In MINC, test messages are multicast from a single source node to multiple destination nodes, which are the receiver nodes under test. In FIG. 2, the single source node is root node 0202, and the destination nodes are leaf nodes node d 210-node g 216. An end user has access to, and control of, source and receiver nodes. Herein, “test messages” may also be referred to as “probe messages”. To simplify the terminology further, “probe messages” may also be referred to as “probes”. In a multicast transmission, a probe is replicated at branch nodes. Separate copies are then forwarded to other branch nodes and to the destination nodes. In an example, packet data network 102 comprises an Internet Protocol (IP) network. A probe comprises one or more packets in which the source IP address is that of source node root node 0202, and the destination IP addresses are those of leaf nodes node d 210-node g 216. Typically, the IP addresses of node d 210-node g 216 are defined elements in a multicast group, which, for example, may be a range of addresses in a multicast subnet.

As shown in FIG. 2, probe i 232, where i is an integer 1, 2, 3 . . . , is transmitted from the source node root node 0202 to branch node node k 204, which then transmits a copy of the probe, shown in the figure as 234, to branch node node b 206. Branch node node k 204 transmits another copy of the probe, shown in the figure as 238, to branch node node c 208. Herein, the term “probe” comprises both the original probe transmitted from the source node root node 0202, and copies of the probe transmitted to branch nodes and destination nodes. In one example, the network parameter under test is loss (within a specified time interval). In FIG. 2, probes which are successfully transmitted are indicated as circles, 232-240. Probes which are lost are indicated as squares, 242 and 244. In this example, the probe reaches branch nodes, node k 204, node b 206, and node c 208. The probe further reaches destination nodes, node d 210 and node f 214; but the probe does not reach destination nodes, node e 212 and node g 216. A series of probes is used to measure the time dependence of network parameters. Note that the time interval between consecutive probes does not need to be constant. Herein, in measurements of loss, a destination node records (also referred to as “observes”) the “arrival” of a probe message. If a probe message does not arrive at a destination node within a user-defined interval, the destination node declares the probe message to be lost.

In embodiments in which the network parameter under test is loss, the multicast process is characterized by node states and link processes. The source node root node 0202 transmits a discrete series of probes probe i, where the index i is an integer 1, 2, 3 . . . . The node state X_l(i) indicates whether probe i has arrived at node l. The value X_l(i)=1 indicates that probe i has arrived at node l. The value X_k(i)=0 indicates that probe i has not arrived at node l, and has therefore been lost. In FIG. 2, probe i successfully arrived at node k 204-node d 210, and node f 214. The probe did not arrive at node e 212 and node g 216. Therefore, X_k(i)=X_b(i)=X_c(i)=X_d(i)=X_f(i)=1; and X_e(i)=X_g(i)=0.

The link process Z_l(i) indicates whether link l is capable of transmission during the interval in which probe i would attempt to reach node l, assuming that probe i were present at the father node f(l). The value Z_l(i)=1 indicates that the link is capable of transmission. The value Z_l(i)=0 indicates that the link is not capable of transmission. If node r is a destination node which is a receiver node under test, then X_r(i) provides loss statistics on the end-to-end path connecting source node root node 0 to node R. The aggregate data collected from a set of receiver nodes {node R} characterizes the set of end-to-end paths from source node root node 0 to each receiver node. A goal of MINC is to use the aggregate data to infer temporal characteristics of loss processes determining the link processes Z_l={Z_l(i)} along individual links internal to the network. Examples of model link loss processes Z_l={Z_l(i)} include Bernoulli, On-Off, and Stationary Ergodic Markov Process of Order r.

An example of MINC, in which the parameter under test is loss, is illustrated in FIG. 3 and described in the corresponding flow chart shown in FIG. 4. FIG. 3 shows a network schematic of a communications system comprising packet data network 302, source root node 0304, and destination nodes node d 306-node g 312. The circle 314 represents a probe. In the simple example described in the flow chart shown in FIG. 4, a series of four probes are multicast from source root node 0304 to destination nodes node d 306-node g 312. In actual test runs, the number of probes is much greater than 4. A test run may comprise 10,000 probes, for example.

In step 402, the probe index i is initialized to 1. In step 404, source root node 0304 multicasts probe i 314 to destination nodes node d 306-node g 312. In step 406, each individual destination node, node d 306-node g 312 collects data from probe i 314. In this instance, the data comprises records (observations) of whether the probe has arrived at a destination node.

The data is collected in a database which may be programmed in source root node 0304, destination nodes node d 306-node g 312, or on a separate host which may communicate with root node 0304 and destination nodes node d 306-node g 312. In step 408, the probe index i is incremented by 1. In step 410, the process returns to step 404, and steps 404-408 are iterated until four probes have been multicast. The process then continues to step 412, in which the probe data is outputted. In step 414, temporal link-loss characteristics of packet data network 302 are inferred from the probe data outputted in step 412. Details of the inference process are discussed below.

An example of data outputted in step 412 is shown in table 412A, which comprises columns (col.) 416-424 and rows 426-434. In row 426, the column headings indicate probe index i col. 416 and destination nodes node d col. 418-node g col. 424. Column 416, rows 428-434, track the probes, probe i,i=1-4. The entries in rows 428-434, col. 418-424, track the set of node states X_l={X_l(i)}, where l=d-g and i=1-4. A node state has the value 1 if the probe arrived (was received), and the value 0 if the probe was lost (was not received).

The process illustrated in the flow chart shown in FIG. 4 is a discrete time loss model with the following loss dependence structure:

• • Spatial. A loss process on one link is independent of the loss process on any other link.
  • Temporal. In previous applications of MINC, a loss process within a single link is independent of time. In advantageous procedures described herein, this constraint is removed, and a larger range of network characteristics may be analyzed. A loss process on any link (except for the link starting from the root node 0)is stationary and ergodic. That is, within a link, packet losses may be correlated, with parameters that in general depend on the link.In examples discussed herein, the temporal characteristics under test comprise measurements of “pass-runs” and “loss-runs” across a link within packet data network B02. Herein, a “pass” refers to a probe which has been successfully transmitted across a link and arrives at a destination node. A “loss” refers to a probe which has been lost during transmission across a link. A “pass-run” refers to a consecutive sequence of passes delimited by a loss before the first pass of the pass-run and a loss after the last pass of the pass-run. Similarly, a “loss-run” refers to a consecutive sequence of losses delimited by a pass before the first loss of the loss-run and a pass after the last loss of the loss-run. As an example, the following sequential data may be collected: a pass-run of 5,000 probes; a loss-run of 2000 probes; a pass-run of 50,000 probes; and a loss-run of 100 probes.

As discussed above, average loss (in a specified time interval) does not provide adequate characterization of links. Examples of more detailed link-loss parameters include the mean length of a pass-run, the mean length of a loss-run, and the probability that the length of a pass-run or loss-run exceeds a specific value. As discussed above, advantageous procedures process the aggregate data recorded (collected) from probes received at the destination nodes to estimate the link-loss parameters of interest for individual links within the network. As discussed above, a “path” is an end-to-end network link connecting one node to another node. A path may comprise multiple links. Herein, “path passage” refers to successful transmission of a probe across a path. Herein, “link passage” refers to successful transmission of a probe across a link. Individual link passage characteristics are inferred from measured path passage characteristics. Below, a system of equations describing path passage characteristics as functions of link passage characteristics is first derived. The path passage characteristics are values which are calculated from the aggregate data. Solutions to the system of equations then yield the link passage characteristics. In some instances, the solutions are approximate, and the approximate solutions yield estimates of the link passage characteristics.

As an example, let P_k be a random variable taking the marginal distribution of a pass-run, then the mean pass-run length is:

$\begin{array}{cc}\begin{array}{c}E\left[P_k\right]=\sum _{j\ge 1}j\mathrm{Pr}\left[P_k=j\right]=\sum _{j\ge 1}\mathrm{Pr}\left[P_k\ge j\right]\\ =\frac{\mathrm{Pr}\left[Z_k\left(1\right)=1\right]}{\mathrm{Pr}\left[Z_k\left(1\right)=1\right]-\mathrm{Pr}\left[Z_k\left(0\right)=1,Z_k\left(1\right)=1\right]}\end{array}\mathrm{where}E[·]\mathrm{means}\mathrm{expected}\mathrm{value}\mathrm{and}\mathrm{Pr}[·]\mathrm{means}\mathrm{probability}\mathrm{of}[·]& \left(\mathrm{Eqn}1\right)\end{array}$

Similarly, values such as mean loss-run length, probability that a pass-run is greater than a specified value, and probability that a loss-run is greater than a specified value may be calculated.

Methods to estimate parameters of interest are described herein. The following parameters and functions are defined herein:

I={i₁, i₂, . . . i_s)} where s is an integer, s≧1,   (Eqn 2)

• • is a set of time indices at which probes are transmitted. It is a generalization of a simple sequence i=1, 2, 3 . . . . That is, the time indices do not need to be equally spaced or even contiguous. For example, data may be collected from probe 7, probe 10, and probe 27.
  • χ_l(I) is a pattern of probes, corresponding to the index set I, which survived to node l.

χ_l(I)={{X_l(i)}: X_l(i₁)=X_l(i₂)= . . . X_l(i_s)=1}  (Eqn 3)

• • _l(I) is a link passage mask, defined such that probes with indices in I, if present at node f(l), will pass to node l, where node f(l) is the father of node l.

_l(I)={{Z_l(i)}: Z_l(i₁)=Z_l(i₂)= . . . Z_l(i_s)=1}  (Eqn 4)

• • Link pattern passage probability is defined by

α_l(I)=Pr[ (I)]=Pr[χ_l(I)|χ_f(l)(I)]  (Eqn 5)

• • where Pr[·] means probability of [·]
  • That is, if a probe has reached the father node node f(l), α_l is the probability that the probe will reach node l across link l.
  • Path pattern probability is defined by

(I)=Pr[χ_l(I)]=α_l(I)(I)   (Eqn 6)

• • That is, is the probability that the probe has successfully reached node l across the full path from root node to node l.In this instance, the path pattern probability is equal to the product of the link pattern probabilities of the individual links from root node 0 to node l:

$\begin{array}{cc}{}_l\left(I\right)=\prod _{w\in a\left(l\right)}{\alpha }_w\left(I\right)\mathrm{where}a\left(l\right)\mathrm{is}\mathrm{the}\mathrm{set}\mathrm{of}\mathrm{ancestors}\mathrm{of}\mathrm{node}l.& \left(\mathrm{Eqn}7\right)\end{array}$

An example is discussed with respect to the tree model previously shown in FIG. 2. The source node is root node 0202. The source node multicasts a sequence of four probes. In Eqn 2, I={i₁=1, i₂=2, i₃=3, i₄=4). The receiver nodes which collect the probes are the destination leaf nodes, node d 210-node g 216. Consider node l=node b 206. In Eqn 4 and Eqn 5, the father node node f(b) of node b 206 is node k 204. In Eqn 7, the set of ancestors of node b 206 is a(b)=(node k 204, node 0202).

Assume that probe 1-probe 4 all arrive at node k 204. Then, in Eqn 3, χ_k(I)={X_k(1)=X_k(2)=X_k(3)=X_k(4)=1}. Further assume that probe 1, probe 2, and probe 3 all arrive at node b 206, but probe 4 is lost. In this instance, in Eqn 4, (I)={Z_b(1)=Z_b(2)=Z_b(3)=1}, and, at node b 206, χ_b(I)={X_b(1)=X_b(2)=X_b(3)=1}. In Eqn 5, the link pattern passage probability is α_b(I)=Pr[(I)=Pr[χ_b(I)|χ_k(I)]. Or, in terms of this simple example, if probe 1 arrives at node k 204, the probability of probe 1 arriving at node b is equal to the probability that the link passage probability across link b 220 is l. A similar analysis applies for the other probes and other nodes.

In Eqn 7, now consider node l=node d 210, one of the receiver nodes which collects data. The set of ancestors of node d, denoted above as α(d) in Eqn 7, comprises {node b 206, node k 204, root node 0202}. For a probe probe i, the probability of path passage (i) from root node 0202 to receiver node d 210 is equal to the product of the probability of link passage across link k 218×the probability of link passage across link b 220×the probability of link passage across link d C22. A similar analysis holds for the other receiver nodes, node e 212-node g 216. A goal of MINC is to use the data collected at receiver nodes, node d 210-node g 216, to estimate the link passage probabilities across links, link k 218-link g 230. In a more generalized example, a goal is to estimate the link pattern passage probability α_l(I) of arbitrary patterns for all internal links l. These can be extracted from Eqn. 6 if (I) is known for all l, l≠0. In general, solving a polynomial equation of order >1 is required.

In an embodiment, path passage probabilities are calculated as a function of link passage probabilities by a process of subtree partitioning, which results in lower order polynomial equations. For example, subtree partitioning may result in a linear equation instead of a quadratic equation. The underlying concept of subtree partitioning is illustrated for a binary tree in the example shown in FIG. 5. In the binary tree T J00, each branch node has two child nodes. Here, probes are multicast from a source node S 502 to branch node node k 504. Tree T 500 is partitioned into two subtrees, subtree T_k,1 506 and subtree T_k,2 508. Branch node node k 504 in tree T 500 is configured as the root node for each of the two subtrees, T_k,1 506 and T_k,2 508. Node k 504 has two child nodes, node d_l,1 510 and node d_l,2 512. Each child node has two child nodes of its own. Node d_l,1 510 is a branch node in subtree T_k,1 506, and node d_l,2 512 is a branch node in subtree T_k,2 508. In turn, branch node node d_l,1 510 has two child nodes, node R_1,1 514 and node R_2,1 516. Similarly, branch node node d_1,2 512 has two child nodes, node R_1,2 518 and node R_2,2 520. Here, node R_1,1 514, node R_2,1 516, node R_1,2 518 and node R_2,2 520 are receiver nodes which receive probes from source node S 502.

In an example for a binary tree with subtree partitioning, path passage probabilities are calculated as follows. The following parameters and functions are defined herein.

$\begin{array}{cc}{Y}_{k,c}\left(i\right)=\underset{j\in R_{k,c}}{}X_j\left(i\right),\mathrm{for}c=\left\{1,2\right\},\mathrm{where}c=1\mathrm{is}\mathrm{subtree}1,& \left(\mathrm{Eqn}8\right)\end{array}$

• • c=2 is subtree 2. Here, V denotes bitwise OR. Y_k,c(i) is a random variable.
  • For c=0, where c=0 refers to the unpartitioned tree,

Y _k,0(I)=Y_k,l(i)VY_k,2(i)   (Eqn 9)

• • Here Y_k,1(i)=1 if there exists a receiver in R_k,1 which receives the i-th probe (else 0). Similarly, Y_k,2(i)=1 if there exists a receiver in R_k,2 which receives the i-th probe (else 0). Y_k,0(I)=1 if at least one of {R_k,1, R_k,2} contains such a receiver (else 0).

γ_k,c(i)=Pr[Y_k,c(i)=1 , for c ε {0, 1, 2}  (Eqn 10)

Then, the values (i) of may be calculated as:

(i)=Pr[χ_k(i)=γ_k,0(i), for k ε R,   (Eqn 11)

• • where R is the set of receiver nodes.Let U be the set of non-root nodes, then U\R is the set of branch nodes (non-root, non-receiver). For k ε U\R, the following value is defined:

$\begin{array}{cc}{\beta }_{k,c}\left(i\right)=\mathrm{Pr}\left[Y_{k,c}\left(i\right)=1|X_k\left(i\right)\right]=\frac{{\gamma }_{k,c}\left(i\right)}{{}_k\left(i\right)}& \left(\mathrm{Eqn}12\right)\end{array}$

Then,

γ_k,0(i)=(i)β_k,0(i)   (Eqn 13)

γ_k,0(i)=(i){1−(1−γ_k,1(i)/(i))(1−γ_k,2(i)/(i))}  (Eqn 14)

Eqn 14 is linear in (i) and can be solved:

$\begin{array}{cc}{}_k\left(i\right)=\frac{{\gamma }_{k,1}\left(i\right){\gamma }_{k,2}\left(i\right)}{{\gamma }_{k,1}\left(i\right)+{\gamma }_{k,2}\left(i\right)-{\gamma }_{k,0}\left(i\right)}& \left(\mathrm{Eqn}15\right)\end{array}$

Summarizing Eqn 11 and Eqn 15:

$\begin{array}{cc}{}_k\left(i\right)={\gamma }_{k,0}\left(i\right),\mathrm{for}k\in R{}_k\left(i\right)=\frac{{\gamma }_{k,1}\left(i\right){\gamma }_{k,2}\left(i\right)}{{\gamma }_{k,1}\left(i\right)+{\gamma }_{k,2}\left(i\right)-{\gamma }_{k,0}\left(i\right)},\mathrm{for}k\in U/R& \left(\mathrm{Eqn}16\right)\end{array}$

If the network comprises an arbitrary tree, in which a branch node may have more than two child nodes, the corresponding equation for (i) is a polynomial equation of order |d_k|−1, where |d_k| is the number of children of node k. In an example, the order of the equation may be reduced (for example, from quadratic to linear) by a more generalized subtree partitioning procedure. An example is shown in FIG. 6, which shows a schematic of an arbitrary tree T 600. As in the binary tree model above, probes are multicast from a source node S 602 to branch node node k 604. In this instance, branch node node k 604 has four child nodes, 610-616. In the subtree partitioning procedure, branch node node k 604 is configured as the root node for two subtrees, denoted T_k,1 606 and T_k,2 608. The union of T_k,1 606 and T_k,2 608 is denoted T_k,0 634, which is the entire subtree under node k 604. Each of the four child nodes 610-616 is then allocated to one of the subtrees. In general, child nodes are allocated to subtrees such that the population of child nodes in each subtree is approximately equal. In some instances, depending on the tree structure, an exactly equal number of child nodes in each subtree may not be achievable (for example, if there are an odd number of child nodes). Herein, “approximately equal” means that the number of child nodes in each subtree are as close as the tree architecture permits in a specific instance. In this instance, child nodes 610 and 612 are allocated to subtree T_k,1 606, and child nodes 614 and 616 are allocated to subtree T_k,2 608. The child nodes are then indexed as node d_1,1 610, node d_2,1 612, node d_1,2 610, and node d_1,2 614. These child nodes then serve as branch nodes for leaf nodes. The leaf nodes in subtree T_k,1 606 are node R_1,1 618, node R_2,1 620, node R_3,1 622, and node R_4,1 624. Similarly, the leaf nodes in subtree T_k,2 608 are node R_1,2 626, node R_2,2 628, node R_3,2 630, and node R_4,2 632. The linear solutions for (i) shown in Eqn 15 hold for arbitrary trees.

As discussed above, Eqn 15 apply for a single probe i. Another parameter of interest is the joint probability of a probe pattern I. In an example in which subtree partitioning is used, this parameter is calculated as follows.

$\begin{array}{cc}{Y}_{k,c}\left(I\right)=\underset{h\in I}{}Y_{k,c}\left(h\right),\mathrm{for}\mathrm{the}\mathrm{two}\mathrm{subtrees}c=\left\{1,2\right\}& \left(\mathrm{Eqn}17\right)\end{array}$

Here denotes bitwise AND.

Y _k,0(I)=Y_k,1(I)Y_k,2(I) , where c=0   (Eqn 18)

• • refers to the unpartitioned tree.Then, corresponding to Eqn 16 for the single probe i

$\begin{array}{cc}{}_k\left(I\right)={\gamma }_{k,0}\left(I\right),\mathrm{for}k\in R{}_k\left(I\right)=\frac{{\gamma }_{k,1}\left(I\right){\gamma }_{k,2}\left(I\right)}{{\gamma }_{k,1}\left(I\right)+{\gamma }_{k,2}\left(I\right)-{\gamma }_{k,0}\left(I\right)}\mathrm{for}k\in U/R& \left(\mathrm{Eqn}19\right)\end{array}$

If subtree partitioning is not used, then the values corresponding to Eqn. 12 are

$\begin{array}{cc}{\beta }_{k,c}\left(I\right)={\gamma }_{k,c}\left(I\right)/{}_k\left(I\right),\mathrm{for}c\in \left\{0,1,\dots \left(d_k-1\right)\right\}& \left(\mathrm{Eqn}20\right)\\ {\beta }_{k,0}\left(I\right)=1-\prod _{c=1}^{d_k-1}\left(1-{\beta }_{k,c}\left(I\right)\right)& \left(\mathrm{Eqn}21\right)\end{array}$

The resulting equation for (I) is not linear, but a polynomial of order |d_k|−1. Subtree partitioning is advantageous because is Eqn 19 linear.

In the subtree partitioning schemes described above, all probes in I passed through node k and reached receivers via nodes all within a single subtree. These schemes do not capture cases in which probes reach receivers for each index in I in a distributed way across the two subtrees, T_k,1 and T_k,2. In a further example of subtree partitioning, this limitation is removed, and (I) may be derived from all trials which imply χ_k(I).

In one example, for I={i,i+1}, and l, m, n, o ε {0,1}:

[l]=Pr[X _k(i)=l]  (Eqn 22)

[lm]=Pr[X _k(i)=l, X_k(i)=(i+1)=m]  (Eqn 23)

$\begin{array}{cc}{Y}_{k,c}\left(i\right)=\underset{j\in R_{k,c}}{}X_j\left(i\right)& \left(\mathrm{Eqn}24\right)\end{array}$ γ_k,c(l)=Pr[Y_k,c=l], for c={1, 2}  (Eqn 25)

γ_k,c[lm]=Pr[Y_k,c(i)=l, Y_k,c(i+1)=m], for c={1, 2}  (Eqn 26)

$\begin{array}{cc}{\gamma }_k\left[\mathrm{lm},\mathrm{no}\right]=\mathrm{Pr}\left[Y_{k,1}\left(i\right)=l,Y_{k,1}\left(i+1\right)=m,Y_{k,2}\left(i\right)=n,Y_{k,2}\left(i+1\right)=o\right]& \left(\mathrm{Eqn}27\right)\end{array}$ β_k,c(l)=Pr[Y_k,c(i)=l|χ(i)], for c={1, 2}  (Eqn 28)

β_k,c(lm)=Pr[Y_k,c(i)=l, Y_k,c(i+1)=m |χ(i)], for c={1, 2}  (Eqn 29)

γ_k[11]=Pr[Y_k(i)=1, Y_k(i+1)=1]  (Eqn 30)

Then trials which imply χ_k(I) are

γ_k[11]=γ_k[10,01]+γ_k[01,10]+γ_k,1[11]+γ_k,2[11]−γ_k[11,11]  (Eqn 31)

where γ_k[10,01] and γ_k[01,10] capture those missed by Y_k,0 From the conditional independence of the two trees,

$\begin{array}{cc}{\gamma }_k\left[11\right]-{\gamma }_{k,1}\left[11\right]-{\gamma }_{k,2}\left[11\right]={\gamma }_{k,}\left[10,01\right]-{\gamma }_{k,}\left[01,10\right]-{\gamma }_{k,}\left[11,11\right]={}_k\left[11\right]\left({\beta }_{k,1}\left[10\right]{\beta }_{k,2}\left[01\right]+{\beta }_{k,1}\left[01\right]{\beta }_{k,2}\left[10\right]+{\beta }_{k,1}\left[11\right]{\beta }_{k,2}\left[11\right]\right)& \left(\mathrm{Eqn}32\right)\end{array}$

As before,

γ_k,c[11]=[11]β_k,1[11], for c={1, 2}  (Eqn 33)

therefore,

γ_k,c[lm]=[11]β_k,c[lm]+([1]-[11])γ_k,c[1]/[1]  (Eqn 34)

• • for [lm]=[01] or [10] and c={1, 2}The β_k,c(lm) can now be eliminated in Eqn 34 and the resulting quadratic equation for [11] may be solved.

For the above tree and subtree partitioning schemes, estimators for parameters of interest may be derived. From n trials, samples of variables Y_k,c(I) are collected for each I of interest. Values of γ_k,c(I) may then be estimated using the empirical frequencies:

$\begin{array}{cc}{\hat{\gamma }}_{k,c}\left(I\right)=\frac{\sum _{i=0}^{n-s-1}Y_{k,c}\left(i+I\right)}{n-s-1}\mathrm{where}s=I& \left(\mathrm{Eqn}35\right)\end{array}$

The values of {circumflex over (γ)}_k,c(I) are then used to define an estimator _k(I) for (i). In the case of subtree partitioning, this is done by substituting into the relevant equation for (i). Otherwise, the unique root in [0,1] of the polynomial is found numerically. Another parameter of interest, the link joint passage probabilities, is estimated by

$\begin{array}{cc}{\hat{\alpha }}_k\left(I\right)=\frac{{\widehat{}}_k\left(I\right)}{{\widehat{}}_{f\left(k\right)}\left(I\right)}& \left(\mathrm{Eqn}36\right)\end{array}$

The analysis above yields three categories of estimators, all of which work on arbitrary trees and arbitrary probe patterns I. These categories are defined herein. “General” , based on Eqn 21, applies to instances in which there is no subtree partitioning and in which k(i) is solved numerically if the tree is non-binary. “Subtree”, based on Eqn 19, applies to instances in which there is subtree partitioning. “Advanced subtree”, based on Eqn 32, yields a quadratic in (I) when using subtree partitioning.

In another embodiment, the parameter of interest is delay. In the examples discussed above, in which the parameter of interest was loss, the multicast process was characterized by node states and link processes. The node state X_l(i) indicated whether probe i had arrived at node l. The link process Z_l(i) indicated whether link l was capable of transmission during the interval in which probe i would have attempted to reach node l, assuming that probe i had been present at the father node f(l). For delay, the multicast process is characterized by two processes. The delay measurement process X_l(i) records the delay along link l. The delay is the difference between the time at which probe i is transmitted from the father node f(l) of node l (assuming probe i has reached f(l)) and the time at which it is received by node l. The link process Z_l(i) is the time delay process which determines the delay encountered by probe i during its transmission from f(l) to node l. In an embodiment, a series of probes, probe i, is transmitted from source node root node 0 to a receiver node node R. At receiver node node R, the total end-to-end path delay from source node root node 0 to node R is recorded. The aggregate data collected from a set of receiver nodes {node R} characterizes the set of end-to-end paths from source node root node 0 to each receiver node. Previous applications of MINC calculated average delays per unit time. See F. Lo Presti et al., “Multicast-based Inference of Network-Internal Delay Distributions,” IEEE/ACM Transactions on Networking, vol. 10(6), pp. 761-775, 2002. An advantageous application of MINC uses the aggregate data to infer temporal characteristics of delay processes determining the link processes Z_l={Z_l(i)} along individual links internal to the network. Examples of link delay processes include Bernoulli Scheme, Stationary Ergodic Semi-Markov Process, and Stationary Ergodic Semi-Markov Process of Order r.

In general, delay values are continuous values from 0 to ∞, (the value ∞ may be used to characterize a lost probe). In one procedure, link delay values are measured as discrete values, which are an integer number of bins with a bin width of q. The set of delay values is then

D={0, q, 2q, . . . , mq, ∞},   (Eqn 37)

• • where mq is a user-definable threshold value. Delays greater than mq are declared to be lost, and the delays are set to ∞.If q is normalized to 1 then the following set is defined:

={0, 1, 2, . . . , m, ∞}  (Eqn 38)

The discrete time discrete state delay process at link k is then {Z_k(i)} and Z_k(i) ε .

An example, in which the parameter under test is delay, is illustrated in the flow chart shown in FIG. 7, which refers to the network schematic previously shown in FIG. 3. A different sequence of probes probe j 316 is now multicast from source root node 0304 to destination nodes node d 306-node g 312. In this example, the probe index j is used to distinguish the delay measurements from the loss measurements (with index i) previously shown in the flowchart of FIG. 3.

In step 702, the probe index j is initialized to 1. In step 704, source root node 0304 multicasts probe j 316 to destination nodes node d 306-node g 312. In step 706, each individual destination node, node d 306-node g 312 collects data from probe j 316. In this instance, data comprises delay values computed from measured arrival times of probe j 316 at each individual destination node, node d 306-node g 312.

As discussed above, the data is collected in a database which may be programmed in source root node 0304, destination nodes node d 306-node g 312, or on a separate host which may communicate with root node 0304 and destination nodes node d 306-node g 312. In step 708, the probe index j is incremented by 1. In step 710, the process returns to step 704, and steps 704-708 are iterated until four probes have been multicast. The process then continues to step 712, in which the probe data is outputted. In step 714, temporal delay characteristics of packet data network 302 are inferred from the probe data outputted in step 712. Details of the inference process are discussed below.

An example of data outputted in step 712 is shown in table 712A, which comprises columns (col.) 716-724 and rows 726-734. In row E26, the column headings indicate probe index j col. 716 and destination nodes node d col. 718-node g col. 724. Column 716, rows 728-734, track the probes, probe j,j=1-4. The entries in rows 728-734, col. 718-724, track the delays between source root node 0304 and destination nodes node d 306-node g 312. In this example, the bin width q is set equal to 1, and the threshold value m is set equal to 150. For j=1, the delay times corresponding to destination nodes node d col. 718-node g col. 724, are, respectively, (1, 4, 6, 2). Similarly, for j=4, the delay times corresponding to destination nodes node d col. 718-node g col. 724, are, respectively, (20, ∞,∞, 150). Here, a value of ∞ indicates that the delay time was >150, and the probe was declared lost.

The delay measurement process at a node k is denoted

{X_k(i)): X_k(i)ε {0, 1, 2, . . ., m ∞},   (Eqn 39)

• • where is the genealogical level with respect to root of node k, =root.For probe i, then,

X _k(i)=Z_k(i)+X_f(k)(i)   (Eqn 40)

which states that the delay between root and node k is equal to the delay between root and f(k) and the incremental delay between f(k) and node k. The total delay from root to node k is then the sum of the delay processes over all the ancestor nodes of node k:

$\begin{array}{cc}{X}_k\left(i\right)=\sum _{j\in a\left(k\right)}Z_j\left(i\right),& \left(\mathrm{Eqn}41\right)\end{array}$

• • where α(k) is the set of ancestor nodes of node k.In addition, the following probability results:

$\begin{array}{cc}\mathrm{Pr}\left[X_k\left(i\right)=p|X_{f\left(k\right)}\left(i\right)=q\right]=\{\begin{array}{cc}0& \mathrm{for}p<q,\\ 1& \mathrm{for}p=q=\infty ,\\ \mathrm{Pr}\left[Z_k\left(i\right)=p-q\right]& \mathrm{otherwise}.\end{array}& \left(\mathrm{Eqn}42\right)\end{array}$

which states that if the delay from root to f(k) is the value q, then the probability that the delay from [root to node k]=p has three outcomes. If p<q, then the probability is obviously 0, otherwise the delay between f(k) and node k is negative (probe i arrives at node k before it arrives at f(k)). If q=∞, then the probability that p=∞ is obviously 1 since if probe i is lost at f(k) it continues to be lost at node k (probe i is not regenerated between f(k) and node k). Otherwise, the probability is the probability that the link delay process Z_k(i) has the value (p−q).

As in the examples described above for a loss process, an embodiment for a delay process is applied to instances with the following dependence structure:

• • Spatial. A delay process on one link is independent of the delay process on any other link.
  • Temporal. In previous applications of MINC, a delay process within a single link is independent of time. In advantageous procedures described herein, this constraint is removed, and a larger range of network characteristics may be analyzed. A delay process on any link is stationary and ergodic. That is, within a link, delays may be correlated, with parameters that in general depend on the link.

Packet delay on a specific network link is equal to the sum of a fixed delay and a variable delay. The fixed delay, for example, may be the minimum delay resulting from processing by network equipment (such as routers and switches) and transmission across physical links (such as fiber or metal cable). The minimum delay is characteristic of networks in which the traffic is low. As traffic increases, a variable delay is introduced. The variable delay, for example, may result from queues in buffers of routers and switches. It may also arise from re-routing of traffic during heavy congestion. In one process, delay is normalized by subtracting the fixed delay from the total delay. For example, for a specific link to a specific receiver, the fixed delay may be set equal to the minimum delay measured over a large number of samples at low traffic. If d_max is the maximum normalized delay measured over the set of receivers, then the threshold m for declaring a packet as lost, may for example, be set to

m=d _max /q, where q is the bin width.   (Eqn 43)

In general, a goal is to estimate the complete family of joint probabilities

Pr[Z _k(i_l)=d₁, Z_k(i₂)=d₂, . . . , Z_k(i_s)=d_s]  (Eqn 44)

• • for any set s≧1 probe indices I={i₁, i₂, . . . i_s} and d₁,d₂, . . . d_s ε Principal values of interest are run distributions and mean run lengths. Let L_k^H denote a random variable which indicates the length of runs of Z_k in a subset H, in which H satisfies a set of user-defined parameters, of full state space . Examples of H are given below. Then, the probability that L_k^H is greater than or equal to a value j is

$\begin{array}{cc}\mathrm{Pr}\left[L_k^H\ge j\right]=\frac{\begin{array}{c}\mathrm{Pr}\left[Z_k\left(j\right)\in H,\dots ,Z_k\left(1\right)\in H\right]-\\ \mathrm{Pr}\left[Z_k\left(j\right)\in H,\dots ,Z_k\left(0\right)\in H\right]\end{array}}{\mathrm{Pr}\left[Z_k\left(1\right)\in H\right]-\mathrm{Pr}\left[Z_k\left(0\right)\in H,Z_k\left(1\right)\in H\right]}& \left(\mathrm{Eqn}45\right)\end{array}$

The mean run length is

$\begin{array}{cc}{\mu }_k^H=E\left[L_k^H\right]=\frac{\mathrm{Pr}\left[Z_k\left(j\right)\in H\right]}{\mathrm{Pr}\left[Z_k\left(1\right)\in H\right]-\mathrm{Pr}\left[Z_k\left(0\right)\in H,Z_k\left(1\right)\in H\right]}& \left(\mathrm{Eqn}46\right)\end{array}$

In Eqn 46, μ_k^H is the ratio of the expected proportion of time spent in runs in the subset H (per unit time index) divided by the expected number of transitions into H (per unit time index). The mean run length of a delay state may be derived if the simplest joint probabilities, with respect to that state may be estimated:

• • for a single probe, Pr[Z_k(i)ε H],
  • for a successive pair of probes, Pr[Z_k(i)ε H, Z_k(i+1)ε H].The tail probabilities of runs in H, Pr[L_k^H≧j], can be obtained from the joint probabilities of the state H for one, two, j, and j+1 probes.

Eqn 45 may be used to partition the link states into two classes. States in subset H are referred to as “bad”. States in \H are referred to as “good”. For example, H may refer to states with a delay greater than a user-defined value d. In which case, μ_k^H is the mean duration of runs in which the delay is at least d.

As in the procedure for estimating temporal loss characteristics, in an embodiment for estimating temporal delay characteristics, the source at the root node multicasts a stream of n probes, and each receiver records the end-to-end delay that it observes. The transmission of probes may then be viewed as an experiment with n trials. The outcome of the i-th trial is the set of discretized source-to-receiver delays

{X_k(i), k ε R}, X_k(i)ε {0, 1, . . . , m, ∞}  (Eqn 47)

To calculate joint probabilities, the following values are defined herein.

I={i₁, i₂, . . . i_s}, as before I is a set of probe indexes,   (Eqn 48)

not necessarily contiguous

_k(I)=[X_k(i₁), X_k(i₂), . . . , X_k(i_s)] is a random vector   (Eqn 49)

_k(I)=[Z_k(i₁), Z_k(i₂), . . . , Z_k(i_s)] is a random vector   (Eqn 50)

, are delay vectors, and ≦ means d_j≦v_j for any j   (Eqn 51)

=[m, m, . . . , m]  (Eqn 52)

=[0, 0, . . . , 0]  (Eqn 53)

Then, the joint link probability is

α_k(I, )=Pr[_k(I)=, for ,≧,   (Eqn 54)

and the joint path passage probability is

$\begin{array}{cc}\begin{array}{c}{}_k\left(I,\right)=\mathrm{Pr}[{}_k\left(I\right)\\ =]\mathrm{for}\le ,\le \\ =\sum _{\le \le }a_k\left(I,\right){}_{f\left(k\right)}\left(-\right)\end{array}& \left(\mathrm{Eqn}55\right)\end{array}$

After the values _k(I, ), for all k ε U, for ≦≦, have been obtained, the following values are recursively deconvolved:

$\begin{array}{cc}\mathrm{For}=,a_k\left(I,\right)=\frac{{}_k\left(I,\right)}{{}_{f\left(k\right)}\left(I,\right)}& \left(\mathrm{Eqn}56\right)\\ \mathrm{For}<\le ,a_k\left(I,\right)=\frac{{}_k\left(I,\right)-\sum _{<\le }{}_{f\left(k\right)}\left(I,\right)a_k\left(I,-\right)}{{}_{f\left(k\right)}\left(I,\right)}& \left(\mathrm{Eq}n57\right)\end{array}$

For the case where ≦ does not hold (that is, at least one element of is ∞), α_k(I, ) is obtained using α_k(I, ), ≦, recursively using the α_k for smaller index sets. For example, for =[d₁=∞, d₂, . . . , d_s], then α_k(I, ) may be re-expressed as follows:

$\begin{array}{cc}{a}_k\left(I,\right)=a_k\left(\left\{i_2,\dots ,i_s\right),\left[d_2,\dots ,d_s\right]\right)-\sum _{v_1\le m}a_k\left(I,\left[v_1,d_2,\dots d_s\right]\right)& \left(\mathrm{Eqn}58\right)\end{array}$

For k ε U, path probabilities (I, ), ≦≦, are estimated by using the principle of subtree partition as follows. Consider branch node k 604 in the tree T 600 (FIG. 6). It is the root of the subtree T_k,0 634, which has receiver nodes R_k,0 where R_k,0 is the combined set of receiver nodes R_1,1 618 to R_4,1 624 and R_1,2 626 to R_4,2 632. The set of child subtrees of node k 604 are divided into two sets, corresponding to two virtual subtrees T_k,1 606 and T_k,2 608. Let j={0, 1, 2} be used to index quantities corresponding to subtrees T_k,0 634, T_k,1 606, and T_k,2 608, respectively. For a set of probe indices I, the following random vectors and probabilities are defined:

$\begin{array}{cc}{Y}_{k,j}\left(i\right)=\underset{r\in R_{k,j}}{\min }X_r\left(i\right),{}_{k,j}\left(I\right)=\left[Y_{k,j}\left(i_I\right),\dots ,Y_{k,j}\left(i_s\right)\right]{\tilde{Y}}_{k,j}\left(i,d\right)=\{\begin{array}{cc}1& \mathrm{if}Y_{k,j}\left(i\right)-X_k\left(i\right)\le d\\ 0& \mathrm{if}Y_{k,j}\left(i\right)-X_k\left(i\right)>d\end{array}& \left(\mathrm{Eqn}59\right)\\ {\widetilde{}}_{k,j}\left(I,\right)=\left[{\tilde{Y}}_{k,j}\left(i_I,d_I\right),\dots ,{\tilde{Y}}_{k,j}\left(i_s,d_s\right)\right]& \left(\mathrm{Eqn}60\right)\\ {\gamma }_{k,j}\left(I,\right)=\mathrm{Pr}\left[{}_{k,j}\left(I\right)\le \right]{\beta }_{k,j}\left(I,,\right)=\mathrm{Pr}\left[{\widetilde{}}_{k,j}\left(I,\right)=\right]& \left(\mathrm{Eqn}61\right)\end{array}$

where ε{0, 1}^|k|. γ_k,j(I, ) is the probability that for each probe index i_l ε I, the minimum delay on any path from source S to receivers in R_k,j, does not exceed d_l ε . On the other hand, β_k,j(I, , ) is the probability that, for each probe index i_l ε I, the minimum delay on any path from node k 604 to receivers in R_k,j is either ≦d_l or >d_l ε depending on whether b_l ε is 1 or 0. Let =[1, . . . ,1]. Then, , β, and γ are related by the following convolution:

$\begin{array}{cc}{\gamma }_{k,j}\left(I,\right)=\sum _{\le \le }{}_k\left(I,\right){\beta }_{k,j}\left(I,-,\right)& \left(\mathrm{Eqn}62\right)\end{array}$

In order to recover (I, )'s from the γ_k(I, )'s which are directly observable from receiver data, the following two properties of β's are used.

• Property 1. This property gives the relationship between β_k,0 and {β_k,1, β_k,2} of the virtual subtrees.

$\begin{array}{cc}{\beta }_{k,0}\left(I,,\right)=1-\prod _{j=1}^2(1-\left(1-{\beta }_{k,j}\left(I,,\right)\right)+1_{I>1}\left(\sum _{\underset{s.t.{}_1{}_2=}{\left\{{}_1\ne ,{}_2\ne \right\}}}\prod _{j=1}^2{\beta }_{k,j}\left(I,,{}_j\right)\right)& \left(\mathrm{Eqn}63\right)\end{array}$

• Property 2. (Recursion over index sets with =) This property allows β_k,j(I, , ≠) to be expressed in terms of β_k,j(I′, , )'s, where I′ I. For instance, if =[b₁=0, b₂, . . . , b_s], and I′={i₂, . . . , i_s}, =[b₂, . . . , b_s], =[d₂, . . . , d_s], then

β_k,j(I, )=β_k,j(I′, , )−β_k,j(I, [1, b₂, . . . , b_s])

which eliminates the 0 at i_l. The above can be applied recursively to eliminate all zeroes, resulting in terms of the form β_k,j(I′, , ), I′ I, |I|−()≦|I′|≦|I|, where () denotes the number of zeroes in . In general

β_k,j(I, , ≠)=(−1)^z(B)β_k,j(I, )+δ_k,j(I, )   (Eqn 64)

where δ_k,j(I, ) is the appropriate summation of β_k,j's for index sets I′ c I. For example, if I={1, 2}, =[0, 1], =[d₁, d₂], then,

β_k,j(I, )=−β_k,j(I, )+β_k,j({2}, [d₂], )

Hence, δ_k,j(I, )=β_k,j({2}, [d₂], ).By using Equation 64 in Equation 63, terms of the type _j≠ can be removed, leaving only terms of type _j=, giving:

$\begin{array}{cc}{\beta }_{k,0}\left(I,,\right)=1-\prod _{j=1}^2\left(1-{\beta }_{k,j}\left(I,,\right)\right)+1_{I>1}\left(\sum _{\underset{s.t.{}_1{}_2=}{\left\{{}_1\ne ,{}_2\ne \right\}}}\prod _{j=1}^2\left\{{\left(-1\right)}^{z\left(B\right)}{\beta }_{k,j}\left(I,,\right)+{\delta }_{k,j}\left(I,,\right)\right\}\right)& \left(\mathrm{Eqn}65\right)\end{array}$

Using Equation 65 and Equation 62, the desired path probabilities for node k 604, _k(I, ), ≦≦, can be computed using the observables γ_k,j(I, ). The recovery of (I, ) from the above equations involves two levels of recursion: (i) over delay vectors, arising due to convolution, (ii) over index sets arising due to summation term involving δ in Equation 65. The δ(I, .,.) only contains terms involving I′ c I and therefore does not contain (I, .). Thus estimation can be performed recursively starting from I={i} when the summation term with δ vanishes and = when the convolution vanishes. Each step of recursion involves solving a quadratic equation in the unknown _k.

The computation of (I, ) for pairs of consecutive probes i.e. I={1, 2}, proceeds as follows (I={1,2} is same as I={i+1}). Due to recursion over index sets, the case of I={1} is considered first.

• Single probes I={1}: The base case of recursion occurs for I={i} and =[0]. To simplify notation, the following are dropped: the index set I, =, and vector notation for delays. For example, β_k,j(I, [d₁], )=β_k,j(d_l). Writing out Equation 65 and Equation 62,

γ_k,j(0)=_k(0)β_k,j(0)

β_k,0(0)=1−(1−β_k,1(0))(1−β_k,2(0))   (Eqn 66)

from which _k(0) is recovered by solving a linear equation as

${}_k\left(0\right)=\frac{{\gamma }_{k,1}\left(0\right){\gamma }_{k,2}\left(0\right)}{{\gamma }_{k,1}\left(0\right)+{\gamma }_{k,2}\left(0\right)-{\gamma }_{k,0}\left(0\right)}$

Substituting back (0) gives the β_k,j(0)'s for use in the next step. Assuming that and β_k,j's have been computed ∀ v₁<d₁, (d₁) is recovered using Equation 65 and Equation 62 as

${\gamma }_{k,j}\left(d_1\right)={}_k\left(0\right){\beta }_{k,j}\left(d_1\right)+{}_k\left(d_1\right){\beta }_{k,j}\left(0\right)+\sum _{0<v_1<d_1}{}_k\left(v_l\right){\beta }_{k,j}\left(d_1-v_1\right)$ $\underset{*}{{\beta }_{k,0}\left(d_1\right)}=1-\underset{*}{\left(1-{\beta }_{k,1}\left(d_1\right)\right)}\underset{*}{\left(1-{\beta }_{k,2}\left(d_1\right)\right)}$

The unknown terms are marked by a “*”. (d₁) is recovered by solving a quadratic equation and substituting back (d₁) gives β_k,j(d₁)'s.

• Pairs of consecutive probes I={1, 2}: Again, to simplify the notation, the following are dropped: the index set I, = and vector notation for delays. For example, β_k,j(I, [d₁, d₂], )=β_k,j(d₁, d₂). The estimation proceeds from delay vector [0, 0] until [m, m]. Assuming that _k and β_j's have been computed for the set

{[v₁, v₂]: v₁≦d₁, v₂≦d₂}\{[d₁, d₂]},

(d₁, d₂) is recovered as follows. Equation 65 and Equation 62 are expanded.

$\begin{array}{cc}{\gamma }_{k,j}\left(d_1,d_2\right)={}_k\left(0,0\right){\beta }_{k,j}\left(d_1,d_2\right)+{}_k\left(d_1,d_2\right){\beta }_{k,j}\left(0,0\right)+\sum _{\begin{array}{c}{v}_1\le d_1\\ \left(v_1,v_2\right)\ne \left(0,0\right),\end{array}}\sum _{\begin{array}{c}{v}_2\le d_2\\ \left(v_1,v_2\right)\ne \left(d_1,d_2\right)\end{array}}{}_k\left(v_1,v_2\right){\beta }_{k,j}\left(d_1-v_1,d_2-v_2\right)\underset{*}{{\beta }_{k,0}\left(d_1,d_2\right)}=1-(1-{\beta }_{k,1}\underset{*}{\left(d_1,d_2\right))}(1-{\beta }_{k,2}\underset{*}{\left(d_1,d_2\right))}+\left(-{\beta }_{k,1}\underset{*}{\left(d_1,d_2\right)}+{\beta }_{k,1}\left(d_2\right)\right)\left(-{\beta }_{k,2}\underset{*}{\left(d_1,d_2\right)}+{\beta }_{k,2}\left(d_1\right)\right)+\left(-{\beta }_{k,1}\underset{*}{\left(d_1,d_2\right)}+{\beta }_{k,1}\left(d_1\right)\right)\left(-{\beta }_{k,1}\underset{*}{\left(d_1,d_2\right)}+{\beta }_{k,1}\left(d_1\right)\right)\left(-{\beta }_{k,2}\underset{*}{\left(d_1,d_2\right)}+{\beta }_{k,2}\left(d_2\right)\right)& \left(\mathrm{Eqn}67\right)\end{array}$

The unknown terms are marked by a “*” and (d₁, d₂) is obtained by solving a quadratic equation.

The parameter γ_k,j(I, ) may be estimated using the empirical frequencies as:

$\begin{array}{cc}{\hat{\gamma }}_{k,j}\left(I,\right)=\frac{\sum _{i=0}^{n-I-1}{{}_{k,j}\left(i+1\right)\le }_{}}{n-I-1}& \left(\mathrm{Eqn}68\right)\end{array}$

The parameter γ_k,j(I, ) is then used to define an estimator (I, ) for (I, ). The parameter {circumflex over (α)}(I, ) is then recursively deconvolved.The mean run length of delay state p ε is estimated using joint probabilities of single and two packet indices as

$\begin{array}{cc}{\hat{\mu }}_k^p=\frac{{\hat{a}}_k\left(p\right)}{{\hat{a}}_k\left(p\right)-{\hat{a}}_k\left(p,p\right)}\mathrm{where}a_k\left(p\right)=a_k\left(\left\{i_1\right\},\left[p\right]\right)\mathrm{and}a_k\left(p,p\right)=a_k\left(\left\{i_1,i_2\right\},\left[p,p\right]\right)& \left(\mathrm{Eqn}69\right)\end{array}$

When delay states are classified into bad H and good G=H states, the mean run length of bad state is estimated using the joint probabilities of single and two packet indices as:

$\begin{array}{cc}{\hat{\mu }}_k^H=\frac{\sum _{p\in H}{\hat{a}}_k\left(p\right)}{\sum _{p\in H}{\hat{a}}_k\left(p\right)-\sum _{p_1\in H}\sum _{p_2\in H}{\hat{a}}_k\left(p_1,p_2\right)}& \left(\mathrm{Eqn}70\right)\end{array}$

A similar expression is used to estimate {circumflex over (μ)}_k^G.

One embodiment of a network characterization system which performs multicast-based inference may be implemented using a computer. As shown in FIG. 8, computer 802 may be any type of well-known computer comprising a central processing unit (CPU) 806, memory 804, data storage 808, and user input/output interface 810. Data storage 808 may comprise a hard drive or non-volatile memory. User input/output interface 810 may comprise a connection to a user input device 822, such as a keyboard or mouse. As is well known, a computer operates under control of computer software which defines the overall operation of the computer and applications. CPU 806 controls the overall operation of the computer and applications by executing computer program instructions which define the overall operation and applications. The computer program instructions may be stored in data storage 808 and loaded into memory 804 when execution of the program instructions is desired. Computer 802 may further comprise a signal interface 812 and a video display interface 816. Signal interface 812 may transform incoming signals, such as from a network analyzer, to signals capable of being processed by CPU 806. Video display interface 816 may transform signals from CPU 806 to signals which may drive video display 820. Computer 802 may further comprise one or more network interfaces. For example, communications network interface 814 may comprise a connection to an Internet Protocol (IP) communications network 826, which may transport user traffic. In one embodiment, the network characterization system further comprises nodes within communications network 826. These nodes may serve as a source node and a set of receiver nodes. As another example, test network interface 818 may comprise a connection to an IP test network 824, which may transport dedicated test traffic. Computer 802 may further comprise multiple communications network interfaces and multiple test network interfaces. In some instances, the communications network 826 and the test network 824 may be the same. Computers are well known in the art and will not be described in detail herein.

The foregoing Detailed Description is to be understood as being in every respect illustrative and exemplary, but not restrictive, and the scope of the invention disclosed herein is not to be determined from the Detailed Description, but rather from the claims as interpreted according to the full breadth permitted by the patent laws. It is to be understood that the embodiments shown and described herein are only illustrative of the principles of the present invention and that various modifications may be implemented by those skilled in the art without departing from the scope and spirit of the invention. Those skilled in the art could implement various other feature combinations without departing from the scope and spirit of the invention.

Claims

1. A method for calculating a temporal delay characteristic of a packet data network comprising a source node, a plurality of receiver nodes, and a plurality of paths, comprising the steps of: recording delays of a plurality of multicast probe messages; and,calculating said temporal delay characteristic from said recorded delays.

2. The method of claim 1 wherein each path in said plurality of paths connects said source node with one of said plurality of receiver nodes and wherein each of said paths comprises at least one link.

3. The method of claim 2 wherein said temporal loss characteristic of said packet data network comprises the temporal loss characteristic of at least one link.

4. The method of claim 1 wherein said multicast probe messages are transmitted from said source node.

5. The method of claim 1 wherein said step of recording delays further comprises the step of recording delays at each of said plurality of receiver nodes.

6. The method of claim 1 wherein said temporal delay characteristic comprises the number of probe messages per unit time having a delay between a first value and a second value.

7. The method of claim 1 wherein said temporal delay characteristic comprises a number of probe messages per unit time having a delay less than a value.

8. The method of claim 1 wherein said temporal delay characteristic comprises a number of probe messages per unit time having a delay greater than a value.

9. The method of claim 1 wherein a probe message from said plurality of probe messages is declared lost if the delay exceeds a threshold value.

10. The method of claim 1 wherein the topology of said packet data network is a binary tree.

11. The method of claim 10 wherein said binary tree is partitioned into two subtrees.

12. The method of claim 1 wherein the topology of said packet data network is an arbitrary tree.

13. The method of claim 12 wherein said arbitrary tree is partitioned into two subtrees.

14. A network characterization system for calculating a temporal delay characteristic of a packet data network comprising a source node, a plurality of receiver nodes, and a plurality of paths wherein each path in said plurality of paths connects said source node with one of said plurality of receiver nodes and wherein each of said paths comprises at least one link, said network characterization system comprising: means for recording delays of a plurality of multicast probe messages; and,means for calculating said delay characteristic from said recorded arrivals.

15. The network characterization system of claim 14 wherein said means for calculating said temporal delay characteristic from said recorded delays further comprises: means for calculating said temporal delay characteristic of at least one link.

16. The network characterization system of claim 14, further comprising: means for multicasting probe messages from said source node.

17. The network characterization system of claim 14 wherein said means for recording delays of a plurality of multicast probe messages further comprises: means for recording delays of a plurality of multicast probe messages at each of said plurality of receiver nodes.

18. The network characterization system of claim 14 wherein said means for calculating said temporal delay characteristic from said recorded delays further comprises means for calculating at least one of an average delay per unit time, a number of probe messages per unit time with delays greater than a first value, a number of probe messages per unit time with delays less than a second value, and a number of probe messages per unit time with delays greater than a third value and less than a fourth value.

19. The network characterization system of claim 14, further comprising: means for partitioning a binary tree into two subtrees.

20. The network characterization system of claim 14, further comprising: means for partitioning an arbitrary tree into two subtrees.

21. A computer readable medium storing computer program instructions for calculating a temporal delay characteristic of a packet data network comprising a source node, a plurality of receiver nodes, and a plurality of paths wherein each path in said plurality of paths connects said source node with one of said plurality of receiver nodes and wherein each of said paths comprises at least one link, said computer program instructions defining the steps of: recording delays of a plurality of multicast probe messages; and,calculating said temporal delay characteristic from said recorded delays.

22. The computer readable medium of claim 21 wherein said computer program instructions defining the step of calculating said delay characteristic from said recorded delays further comprise computer program instructions defining the step of: calculating said temporal delay characteristic of at least one link.

23. The computer readable medium of claim 21 wherein said computer program instructions further comprise computer program instructions defining the step of: transmitting multicast probe messages from said source node.

24. The computer readable medium of claim 21 wherein said computer program instructions defining the step of recording delays of a plurality of multicast probe messages further comprise computer program instructions defining the step of: recording delays of a plurality of multicast probe messages at each of said plurality of receiver nodes.

25. The computer readable medium of claim 21 wherein said computer program instructions defining the step of calculating said temporal delay characteristic from said recorded delays further comprise computer program instructions defining the step of: calculating at least one of an average delay per unit time, a number of probe messages per unit time with delays greater than a first value, a number of probe messages per unit time with delays less than a second value, and a number of probe messages per unit time with delay values between a third value and a fourth value.