Abstract
A wave construction method is described that can be used to generate a new and different wave front, which is not parallel to the natural expanding wave front from the emitted signal generated from the source array. Restated, a wave construction method is described that can be used to generate a new and different wave front which is not perpendicular to the direction from the source array to the center of the wave front. This method can also be used to shape a wave in the far field or near field, by changing the locations of the computed points, and generating a linear or non-linear shaped wave (front). This disclosure allows for the controlling, rotating, or shaping radio frequency or acoustic waves in free space or in a fluid.
Claims
1. A method of wave construction for controlling, rotating, or shaping radio frequency or acoustic waves comprising: determining a multiplicity M of points in the far field of a radio frequency (RF) antenna or acoustic transducer array comprising M antennas or transducers, determining electric or acoustic field voltages and phases at said M points, wherein said determined voltages and phases emulate the same equipotential voltage and phase characteristics as a natural expanding wave from the array, but rotated or wrapped onto a different virtual surface, computing a set of M complex array weights h for the M antennas or transducers based on the determined points and the corresponding determined voltages and phases, wherein the weights h are computed according to
h=k.sub.oR.sub.xx.sup.1V where h is a vector of the weights h for each antenna, k.sub.o is a constant, R.sub.xx is an MM matrix comprising elements 1/r.sub.ij exp(jkr.sub.ij), where r.sub.ij is the distance between the ith antenna or transducer and the jth point, and V is a vector of the determined voltages and phases, generating a signal for each of the M antennas or transducers, wherein each said signal has a same message signal content of a transmission from a single source, and transmitting, from each antenna or transducer, a complex signal formed from the multiplication of a respective generated signal and a respective computed complex array weight, wherein the transmitted complex signals combine at the selected far field points to produce a wavefront angle that is not perpendicular to the direction or location of the array antennas or transducers.
2. The method of claim 1, wherein the wave is shaped in the far field by changing locations of the multiplicity of points and generating a linear or non-linear shaped wave.
3. The method of claim 1, wherein the wavefront angle is rotated with respect to the natural expanding wave by a predetermined angle.
4. The method of claim 1, wherein the rotated wavefront is received by a passive direction finding system or radar with an estimated angle that is not perpendicular to the direction or location of the array antennas or transducers.
5. The method of claim 1, wherein the multiplication of the complex signals by the weights produces constructive and destructive interference at the multiplicity of determined points in the far field.
6. The method of claim 1, wherein the multiplicity of points are determined along a virtual line that is rotated from the natural expanding wave by a predetermined angular value.
7. The method of claim 1, wherein the multiplicity of points are determined along a virtual curve that is rotated from the natural expanding wave by a predetermined angular value.
8. The method of claim 1, wherein the multiplicity of points are determined along a virtual line that is rotated from the natural expanding wave by a predetermined angular value, and wherein the transmitted complex signals combine constructively and destructively at said points.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) For a more complete understanding of the present invention, reference is now made to the following descriptions taken in conjunction with the accompanying drawings, in which:
(2) FIG. 1A-FIG. 1A illustrates the RF beam generated by the system of FIG. 2A, that radiates (or travels) along the direction of propagation of the wave. For this illustration an (RF) Dish or directive antenna is shown. Another embodiment of this technique would be an acoustic wave generated by the system of FIG. 2B. For acoustics, the wave emitter/radiator would be an acoustic transducer.
(3) FIG. 1B-FIG. 1B illustrates the RF beam that radiates (or travels) along the direction of propagation of the wave, generated from a multiplicity of RF antennas. For this illustration an (RF) Array of antennas is shown. Another embodiment of this technique would be an acoustic wave. For acoustics, the wave generation array would be an acoustic transducer array.
(4) FIG. 1C-FIG. 1C illustrates multiple RF beams radiates from a single antenna array. Each beam has its own unique direction of propagation of its wave. The waves are enabled by the multiplication of a combined set of beamformed weights h=h.sub.1+h.sub.2+ . . . +h.sub.N, where each h.sub.i (i=1, 2, . . . , N) is the individual complex multiplicative weight for each i.sup.th beam.
(5) FIG. 1D-FIG. 1D illustrates the outward motion of originally curved wave into planar waves in the far field. As a wave front gets further away from its source of generation, the expanding sphere of the wave front becomes planar in a finite local region, L. The local region is defined where L<<r. Furthermore, the planar wave front is perpendicular to the direction of propagation of the wave.
(6) FIG. 1E-FIG. 1E illustrates one embodiment resultant, that of a rotated wave. A rotated wave is generated from the Wave Mechanics system and mathematics. If received from a passive DF system or radar receiving antenna or array, the receiving system will estimate or compute that the wave originated from a direction that is not in the direction (of wave propagation) of the actual source array, shown at the bottom of the diagram.
(7) FIG. 2A-FIG. 2A shows a block diagram of one embodiment of the Wave Mechanics transmission system. This particular embodiment of the invention includes, but is not limited to, baseband signal conversion to RF for a multiplicity (array) of antennas. The source signal generator, 104, produces a digital signal that is processed by the DSP processor block, 103, which also multiplies the signal, s(t), by the weight vector h, and forwards each antenna signal to the multi-channel Digital to RF converter, 102. As one embodiment of this invention, the processed signal is first converted to analog, within a spatial channel, via a Digital-to-Analog Converter (DAC), and then RF mixed or upconverted to RF. This conversion process is usually performed coherently, however, there are other embodiments that would not require coherent upconversion.
(8) FIG. 2B-FIG. 2B shows a block diagram of another embodiment of the Wave Mechanics transmission system. This particular embodiment of the invention includes, but is not limited to, baseband to acoustic signal conversion for a multiplicity (array) of acoustic transducers. The source signal generator, 114, produces a digital signal that is processed by the DSP processor block, 113. which also multiplies the signal, s(t), by the weight vector, h, and forwards each antenna signal to the Digital to acoustic signal converter, 112. As one embodiment of this invention, the processed signal is first converted to analog, within each spatial channel, via a Digital-to-Analog Converter (DAC), and then converted to an acoustic signal. This conversion process is usually performed coherently, however, there are other embodiments that would not require coherent conversion.
(9) Comment: FIGS. 3A through 6A refer to RF techniques and mathematics, but are all equivalent to acoustic mathematics and applications.
(10) FIG. 3A-FIG. 3A illustrates a multiplicity of RF antennas, 101a through 101M, as an array, and the far field point, 105a, at which the Wave Mechanics process is to be focused to. It should be noted that the actual distance from the multiplicity of antennas, 101a through 101M, to the far field point, 105a, would be larger than the physical size of the array of antennas.
(11) FIG. 3B-FIG. 3B illustrates the ranges from each antenna in the array to the far field point, 105a. For example, r.sub.1 is the range (distance) from antenna #1 to the far field point, 105a. It should be noted that the actual distance from the multiplicity of antennas, 101a through 101M, to the far field point, 105a, would be larger than the physical size of the array of antennas.
(12) FIG. 3C-FIG. 3C illustrates a multiplicity of far field points, 105a through 105M, and the corresponding multiplicity of antennas, 101a through 101M, in the array. Again it should be noted that the actual distance from the multiplicity of antennas, 101a through 101M, to the multiplicity of far field points, 105a through 105M, would be larger than the physical size of the array of antennas.
(13) FIG. 3D-FIG. 3D illustrates the nomenclature for the range r.sub.ij, where the i.sup.th component specifies the source antenna number (or transducer reference number, for acoustics) and the j.sup.th component specifies the far field point reference number. For example, range r.sub.23 is the range between antenna #2 and far field point #3, 105c.
(14) FIG. 3E-FIG. 3E illustrates the set of ranges r.sub.ij from each antenna, to far field point #1, 105a. The total field at far field point #1, 105a, is the summation of the fields generated from the multiplicity of antennas, 101a through 101M, with respective ranges r.sub.11, r.sub.21, . . . , r.sub.M1. This collection of ranges to a single point can be denoted as a butterfly. It is similar in concept to the butterfly used in generating Fast Fourier Transforms (FFT) in Digital Signal Processing. This can be denoted the butterfly for far field point #1, 105a.
(15) FIG. 3F-FIG. 3F illustrates the set of ranges r.sub.ij from each antenna, to far field point #2, 105b. The total field at far field point #2, 105b, is the summation of all the fields generated from the multiplicity of antennas, 101a through 101M. The total field at far field point #2, 105b, is the summation of all the fields generated from the multiplicity of antennas, 101a through 101M, with respective ranges r.sub.12, r.sub.22, . . . r.sub.M2. This is therefore the butterfly for far field point #2, 105b.
(16) FIG. 3G-FIG. 3G illustrates the collection of butterflies, from all antennas, to the respective far field points, 105a through 105M.
(17) FIG. 4A-FIG. 4A illustrates the growing set of spherical wave fronts, radiated from the source antenna. As the wave front travels outward from the source antenna, the curvature of each (spherical) wave front is reduced; across a given fixed distance, L. At very far distances, small sections of the wave front, L, are effectively planar (flat). Additionally, the planar section is perpendicular to the line of direction or direction of wave propagation.
(18) FIG. 4B-FIG. 4B illustrates one embodiment of the Wave Mechanics technique, that of denoting a multiplicity of far field points in space, centered at a given (arbitrary) point in space, to generate a rotated wave. The multiplicity of far field points are spaced along a line that has been rotated from the original planar wavefront, by an angular value of . This rotation angle can be at any determined angular value. Therefore, the specific location of the far field points, 105a through 105M, have been computed from the location of the Antenna Array points, 101a through 101M, and a specified point along the rotated wave line.
(19) FIG. 4C-FIG. 4C illustrates the wave rotation embodiment of the Wave Mechanics technique, that of denoting a multiplicity of far field points in space, centered at a given point in space, to generate a rotated wave. A receiver, such as a passive Direction Finding System, or radar system, at this location would measure the energy, at a given frequency, and estimate that the wave front bearing angle to a virtual or non-existent antenna or array, as shown in the diagram. This mechanism would thus fool or spoof such a passive or active receiving array.
(20) FIG. 4D-FIG. 4D illustrates another embodiment of the Wave Mechanics technique, that of shaping an arbitrary wave in the far field. The shaping of this wave could be multiple configurations, specified by the location of far field points 105a through 105M.
(21) FIG. 4E-FIG. 4E illustrates an effect of the wave rotation mechanism. While a finite number of far field points, M, 105a through 105M, have been specified, the effect of the wave rotation mechanism is to generate similar phase and amplitude to the M points, but at also an an infinite collection of points along the rotated wave. Note, this rotated wave will have a finite length, L, determined by many factors, such as the number of elements in the transmit array, the operational frequency, and other physical factors.
(22) FIG. 5A-FIG. 5A illustrates another potential effect of the wave rotation mechanism. The rotated wave will have a finite width and possibly length, determined by many factors, such as the number of elements in the transmit array, the operational frequency, and other physical factors. Therefore, the rotated wave would appear to travel in a rotated corridor, different in propagation direction than the original wave propagation angle. This corridor could have finite length or infinite length.
(23) FIG. 5B-FIG. 5B illustrates another potential effect of the wave rotation mechanism, for the wave traveling in a corridor This corridor could have finite length, L.
(24) FIG. 6A-FIG. 6A shows a comparison of a conventional beamformer, of M antennas, to the Wave Mechanics principle. The traditional (conventional) Beamformer computes a set of complex digital weights, h, to generate RF constructive interference at a single point in space. Due to geometry, where r.sub.ij>>>, the single far field point in space can be any far field point along a line from the transmit antenna out to infinity. In contrast, the Wave Mechanics technology computes h to generate constructive interference, simultaneously, at a multiplicity of points in space (e.g. the Far Field), and generates the wave rotation or shaping within a finite width corridor in the Far Field.
DETAILED DESCRIPTION OF THE INVENTION
(25) FIG. 2A shows the block diagram of one embodiment of the Wave Mechanics transmission system; for an antenna/array (RF) system. This is comprised of a multiplicity of M antenna elements, 101a through 101M, each fed by a coherent (in phase) RF converted signal. This embodiment of the present invention includes, but is not limited to, baseband signal conversion to RF, for a multiplicity (array) of antennas. The source signal generator, 104, produces a digital signal that is processed by the DSP processing block, 103, which also multiplies the signal, s(t), by the weight vector, h, and forwards each antenna signal to the Digital to RF converter block, 102.
(26) As one embodiment of this invention, the processed signal is first converted to analog, within a spatial channel, via a Digital-to-Analog Converter (DAC), and then RF mixed or upconverted to RF. This conversion process is usually performed coherently, however, there are other embodiments that would not require coherent upconversion. The Digital to RF Conversion can be accomplished via an RF Exciter, a Digital to Analog Converter (DAC), an RF Transceiver, or a split RF synthesizer. The Signal Source could be any type of baseband source that is then digitized; prior to the Digital to RF Conversion stage. The Digital Signal Processing (DSP) processor both computes the optimal weight vector, h, as well as performs the real time multiplication at the baseband ample rate, of:
output=h(t).Math.s(t)(1
(27) The process inputs the digitized signal, s(t), copies the signal M times, and multiplies each sample; for the same time instant, by h.sub.i, where i=1, 2, . . . , M is the antenna reference number.
(28) Note, this is similar to the conventional beamformer process, however only in the generation of the output (s(t).Math.h). For the Wave Mechanics technology, h is generated completely different from the method the conventional beamformer uses to compute h.
(29) FIG. 2B shows the general embodiment of the same system, however, for acoustics applications and frequencies.
(30) For a tone that is transmitted (radiated) from an Antenna #1, at a distance, r, in the far field from the antenna, the field voltage for any far field distance (r) frequency (f) and time (t) can be represented as a traveling wave:
(31)
(32) In volts per meter. Where:
(33) r=displacement (distance) from antenna #1 to distance of a given point
(34) f=frequency of the wave
(35) t=time, and
(36) Where >>: the far field definition. The wavelength can be written in terms of the speed of light, c, and frequency, f, as:
(37)
for acoustics, c would be the speed of sound in the fluid or air.
(38) FIG. 3A illustrates a multiplicity of RF antennas, 101a through 101M, as an array, and the far field point, 105a, at which the Wave Mechanics process is to be applied to. As shown in FIG. 3B, the lengths can be denoted as r.sub.1, r.sub.2, . . . , r.sub.M and resulting voltage at the far field point as:
(39)
(40) This is the voltage sum, from the M antennas, each with a different distance and wave reference time, t.sub.i. Assume that the transmitted signals from each antenna are now coherent (e.g. synchronized in time), then
t=t.sub.1=t.sub.2= . . . t.sub.M(5
(41) Relationship (4) can therefore be expressed as:
(42)
(43) By weighing each signal, transmitted from each antenna, with vector h=[h.sub.1, h.sub.2, . . . , h.sub.M], the weighted sum for (s) can be expressed as:
(44)
(45) This can be expressed in vector form as:
(46)
(47) Or in compact form:
V.sub.W(f,t)=h.sup.T.Math.V(f,t)(9
(48) The scalar V.sub.W(f,t) is a maximum when h=conjugate[V(f,t)]. This is an example of simple (conventional) RF beamforming. Without loss of generality, this derivation and expression also applies to acoustic beamforming.
(49) Wave Mechanics operates by generating a collection of points in the far field. This is shown by FIG. 3C. Note that the drawing shows the Fair Field points close to the array antennas. For sake of argument, and not requiring a very large drawing, it should be noted that the actual distance from the multiplicity of antennas, 101a through 101M, to the multiplicity of far field points, 105a through 105M, would be much larger than the physical size of the array of antennas.
(50) For consistency, the two dimensional displacement vector will be denoted by using the convention of:
r.sub.ij=r.sub.from,to
(51) That is, the first component in the subscript is the value of the antenna reference number in the array, and the second component in the subscript is the referenced far field point. Thus the from value is the antenna (number), and the to value is the reference point in the far field. This is shown more clearly in FIG. 3D.
(52) The weighted voltage at Far Field Point #1, shown by FIG. 3E, is expressed as:
(53)
(54) FIG. 3E illustrates the set of ranges from each antenna, to far field point #1, 105a. The total field at far field point #1, 105a, is the summation of the fields generated from the multiplicity of antennas, 101a through 101m, with respective ranges r.sub.11, r.sub.21, . . . r.sub.M1. This collection of ranges to a single point can be denoted as a butterfly. It is similar to the butterfly used in generating Fast Fourier Transforms (FFT) in Digital Signal Processing.
(55) A finite bandwidth signal, s(t), can be coherently injected into each antenna. Therefore (10) can be expressed, with s(t), as:
(56)
(57) The same weights, h=[h.sub.1, h.sub.2, . . . , h.sub.M], can be used to force a voltage at the second (arbitrarily chosen) point, shown by FIG. 3F, in the far field with:
(58)
(59) FIG. 3F illustrates the set of ranges r.sub.ij from all antennas, to far field point #2, 105b. The total field at far field point #2, 105b, is the summation of the fields generated from the multiplicity of antennas, 101a through 101m. The total field at far field point #2, 105b, is the summation of the fields generated from the multiplicity of antennas, 101a through 101m, with respective ranges r.sub.12, r.sub.22, . . . r.sub.M2. This is another butterfly, with ranges from all antennas, yet to different far field point, 105b.
(60) This can be continued, to the M.sup.th far field point, as:
(61) 0
(62) The relationships in (11), (12), through (13) can be expressed in matrix form as:
(63)
(64) This embodiment can be represented by the drawing in FIG. 3G. FIG. 3G illustrates the collection of butterflies, from all antennas, to the respective far field points, 105a through 105m.
(65) It should be noted that the message signal, s(t), can be literally any (modulated) signal with finite bandwidth. This can include a Digital Radio Frequency Memory (DRFM) signal.
(66) Notice that since all signals are coherently RF converted, with synchronized initial phases, then the time dependence is the same for all components. This time dependence can be removed from all matrix values, to a constant multiplied by the matrix, expressed as:
(67)
(68) The expression in (15) can be rewritten in compact form as:
s(t)e.sup.jtR.sub.xxh=V(16
or
(69)
(70) Solving for h:
(71)
(72) Note that:
(73)
(74) can be treated as a constant. Therefore:
h=k.sub.oR.sub.xx.sup.1V(20
(75) Which is a more compact form for the expression.
(76) It should be noted, that in the very far field, where r.sub.ij>>>, then
r.sub.11r.sub.12r.sub.21 . . . r.sub.MMr.sub.o(21
(77) Where r.sub.o is the approximate distance from the center of the array, to the center of the collection of Far Field points. The expression in (15) can be rewritten as:
(78)
(79) Or in compact form:
(80)
(81) Again, the value in the parenthesis can be treated as a constant:
(82)
(83) Therefore:
h=k.sub.oR.sub.xx.sup.1V(25
(84) This is a more compact form for the expression. This value of h, for Wave Mechanics, operates on a multiplicity of far field points. This is in comparison to h in (9), the beamformer case, which only operates on a single far field point.
(85) Applications:
(86) This relationship (25) can be used in a number of ways.
(87) Plane Wave Rotation:
(88) Wave mechanics can be used to generate a line of points, with equal amplitude and phase, yet not from an angle perpendicular to the pointing vector direction between the center of the array to the center of the collection of points. For this case, all the values in V from (25), would be the same; both in phase and amplitude. However, the Far Field points are not in a line parallel with the original plane wave wave-front. This is shown more clearly in FIGS. 4A and 4B. In FIG. 4A, a single antenna, or set of beamformed antennas (in an array) generates an expanding wavefront, for a given signal generated at a specific point in time. In the far field, due to the expanding size of the spherical wave, which denotes the signal information at a very small fraction of time, the wave front is nearly planar within a small region of the local space.
(89) In FIG. 4B, the exact locations are computed for a set of M points (show by red dots) in the Far Field, arranged along in a line. This line, of points, has been rotated from the original/natural plane wave by an angle . Using this angle , the line for the rotation wave has been computed, and points have been arranged along the line. Location and separation distance of these points are based on desired characteristics of the rotated wave. The new direction of propagation, of this rotated wave, is shown by the red array. The newly generated wave front has wave characteristics identical to that of a natural wave, but launched by a virtual (non-existent) antenna or array that is far to the left of the actual array, shown as a collection of M antennas in the drawing. This is shown more clearly in FIG. 4C. Applications include, but are not limited, to spoofing or fooling (RF) Surface to Air Missile systems, incoming missiles, and (Acoustic) decoys to fool torpedo's or submarine acoustic detection and tracking systems.
(90) Wave Shaping with Similar Phase and Amplitude:
(91) Wave mechanics can be used to shape the wave, at some location in the Far Field. FIG. 4D shows a collection of wave points, shaped in a half circle. The apparent line of wave propagation appears as coming from multiple angles, for a potential receiver at the center of the circle. For this case, again, all components of V are the same value; in both amplitude and phase.
(92) Wave Shaping with Phase and/or Amplitude
(93) An application could be constructed with a desired shaped or rotated wave, but with dis-similar phase and or amplitude values in V.
(94) Comments:
(95) For a given carrier frequency, the matrix Rxx is based solely on geometry; that is displacements between each antenna and each Far Field location point.
(96) The Wave Mechanics resultant effectively generates a greater multiplicity of N beamformed points, along the rotation of shaped where, where N>>>M. This is more clearly shown in FIG. 4E.
(97) The Rotated Plane Wave line, shown in FIG. 4E consists of the computed and described M Far Field point (shown in solid red) as well as the near-infinite non-computed points (shown in non-solid red) along the rotated wave line. These non-computed points, as well as all computed (solid red) lines will have similar beamformed power levels. This is analogous to a beamformed wall rather than a single beamformed point. Therefore, the rotated wave will have similar power level across the full wave corridor, similar to a section of the original (non-rotated) wave.
(98) In certain embodiments, the rotated wave(s) could travel in a finite width and finite length corridor. For certain embodiments, unlike conventional RF Beamforming, the rotated waves might not propagate to infinity. Thus, the wave rotation effect could be limited in range, from the center of the collection of Far Field computed points. This is shown more clearly in FIGS. 5A and 5B. The Corridor width, W, and length, L, are determined by the number of antennas in the array, M, the frequency, and the separation distance between the Far Field (computed) points.
(99) Wave Mechanics is very different from conventional RF (or acoustic) beamforming. In a conventional Beamformer, h only operates on the first (or a single row) of Rxx. That is, beamforming only computes h based on a single far field point; which is expressed as a single row of Rxx. In the Wave Mechanics technique, h includes the contribution all rows of Rxx, and thus simultaneously operates on all rows of Rxx.
(100) As shown in FIG. 6A, the traditional (conventional) Beamformer computes h to generate RF constructive interference at single point in space. Due to geometry, where R.sub.ij>>>, this generates literally an infinite corridor or beam in the Far Field along the natural wave of propagation. In contrast, the Wave Mechanics technology computes h to generates constructive interference, simultaneously, at a multiplicity of points in space (e.g. the Far Field), and generates the wave rotation or shaping within a finite width corridor in the Far Field.