Intermediate Frequency Spectral Placement of Bandpass Sampled Signals

Information

  • Patent Application
  • 20160043749
  • Publication Number
    20160043749
  • Date Filed
    August 05, 2014
    10 years ago
  • Date Published
    February 11, 2016
    8 years ago
Abstract
Systems and methods for bandpass sampling of a continuous analog signal according to several embodiments can use a linear relationship between sampling frequency and intermediate frequency
Description
FIELD OF THE INVENTION

The present invention pertains generally to methods for designing radios. More specifically, the present invention according to several embodiments covers a computationally efficient procedure for designing radios that use bandpass sampling to accomplish analog to digital conversion of signals. The present invention is particularly, but not exclusively, useful for designing radios that accomplish bandpass sampling using a linear relationship between intermediate frequency and sampling frequency design parameters, in order to simplify the radio design process.


BACKGROUND OF THE INVENTION

An analog signal is continuous in time, and it can often be necessary to convert the analog signal to a flow of digital values (i.e, a digital signal), which can typically be further processed more easily than the original analog signal. An analog to digital converter (ADC) can be used to accomplish this conversion. An ADC can operate by sampling the value of the input at discrete intervals in time. Provided that the input is sampled at a rate above the Nyquist rate (defined as twice the highest frequency component of the signal), then the digital signal is an accurate representation of the analog signal, and the analog signal can be reconstructed from the digital signal.


However, if the signal is sampled at a rate less than the Nyquist rate, the signal can be incorrectly converted to a signal that is not really representative of the original analog signal. This is commonly referred to as aliasing. Aliasing can occur because instantaneously sampling a function at two or fewer times per cycle can result in missed cycles, and therefore the appearance of an incorrect lower converted frequency. For example, a 2 kHz sine wave being sampled at 1.5 kHz could be reconstructed as a 500 Hz sine wave. It is therefore important to define the rate at which new digital values can be sampled from the analog signal. The rate at which new values are sampled can be called the sampling rate or sampling frequency of the converter.


Bandpass sampling is a widely used approach to sample a continuous-time bandpass signal, and in so doing, downconvert the signal to a lower intermediate frequency (IF). When used in software-defined radios (SDR), bandpass sampling can allow the ADC to be placed closer to the antenna and can eliminate the need for certain analog circuitry components such as mixers, along with their associated nonlinear characteristics. To take advantage of these attributes, however, an appropriate designed IF and sampling frequency must be selected.


Selecting the “right” sampling frequency (Fs) and intermediate frequency (IF or fIF) can be hard to do, because to date, only a nonlinear relationship between the intermediate frequency (IF) and sampling frequency (or rate) has been presented when bandpass sampling is used. Cognitive radios (CR) and SDR's have the ability to vary the sampling rate used by their ADC, which changes the IF; but due to the nonlinearity of the aforementioned relationship, it can be challenging (and computationally intensive) to first choose a desired IF and then select the corresponding sampling rate without delving into an iterative, inefficient process.


In view of the above, it is an object of the present invention to provide a bandpass sampling method and device employing such design methods, which use a linear relationship between sampling frequency and IF. Another object of the present invention is to provide a bandpass sampling method and devices that can be used in place of the previously known nonlinear one, and which can find a sampling frequency without iteration in applications where the IF is important. Still another object of the present invention is to provide a bandpass sampling method and device wherein IF is chosen first as the primary design parameter and sampling is then chosen in a deterministic manner, without iteration. Yet another object of the present invention is to provide a bandpass sampling method and devices, which can reduce the complexity of the receiver. Another object of the present invention is to provide a bandpass sampling method and devices, which allow for more efficient allocation of bandwidth resources when multiple signals are being bandpass sampled. Yet another object of the present invention is to provide a bandpass sampling method and devices, which allows for more efficient carrier aggregation. Still another object of the present invention to provide a bandpass sampling method and devices, which can be easy to implement and deploy in a cost-efficient manner.


SUMMARY OF THE INVENTION

The methods for bandpass sampling of a continuous analog signal using a linear relationship between sampling frequency and intermediate frequency, and devices accomplishing such methods (such as analog to digital converters, ADC's, by way of non-limiting example) according to several embodiments can use the linear relationship







f
IF

=



(

-
1

)


n
+
1




(



-



n
2






F
S


+

f
C


)






where fIF is the intermediate frequency, Fs is the sampling frequency, fc is the carrier frequency for the signal and n is an integer such that







1

n





f
U

B




,




where B is the bandwidth and fU is the upper frequency (i.e., highest frequency component) of the signal and └ ┘ denotes the mathematical floor function. In cases where there is one signal at a single carrier frequency, the sampling frequency Fs can be chosen first and the IF can be calculated without iteration using the relationship above. Or, in instances where IF is deemed an important design parameter, the IF can be selected first, and sampling frequencies which will not cause aliasing can be solved by calculating Fs given IF using the equation above for all values of n so that







1

n





f
U

B




,




and then selecting a particular Fs. In some embodiments, the IF can be chosen to be Fs/4. In cases where there are multiple signals at multiple carrier frequencies, the above linear relationship can be used to select a particular Fs that simplifies processing when used along with a polyphase downconverter channelizer.





BRIEF DESCRIPTION OF THE DRAWINGS

The novel features of the present invention will be best understood from the accompanying drawings, taken in conjunction with the accompanying description, in which similarly-referenced characters refer to similarly-referenced parts, and in which:



FIG. 1A is a prior art graph in the X(f) frequency domain of intermediate frequency (IF) versus sampling frequency (Fs) for the original spectrum of a representative signal;



FIG. 1B is a prior art graph in the X(f) frequency domain of within the first Nyquist zone (±Fs/2) of intermediate frequency (IF) versus sampling frequency (Fs) for a normal spectral placement;



FIG. 1C is a prior art graph in the X(f) frequency domain of within the first Nyquist zone (±Fs/2) of intermediate frequency (IF) versus sampling frequency (Fs) for an inverted spectral placement;



FIG. 2 is a graphs k versus n, which is used to demonstrate the proof of the linear relationship between IF and Fs for the systems and methods of the present invention;



FIG. 3 is a graph of sampling frequency versus intermediate frequency for a signal having a carrier frequency of 100 MHz and a bandwidth of 20 MHz which compares the results of the linear methods of the present invention to the nonlinear methods of the prior art; and,



FIG. 4 is a diagram, which illustrates steps that can be taken to practice the methods of the present invention according to several embodiments;



FIG. 5A is a block diagram of a polyphase downconverter channelizer, which can incorporate some of the methods of the present invention according to several embodiments;



FIG. 5B is a graph in the X(f) frequency domain of intermediate frequency (IF) versus sampling frequency (Fs), which illustrates how an IF signal can fall within two separate channels; and,



FIG. 5C is a graph in the X(f) frequency domain of intermediate frequency (IF) versus sampling frequency (Fs), which illustrates how each IF signal can fall within its own channel, if Fs is chosen correctly according to the methods of the present invention according to several embodiments.





DETAILED DESCRIPTION OF THE EMBODIMENTS
1. Prior Art Bandpass Sampling Relationships

In the prior art, for a real, continuous-time bandpass signal x(t) with carrier frequency fc and bandwidth B (i.e. X(f) is only nonzero on the intervals (−fc−B/2, −fc+B/2) and (fc−B/2, fc+B/2)), the valid sampling frequencies that will not result in aliasing can be given by:











2






f
U


n



F
S




2






f
L



n
-
1






Eq
.





(
1
)








where fL=fc−B/2, fU=fc+B/2, and n is an integer such that









1

n





f
U

B







Eq
.





(
2
)








Instead of using the above equations, alternative equations in the prior art have been presented for finding valid sampling frequencies that included determining the IF. For example, the IF of the sampled signal in the prior art can be given by:










f
IF

=

{




rem
(


f
C

,

F
S


)








f
C



F
S

2









even







F
S

-

rem


(


f
C

,

F
S


)










f
C



F
S

2









odd









Eq
.





(
3
)








where rem is the remainder after division. According to Eq. (3), the IF is a nonlinear function of the sampling frequency. To make use of Eq. (3), the person of skill in the relevant art could: 1) Use Eq. (1) to select a valid sampling frequency FS; and, 2) Use Eq. (3) to determine the IF. Note that due to the nonlinearity of Eq. (3), FS needs to be selected first and then the IF can be calculated. The opposite order would be challenging to do, i.e. selecting the IF first and then calculating FS. Thus, Eq. (3) can be difficult to use in applications where the IF is an important design parameter.


2. Derivation of Linear Relationship Between if and FS

In order to derive a linear relationship between the IF and sampling frequency, it can be necessary to briefly review some of the basics of sampling theory. Using the aforementioned assumption that x(t) is a real, continuous-time bandpass signal with carrier frequency fc and bandwidth B, sampling x(t) at a sampling frequency FS can be equivalent to multiplying x(t) by a continuous-time impulse train with period T=1/FS and passing the resulting signal through an ideal lowpass filter with cutoff frequency FS/2. Since multiplication in the time domain is equivalent to convolution in the frequency domain, images of X(f) are centered at integer multiples of FS (since the Fourier transform of an impulse train with period T is also an impulse train with impulses at kFS where k=0, ±1, ±2, . . . ).


Assume also that FS is chosen such that it satisfies Eq. (1). This assumption applies for the rest of the derivation. Then, the IF signal is the image of the positive and negative spectrums of X(f) that fall within the first Nyquist zone (±Fs/2). This is shown in graphs 10, 12 and 14 of FIGS. 1A, 1B and 1C, respectively, which can represent the original spectrum, the normal spectral placement and the inverted spectral placement, respectively, of X(f). The graphical representations can be expressed quantitatively as the IF signal is the k0 image of the positive and negative spectrums of X(f) that satisfies the following:










k
0

=




min

k



0





k






s
.
t
.





kF
S




-

f
C




-


F
S

2







Eq
.





(
4
)








where custom-character0 denotes the set of nonnegative integers. The resulting IF can be given by the following:










f
IF

=

{






-

k
0




F
S


+

f
C





n





odd








k
0



F
S


-

f
C





n





even









Eq
.





(
5
)








For the IF and sampling frequency to be linearly related, k0 needs to be deterministic, which will now be proved, i.e., Eq. (4) needs to be solved. To solve Eq. (4), the constraint k in Eq. (4) can be rewritten as









k




f
C


F
S


-

1
2






Eq
.





(
6
)








For a fixed value of n, the smallest and largest values of FS in Eq. (6) can be given by Eq. (1). Substituting these two values of FS into Eq. (6) produces two equations for the lower bound of k. These two lower bounds are evaluated in the following subsections.


2.1. Smallest Lower Bound for k in Equation Eq. (6)


The smallest the right-hand side (RHS) of Eq. (6) can be is when the largest value of FS from Eq. (1) is substituted into Eq. (6). Eq. (6) can then be simplified to









k



n
-
2
+

B

2






f
C





2
-

B

f
C








Eq
.





(
7
)








In general, the bandwidth to carrier frequency ratio B/fc of a bandpass signal is greater than 0. Also note that since n is fixed, there is another implicit constraint imposed on B/fc from Eq. (2). More specifically, the possible values of B/fc can be









0
<

B

f
C




1

n
-

1
2







Eq
.





(
8
)








Looking at Eq. (7), the RHS is smallest when B/fc=c, where c is a very small positive number close to 0. Substituting this into Eq. (7) results in









k



n
-
2
+

ɛ
2



2
-
ɛ






Eq
.





(
9
)








The RHS of Eq. (7) is largest when B/fc is largest. Equation Eq. (7) then becomes









k




n
2

-


5
2


n

+

3
2




2





n

-
2






Eq
.





(
10
)








Referring now to FIG. 2, the top graph 22 of FIG. 2 illustrates the smallest integer value of k that satisfies Eq. (10). As seen in the figure, the smallest integer values of k that satisfy Eq. (9) and Eq. (10) are both equivalent to └n/2┘ This can also be verified by substituting k=n/2 for n even and k=(n−1)/2 for n odd into Eq. (9) and Eq. (10) and seeing that both inequalities can be still satisfied, while k=(n−2)/2 and k=(n−3)/2 do not satisfy them. Thus, the smallest integer value of k that satisfies Eq. (7) is the same for all values of B/fc in Eq. (8), and Eq. (7) simplifies to









k
>



n
2







Eq
.





(
11
)








2.2. Largest Lower Bound for k in Equation Eq. (6)


For the present invention according to several embodiments, the proof in this section is analogous to the previous section. The largest the RHS of Eq. (6) is when the smallest value of FS from Eq. (1) is substituted into Eq. (6). Eq. (6) can then become









k



n
-
1
-

B

2






f
C





2
+

B

f
C








Eq
.





(
12
)








The RHS of Eq. (12) is smallest when B/fc is largest from Eq. (8). Eq. (12) can then become









k



n
-

3
2


2





Eq
.





(
13
)








The RHS of Eq. (12) is largest when B/fc=ε (where ε is defined above). This makes Eq. (12) simplify to









k



n
-
1
-

ɛ
2



2
+
ɛ






Eq
.





(
14
)








Referring again to FIG. 2, the plots of equations Eq. (13) and Eq. (14) can be seen at the bottom graph 24 of FIG. 2 for the smallest integer value of k and ε=10−12. As seen in the figure and can be verified just like the previous section, the smallest integer values of k that satisfy Eq. (13) and Eq. (14) are both equivalent to └n/2┘. Thus, Eq. (12) also simplifies to Eq. (11), i.e.






k





n
2



.





2.3. Linear Relationship


Since the two lower bounds of k are always the same for a fixed value of n, the solution to Eq. (4) is










k
0

=



n
2







Eq
.





(
15
)








where n is from Eq. (2). The IF in Eq. (5) can then be rewritten as










f
IF

=



(

-
1

)


n
+
1




(



-



n
2






F
S


+

f
C


)






Eq
.





(
16
)








Equation (16) illustrates that the IF and FS can be expressed as a linear relationship, conditioned on n.


Referring now to FIG. 3, FIG. 3 illustrates a visual comparison of the nonlinear Eq. (3) of the prior art, and the linear Eq. (16), which is used in the systems and methods of the present invention. More specifically, FIG. 3 graphs Eq. (3) and Eq. (16) for the case when fc=100 MHz and B=20 MHz. As can be seen in FIG. 3, the linear relationship in Eq. (16) agrees well with non-linear relationship in the prior art, as described in Eq. (3). Due to the linearity of Eq. (16), it can be used in any case that the nonlinear Eq. (3) is used, plus additional cases. If the sampling frequency is the most important design parameter and the IF does not matter, then one can: Use Eq. (1) to select a valid FS, and in so doing, determine n at the same time; and, Use Eq. (16) to determine the IF.


However, if the IF is an important design parameter, and referring now to method 40 in FIG. 4, if there is flexibility in selecting the sampling frequency, then one can: 1) Select the fIF to be used in Eq. (16), as shown by step 42, and then solve for FS for all values of n in Eq. (2), as shown by box 44; 2) Check which values of FS can be valid using Eq. (1), box 46; and, 3) Select a particular FS that meets other design parameters, such as ADC sampling capabilities, as illustrated by box 48 in FIG. 4. Some applications where the IF can be an important design parameter and Eq. (16) would be useful can be discussed below.


3. Applications

3.1 Digital Down Conversion of a Signal


For the methods of the present invention, an application where the IF can be important (the primary design parameter) when using bandpass sampling can be digital down conversion (DDC). This can be because the complexity of the DDC design can be reduced if both the IF and sampling frequency can be chosen carefully. In the prior art, a digital receiver has been presented where the IF was chosen to be FS/4. By doing so, the discrete-time cosine and sine functions that can be used to shift the real IF signal to complex baseband become [1, 0, −1, 0, . . . ] and [0, −1, 0, 1, . . . ], respectively, i.e., integer numbers. Since Eq. (16) is linear, it can easily be used to find sampling frequencies that place the IF at FS/4 by solving the following equations and checking the validity with Eq. (1):












F
S

4

=



(

-
1

)


n
+
1




(



-



n
2






F
S


+

f
C


)



,






where





n

=
1

,
2
,





,





f
U

B



.





Eq
.





(
17
)








Also, Eq. (17) can be used in conjunction with other linear constraints to determine FS, such as requiring FS to be an integer multiple of the symbol rate (the number of symbol changes to the transmission medium per second using a digitally modulated signal) in digital communications.


3.2. Multiple Signal Reception


The IF can also be of primary importance in applications where bandpass sampling is used to receive multiple signals at different carrier frequencies simultaneously. One such application is carrier aggregation. Carrier aggregation has been proposed in the prior art for LTE-Advanced to increase the data rates of cellular networks, and it allows for the use of up to 5 carriers (each with bandwidths up to 20 MHz) that can be spaced apart in frequency. Bandpass sampling can be used to sample the signals on all (up to 5) carriers simultaneously.


The systems and methods of the present invention according to several embodiments can accomplish carrier aggregation. To do this, it can be necessary to make sure the resulting IF signals do not overlap. For example, assume two carriers need to be received at carrier frequencies fC1 and fC2, with bandwidths B1 and B2, respectively. Further assume that two potential ranges of sampling frequencies, denoted as FS1 and FS2, can be found using Eq. (1) and Eq. (2), where the frequencies in FS1 can sample the signal at fC1 with no aliasing and the frequencies in FS2 can sample the signal at fC2 with no aliasing. Denote the corresponding values of n from Eq. (2) as n1 and n2. If the intersection of FS1 and FS2 is not empty, there may be a subset of the intersection that will not result in the 2 IF signals overlapping. Equation Eq. (16) can be used to find the subset of sampling frequencies corresponding to this subset (if it exists) by solving either Eq. (18a) or Eq. (18b):











f

IF





1


+


B
1

2





f

IF





2


-


B
2

2






(

18

a

)








f

IF





2


+


B
2

2





f

IF





1


-


B
1

2






(

18

b

)







where Eq. (16) can be substituted for fIF1 and fIF2 with n replaced with n1 and n2, respectively, and FS is constrained to be in the intersection of FS1 and FS2. If a solution is found to (18a) (or (18b)), then the IF signal corresponding to fC1 is located at a lower (higher) frequency than the IF signal corresponding to fC2. If neither (18a) nor (18b) can be satisfied, then a different range of sampling frequencies needs to be checked. Note that the linearity of Eq. (16) of the methods of the present invention makes (18a) or (18b) straightforward to solve, whereas the nonlinearity of Eq. (3) of the prior art makes these difficult to solve.


The above process described the process for simultaneously sampling two signals with arbitrary carrier frequencies and bandwidths. It should be appreciated, however, that the process could be extended to the carrier aggregation for the aforementioned 5 carriers, or even any number of carriers (each carrier with an arbitrary bandwidth), by comparing each carrier's IF signal to every other carrier's IF signal to ensure that there is no IF overlap (and thus, no aliasing). In this manner, carrier aggregation for N sub-bandwidths can be accomplished more efficiently and with less computational resources than the devices and methods of the prior art.


The above procedure can be applicable to situations where the carrier frequencies of the signals are either static or dynamic. Two commercial technologies where it is applicable are global navigation satellite systems (GNSS) and Long Term Evolution-Advanced (LTE-Advanced). With respect to GNSS, currently there are four GNSS signals: GPS, GLONASS, GALILEO, and COMPASS. GPS and GLONASS are operational, while GALILEO and COMPASS are still under development. Each GNSS can use multiple carrier frequencies, and it can be desirable for a user to receive signals simultaneously from multiple GNSS, thus improving positioning performance since more satellites are available, when compared to a single GNSS. Therefore, the algorithm described for multiple signal reception can be used to sample signals from some or all four types of GNSS simultaneously, despite each type having a different carrier frequency. Stated differently, since positioning signals i.e., global navigation satellite system (GNSS) signals and signals of opportunity (SoOP) which are signals that can be used for position estimation even though these signals were not originally intended for positioning applications. Examples of SoOP can include cellular signals like LTE, WiFi, and an Advanced Television Systems Committee (ATSC)-standard digital television signal. occupy different frequencies, the systems and methods of the present invention can allow a receiver to calculate which sampling frequency to use to sample all of these signals simultaneously with one radio frequency (RF) front end. From this one sampled signal, the receiver can utilize multiple positioning signals to more accurately determine its position.


3.3. Polyphase Downconverters


Just like Eq. (16) can be used to simplify DDC when one signal is received, Eq. (16) can also be used to simplify the DDC of multiple signals. The standard M-path polyphase downconverter channelizer as known in the prior art can be commonly used to efficiently downconvert multiple signals when their carrier frequencies are equally spaced apart and they have equal bandwidths. If the carrier frequencies of the signals can be equal to integer multiples of FS/M (i.e. the output sampling rate of the M channels of the polyphase channelizer), then all of the signals will be shifted to baseband by the channelizer. If the signals have arbitrary carrier frequencies and bandwidths, a receiver design can be used, which consists of a modified version of the standard M-path polyphase downconverter channelizer followed by additional processing. When bandpass sampling is used in a receiver followed by a polyphase downconverter, the carrier frequencies of the input to the downconverter can be the IFs of the sampled signals.



FIG. 5A depicts an example of such a downconverter 50. Graph 52 in FIG. 5B depicts a scenario where the IFs are equal to integer multiples of FS/M and one of the IF signals falls within two input channels of the downconverter 50 that is shown in FIG. 5A. This scenario can require additional processing after the polyphase downconverter discussed to shift the signals to baseband. However (depending on the particular carrier frequencies and bandwidths of the signals before being sampled by the ADC), it may be possible to select a sampling frequency Fs using the methods of the present invention such that the IF's of the sampled signals can be equal to integer multiples of FS/M, which can result in the scenario depicted in graph 54 as depicted in FIG. 5C. Then, all of the signals will be at baseband at the output of the downconverter 50 and the additional processing that is required in the prior art can be avoided. The linearity of Eq. (16) of the methods of the present invention can be used to accomplish this. More specifically, and as discussed in the previous section, Eq. (16) and Eq. (18) can be used to help find a range of sampling frequencies that do not result in the overlap of the IF signals. Then, similar to Eq. (17), Eq. (16) can be used to select a particular sampling frequency within the range that places some or all of the IF signals either at integer multiples of FS/M or such that each IF signal falls entirely within an input channel of the downconverter as shown in FIG. 5C.


3.4. Cognitive Radio


With the ongoing research in SDR, cognitive radio, and dynamic spectrum access, a radio may have to receive a signal whose carrier frequency changes over time using bandpass sampling (assuming the radio has a tunable bandpass filter at its RF front end). This may be due to either changes in the environment or due to action taken by the operator. But, no matter what the reasons for a change in carrier frequency, the systems and methods of the present invention can be used to quickly and dynamically determine a sampling frequency that would place the IF at a desired location, such as FS/4, to simplify the DDC process.


The use of the terms “a” and “an” and “the” and similar references in the context of describing the invention (especially in the context of the following claims) is to be construed to cover both the singular and the plural, unless otherwise indicated herein or clearly contradicted by context. The terms “comprising,” “having,” “including,” and “containing” can be to be construed as open-ended terms (i.e., meaning “including, but not limited to,”) unless otherwise noted. Recitation of ranges of values herein can be merely intended to serve as a shorthand method of referring individually to each separate value falling within the range, unless otherwise indicated herein, and each separate value is incorporated into the specification as if it were individually recited herein. All methods described herein can be performed in any suitable order unless otherwise indicated herein or otherwise clearly contradicted by context. The use of any and all examples, or exemplary language (e.g., “such as”) provided herein, is intended merely to better illuminate the invention and does not pose a limitation on the scope of the invention unless otherwise claimed. No language in the specification should be construed as indicating any non-claimed element as essential to the practice of the invention.


Preferred embodiments of this invention can be described herein, including the best mode known to the inventors for carrying out the invention. Variations of those preferred embodiments may become apparent to those of ordinary skill in the art upon reading the foregoing description. The inventors expect skilled artisans to employ such variations as appropriate, and the inventors intend for the invention to be practiced otherwise than as specifically described herein. Accordingly, this invention includes all modifications and equivalents of the subject matter recited in the claims appended hereto as permitted by applicable law. Moreover, any combination of the above-described elements in all possible variations thereof is encompassed by the invention unless otherwise indicated herein or otherwise clearly contradicted by context.

Claims
  • 1. A method for selecting a sampling frequency Fs when using bandpass sampling to receive a continuous analog signal, said method comprising the steps of: A) choosing an intermediate frequency fIF B) solving for
  • 2.-5. (canceled)
  • 6. The method of claim 1, wherein said fIF from said step A) is Fs/4.
  • 7. The method of claim 1, wherein said signal is a digital communication signal and said Fs from said step D) is an integer multiple of the signal's symbol rate.
  • 8.-19. (canceled)
  • 20. An improved analog to digital converter (ADC): said ADC using a predetermined intermediate frequency fIF;said ADC accomplishing bandpass sampling of a continuous signal at a sampling frequency Fs without iteration of said fIF, based on the relationship
FEDERALLY-SPONSORED RESEARCH AND DEVELOPMENT

The United States Government has ownership rights in this invention. Licensing inquiries may be directed to Office of Research and Technical Applications, Space and Naval Warfare Systems Center, Pacific, Code 72120, San Diego, Calif., 92152; telephone (619) 553-5118; email: ssc_pac_t2@navy.mil, referencing NC 102470.