nal transmission among CubeSats. Then, the signals to each
CV are detected by signal separation, and the directions of
the propagation paths (i.e., AoDs and AoAs) are estimated.
Finally the bi-directional LOS vectors between each CV
and DV with respect to the corresponding frame bodies are
deduced.
A. Array Signal Transmission and Detection
First, all of the DVs are allowed to transmit asyn-
chronously, and thus, the complicated synchronization and
coordination among the CubeSats are avoided. The multi-
CubeSat detection (i.e., separation of the signals from dif-
ferent DVs) is realized by the CDMA technology. Each DV
is assigned a predefined spread spectrum sequence, i.e., pe-
sudonoise pseudo-noise (PN) sequence with the length of
L bits. Readers are referred to [39] and [40] for a detailed
explanation of the characteristics and generation of PN-
sequences. The spreading sequence assigned to the kth DV
is denoted by a
k
(τ ) and expressed as
a
k
(τ ) =
L
l=1
b
k,l
rect
b
(
τ − lT
b
)
,b
k,l
∈
{
+1, −1
}
(1)
where rect
b
(
τ
)
denotes a rectangular pulse with the dura-
tion of T
b
and power of P
b
=
T
b
0
rect
2
b
(τ )dτ. The duration
of the spreading sequence is T
s
= LT
b
. Please note that
the spread spectrum sequences of the K DVs are selected
such that their autocorrelation functions are pulse functions
and they are orthogonal to each other. The cross-correlation
functions of the spreading sequences are
a
k
1
(τ ),a
k
2
(τ )=a
k
1
(τ ) ∗ a
∗
k
2
(−τ ) =
ε
k
1
,k
2
,k
1
= k
2
LP
b
δ(τ ),k
1
= k
2
(2)
where < ·, · > denotes the sliding inner-production,
ε
k
1
,k
2
≈ 0 is a very small pesudorandom number, and δ(τ )
is the delta function.
T
S
0
a
2
k
1
(τ )dτ = LP
b
is the peak value
of the autocorrelation f unction of a
k
1
(τ ).
A DV (transmitter) modulates a carrier with its spread-
ing sequence circularly (e.g., via BPSK). The modulated
signal is transmitted on its M antennas in a time division
multiplexing manner. By using a microwave switch, the an-
tennas transmit the modulated signal one-by-one, and each
transmits at least 2N spreading sequences. Meanwhile, con-
trolled by a microwave switch, the N antennas of the array
on a CV (receiver) receive the signals from the DVs se-
quentially. Each antenna keeps receiving for the duration of
2T
s
such that it can capture at least one complete spreading
sequence from a DV. After all the antennas on the DVs have
transmitted their spreading sequences, every pair of the DV
and CV antennas have transferred a complete spreading
sequence. The MIMO links between a DV and CV are il-
lustrated in Fig. 1.
Determined by the relative positions and attitudes of
the CubeSats, the signal over the link from the mth (m =
1, 2,...,M) antenna on DV
k
to the nth (n = 1, 2,...,N)
antenna on CV
i
is
y
k,m,i,n
(τ ) = g
k,i
t
k,m
(φ
k,i
,θ
k,i
)r
i,n
(φ
k,i
,θ
k,i
)a
k
(τ − τ
k,i
)
(3)
where
g
k,i
is the free-space path loss depending on the dis-
tance between DV
k
and CV
i
and a
k
(τ − τ
k,i
) is the circular
shift of a
k
(τ )byτ
k,i
, the propagation delay from DV
k
to
CV
i
.In(3),t
k,m
(φ
k,i
,θ
k,i
) is the complex gain of the mth
Tx antenna, including the antenna radiation pattern at the
AoD of (φ
k,i
,θ
k,i
) and the phase shift with respect to the
reference antenna for its position. Similarly, r
i,n
(φ
k,i
,θ
k,i
)
is the complex gain of the nth Rx antenna, including the
antenna pattern at the AoA of (φ
k,i
,θ
k,i
) and the phase shift
with respect to the reference antenna.
As all of the K DVs are broadcasting simultaneously,
the signals received on the CVs’ antennas are the superpo-
sition of all the DVs’ signals. The received signal on the
nth antenna of CV
i
is
y
i,n
(τ ) =
K
k=1
y
k,m,i,n
(τ ) + w
i,n
(τ )(4)
where
y
k,m,i,n
(τ ) is the signal from the mth antenna on
DV
k
as given in (3) and w
i,n
(τ ) is the noise through the nth
receiving chain on CV
i
.
We can then separate and detect each DV’s signal based
on the orthogonality among the spreading sequences. Each
CV performs sliding correlation between the received
signal
y
i,n
(τ ) and the spreading sequence a
k
(τ ) in the delay
domain to extract the signal from DV
k
. The result is given
in (5) shown at the bottom of the next page. By plugging (2)
into (5), we can obtain (6) shown at the bottom of the next
page, where A
k,m,i,n
LP
b
g
k,i
t
k,m
(φ
k,i
,θ
k,i
)r
i,n
(φ
k,i
,θ
k,i
)
is the complex amplitude of the delta func-
tion in the cross-correlation function and ε
i,n
K
k
=1,k
=k
g
k
,i
t
k
,m
(φ
k
,i
,θ
k
,i
)r
i,n
(φ
k
,i
,θ
k
,i
)ε
k,k
≈ 0isthe
inter-CubeSat interference received on the nth antenna
of CV
i
. Because of the orthogonality property of the
spreading sequences, the cross-correlation result in (6)
generates a peak at the excess delay of τ
k,i
with the
complex amplitude of A
k,m,i,n
, which is called the link
impulse response (LIR). Please note that the amplitude of
the LIR is amplified by LP
b
and the interference from
other DVs is removed by the sliding correlation, which is
the spread spectrum gain. Therefore, employing spread
spectrum sequences provides the advantages of mitigating
inter-CubeSat interference and high signal-to-interference
ratio (SIR).
Although the following derivation is based on an ideal
signal model with the Gaussian noise ignored for easy pre-
sentation, important insights can be obtained. As shown in
Fig. 1, the spreading sequence from a DV is transmitted
over the MIMO links to both CVs. CV
i
uses the spread-
ing sequence for DV
k
, a
k
(τ ), to perform sliding correlation
with the received signals on its N antennas,
y
i,n
(τ )for
n = 1, 2,...,N. Thus, it obtains the LIR between each
pair of antennas on DV
k
and CV
i
, A
k,m,i,n
. The LIRs are
organized in an M × N matrix as given in (7) shown at the
WANG ET AL.: MULTI-CUBESAT RELATIVE POSITION AND ATTITUDE DETERMINATION 3381