US20160280344A1 - Underwater Acoustic Array, Communication and Location System - Google Patents
Underwater Acoustic Array, Communication and Location System Download PDFInfo
- Publication number
- US20160280344A1 US20160280344A1 US14/483,631 US201414483631A US2016280344A1 US 20160280344 A1 US20160280344 A1 US 20160280344A1 US 201414483631 A US201414483631 A US 201414483631A US 2016280344 A1 US2016280344 A1 US 2016280344A1
- Authority
- US
- United States
- Prior art keywords
- data
- cross
- generate
- sync
- location
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004891 communication Methods 0.000 title claims abstract description 53
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 4
- 230000001934 delay Effects 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 3
- 238000000034 method Methods 0.000 description 40
- 230000008569 process Effects 0.000 description 35
- 238000012937 correction Methods 0.000 description 20
- 230000006870 function Effects 0.000 description 19
- 238000010586 diagram Methods 0.000 description 17
- 230000009189 diving Effects 0.000 description 10
- 238000011084 recovery Methods 0.000 description 7
- 230000005540 biological transmission Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 239000003999 initiator Substances 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000003491 array Methods 0.000 description 3
- 239000000872 buffer Substances 0.000 description 3
- 239000004020 conductor Substances 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- 101001122448 Rattus norvegicus Nociceptin receptor Proteins 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005266 casting Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 238000005538 encapsulation Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 235000014653 Carica parviflora Nutrition 0.000 description 1
- 241000243321 Cnidaria Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 208000013407 communication difficulty Diseases 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 230000006837 decompression Effects 0.000 description 1
- 230000009429 distress Effects 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000003278 mimic effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000010079 rubber tapping Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000010454 slate Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B63—SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
- B63C—LAUNCHING, HAULING-OUT, OR DRY-DOCKING OF VESSELS; LIFE-SAVING IN WATER; EQUIPMENT FOR DWELLING OR WORKING UNDER WATER; MEANS FOR SALVAGING OR SEARCHING FOR UNDERWATER OBJECTS
- B63C11/00—Equipment for dwelling or working underwater; Means for searching for underwater objects
- B63C11/02—Divers' equipment
- B63C11/26—Communication means, e.g. means for signalling the presence of divers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S1/00—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith
- G01S1/72—Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith using ultrasonic, sonic or infrasonic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S11/00—Systems for determining distance or velocity not using reflection or reradiation
- G01S11/14—Systems for determining distance or velocity not using reflection or reradiation using ultrasonic, sonic, or infrasonic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
- G01S3/802—Systems for determining direction or deviation from predetermined direction
- G01S3/808—Systems for determining direction or deviation from predetermined direction using transducers spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
- G01S3/802—Systems for determining direction or deviation from predetermined direction
- G01S3/808—Systems for determining direction or deviation from predetermined direction using transducers spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
- G01S3/8083—Systems for determining direction or deviation from predetermined direction using transducers spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems determining direction of source
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/02—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
- G01S5/14—Determining absolute distances from a plurality of spaced points of known location
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B11/00—Transmission systems employing sonic, ultrasonic or infrasonic waves
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B13/00—Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
- H04B13/02—Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy
Definitions
- aspects of the present invention relate to locating underwater divers using acoustic signals.
- aspects of the present invention relate to the use of acoustic transducers and associated signal processing techniques to transmit and receive signals underwater in order to locate a diver or other object in three dimensions.
- Scuba diving is a unique and enjoyable recreational experience. It is estimated that, every day, over 75,000 persons participate in the sport at thousands of diving resorts and operations worldwide. In addition, various commercial and military operations utilize scuba divers to perform activities as search and rescue, salvage, underwater construction and repair activities, and military reconnaissance.
- one aspect of the uniqueness of diving is that underwater communication is extremely limited and communication between divers and the surface is almost non-existent. Visual interaction between divers is more often than not impossible, particularly when distances between divers increase and when divers are at different depths. Divers within a few feet of each other, but at different depth planes, may not have any visibility of the other diver. While the diving industry has developed several techniques to facilitate underwater communication (e.g. hand signals, writing on a slate, tapping on one's tank, and some electronic communication devices), most of them require close proximity for the communicating divers and, in the case of electronic means, are prohibitively expensive for all but industrial divers and military operations.
- underwater communication presents complicated problems due to the tendency of electromagnetic and ultrasonic waves to attenuate or otherwise deteriorate when travelling or propagating underwater.
- location information as well as the need to communicate text-based messages, provides another hurdle that cannot easily be solved with known underwater communication systems.
- the need to add security or otherwise covert communication features to these devices adds yet another level of complication for acoustic-based communication systems.
- An underwater communication and location device includes an array of omni-directional transducers and a transmitter.
- the transmitter is adapted to generate a code element; generate a data element; generate a modulation element; combine the code element, the data element and the modulation element into an analog waveform; and transmit the analog waveform through at least one of the transducers.
- a communication system includes a transmitter, a receiver and a processor.
- the processor is adapted to receive incoming data, generate cross-correlation data, pass the cross-correlation data into a block delay calculator and calculate the average phase delay between the cross-correlation data.
- a communication system including a transmitter, a receiver and a processor.
- the processor is adapted to, align a target synchronization pattern with incoming data by cross-correlating segments of incoming data with corresponding segments of the target synchronization pattern to form block cross-correlation data.
- An underwater communication system includes a plurality of location and messaging units for sending and receiving ultrasonic energy through water and a transceiver.
- Each of the location and messaging units has at least one ultrasonic transducer.
- the transceiver is adapted to generate a code element; generate a data element; generate a modulation element; combine the code element, the data element and the modulation element into an analog wave form; transmit the analog wave form through at least one of the ultrasonic transducers of at least one of the messaging units; and receive the analog waveform generated by at least one of the location and messaging units.
- FIG. 1 is a schematic overview of a typical diving scenario
- FIGS. 2A and 2B are of a form used for casting a negative-image form that is used to accurately position the transducers at their proper locations prior to encapsulation used in accordance with aspects of the present invention
- FIG. 3 shows a transducer array prior to encapsulation used in accordance with aspects of the present invention
- FIGS. 4A and 4B show a cross section of an encapsulated transducer array used in accordance with aspects of the present invention
- FIG. 5 shows a constellation map of a modulation scheme used in a communication channel example as used in connection with aspects of the present invention
- FIG. 6A shows a transmitter block diagram used in connection with aspects of the present invention
- FIG. 6B shows a detailed transmitter block diagram used in connection with aspects of the present invention.
- FIG. 7 shows a receiver block diagram used in connection with aspects of the present invention.
- FIG. 8 shows a find/align sync block diagram used in connection with aspects of the present invention
- FIG. 9 shows a find sync block diagram used in connection with aspects of the present invention.
- FIG. 10 shows an align sync block diagram used in connection with aspects of the present invention
- FIG. 11 shows a cross-correlation block diagram used in connection with aspects of the present invention.
- FIG. 12 shows a multi-channel align sync section block diagram used in connection with aspects of the present invention
- FIG. 13 shows a phase correction block diagram used in connection with aspects of the present invention
- FIG. 14 shows an accumulate block diagram used in connection with aspects of the present invention
- FIG. 15 shows a Costas loop block diagram used in connection with aspects of the present invention.
- FIG. 16 shows a chart showing the magnitudes of several align sync correlations in accordance with one embodiment of the present invention
- FIG. 17 shows a chart showing the cross-correlations of the magnitudes of the align sync correlations in accordance with one embodiment of the present invention
- FIG. 18 is a graphic representation of a transducer array in accordance with aspects of the present invention and how a source signal is received and the several individual transducers in the array;
- FIG. 19 shows the Process Data Section used in connection with aspects of the present invention.
- FIG. 1 a generalized schematic overview 100 is shown of an underwater diving scenario involving two or more divers (three in the example of FIG. 3 ) and a surface vessel.
- each of the surface vessel 102 , the surfaced diver 106 , and each of the submerged divers 110 and 114 carry or are otherwise equipped with a device that can both generate and receive an acoustic signal, such as an array of acoustic transducers. While not shown in FIG. 1 , each of the acoustic transducer arrays produces and receives an omni-directional acoustic signal.
- FIG. 1 a generalized schematic overview 100 is shown of an underwater diving scenario involving two or more divers (three in the example of FIG. 3 ) and a surface vessel.
- each of the surface vessel 102 , the surfaced diver 106 , and each of the submerged divers 110 and 114 carry or are otherwise equipped with a device that can both generate and receive an acoustic signal, such as an array of acoustic transducers.
- each of the signal devices such as the acoustic transducer array described above, are adapted to both generate and receive an acoustic signal and coincidentally determine the location of the device that generated that signal, both in terms of angular position and linear distance.
- an array 300 in this example, consists of spherical piezo transducers 304 , 306 , 308 and 310 held in a tetrahedral array configuration by an inner array structure 302 .
- the inner array structure 302 is generated in a casting process as a negative-image form by using the tree-like structure 200 , found in FIGS. 2A and 2B .
- the locations of the piezo transducers 304 , 306 , 308 and 310 are defined by the positions of the balls 202 , 204 , 206 and 208 on the ends of the tree-like structure 200 , therefore it is important that the balls 202 , 204 , 206 and 208 are accurately positioned.
- the transducer array, 300 containing the actual piezo transducers 304 , 306 , 308 and 310 , is encapsulated in a protective material
- FIGS. 2A and 2B provide details about the geometric orientation of the plurality of individual transducers that comprise the array 300 . As shown in FIGS.
- transducers 304 , 306 , 308 and 310 whose locations are defined by the balls 202 , 204 , 206 and 208 are arranged in a generally tetrahedral configuration and are maintained in this configuration though the use of a support structure 210 . While specific dimensional data is provided with the example in FIGS. 2A and 2B , it should be understood that the nominal dimensions between the individual transducers may range in various embodiments, where a minimum of four non-coplanar transducers are required for purposes of locating another similarly equipped diver or vehicle in three dimensional space, resulting in a unique geometric solution for said diver or vehicle. The dimensioning shown in FIGS.
- FIG. 2A and 2B is relative and can be applied to any scale, as long as the relative distances between the individual transducers remains in alignment with the shown dimensions.
- these relative dimensions are only one embodiment.
- FIG. 2B shows a 120° angle between the transducers and where the transducers are all located the same distance from the center of the array, in other embodiments, that distance from the center may range from less than a centimeter to over a meter, depending on the frequencies that are used for communication.
- Other embodiments may include other non-tetrahedral configurations, as long as there are a minimum of four transducers that are not all contained within the same plane. Additional transducers; exceeding four may be used, and may offer signal processing benefits over just four transducers.
- an encapsulated transducer array 300 is shown that includes the four transducers 304 , 306 , 308 , and 310 held in geometric orientation by a pliable but sturdy inner-array form 302 preferably made out of an acoustically neutral material such as Rho-C.
- Rho-C is a polymeric material with acoustic properties that very closely mimic those of water.
- the formed transducer array 300 is then encapsulated and incorporated into a more complete form with a housing that contains electronics and cabling for use in connection with an operational device.
- FIGS. 4A and 4B show one embodiment of a completed transducer array as an assembly 400 that can be connected to the associated electronics and control mechanisms described below.
- Transducers 304 , 306 and 308 are connected via conductors 402 through a channel 406 , mechanical mounting assemblies 410 and 412 where they emerge in a conductor cluster 415 .
- the completed assembly may be mounted or otherwise engaged on a dive computer, board or other mechanism that is ported by an underwater diver.
- conductor assembly 415 is coupled to an electronics package that performs one or more of the functions described in conjunction with FIGS. 6-15 below.
- the outgoing data needs to be processed into a form ready for transmission, and the incoming data needs to be decoded and translated into information that indicates the distance and vector of the source as well as to identify the source (e.g. a particular diver from among several on a dive or underwater mission).
- a communications channel and associated software and algorithms may be utilized to package the data for transmission and to interpret the data coming from the transducers.
- the transmit function 600 represents the first step in a communications channel protocol implemented by a diver location system constructed in accordance with aspects of the present invention and begins the process of identifying a particular diver and establishing a communications link for both messaging and location identification.
- an ASCII text is interleaved at 620 in order to distribute burst errors over the entire transmission, and thereby turn them into shorter errors that are correctable by a forward error correcting code encoder 630 .
- Barker code 640 is pre-pended to the data, then sync data is pre-pended to this array of data in the form of a stream of 1's (or 0's), after which the full array of data is combined with (mixed or added with) a pseudo random pattern generated at 650 .
- the resulting pseudo-random data is translated into hexadecimal format at 655 in order to choose the appropriate modulated wavelet at 660 for a given pseudo-random data hexadecimal symbol value.
- Establishing the pseudo-random communication code 650 is a function that is performed during an initialization procedure where communication codes are established in the form of a lookup table.
- the wavelet modulation setup block 660 also performed during initialization, develops 16 modulation wavelets, one for each symbol type that can be transmitted. In this embodiment, these wavelets are quantized and scaled to provide four transmit amplitude ranges that vary by factors of 2 ⁇ in voltage (4 ⁇ for power) to be selected as needed to either conserve power or increase the transmit distance, using the Read-Only Memory (ROM) tables described below.
- ROM Read-Only Memory
- FIG. 6B shows a more detailed version of the transmit block diagram 600 initially shown in FIG. 6A , expanding on the details of the pseudo-random communication code 650 and the modulation wavelet setup function 660 .
- Pseudo-random communication code 650 further includes a pattern generator 651 and a hex convertor coupled with a memory table 653 .
- pseudo-random communication code 650 generates a preset pattern for all users of the system to communicate to/with.
- Output of the pseudo-random communication code 650 is added to the data code coming from the data transmit function 605 .
- a data bit is added to each pseudo-random hex bit coming from the pseudo-random communications code 650 .
- the modulation wavelet setup block 660 is used to pre-construct modulated wavelets, where in this embodiment constitute a form of 16-Quadrature Amplitude Modulation (16-QAM), according to the constellation map of complex in-phase and quadrature components in FIG. 5 , showing how each hexadecimal value is modulated, in terms of amplitude and phase.
- the wavelet selector 669 is used to select modulated wavelets based on the values of the hex-formatted pseudo-random data 652 .
- the modulated pseudo-random data now goes into a digital-to-analog converter 670 and then transmit driver 675 before being passed to the transmit piezo transducer at 680 .
- Modulator block 660 includes a constellation map 661 , carrier 665 , modulator 662 , quantizer 663 and a series of ROM tables 664 , one for each transmit amplitude for an embodiment where various amplitudes are used, mostly in the interest of conserving power while only transmitting short distances.
- modulator block 660 constructs different wavelets (16 in this example according to the matrix shown in FIG. 5 ) that each represents a potential wavelet signal for transmitting by the piezo transducer.
- four communication channels are used as an example, but it is contemplated that the scope of the invention can be expanded for other numbers of communication channels.
- Range selector 668 takes the output of the modulation wavelet setup block (for the desired transmit amplitude) and feeds that data to a wavelet selector 669 .
- Wavelet selector 669 uses information comprising the pseudo-random hex code from 652 for selecting the appropriate modulation wavelet from 668 and the resulting information is passed through a digital to analog converter 670 , then through a driver 675 and on to one or more piezo transducers at 680 as a transmitted analog waveform.
- a text message sent to a diver such as “jump in” is converted to hexadecimal ASCII format at 610 and is then interleaved, bit by bit at 620 .
- Interleaving mixes up the data to provide robustness against large acoustical spikes that might otherwise overwhelm the data recovery channel.
- a forward error correcting code FEC
- the FEC routine creates a 1 bit delay shift register of the data and a 2 bit delay shift register. From the three streams of bits (including un-delayed), two XOR arrays are created.
- XOR1 performs an exclusive-or operation on the undelayed data array, the 1 bit delay and the 2-bit delay, resulting in the XOR1 array.
- the XOR2 is produced by an exclusive-or of the undelayed data array and the 2-bit delay array.
- the two XOR arrays are then interleaved, bit by bit, to produce the FEC encoded data coming out of 630
- the Barker code is added at 640 , then the sampling rate for the Barker and FEC data is increased by a factor of 4 in this case (where a “0” translates to “0000” and “1” translates to “1111”).
- the sync data is pre-pended to this array in the form of a stream of 1's (or 0's), and now the entire data stream is ready to be combined with (modulo-2 added or mixed with) the pseudo-random number (PRN) from 650 .
- PRN pseudo-random number
- Data from 605 is stepped through one bit at a time, and the pseudo-random pattern (PRN) one nibble (4 bits) at a time. If the data bit is 0, the next nibble of the PRN is added to the output array. If the data bit is 1, the next nibble of the PRN is inverted and it is added to the output. In this example the PRN is mixed with the data stream.
- the final step of this process which produces data that can then be transmitted is to step through the data array, nibble by nibble, and to use these hex values to index an array of mod-wavelets (containing segments of modulated waveforms).
- the mod-wavelet array is in one example 16 columns by 72 rows of samples, so for every nibble (or hex value) one of 16 columns is selected. So the algorithm takes a nibble of data at a time and writes the 72 sample values corresponding to the nibble column. This process of compiling a collection of mod-wavelets is continued for the entire length of data transmitted.
- a receive block diagram 700 is shown that illustrates an embodiment of how four-channels of piezo transducers interpret source diver coordinates 754 and receive user data 756 from another diver or location system.
- Piezo array 702 comprises four separate piezo transducers and thus generates a four-channel signal represented by communication channels 704 a - d .
- DASS, DSP and micro-Blaze (or another soft core processor) portions of the receive function are separated in FIG. 7 but can all be implemented in a single board, FPGA or functioning ASIC processor.
- Analog-to-digital converters 706 a - d pass the received four channel signal to a data acquisition module 712 .
- SRAM 708 and FIFO unit 710 comprise a diagnostics unit for analyzing communications channel data coming from the piezo array 702 .
- Communication channels exit the data acquisition module 712 at 713 and are processed by multipliers 726 , root mean squared calculation module 714 , a hardware gain select module 778 and output to a four channel phase lock loop 724 , such as a Costas loop.
- a module 720 for selecting the maximum RMS value from 714 which is then passed to a module 722 for calculating the soft gain adjustment that is looped back to the multipliers 726 , in order to maintain a fairly constant signal amplitude.
- Output from the phase-lock loop 724 (detailed in FIG. 15 ) is demodulated and filtered data, also called base-band filtered data, or BB-filtered cycle data 728 , which is then sent to the find/align sync module 800 (See FIG. 8 ).
- frequency correction data 730 is fed back to phase lock loop 724 as needed.
- a signal indicating that start of data was found, start data accum 734 is used to start processes in the phase correction block 748 (detailed in FIG. 13 ) and the accumulate data block 746 (detailed in FIG. 14 ).
- Equalization filter values 732 is also passed to both the phase correction block 748 and the accumulate data block 734 .
- the number corresponding to the selected primary channel used for data recovery (channel select index) 740 is passed to selector 742 for the purpose of selecting that channel of BBF data 744 (the one with the best signal) to be passed to phase correction block 748 and data accumulator 746 .
- the standard deviation (or RMS value) of BBF data 738 is used for proper data scaling within the accumulate data process 746 .
- the find/align sync module 800 passes the align-sync cross-correlations 736 to the triangulation block 750 to calculate the source diver coordinates 754 .
- Text Messaging Data is processed at 752 from data collection module 746 before being sent out as user text data 756 .
- the messaging component of the communication channel shown as “User data” is based on data collection from just one of the four channels shown in the example of FIG. 8 .
- Data accumulator 746 uses the messaging data (selected BBF data 744 ) from that data stream.
- the Phase-Lock Loop/Demodulator (in the form of a Costas Loop) 1500 is shown in more detail.
- the upper portion of the loop demodulates the inphase (real) component of raw receive data.
- the lower portion of the loop demodulates the quadrature (imaginary) component of raw receive data.
- Each component is multiplied by either a cosine or a sine function that comes from a numerically-controlled oscillator (NCO), which is basically a process consisting of cosine and sine functions that have phase and frequency parameters incorporated in them.
- NCO numerically-controlled oscillator
- the data is then low-pass filtered to preserve only the demodulated waveform, after which consecutive groups of four samples each (that make up a carrier cycle for this embodiment) are averaged resulting in a stream of complex data, now running at approximately the carrier cycle rate, one value per cycle.
- This complex data is then base-band filtered to improve symmetry and to reject noise, where the result is referred to as BB-filtered cycle data 728 ; these are complex values where both magnitude and phase are still preserved.
- This data is used extensively throughout find sync, align sync and data recovery. After sync has been found and aligned, frequency correction and the initial phase correction to the carrier is adjusted according to phase and frequency errors measured and averaged over the recovered Barker code data. After Barker code data, continual phase corrections are made based on incoming user data.
- find/align sync 800 block diagram is shown in more detail.
- Find/align sync is one component in the receive functionality of the present embodiment.
- Various components from find/align sync module 800 that pass back to or are incoming from receive module 700 are shown within FIG. 8 , such as start data accumulator 734 , BB-Filtered cycle data 728 , align-sync cross-correlations 736 , channel selector 740 , frequency corrector 730 , and equalization filter 732 .
- find/align sync module 800 the four communication channels of BB-filtered cycle data 728 enter from 724 into a series of four find sync/align sync modules shown as 802 , 804 , 806 and 808 .
- Each of the modules 802 - 808 include a channel find block (‘a’ in each block) and a channel align block (‘b’ in each block) working in concert with each other to calculate and output four functions that enable a sync pattern to be found and aligned.
- Outputs from each of the find sync/align sync blocks include a sync found signal 810 a - 810 d , an align complete signal 812 a - 812 d , a block delay value 814 a - 814 d , and align sync cross correlations 816 a - 816 d .
- Each of the sync found signals are passed to a sync found logic element 822 .
- Each of the align complete signals are passed to an align complete logic element 824 . If sync pattern is found at 822 , the resulting signal is passed through the align sync blocks until sync alignment is complete at 824 . If sync is not found at 822 (find sync), the process is repeated until sync is found.
- align sync cross-correlation waveforms 736 are passed to the triangulation section 750 to determine where the source diver is located.
- the align sync cross-correlation waveforms 736 are also used in Multi-Channel Align Sync 830 to identify the primary channel used for text data recovery and to properly align that data.
- Data read from the align sync modules 802 b - 808 b are also fed to a multi-channel align sync 830 which exports signals such as selected channel number 740 , frequency correction 730 and equalization filter values 732 .
- Align sync cross-correlations 736 is also output from the align sync block in order to triangulate the source diver location at 750 ( FIG. 7 ).
- a detailed embodiment of a find-sync process 900 is shown.
- the functions within dashed-line labeled 802 a correspond to any of find sync blocks 802 a - 808 a shown previously in FIG. 8 .
- BB-filtered cycle data (base-band filtered data) is received from 728 into a buffer 902 .
- One or more sync cross-correlation blocks 904 a - 904 n pass block cross-correlation data (the result of cross-correlating a segment of the complex target sync pattern with a segment of incoming complex base-band filtered data 728 , as detailed in FIG. 11 ) into the block delay calculation section 908 .
- functions such as multipliers, normalizers, summation and conjugation are applied to the complex block cross-correlations from 904 a - 904 n to calculate the average phase delay between adjacent block cross-correlations (also called the average block delay).
- This average block delay “A” is used to generate an array of complex unity-length rotational vectors ( ⁇ 1 through ⁇ 11 ) that are multiplied with each of the block cross-correlations 904 a - 904 n within the block delay compensation section 910 to individually correct for the delays inherent in each of the block cross-correlations 904 a - 904 n , due to Doppler shift and other sources of frequency offset.
- sync found signal 914 is eventually generated, as shown in FIGS. 8 and 9 as 810 a - 810 d , if and when the magnitude of the summed cross-correlations exceed a predefined threshold (shown as an input to a comparator 912 ).
- Align sync is a simplified version of find sync that accurately aligns the sync pattern with the complex BB-filtered cycle data 728 .
- the functions within the dashed-line labeled as 802 b - 808 b correspond to any of align sync blocks 802 b - 808 b shown previously in FIG. 8 .
- BB-filtered cycle data is received from 724 into a buffer 1002 .
- One or more sync cross-correlation blocks 1004 a - 1004 n pass the complex block cross-correlation data into the block delay compensation process 1008 .
- complex functions such as multipliers, normalizers and summation are used to process the data, where the complex rotational vectors ( ⁇ 1 through ⁇ 11 ) generated earlier during find sync 900 are multiplied with each of the block cross-correlations 1004 a - 1004 n within the block delay compensation section 1008 to individually correct for the delays inherent in each of the complex block cross-correlations 1004 a - 1004 n , due to Doppler shift and other sources of frequency offset. At this point the individually compensated block cross-correlations can be summed and yield outputs 816 a - d .
- FIG. 11 shows the details of an embodiment of the Sync Cross-Correlation Blocks, such as block 1004 a shown in FIG. 10 , as they are used in one channel of the find sync and align sync processes, to perform complex block cross-correlations between the target sync pattern and BB-filtered cycle data.
- FIG. 12 shows the details of an embodiment of the multi-channel align sync section of the receive block diagram.
- the multi-channel align sync process 1200 described in FIG. 12 identifies the absolute maximum value of each of the align sync cross-correlations 816 a - d from each channel.
- the dominant channel is the one with the largest absolute maximum cross-correlation value, identified as a selected channel index 1222 after a max selection process.
- Another method can be used by taking the largest amplitude channel identified during the gain control process. This selected channel index is used in various selection processes for selecting the various types of data corresponding to that channel.
- the align sync cross-correlations 816 a - 816 d of the four channels are each held in separate cross-correlation registers 1202 a - 1202 d , respectively.
- these registers are 48 values in length over four channels and are processed through a function that calculates the absolute peak value of each channel at 1204 a - 1204 d , after which the maximum of the values coming from 1204 a - d is determined, resulting in 1222 .
- a selection process 1206 such as a multiplexer switch uses the selected channel index 1222 , and finds the maximum signal and passes it through a process 1220 that subsequently determines the 12 largest-magnitude contiguous values that define a particular symbol.
- the indexing for these contiguous values is passed to selection process 1208 , an equalization filter 1209 , which takes the complex conjugates of the 12 contiguous values and normalizes them to determine the complex values of the equalization filter.
- This equalization filter is applied later to BB-filtered data for data recovery.
- Channel block delays 814 a - 814 d are passed through a selection multiplexer 1214 (using selected channel index 1222 or 740 in FIG. 7 ) and into a frequency correction process 1218 .
- frequency correction process 1218 takes the angle of the selected block delay in the complex plane, then multiplies that number by a constant to obtain a normalized frequency correction factor.
- FIG. 13 shows one embodiment of a phase correction process 1300 , such as the phase correction 748 shown previously in FIG. 7 .
- the BB-filtered cycle data 728 of the selected channel, as selected by the channel index 1222 is passed to this correction process 1300 , where the phase angle is determined in the angle process 1302 by taking the four-quadrant arctangent of each value in the complex plane of BB-filtered data; this is considered the current phase angle of the incoming data.
- the alignment of the incoming data relative to the target sync pattern, consisting of pseudo-random data (PRN data) that was used to generate the transmitted waveform is known, which allows the current phase angle of the incoming data to be paired with the ideal phase angle associated with the modulated PRN data.
- PRN data pseudo-random data
- phase correction is set to zero.
- a known Barker code that is not mixed with the variable user data.
- the Barker Code values that were known to be transmitted are used to set the initial value for phase correction.
- Each value of the Barker determines the switch position 1316 to be selected, either the error associated with the non-inverted symbol or that of the inverted symbol. For example, a Barker value of “0” means use that the transmitted symbol was originally uninverted, therefore switch to the uninverted error. In a similar manner, a value of “1” means use the inverted error.
- Barker is the remaining unknown incoming data, where it is not known whether to use the uninverted or the inverted error; for this, use the switch position pointing to the smallest error.
- the resulting phase error values are added to a moving average value from 1318 , to generate the phase correction value 1320 that is passed to the Costas Loop of the selected channel to continue with phase-locked data recovery.
- FIG. 14 shows one embodiment of an accumulate data process 1400 such as the data accumulation process 746 shown in FIG. 7 .
- Selected BB-filtered data 744 running at the cycle-rate is multiplied by the equalization filter 732 .
- This multiplication involves multiplying a pair of complex numbers for every cycle over the entire symbol being evaluated to generate equalized cycle-rate data 1402 .
- All of the equalized cycle values for an entire symbol are averaged together to generate equalized symbol-rate data, which is then normalized by the ratio of the root-mean-square of the sync symbols over a statistically large sample size divided by the root-mean-square of the BB-filtered values over a statistically large sample size, to obtain normalized and equalized symbol-rate data 1412 .
- the process data function 1900 takes the PRN-removed symbol-rate data 1410 and reduces it to user data.
- the data 1410 is separated into real and imaginary components. Values for components for all of the values that constitute an encoded bit are averaged together to form encoded data 1902 .
- the data is further processed using a Viterbi decoder to generate decoded data 1904 .
- user data 1906 is recovered after a de-interleaving process.
- FIG. 16 shows the align sync cross-correlated data for one embodiment, referenced earlier as 736 in FIG. 7 .
- This data in its complex form was generated by the align sync process, through which the target sync pattern was cross-correlated with bb-filtered cycle data 728 .
- the data 736 is shown here with the phase information removed, as only absolute values. Where the peaks occur in this data 736 , in FIG. 16 , is where bb-filtered cycle data 728 for each channel is maximally aligned with the sync pattern.
- the time delay between these waveforms can be used through a geometric triangulation process, and by knowing the relative positions of the transducer array elements in three dimensions, and by knowing the speed of sound in the medium, a determination can be made as to the direction of the source of this transmission.
- using this align sync cross-correlated data 736 in its raw form can be subject to unnecessary errors attributed to the asymmetries and irregularities in these waveforms.
- Additional accuracy can be obtained by employing second degree polynomial fitting to accurately locate the peaks.
- Another method can be used by measuring the relative locations of the centroids of each waveform. Regardless of the specific means to determine the relative delays between channels, that relative delay data is applied to a geometric model of the transducer array to determine the direction to the source diver responsible for the transmission being received. The distance to the source diver is determined by measuring the round-trip time-of-flight of transmissions passed between the two divers and by knowing the propagation speed of sound in the medium.
- FIG. 18 shows an embodiment of the transducer array 1800 with four individual transducers 1812 a , 1814 a , 1816 a and 1818 a on or near the initiating diver or vessel (an initiator), along with a point source 1810 that represents a signal originating from a transmitter located on another diver or a surface vessel (a responder). As shown by distance lines 1812 b , 1814 b , 1816 b and 1818 b , the linear distance between the responder point source 1810 and each of the transducers (on the initiator) is different.
- the direction and distance between the transducer array 1800 and the point source 1810 can be calculated.
- Embodiments of a communication system described herein may be incorporated into one of many different types of underwater dive systems that may be utilized by recreational, commercial, industrial and military industries. Aspects of a communication system constructed in accordance with the present invention may be incorporated into a handheld dive computer underwater dive platform, or directly onto the BCD (buoyancy control device) of a diver.
- BCD buoyancy control device
- Visual display means may also be incorporated into such a system to enable a diver to see and send messages to one or more other divers while utilizing aspects of the communication system described herein.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Mechanical Engineering (AREA)
- Ocean & Marine Engineering (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
Description
- This application is a continuation of U.S. patent application Ser. No. 13/445,108 filed Apr. 12, 2012 and entitled “UNDERWATER ACOUSTIC ARRAY, COMMUNICATION AND LOCATION SYSTEM,” the entire disclosure of which is hereby incorporated by reference for all purposes, as if fully set forth herein.
- This material is based upon work supported by the United States Special Operation Command under Contract No. H92222-10-C-0018.
- Aspects of the present invention relate to locating underwater divers using acoustic signals. In particular, aspects of the present invention relate to the use of acoustic transducers and associated signal processing techniques to transmit and receive signals underwater in order to locate a diver or other object in three dimensions.
- Scuba diving is a unique and enjoyable recreational experience. It is estimated that, every day, over 75,000 persons participate in the sport at thousands of diving resorts and operations worldwide. In addition, various commercial and military operations utilize scuba divers to perform activities as search and rescue, salvage, underwater construction and repair activities, and military reconnaissance.
- Operationally, one aspect of the uniqueness of diving is that underwater communication is extremely limited and communication between divers and the surface is almost non-existent. Visual interaction between divers is more often than not impossible, particularly when distances between divers increase and when divers are at different depths. Divers within a few feet of each other, but at different depth planes, may not have any visibility of the other diver. While the diving industry has developed several techniques to facilitate underwater communication (e.g. hand signals, writing on a slate, tapping on one's tank, and some electronic communication devices), most of them require close proximity for the communicating divers and, in the case of electronic means, are prohibitively expensive for all but industrial divers and military operations.
- Limited underwater visibility exacerbates communication difficulty. The best recreational ocean diving sites around the world may have 150-200 feet of underwater visibility. Most have 60 feet or less of visibility. With few exceptions, visibility at inland sites such as lakes, rivers and quarries drops below 20 feet. Considering that ocean currents in many dive areas may have a velocity of one to two knots, a diver in perfect visibility conditions can drift out of sight in less than 60 seconds.
- Visibility is further reduced by underwater topography which may include coral or rock formations. Night-diving conditions obviously limit communication even further and complicates the divemaster's supervisory responsibilities.
- Major scuba-diving-certifying organizations have attempted to mitigate these risks by establishing well-accepted rules: always dive with and stay close to a “buddy”; evaluate conditions carefully and seek orientation with a local dive shop before diving; plan the dive carefully, follow the plan once underwater, surface when one becomes separated from the group; etc. The fact remains, though, that the communication options available to the average recreational diver when in distress or when separated from the group are extremely limited. In contrast with many land-based activities, divers do not have the option of carrying emergency rescue beacons or other long-range communication options. At best, some recreational divers carry only a simple whistle or inflatable tube for use at the surface.
- Simple and unsophisticated antenna-based location systems have been developed in the past, such as those found in U.S. Pat. No. 7,388,512, but fail to provide reliable and accurate directional and distance location information that is useful in the low-visibility environments described above.
- Furthermore, underwater communication presents complicated problems due to the tendency of electromagnetic and ultrasonic waves to attenuate or otherwise deteriorate when travelling or propagating underwater. The need for location information, as well as the need to communicate text-based messages, provides another hurdle that cannot easily be solved with known underwater communication systems. Finally, the need to add security or otherwise covert communication features to these devices adds yet another level of complication for acoustic-based communication systems.
- An underwater communication and location device includes an array of omni-directional transducers and a transmitter. The transmitter is adapted to generate a code element; generate a data element; generate a modulation element; combine the code element, the data element and the modulation element into an analog waveform; and transmit the analog waveform through at least one of the transducers.
- A communication system includes a transmitter, a receiver and a processor. The processor is adapted to receive incoming data, generate cross-correlation data, pass the cross-correlation data into a block delay calculator and calculate the average phase delay between the cross-correlation data.
- A communication system including a transmitter, a receiver and a processor. The processor is adapted to, align a target synchronization pattern with incoming data by cross-correlating segments of incoming data with corresponding segments of the target synchronization pattern to form block cross-correlation data.
- An underwater communication system includes a plurality of location and messaging units for sending and receiving ultrasonic energy through water and a transceiver. Each of the location and messaging units has at least one ultrasonic transducer. The transceiver is adapted to generate a code element; generate a data element; generate a modulation element; combine the code element, the data element and the modulation element into an analog wave form; transmit the analog wave form through at least one of the ultrasonic transducers of at least one of the messaging units; and receive the analog waveform generated by at least one of the location and messaging units.
- Other embodiments will become known to one of skill in art after reading the following specification in conjunction with the figures and claims.
-
FIG. 1 is a schematic overview of a typical diving scenario; -
FIGS. 2A and 2B are of a form used for casting a negative-image form that is used to accurately position the transducers at their proper locations prior to encapsulation used in accordance with aspects of the present invention; -
FIG. 3 shows a transducer array prior to encapsulation used in accordance with aspects of the present invention; -
FIGS. 4A and 4B show a cross section of an encapsulated transducer array used in accordance with aspects of the present invention; -
FIG. 5 shows a constellation map of a modulation scheme used in a communication channel example as used in connection with aspects of the present invention; -
FIG. 6A shows a transmitter block diagram used in connection with aspects of the present invention; -
FIG. 6B shows a detailed transmitter block diagram used in connection with aspects of the present invention; -
FIG. 7 shows a receiver block diagram used in connection with aspects of the present invention; -
FIG. 8 shows a find/align sync block diagram used in connection with aspects of the present invention; -
FIG. 9 shows a find sync block diagram used in connection with aspects of the present invention; -
FIG. 10 shows an align sync block diagram used in connection with aspects of the present invention; -
FIG. 11 shows a cross-correlation block diagram used in connection with aspects of the present invention; -
FIG. 12 shows a multi-channel align sync section block diagram used in connection with aspects of the present invention; -
FIG. 13 shows a phase correction block diagram used in connection with aspects of the present invention; -
FIG. 14 shows an accumulate block diagram used in connection with aspects of the present invention; -
FIG. 15 shows a Costas loop block diagram used in connection with aspects of the present invention; -
FIG. 16 shows a chart showing the magnitudes of several align sync correlations in accordance with one embodiment of the present invention; -
FIG. 17 shows a chart showing the cross-correlations of the magnitudes of the align sync correlations in accordance with one embodiment of the present invention; -
FIG. 18 is a graphic representation of a transducer array in accordance with aspects of the present invention and how a source signal is received and the several individual transducers in the array; and -
FIG. 19 shows the Process Data Section used in connection with aspects of the present invention. - With reference to
FIG. 1 , a generalizedschematic overview 100 is shown of an underwater diving scenario involving two or more divers (three in the example ofFIG. 3 ) and a surface vessel. In the diving scenario ofFIG. 1 , each of thesurface vessel 102, the surfaceddiver 106, and each of the submergeddivers FIG. 1 , each of the acoustic transducer arrays produces and receives an omni-directional acoustic signal. InFIG. 1 ,surface vessel 102 producessignal 104, surfaceddiver 106 produces signal 108 and submergeddivers produce signals - With reference to
FIGS. 2A, 2B and 3 , anarray 300, in this example, consists of sphericalpiezo transducers inner array structure 302. Theinner array structure 302 is generated in a casting process as a negative-image form by using the tree-like structure 200, found inFIGS. 2A and 2B . The locations of thepiezo transducers balls like structure 200, therefore it is important that theballs piezo transducers FIGS. 2A and 2B provide details about the geometric orientation of the plurality of individual transducers that comprise thearray 300. As shown inFIGS. 2A, 2B and 3 ,transducers balls support structure 210. While specific dimensional data is provided with the example inFIGS. 2A and 2B , it should be understood that the nominal dimensions between the individual transducers may range in various embodiments, where a minimum of four non-coplanar transducers are required for purposes of locating another similarly equipped diver or vehicle in three dimensional space, resulting in a unique geometric solution for said diver or vehicle. The dimensioning shown inFIGS. 2A and 2B is relative and can be applied to any scale, as long as the relative distances between the individual transducers remains in alignment with the shown dimensions. However, these relative dimensions are only one embodiment. For example, whileFIG. 2B shows a 120° angle between the transducers and where the transducers are all located the same distance from the center of the array, in other embodiments, that distance from the center may range from less than a centimeter to over a meter, depending on the frequencies that are used for communication. Other embodiments may include other non-tetrahedral configurations, as long as there are a minimum of four transducers that are not all contained within the same plane. Additional transducers; exceeding four may be used, and may offer signal processing benefits over just four transducers. - With reference to
FIG. 3 , an encapsulatedtransducer array 300 is shown that includes the fourtransducers array form 302 preferably made out of an acoustically neutral material such as Rho-C. Rho-C is a polymeric material with acoustic properties that very closely mimic those of water. The formedtransducer array 300 is then encapsulated and incorporated into a more complete form with a housing that contains electronics and cabling for use in connection with an operational device.FIGS. 4A and 4B show one embodiment of a completed transducer array as anassembly 400 that can be connected to the associated electronics and control mechanisms described below.Transducers transducer 310 is hidden in the view ofFIG. 4A ) are connected viaconductors 402 through achannel 406, mechanical mountingassemblies conductor cluster 415. The completed assembly may be mounted or otherwise engaged on a dive computer, board or other mechanism that is ported by an underwater diver. In one embodiment,conductor assembly 415 is coupled to an electronics package that performs one or more of the functions described in conjunction withFIGS. 6-15 below. - While the individual transducers integrated into the transducer array are important for transmitting outgoing signals and capturing the incoming signals from the source diver or other transmitter, the outgoing data needs to be processed into a form ready for transmission, and the incoming data needs to be decoded and translated into information that indicates the distance and vector of the source as well as to identify the source (e.g. a particular diver from among several on a dive or underwater mission). Described below are various embodiments of a communications channel and associated software and algorithms that may be utilized to package the data for transmission and to interpret the data coming from the transducers.
- With reference to
FIG. 6A , a high-level block diagram of a transmitfunction 600 is shown. The transmitfunction 600 represents the first step in a communications channel protocol implemented by a diver location system constructed in accordance with aspects of the present invention and begins the process of identifying a particular diver and establishing a communications link for both messaging and location identification. At 610 an ASCII text is interleaved at 620 in order to distribute burst errors over the entire transmission, and thereby turn them into shorter errors that are correctable by a forward error correctingcode encoder 630.Barker code 640 is pre-pended to the data, then sync data is pre-pended to this array of data in the form of a stream of 1's (or 0's), after which the full array of data is combined with (mixed or added with) a pseudo random pattern generated at 650. The resulting pseudo-random data is translated into hexadecimal format at 655 in order to choose the appropriate modulated wavelet at 660 for a given pseudo-random data hexadecimal symbol value. - Establishing the
pseudo-random communication code 650 is a function that is performed during an initialization procedure where communication codes are established in the form of a lookup table. The waveletmodulation setup block 660, also performed during initialization, develops 16 modulation wavelets, one for each symbol type that can be transmitted. In this embodiment, these wavelets are quantized and scaled to provide four transmit amplitude ranges that vary by factors of 2× in voltage (4× for power) to be selected as needed to either conserve power or increase the transmit distance, using the Read-Only Memory (ROM) tables described below. -
FIG. 6B shows a more detailed version of the transmit block diagram 600 initially shown inFIG. 6A , expanding on the details of thepseudo-random communication code 650 and the modulationwavelet setup function 660. However, in general, like references fromFIGS. 6A and 6B describe similar elements within the overall transmitfunction 600.Pseudo-random communication code 650 further includes apattern generator 651 and a hex convertor coupled with a memory table 653. In one embodiment,pseudo-random communication code 650 generates a preset pattern for all users of the system to communicate to/with. Output of thepseudo-random communication code 650 is added to the data code coming from the data transmitfunction 605. In one embodiment, a data bit is added to each pseudo-random hex bit coming from thepseudo-random communications code 650. - The modulation
wavelet setup block 660 is used to pre-construct modulated wavelets, where in this embodiment constitute a form of 16-Quadrature Amplitude Modulation (16-QAM), according to the constellation map of complex in-phase and quadrature components inFIG. 5 , showing how each hexadecimal value is modulated, in terms of amplitude and phase. Thewavelet selector 669 is used to select modulated wavelets based on the values of the hex-formattedpseudo-random data 652. The modulated pseudo-random data now goes into a digital-to-analog converter 670 and then transmitdriver 675 before being passed to the transmit piezo transducer at 680. -
Modulator block 660 includes aconstellation map 661,carrier 665,modulator 662,quantizer 663 and a series of ROM tables 664, one for each transmit amplitude for an embodiment where various amplitudes are used, mostly in the interest of conserving power while only transmitting short distances. In one embodiment,modulator block 660 constructs different wavelets (16 in this example according to the matrix shown inFIG. 5 ) that each represents a potential wavelet signal for transmitting by the piezo transducer. In describing aspects of this invention, four communication channels are used as an example, but it is contemplated that the scope of the invention can be expanded for other numbers of communication channels. -
Range selector 668 takes the output of the modulation wavelet setup block (for the desired transmit amplitude) and feeds that data to awavelet selector 669.Wavelet selector 669 uses information comprising the pseudo-random hex code from 652 for selecting the appropriate modulation wavelet from 668 and the resulting information is passed through a digital toanalog converter 670, then through adriver 675 and on to one or more piezo transducers at 680 as a transmitted analog waveform. - With continuing attention to the transmit block diagram 600 shown in
FIG. 6B , an example of the data interleaving is described. For example, a text message sent to a diver such as “jump in” is converted to hexadecimal ASCII format at 610 and is then interleaved, bit by bit at 620. Interleaving mixes up the data to provide robustness against large acoustical spikes that might otherwise overwhelm the data recovery channel. Once the data is interleaved, a forward error correcting code (FEC) is added at 630. In one example, the FEC routine creates a 1 bit delay shift register of the data and a 2 bit delay shift register. From the three streams of bits (including un-delayed), two XOR arrays are created. XOR1 performs an exclusive-or operation on the undelayed data array, the 1 bit delay and the 2-bit delay, resulting in the XOR1 array. The XOR2 is produced by an exclusive-or of the undelayed data array and the 2-bit delay array. The two XOR arrays are then interleaved, bit by bit, to produce the FEC encoded data coming out of 630 - The Barker code is added at 640, then the sampling rate for the Barker and FEC data is increased by a factor of 4 in this case (where a “0” translates to “0000” and “1” translates to “1111”). The sync data is pre-pended to this array in the form of a stream of 1's (or 0's), and now the entire data stream is ready to be combined with (modulo-2 added or mixed with) the pseudo-random number (PRN) from 650. Data from 605 is stepped through one bit at a time, and the pseudo-random pattern (PRN) one nibble (4 bits) at a time. If the data bit is 0, the next nibble of the PRN is added to the output array. If the data bit is 1, the next nibble of the PRN is inverted and it is added to the output. In this example the PRN is mixed with the data stream.
- In this embodiment, the final step of this process which produces data that can then be transmitted is to step through the data array, nibble by nibble, and to use these hex values to index an array of mod-wavelets (containing segments of modulated waveforms). The mod-wavelet array is in one example 16 columns by 72 rows of samples, so for every nibble (or hex value) one of 16 columns is selected. So the algorithm takes a nibble of data at a time and writes the 72 sample values corresponding to the nibble column. This process of compiling a collection of mod-wavelets is continued for the entire length of data transmitted.
- With reference to
FIG. 7 a receive block diagram 700 is shown that illustrates an embodiment of how four-channels of piezo transducers interpret source diver coordinates 754 and receiveuser data 756 from another diver or location system.Piezo array 702 comprises four separate piezo transducers and thus generates a four-channel signal represented bycommunication channels 704 a-d. DASS, DSP and micro-Blaze (or another soft core processor) portions of the receive function are separated inFIG. 7 but can all be implemented in a single board, FPGA or functioning ASIC processor. Analog-to-digital converters 706 a-d pass the received four channel signal to adata acquisition module 712.SRAM 708 andFIFO unit 710 comprise a diagnostics unit for analyzing communications channel data coming from thepiezo array 702. Communication channels exit thedata acquisition module 712 at 713 and are processed bymultipliers 726, root mean squaredcalculation module 714, a hardware gainselect module 778 and output to a four channelphase lock loop 724, such as a Costas loop. Also present in the hardware is amodule 720 for selecting the maximum RMS value from 714 which is then passed to amodule 722 for calculating the soft gain adjustment that is looped back to themultipliers 726, in order to maintain a fairly constant signal amplitude. Output from the phase-lock loop 724 (detailed inFIG. 15 ) is demodulated and filtered data, also called base-band filtered data, or BB-filteredcycle data 728, which is then sent to the find/align sync module 800 (SeeFIG. 8 ). - From the find/align
sync module 800,frequency correction data 730 is fed back tophase lock loop 724 as needed. A signal indicating that start of data was found, start data accum 734, is used to start processes in the phase correction block 748 (detailed inFIG. 13 ) and the accumulate data block 746 (detailed inFIG. 14 ). Equalization filter values 732 is also passed to both thephase correction block 748 and the accumulatedata block 734. From find/alignsync module 800, the number corresponding to the selected primary channel used for data recovery (channel select index) 740 is passed toselector 742 for the purpose of selecting that channel of BBF data 744 (the one with the best signal) to be passed to phasecorrection block 748 anddata accumulator 746. The standard deviation (or RMS value) ofBBF data 738 is used for proper data scaling within the accumulatedata process 746. The find/alignsync module 800 passes the align-sync cross-correlations 736 to thetriangulation block 750 to calculate the source diver coordinates 754. Text Messaging Data is processed at 752 fromdata collection module 746 before being sent out asuser text data 756. The messaging component of the communication channel shown as “User data” is based on data collection from just one of the four channels shown in the example ofFIG. 8 .Data accumulator 746 uses the messaging data (selected BBF data 744) from that data stream. - Once data has been properly gain adjusted, at the outputs of the
multipliers 726, it is ready to be phase-locked and demodulated. With reference toFIG. 15 , the Phase-Lock Loop/Demodulator (in the form of a Costas Loop) 1500 is shown in more detail. The upper portion of the loop demodulates the inphase (real) component of raw receive data. The lower portion of the loop demodulates the quadrature (imaginary) component of raw receive data. Each component is multiplied by either a cosine or a sine function that comes from a numerically-controlled oscillator (NCO), which is basically a process consisting of cosine and sine functions that have phase and frequency parameters incorporated in them. The data is then low-pass filtered to preserve only the demodulated waveform, after which consecutive groups of four samples each (that make up a carrier cycle for this embodiment) are averaged resulting in a stream of complex data, now running at approximately the carrier cycle rate, one value per cycle. This complex data is then base-band filtered to improve symmetry and to reject noise, where the result is referred to as BB-filteredcycle data 728; these are complex values where both magnitude and phase are still preserved. This data is used extensively throughout find sync, align sync and data recovery. After sync has been found and aligned, frequency correction and the initial phase correction to the carrier is adjusted according to phase and frequency errors measured and averaged over the recovered Barker code data. After Barker code data, continual phase corrections are made based on incoming user data. - With reference to
FIG. 8 , the find/alignsync 800 block diagram is shown in more detail. Find/align sync is one component in the receive functionality of the present embodiment. Various components from find/alignsync module 800 that pass back to or are incoming from receivemodule 700 are shown withinFIG. 8 , such asstart data accumulator 734, BB-Filteredcycle data 728, align-sync cross-correlations 736,channel selector 740,frequency corrector 730, andequalization filter 732. - Within find/align
sync module 800, and for this embodiment, the four communication channels of BB-filteredcycle data 728 enter from 724 into a series of four find sync/align sync modules shown as 802, 804, 806 and 808. Each of the modules 802-808 include a channel find block (‘a’ in each block) and a channel align block (‘b’ in each block) working in concert with each other to calculate and output four functions that enable a sync pattern to be found and aligned. Outputs from each of the find sync/align sync blocks include a sync found signal 810 a-810 d, an align complete signal 812 a-812 d, a block delay value 814 a-814 d, and align sync cross correlations 816 a-816 d. Each of the sync found signals are passed to a sync foundlogic element 822. Each of the align complete signals are passed to an aligncomplete logic element 824. If sync pattern is found at 822, the resulting signal is passed through the align sync blocks until sync alignment is complete at 824. If sync is not found at 822 (find sync), the process is repeated until sync is found. Similarly, if align is complete at 824 (after sync was found), the resulting alignsync cross-correlation waveforms 736 are passed to thetriangulation section 750 to determine where the source diver is located. The alignsync cross-correlation waveforms 736 are also used inMulti-Channel Align Sync 830 to identify the primary channel used for text data recovery and to properly align that data. Data read from thealign sync modules 802 b-808 b are also fed to amulti-channel align sync 830 which exports signals such as selectedchannel number 740,frequency correction 730 and equalization filter values 732. Align sync cross-correlations 736 is also output from the align sync block in order to triangulate the source diver location at 750 (FIG. 7 ). - With reference to
FIG. 9 , a detailed embodiment of a find-sync process 900 is shown. As represented inFIG. 9 , the functions within dashed-line labeled 802 a correspond to any of find sync blocks 802 a-808 a shown previously inFIG. 8 . BB-filtered cycle data (base-band filtered data) is received from 728 into abuffer 902. One or more sync cross-correlation blocks 904 a-904 n pass block cross-correlation data (the result of cross-correlating a segment of the complex target sync pattern with a segment of incoming complex base-band filtereddata 728, as detailed inFIG. 11 ) into the blockdelay calculation section 908. As shown in 908, functions such as multipliers, normalizers, summation and conjugation are applied to the complex block cross-correlations from 904 a-904 n to calculate the average phase delay between adjacent block cross-correlations (also called the average block delay). This average block delay “A” is used to generate an array of complex unity-length rotational vectors (Â1 through Â11) that are multiplied with each of the block cross-correlations 904 a-904 n within the blockdelay compensation section 910 to individually correct for the delays inherent in each of the block cross-correlations 904 a-904 n, due to Doppler shift and other sources of frequency offset. From thedelay compensation process 910, sync foundsignal 914 is eventually generated, as shown inFIGS. 8 and 9 as 810 a-810 d, if and when the magnitude of the summed cross-correlations exceed a predefined threshold (shown as an input to a comparator 912). - With reference to
FIG. 10 , a detailed embodiment of an align-sync process 1000 is shown. Align sync is a simplified version of find sync that accurately aligns the sync pattern with the complex BB-filteredcycle data 728. As represented inFIG. 10 , the functions within the dashed-line labeled as 802 b-808 b correspond to any ofalign sync blocks 802 b-808 b shown previously inFIG. 8 . BB-filtered cycle data is received from 724 into abuffer 1002. One or more sync cross-correlation blocks 1004 a-1004 n pass the complex block cross-correlation data into the blockdelay compensation process 1008. As shown in 1008, complex functions such as multipliers, normalizers and summation are used to process the data, where the complex rotational vectors (Â1 through  11) generated earlier duringfind sync 900 are multiplied with each of the block cross-correlations 1004 a-1004 n within the blockdelay compensation section 1008 to individually correct for the delays inherent in each of the complex block cross-correlations 1004 a-1004 n, due to Doppler shift and other sources of frequency offset. At this point the individually compensated block cross-correlations can be summed and yield outputs 816 a-d. These align sync processes are repeated with each single-step shift of BB-filtered cycle data (rather than the more coarse shifting used to just find sync) in the vicinity of where sync was found. Having this high-resolution cross-correlation data allows for the incoming data to be accurately aligned for data recovery, and provides the cross-correlation waveforms needed for triangulation of the source diver location. Fromdelay compensation process 1008, block delays 814 a-814 d and the complex align sync delay-compensated cross-correlations 816 a-816 d are generated.Buffers -
FIG. 11 shows the details of an embodiment of the Sync Cross-Correlation Blocks, such asblock 1004 a shown inFIG. 10 , as they are used in one channel of the find sync and align sync processes, to perform complex block cross-correlations between the target sync pattern and BB-filtered cycle data. -
FIG. 12 shows the details of an embodiment of the multi-channel align sync section of the receive block diagram. The multi-channelalign sync process 1200 described inFIG. 12 identifies the absolute maximum value of each of the align sync cross-correlations 816 a-d from each channel. In this example, the dominant channel is the one with the largest absolute maximum cross-correlation value, identified as a selectedchannel index 1222 after a max selection process. Another method can be used by taking the largest amplitude channel identified during the gain control process. This selected channel index is used in various selection processes for selecting the various types of data corresponding to that channel. The align sync cross-correlations 816 a-816 d of the four channels are each held in separate cross-correlation registers 1202 a-1202 d, respectively. In one embodiment these registers are 48 values in length over four channels and are processed through a function that calculates the absolute peak value of each channel at 1204 a-1204 d, after which the maximum of the values coming from 1204 a-d is determined, resulting in 1222. Aselection process 1206, such as a multiplexer switch uses the selectedchannel index 1222, and finds the maximum signal and passes it through aprocess 1220 that subsequently determines the 12 largest-magnitude contiguous values that define a particular symbol. Once this sequence is found, the indexing for these contiguous values is passed toselection process 1208, anequalization filter 1209, which takes the complex conjugates of the 12 contiguous values and normalizes them to determine the complex values of the equalization filter. This equalization filter is applied later to BB-filtered data for data recovery. - Channel block delays 814 a-814 d (See
FIGS. 10 and 12 ) are passed through a selection multiplexer 1214 (using selectedchannel index FIG. 7 ) and into afrequency correction process 1218. In one embodiment,frequency correction process 1218 takes the angle of the selected block delay in the complex plane, then multiplies that number by a constant to obtain a normalized frequency correction factor. -
FIG. 13 shows one embodiment of aphase correction process 1300, such as thephase correction 748 shown previously inFIG. 7 . The BB-filteredcycle data 728 of the selected channel, as selected by thechannel index 1222 is passed to thiscorrection process 1300, where the phase angle is determined in theangle process 1302 by taking the four-quadrant arctangent of each value in the complex plane of BB-filtered data; this is considered the current phase angle of the incoming data. After align complete, the alignment of the incoming data relative to the target sync pattern, consisting of pseudo-random data (PRN data), that was used to generate the transmitted waveform is known, which allows the current phase angle of the incoming data to be paired with the ideal phase angle associated with the modulated PRN data. Initially, phase correction is set to zero. Directly following the sync pattern used for alignment, is a known Barker code that is not mixed with the variable user data. When phase correction is first calculated, the Barker Code values that were known to be transmitted are used to set the initial value for phase correction. Each value of the Barker determines theswitch position 1316 to be selected, either the error associated with the non-inverted symbol or that of the inverted symbol. For example, a Barker value of “0” means use that the transmitted symbol was originally uninverted, therefore switch to the uninverted error. In a similar manner, a value of “1” means use the inverted error. After Barker, is the remaining unknown incoming data, where it is not known whether to use the uninverted or the inverted error; for this, use the switch position pointing to the smallest error. The resulting phase error values are added to a moving average value from 1318, to generate thephase correction value 1320 that is passed to the Costas Loop of the selected channel to continue with phase-locked data recovery. -
FIG. 14 shows one embodiment of an accumulatedata process 1400 such as thedata accumulation process 746 shown inFIG. 7 . Selected BB-filtereddata 744 running at the cycle-rate is multiplied by theequalization filter 732. This multiplication involves multiplying a pair of complex numbers for every cycle over the entire symbol being evaluated to generate equalized cycle-rate data 1402. All of the equalized cycle values for an entire symbol are averaged together to generate equalized symbol-rate data, which is then normalized by the ratio of the root-mean-square of the sync symbols over a statistically large sample size divided by the root-mean-square of the BB-filtered values over a statistically large sample size, to obtain normalized and equalized symbol-rate data 1412. Thedata 1412 when combined with thePRN Symbols 1414 using the a combination of various functions such as multipliers and adders that are applied separately to the real and imaginary parts of each data set, as shown in the dashedregion 1408, produce PRN-removed symbol-rate data 1410. - In reference to
FIG. 19 , the process data function 1900 takes the PRN-removed symbol-rate data 1410 and reduces it to user data. First, thedata 1410 is separated into real and imaginary components. Values for components for all of the values that constitute an encoded bit are averaged together to form encodeddata 1902. The data is further processed using a Viterbi decoder to generate decodeddata 1904. And finally, user data 1906 is recovered after a de-interleaving process. -
FIG. 16 shows the align sync cross-correlated data for one embodiment, referenced earlier as 736 inFIG. 7 . This data in its complex form was generated by the align sync process, through which the target sync pattern was cross-correlated with bb-filteredcycle data 728. Thedata 736 is shown here with the phase information removed, as only absolute values. Where the peaks occur in thisdata 736, inFIG. 16 , is where bb-filteredcycle data 728 for each channel is maximally aligned with the sync pattern. Since all channel waveforms have the same reference (the same sync pattern), the time delay between these waveforms can be used through a geometric triangulation process, and by knowing the relative positions of the transducer array elements in three dimensions, and by knowing the speed of sound in the medium, a determination can be made as to the direction of the source of this transmission. However, using this align synccross-correlated data 736 in its raw form can be subject to unnecessary errors attributed to the asymmetries and irregularities in these waveforms. Where such undesired artifacts can be attributed to the time-domain impulse response of the entire system and possibly other sources, most of which are in common to all channels, much of these common effects can be removed by performing additional cross-correlations among these four waveforms to arrive at the well-behaved waveforms shown inFIG. 17 . Once this second round of cross-correlations are performed, the data is normalized and truncated below the 50% threshold to remove baseline noise and inter-symbol effects. The results are very clean cross-correlation waveforms, as shown inFIG. 17 , that still have the relative time delays between the channels preserved. At this point, the relative delays can be arrived at by measuring the relative peak locations. Additional accuracy can be obtained by employing second degree polynomial fitting to accurately locate the peaks. Another method can be used by measuring the relative locations of the centroids of each waveform. Regardless of the specific means to determine the relative delays between channels, that relative delay data is applied to a geometric model of the transducer array to determine the direction to the source diver responsible for the transmission being received. The distance to the source diver is determined by measuring the round-trip time-of-flight of transmissions passed between the two divers and by knowing the propagation speed of sound in the medium. -
FIG. 18 shows an embodiment of thetransducer array 1800 with fourindividual transducers point source 1810 that represents a signal originating from a transmitter located on another diver or a surface vessel (a responder). As shown bydistance lines responder point source 1810 and each of the transducers (on the initiator) is different. Using the arrival times of the signal generated by thepoint source 1810 at the transducers (as described above), and by measuring the round-trip time-of-flight between the two divers (from initiator to responder, and from responder back to initiator), the direction and distance between thetransducer array 1800 and thepoint source 1810 can be calculated. - Embodiments of a communication system described herein may be incorporated into one of many different types of underwater dive systems that may be utilized by recreational, commercial, industrial and military industries. Aspects of a communication system constructed in accordance with the present invention may be incorporated into a handheld dive computer underwater dive platform, or directly onto the BCD (buoyancy control device) of a diver.
- Visual display means may also be incorporated into such a system to enable a diver to see and send messages to one or more other divers while utilizing aspects of the communication system described herein.
- Aspects of the communication systems described herein may be combined with other known dive computer functions such as depth, dive time, tank pressure, decompression time as well as other known aspects of diving operations.
- Those skilled in the art can readily recognize that numerous variations and substitutions may be made in the invention, its use, and its configuration to achieve substantially the same results as achieved by the embodiments described herein. Accordingly, there is no intention to limit the invention to the disclosed exemplary forms. Many variations, modifications, and alternative constructions fall within the scope and spirit of the disclosed invention.
Claims (11)
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/483,631 US9444556B1 (en) | 2012-04-12 | 2014-09-11 | Underwater acoustic array, communication and location system |
US14/682,731 US9503202B2 (en) | 2012-04-12 | 2015-04-09 | Underwater acoustic array, communication and location system |
US15/166,812 US20160315717A1 (en) | 2012-04-12 | 2016-05-27 | Underwater acoustic array, communication and location system |
US15/289,482 US20170026135A1 (en) | 2012-04-12 | 2016-10-10 | Underwater acoustic array, communication and location system |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/445,108 US8842498B2 (en) | 2012-04-12 | 2012-04-12 | Underwater acoustic array, communication and location system |
US14/483,631 US9444556B1 (en) | 2012-04-12 | 2014-09-11 | Underwater acoustic array, communication and location system |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/445,108 Continuation US8842498B2 (en) | 2012-04-12 | 2012-04-12 | Underwater acoustic array, communication and location system |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/682,731 Continuation-In-Part US9503202B2 (en) | 2012-04-12 | 2015-04-09 | Underwater acoustic array, communication and location system |
Publications (2)
Publication Number | Publication Date |
---|---|
US9444556B1 US9444556B1 (en) | 2016-09-13 |
US20160280344A1 true US20160280344A1 (en) | 2016-09-29 |
Family
ID=49324960
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/445,108 Active 2033-03-28 US8842498B2 (en) | 2012-04-12 | 2012-04-12 | Underwater acoustic array, communication and location system |
US14/483,631 Active 2032-08-12 US9444556B1 (en) | 2012-04-12 | 2014-09-11 | Underwater acoustic array, communication and location system |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/445,108 Active 2033-03-28 US8842498B2 (en) | 2012-04-12 | 2012-04-12 | Underwater acoustic array, communication and location system |
Country Status (3)
Country | Link |
---|---|
US (2) | US8842498B2 (en) |
EP (1) | EP2837208B1 (en) |
WO (1) | WO2013154862A1 (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110244288A (en) * | 2019-07-01 | 2019-09-17 | 中国人民解放军军事科学院国防科技创新研究院 | A kind of sonar array signal processing method based on focusing playback principle |
WO2021067919A1 (en) * | 2019-10-04 | 2021-04-08 | Woods Hole Oceanographic Institution | Doppler shift navigation system and method of using same |
WO2022060720A1 (en) * | 2020-09-16 | 2022-03-24 | Woods Hole Oceanographic Institution | Single-receiver doppler-based sound source localization to track underwater target |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9154248B2 (en) * | 2011-02-14 | 2015-10-06 | Sirius Xm Radio Inc. | Method and apparatus for enhanced playback of content while switching among channels of broadcast or streamed content while being received |
US9503202B2 (en) * | 2012-04-12 | 2016-11-22 | Ceebus Technologies, Llc | Underwater acoustic array, communication and location system |
US8842498B2 (en) * | 2012-04-12 | 2014-09-23 | Ceebus Technologies Llc | Underwater acoustic array, communication and location system |
MX2014015006A (en) * | 2012-06-07 | 2015-05-11 | California Inst Of Techn | Communication in pipes using acoustic modems that provide minimal obstruction to fluid flow. |
GB2523131A (en) * | 2014-02-13 | 2015-08-19 | Amir Emad Fakhry Gerges | Managing diving activities by wireless data communication |
US9995824B2 (en) * | 2014-04-09 | 2018-06-12 | Thomas Danaher Harvey | Methods and system to assist search for lost and submerged objects |
US10871566B2 (en) * | 2014-04-09 | 2020-12-22 | Thomas Danaher Harvey | Methods and system to assist search and interception of lost objects |
US9794737B2 (en) * | 2014-09-08 | 2017-10-17 | The United States Of America, As Represented By The Secretary Of The Army | Underwater signal conversion |
US9667349B1 (en) | 2015-04-15 | 2017-05-30 | The United States Of America As Represented By The Secretary Of The Navy | Dynamic range extension of heterodyne fiber-optic interferometers via instantaneous carrier measurement |
CN104836730B (en) * | 2015-05-08 | 2018-07-24 | 东南大学 | It is a kind of to be applied to clustering route protocol method of the three-dimensional underwater acoustic sensor based on layering |
RU2623423C2 (en) * | 2015-07-14 | 2017-06-26 | Дмитрий Павлович Фролов | Personal underwater resque and navigation system |
US9520910B1 (en) * | 2015-09-24 | 2016-12-13 | Nxp B.V. | Receiver component and method for enhancing a detection range of a time-tracking process in a receiver |
CN106741658B (en) * | 2016-11-09 | 2018-03-13 | 哈尔滨工程大学 | A kind of acoustic marker automatic distributing and discharging structure |
WO2018101864A1 (en) * | 2016-11-30 | 2018-06-07 | Saab Ab | Sonar device with holder |
WO2018101865A1 (en) * | 2016-11-30 | 2018-06-07 | Saab Ab | Sonar with damping structure |
US11990948B2 (en) * | 2018-08-09 | 2024-05-21 | Onesubsea Ip Uk Limited | Subsea telecommunication system |
US10611445B1 (en) | 2018-09-19 | 2020-04-07 | Garmin Switzerland Gmbh | Wearable electronic device for detecting diver respiration |
US10358199B1 (en) | 2018-09-19 | 2019-07-23 | Garmin Switzerland Gmbh | Scuba tank air pressure monitor system |
US11472524B2 (en) | 2020-03-30 | 2022-10-18 | Raytheon Company | Integrated acoustic localization and communications array |
CN112147578B (en) * | 2020-08-21 | 2024-03-19 | 中国船舶重工集团公司第七一五研究所 | High-precision deep water transmitting array and multi-element vertical receiving array element positioning system and method |
CN113911298B (en) * | 2021-10-14 | 2022-07-22 | 厦门理工学院 | Intelligent device for assisting underwater search positioning |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3944967A (en) * | 1974-12-23 | 1976-03-16 | The United States Of America As Represented By The Secretary Of The Navy | Updated diver navigation device |
US4787252A (en) | 1987-09-30 | 1988-11-29 | Panametrics, Inc. | Differential correlation analyzer |
US4951263A (en) | 1989-07-20 | 1990-08-21 | Sandia Research Associates, Inc. | Spread spectrum underwater location beacon system |
IL121561A (en) | 1997-08-18 | 2000-10-31 | Divecom Ltd | Underwater communication apparatus and communication network |
US6130859A (en) | 1997-12-01 | 2000-10-10 | Divecom Ltd. | Method and apparatus for carrying out high data rate and voice underwater communication |
US6272072B1 (en) | 1998-08-14 | 2001-08-07 | Wulich Wave Ltd. | Underwater communication method, device, and system utilizing a doppler frequency shift |
US6674687B2 (en) | 2002-01-25 | 2004-01-06 | Navcom Technology, Inc. | System and method for navigation using two-way ultrasonic positioning |
US7146014B2 (en) | 2002-06-11 | 2006-12-05 | Intel Corporation | MEMS directional sensor system |
EP1615203A1 (en) | 2004-07-07 | 2006-01-11 | Nederlandse Organisatie Voor Toegepast-Natuurwetenschappelijk Onderzoek Tno | Ultrasonic transducer system |
US7388512B1 (en) | 2004-09-03 | 2008-06-17 | Daniel F. Moorer, Jr. | Diver locating method and apparatus |
IL164619A0 (en) | 2004-10-14 | 2005-12-18 | Arrangement for charging a battery of an underwater portable apparatus, which prevents discharge of the battery when the apparatus is underwater | |
EP1932145A2 (en) | 2005-08-29 | 2008-06-18 | VAN DER MERWE, Jacobus Johannes | Method of amplitude modulating a message signal in the audible frequency range onto a carrier signal in the ultrasonic frequency range |
US8654610B2 (en) | 2009-12-16 | 2014-02-18 | Shb Instruments, Inc. | Underwater acoustic navigation systems and methods |
US9503202B2 (en) * | 2012-04-12 | 2016-11-22 | Ceebus Technologies, Llc | Underwater acoustic array, communication and location system |
US8842498B2 (en) * | 2012-04-12 | 2014-09-23 | Ceebus Technologies Llc | Underwater acoustic array, communication and location system |
FI125009B (en) | 2013-09-10 | 2015-04-30 | Suunto Oy | Underwater communication systems and associated communication methods and devices |
-
2012
- 2012-04-12 US US13/445,108 patent/US8842498B2/en active Active
-
2013
- 2013-04-02 EP EP13775859.5A patent/EP2837208B1/en active Active
- 2013-04-02 WO PCT/US2013/034898 patent/WO2013154862A1/en active Application Filing
-
2014
- 2014-09-11 US US14/483,631 patent/US9444556B1/en active Active
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110244288A (en) * | 2019-07-01 | 2019-09-17 | 中国人民解放军军事科学院国防科技创新研究院 | A kind of sonar array signal processing method based on focusing playback principle |
WO2021067919A1 (en) * | 2019-10-04 | 2021-04-08 | Woods Hole Oceanographic Institution | Doppler shift navigation system and method of using same |
WO2022060720A1 (en) * | 2020-09-16 | 2022-03-24 | Woods Hole Oceanographic Institution | Single-receiver doppler-based sound source localization to track underwater target |
Also Published As
Publication number | Publication date |
---|---|
EP2837208A1 (en) | 2015-02-18 |
WO2013154862A1 (en) | 2013-10-17 |
US9444556B1 (en) | 2016-09-13 |
EP2837208B1 (en) | 2019-05-08 |
EP2837208A4 (en) | 2016-08-03 |
US20130272098A1 (en) | 2013-10-17 |
US8842498B2 (en) | 2014-09-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9444556B1 (en) | Underwater acoustic array, communication and location system | |
US9503202B2 (en) | Underwater acoustic array, communication and location system | |
US20170026135A1 (en) | Underwater acoustic array, communication and location system | |
US8139442B2 (en) | Robust underwater communication system | |
EP1378079B1 (en) | Communication system for underwater use | |
US6980487B2 (en) | Underwater communication system | |
US20050249036A1 (en) | Underwater location apparatus | |
CN103823205A (en) | Underwater locating navigation system and method | |
CN114679199B (en) | Underwater sound positioning and communication integrated signal design method | |
Baggeroer et al. | DATS-a digital acoustic telemetry system for underwater communications | |
WO2009027639A1 (en) | Underwater communications | |
Muller et al. | Underwater GPS System for Autonomous Underwater Wireless Drone Control [J] | |
Pennec et al. | On joint acoustic communication and positioning through MFSK-modulated signals and Costas arrays | |
Rodionov et al. | OFDM-Based Underwater Acoustic Communication System Designing for Under-Ice and Cold-Water Applications | |
Shustov et al. | COMPACT POSITIONING, DATE TRANSMISSION AND VOICE COMMUNICATION SYSTEM FOR UNDERWATER APPLICATIONS | |
Kurochkin et al. | Development of a Hydroacoustic Modem as a Research Project of the Department of Marine Electronics (SMTU) | |
Chase | A Tracking and Telemetry System for Severe Multipath Acoustic Channels | |
JP2021108432A (en) | Underwater communication system | |
Pajovic | TURBO-CODED FREQUENCY DIVISION |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: CEEBUS TECHNOLOGIES, LLC, COLORADO Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:CAHALAN, DENNIS;NOONAN, TED;HENDRICKS, MICHEL;AND OTHERS;REEL/FRAME:033723/0217 Effective date: 20120717 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
FEPP | Fee payment procedure |
Free format text: MAINTENANCE FEE REMINDER MAILED (ORIGINAL EVENT CODE: REM.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
FEPP | Fee payment procedure |
Free format text: SURCHARGE FOR LATE PAYMENT, SMALL ENTITY (ORIGINAL EVENT CODE: M2554); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YR, SMALL ENTITY (ORIGINAL EVENT CODE: M2551); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY Year of fee payment: 4 |
|
FEPP | Fee payment procedure |
Free format text: MAINTENANCE FEE REMINDER MAILED (ORIGINAL EVENT CODE: REM.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
FEPP | Fee payment procedure |
Free format text: 7.5 YR SURCHARGE - LATE PMT W/IN 6 MO, SMALL ENTITY (ORIGINAL EVENT CODE: M2555); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YR, SMALL ENTITY (ORIGINAL EVENT CODE: M2552); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY Year of fee payment: 8 |