NO336331B1 - Fremgangsmåte for prosessering av ekkoer fra utsendte signaler. - Google Patents

Fremgangsmåte for prosessering av ekkoer fra utsendte signaler.

Info

Publication number
NO336331B1
NO336331B1 NO20053337A NO20053337A NO336331B1 NO 336331 B1 NO336331 B1 NO 336331B1 NO 20053337 A NO20053337 A NO 20053337A NO 20053337 A NO20053337 A NO 20053337A NO 336331 B1 NO336331 B1 NO 336331B1
Authority
NO
Norway
Prior art keywords
spectrum
frequency
doppler
signal
copy
Prior art date
Application number
NO20053337A
Other languages
English (en)
Other versions
NO20053337D0 (no
NO20053337L (no
Inventor
Yves Doisy
Laurent Deruaz
Pierre Alinat
Marc Masoni
Original Assignee
Thales Sa
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Thales Sa filed Critical Thales Sa
Publication of NO20053337D0 publication Critical patent/NO20053337D0/no
Publication of NO20053337L publication Critical patent/NO20053337L/no
Publication of NO336331B1 publication Critical patent/NO336331B1/no

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/02Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
    • G01S15/50Systems of measurement, based on relative movement of the target
    • G01S15/52Discriminating between fixed and moving objects or between objects moving at different speeds
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/02Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
    • G01S15/50Systems of measurement, based on relative movement of the target
    • G01S15/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S15/582Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse-modulated waves and based upon the Doppler effect resulting from movement of targets
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52001Auxiliary means for detecting or identifying sonar signals or the like, e.g. sonar jamming signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/523Details of pulse systems
    • G01S7/526Receivers
    • G01S7/527Extracting wanted echo signals
    • G01S7/5273Extracting wanted echo signals using digital techniques
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2626Arrangements specific to the transmitter only
    • H04L27/2627Modulators
    • H04L27/2628Inverse Fourier transform modulators, e.g. inverse fast Fourier transform [IFFT] or inverse discrete Fourier transform [IDFT] modulators

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Discrete Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Signal Processing (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Holo Graphy (AREA)
  • Dc Digital Transmission (AREA)

Description

Oppfinnelsens bakgrunn
Oppfinnelsens område er fremgangsmåter for prosessering av sonarsignaler. Oppfinnelsen vedrører mer spesifikt sonar ekkoer som kommer fra objekter som er bestrålt av sonarsendinger med et "kam"-linjespektrum sammensatt av spektrallinjer.
Oppsummering av oppfinnelsen
Én fremgangsmåte for utnyttelse av sonarsignaler, kjent fra tidligere teknikk, omfatter sending av et kjent signal inn i mediet og prosessering av de mottatte signalene ved hjelp av korrelasjon med signaler som representerer kopier av det sendte signalet når det reflekteres fra et bevegelig objekt og mottas av sonaren. Disse kopiene, eller dopplerkopiene, korreleres med signalet som faktisk mottas for å skille signalet som kommer fra objektet fra bakgrunnsstøyen og fra etterklangen.
Utsending av signaler som har et linjespektrum, med andre ord et spektrum sammensatt av distinkte linjer, er videre av stor interesse for detektering av hurtiggå-ende eller saktegående undervannsobjekter i et miljø som er begrenset av etter-klang. Riktignok, så lenge dopplerforskyvningen av objektene er større enn bredden av de elementære linjene i spektret, kan etterklangsstøyen reduseres ved hjelp av dopplerfiltrering. En slik dopplerfiltreringsprosess er beskrevet i europeisk patent EP1078280 som vedrører en fremgangsmåte for detektering av bevegelige objekter ved hjelp av en aktiv sonar.
Denne type dopplerfiltrering omfatter gjennomføring av en flerkopikorrelasjon av det mottatte signalet. Flerkopikorrelasjon omfatter gjennomføring av, på en kjent måte, korrelasjonsoperasjoner på det mottatte signalet med ulike kopier, eller dopplerkopier, av det sendte signalet. Hver dopplerkopi viser det samme spektret som signalet reflektert fra et objekt med en gitt hastighet, hvilket forklarer hvorfor de refereres til som dopplerkopier av det sendte signalet.
Korrelasjonsoperasjonene gjennomføres samtidig, eller i samme tidsperiode, slik at utnyttelsen av, i aktiv sonarmodus, utsendinger av linjespektrumsignaler medfører at det utføres like mange samtidige korrelasjonsoperasjoner som det er brukt dopplerfiltre. Følgelig, når antall brukte dopplerkopier er stort, kan prosesseringen av flerkopikorrelasjon kreve stor behandlingskapasitet.
Videre er det kjent at, for signaler som har et stort BT-produkt, for eksempel større enn 50, kan dopplerfiltreringsoperasjonen ved korrelasjon med fordel utføres i spektraldomenet. Bokstavene B og T i BT-produktet benevner henholdsvis det sendte signalets bånd B og varighet T. Faktisk viser det seg at overføring til spektraldomenet, for eksempel ved hjelp av en FFT (Fast Fourier Transform), er mer fordelaktig med hensyn til beregningsbelastning, enn en direkte korrelasjon i tidsdomenet.
Likevel, når korrelasjonen utføres med flere dopplerkopier, medfører denne utførel-sesformen for dopplerfiltrering en beregningsbelastning som fremdeles er altfor stor. Med dette i minne erindres det at økningen i beregningsbelastning påført av flerkopiprosesseringen er stort sett proporsjonal med antall brukte kopier. I praksis er antall kopier i størrelsesorden flere titalls, eller til og med flere hundre. Fransk patent FR 2853417 nevner i dette henseende et antall dopplerkopier større enn 200 nødvendig for dopplerfiltreringen.
US 6314053 viser en metode for Dopplerprosessering av ekko av sendte signaler, som har et linjespektrum, med en flerkopikorrelasjon av det mottatte signalet.
US 5557637 viser prosessering av signaler som er påvirket av Dopplereffekten fra bevegelige kjøretøy eller satellitter, hvor det benyttes Fast Fourier Transform, fol-ding og reindeksering for å redusere beregningsbehovet av de prosesserte ra-dar/sonar-signalene.
Andre metoder for prosessering av sonar-signaler er kjent fra:
CHAN F. et al., «A non-linear phase-only algorithm for active sonar signal processing», OCEANS '97, MTS/IEEE Conference Proceedings, Halifax, NS, Canada 6-9 Oet. 1997, side 506-511.
HERMAND J-P, «Model-based matched filter processing for delay-Doppler meas-urement in a multipath dispersive ocean channel», OCEANS '93, Engineering in Harmony with Ocean, Proceedings, Victoria, BC, Canada 18-21 Oet. 1993, side 306-311.
MIO K. et al., "Space time adaptive processing for low frequency sonar, OCEANS 2000, MTS/IEEE conference and Exhibition, Sept. 11-14, 2000, side 1315-1319.
Den foreliggende oppfinnelsen består av en fremgangsmåte som tillater at beregningsbelastningen forårsaket av flerkopikorrelasjonen reduseres. Med denne hen sikt utfører fremgangsmåten i henhold til oppfinnelsen en flerkopikorrelasjonsope-rasjon som omfatter minst følgende trinn: - et trinn for tid-frekvens-transformasjonen av det mottatte signalet; - et trinn for dannelsen av syntetiske spektra, hvor hvert genererte spektrum dannes ved at det velges bestemte spektralkomponenter fra spektret til det mottatte signalet og at disse spektralkomponentene sammenkjedes («concatenation»); - et trinn for beregning av interspektralproduktet, som er produktet av hvert syntetisk spektrum og det sammenkjedede konjugerte spektret til den korresponderende dopplerkopien; - et trinn for utførelse av flerkopikorrelasjonsoperasjonen hvor produktet av det sammenkjedede konjugerte spektret til hver dopplerkopi og det syntetiske spektret til det korresponderende mottatte signalet beregnes; - et trinn for gjennomføring av frekvens-tid-transformasjonen av de fremskaffede interspektralproduktene, hvilket gjør det mulig for systemet å returnere til tidsdomenet.
Fremgangsmåten i henhold til oppfinnelsen har den fordelen at flerkopikorrelasjonen utføres kun på de nyttige delene av spektret. Ved å eliminere unødvendige beregninger, medfører den dermed at antall beregninger som utføres reduseres betraktelig, spesielt under trinnene for interspektralberegning og frekvens-tid-transformasjon. Den medfører dermed at totalt antall beregninger som gjennomfø-res reduseres eller, uttrykt på en annen måte, at det gjennomføres, for samme antall beregninger, en dopplerfiltreringsoperasjon med høyere oppløsning.
Fremgangsmåten i henhold til oppfinnelsen har den ekstra fordelen at den er enkel å implementere. Den er også spesielt tilpasset systemer som sender ut bølger med spektrum som er sammensatt av periodiske linjer, hvor hver linje selv har et ele-mentært spektrum.
Kort beskrivelse av tegningene
Andre trekk og fordeler vil klargjøres ved hjelp av beskrivelsen illustrert i de ved-lagte figurene som viser:
- figur 1, en tidsavhengig fremstilling av et eksempel på et signal som har et linjespektrum; - figur 2, den skjematiske spektralfremstillingen av signalet illustrert i figur 1; - figur 3, en fremstilling av hovedtrinnene i en konvensjonell fremgangsmåte for flerkopikorrelasjon, hvor korrelasjonen gjennomføres i spektraldomenet; - figur 4 viser også et begrenset snitt av spektret i figur 2 og av gruppen med dopplerkanaler definert av frekvensanalyseoperasjonen gjennomført på signalet; - figur 5, den sammenkjedede spektralrepresentasjonen til signalet illustrert i figur i; - figur 6, flytdiagrammet for hovedtrinnene i fremgangsmåten i henhold til oppfinnelsen.
Detaljert beskrivelse av oppfinnelsen
Som nevnt tidligere kan fremgangsmåten i henhold til oppfinnelsen anvendes spesielt på aktive sonarsystemer som sender ut bølger med et spektrum formet som en vanlig rekke, langs frekvensaksen, med elementære spektra eller linjer. For å tydeliggjøre beskrivelsen, er oppfinnelsen beskrevet i dokumentet ved hjelp av en foretrukket eksemplifiserende utførelsesform hvor den utsendte bølgen består av pulser dannet av den veide summen av N suksessive elementære pulser modulert i frekvens i henhold til den samme loven, for eksempel lineær eller hyperbolsk. En illustrasjon i tidsdomenet til en puls av denne typen vises i figur 1. I denne figuren kan de N elementære pulsene 11 som har amplitude bestemt av anvendt veielov sees. Denne veiingen definerer en modulasjonskurve a(t) 12 og bestemmer formen på spektret til pulsen dermed generert. Det dannede signalet, med varighet T, er dermed sammensatt av N pulser med varighet T/N. Hensikten med veiingen anvendt på hele signalet er å kontrollere spektralbredden til det utsendte signalet og, mer spesifikt, å redusere nivået på de sekundære spektrallinjene. Følgelig kan det sendte signalet uttrykkes ved følgende ligning: hvor p(t) er uttrykket for en elementær puls, og p(t) er lik 0 når t er utenfor intervallet [0 T/N]. Signalet a(t) representerer anvendt veiing, som refereres til som modulasjonskurven til signalet. Når det gjelder spektret til det sendte signalet, kan dette uttrykkes ved følgende ligning:
I ligning [2] representerer P(f) spektret til den elementære pulsen og kan uttrykkes ved:
Tilsvarende kan spektret til modulasjonskurven a(t) til det sendte signalet uttrykkes ved:
Spektret E(f) vises skjematisk i figur 2. Dens totale bredde er båndet B til den elementære pulsen. Den er sammensatt av spektrallinjer 21 posisjonert hver n<*>N/T langs frekvensaksen.
Amplitudeveiingen a(t) har den effekten at den begrenser bredden til hver linjes spektrum til et intervall p/T rundt midtfrekvensen. Faktoren p representerer en ko-effisient som er avhengig av valgt lov for amplitudeveiing a(t). Dersom, for eksempel, det bestemmes at, i henhold til en foretrukket utførelsesform, a(t) = sin2(7tt/T), vil koeffisienten p få verdien p = 4.
Fremstillingen i figur 2 demonstrerer at spektret til det sendte signalet omfatter BT/N frekvensområder med bredde p/T, tilsvarende spektrallinjene 21, hvor ampli-tuden til spektret er signifikant, atskilt av frekvensområder 22 med bredde stort sett lik N/T, og hvor amplitudene til spektralkomponentene er svært lave. Det kan også observeres at disse områdene med svært lave amplituder representerer en signifikant del av båndet til det sendte signalet. Denne markante distribusjonen av spektralamplituden er fordelaktig benyttet av fremgangsmåten i henhold til oppfinnelsen.
Generelt sett omfatter en kjent fremgangsmåte for gjennomføring av doppleranalysen av det mottatte signalet i området sonarekkoprosessering gjennomføring av korrelasjonen av det mottatte signalet med signaler som representerer kopier av det sendte signalet men med en frekvensforskyvning eller fasedreiing tilsvarende dreiingen som det sendte signalet ville ha gjennomgått under dets refleksjon fra et bevegelig objekt. I denne type prosessering med flerkopikorrelasjon, er antall an-vendte dopplerkopier avhengig av frekvensområdet som skal analyseres og av ko-deoppløsningen. Dette antallet er også avhengig av beregningskapasiteten til ma-skinen som må gjennomføre de korresponderende aritmetiske operasjonene. Denne tvillingavhengigheten medfører, i praksis, et vanskelig kompromiss mellom den ønskede kvaliteten på doppleranalysen og beregningstiden nødvendig for å oppnå resultatet.
Figur 3 viser det skjematiske flytdiagrammet for en konvensjonell fremgangsmåte for prosessering ved hjelp av flerkopikorrelasjon. Med hensikt å begrense antall beregninger som må gjennomføres, gjennomføres denne operasjonen i frekvensdomenet. Av denne grunn gjennomgår det digitaliserte mottatte signalet 31, for eksempel, en tid-frekvens-transformasjon 32. I en foretrukket utførelsesform, gjennomføres tid-frekvens-transformasjonen ved hjelp av en FFT (Fast Fourier Transform) utført i en skyvevindumodus ("sliding window") over et antall punkter NI bestemt av ønsket antall dopplerkanaler. FFT-operasjonen kan for eksempel gjennomføres i en skyvevindumodus over et antall punkter NI = 2<*>Fe<*>T, med et vindu lik Fe<*>T punkter fra én FFT-operasjon til den neste. I dette tilfellet, representerer Fe punktfrekvensen og T varigheten av det mottatte signalet.
Frekvenssignalet 33 som dermed oppnås utsettes deretter for flere samtidige korrelasjonsoperasjoner 34. Hver operasjon omfatter beregning av, på en kjent måte, interspektralresultanten av produktet til spektret til det mottatte signalet og spektret til et lagret signal 35, som kalles dopplerkopi. Spektret til hver dopplerkopi korresponderer med det til det utsendte signalet som har gjennomgått en homotetisk frekvenstransformasjon. Den homotetiske frekvenstransformasjonen korresponderer til én av kanalene i det analyserte dopplerbåndet.
Hver dopplerkopi kan karakteriseres av dopplerparameteren 8 definert ved uttrykket 8 = 1 + 2V/c hvor V representerer radialhastigheten til objektet i forhold til far-tøyet som bærer sonaren og c lydhastigheten i mediet. Den kan da uttrykkes som: og dens spektrum følger uttrykket: Utrykket for C(f) er analogt til utrykket for E(f), med unntak av at det er sammensatt av linjer sentrert på frekvensene 5<*>n<*>N/T som representerer en enkel nomote-tisk transformasjon av spektret.
Resultatet 36 av hver korrelasjonsoperasjon transponeres deretter til tidsdomenet for å kunne anvendes. Transponeringen gjennomføres fortrinnsvis av inverstransformasjonen 37 som brukes for overføring til frekvensdomenet, av en frekvens-tid-transformasjon av type invers DFT eller invers FFT eller annet, gjennomført over Ni utvalg eller punkter.
Når det gjelder det sendte signalet, observeres det at spektrene til dopplerkopiene er ulik null over BT/N frekvensområder som korresponderer med linjene til spektret til det sendte signalet oversatt til frekvens. Når det gjelder det sendte signalet, er bredden til hver spektrallinje lik p/T, og linjene atskilles av 5N/T på frekvensaksen.
Utførelse av flerkopikorrelasjonsberegningen er ensbetydende med utførelse av Ndsamtidige korrelasjonsoperasjoner, med andre ord, i spektraldomenet,
Ndinterspektral-multipliseringsoperasjoner. Overføring til spektraldomenet tillater fordelaktig å erstatte tidskorrelasjonoperasjonene, som involverer på en kjent måte summene av produktene, med enkle multiplikasjonsoperasjoner som er mindre komplekse og dermed mindre kostbare når det gjelder beregningsbelastning. Likevel, da det involverer Ndsamtidige operasjoner, forblir belastningen høy. Det min-nes om at, i denne sammenheng, størrelsen på beregningsbelastningen er avhengig av antall Nddopplerkanaler som skal prosesseres og også av produktet BT, hvor B representerer båndbredden til den elementære pulsen 11 og T varigheten til de N elementære pulsene som danner det sendte signalet.
Flytdiagrammet i figur 3 gjør det mulig å håndtere den relative veiingen i den totale beregningsbelastningen for hvert trinn i fremgangsmåten som fremheves. Det vil dermed observeres at trinnene som korresponderer med operasjonene 34 og 37 har en veiing som er svært følsom i beregningsbelastningen, både når det gjelder antall beregninger nødvendig for hver operasjon og når det gjelder samtidigheten av deres utførelse for hver av Nddopplerkanalene. Til sammenligning har trinnet som korresponderer med operasjonen 32 for FFT frekvensanalysen av det mottatte signalet kun en lav relativ vekt.
I tilfellet hvor den akseptable beregningsbelastningen ikke gjør det mulig å gjen-nomføre flerkopikorrelasjonsoperasjonen som den beskrevet, vil den konvensjonelle fremgangsmåten for prosessering illustrert i figur 3 ved første blikk bare tilby to løsningstyper. En første løsning omfatter derfor minsking av antall dopplerkanaler som brukes samtidig, hvilket gjør det mulig å redusere antall Ndkorrelasjonsoperasjoner som gjennomføres samtidig. Den andre løsningen omfatter begrensing av båndbredden okkupert av den elementære pulsen noe som kan være uforenelig med kravene forbundet med signalprosesseringsfunksjonene som anvendes andre steder.
Konfrontert med disse restriksjonene, foreslår fremgangsmåten i henhold til oppfinnelsen en løsning som er et alternativ til de tidligere nevnte løsningene. Med denne hensikt, består prinsippet til fremgangsmåten i henhold til oppfinnelsen hovedsakelig av å betrakte spektret til det sendte signalet som en samling av elementære spektra formet av spektrallinjene 21 som danner det globale spektret. Disse elementære spektrene atskilles av frekvensområder som er overflødige i dopplerpro-sesseringen og hvor signalnivået er svært lavt. Tilsvarende, da det mottatte signalet korreleres mot dopplerforskjøvede kopier av det sendte signalet, omfatter det oppnådde interspektret også frekvensområder hvor frekvenskomponentene er vel-dig svake. Denne betraktningen fører til introduksjonen av begrepet effektiv båndbredde, hvilket kan defineres av følgende ligning:
I ligningen ovenfor, defineres det effektive båndet Befrsom produktet BT/N av antall linjer som danner signalspektret og spektralbredden pB/N til en linje.
Begrepet effektiv båndbredde, som definert her, tillater det faktumet at alle punktene som er nyttige for gjennomføring av en doppleranalyse av det mottatte signalet, ved hjelp av digital korrelasjon i spektraldomenet, ikke er i B, men i Befrfor at de skal tas med i beregningen.
Fremgangsmåten i henhold til oppfinnelsen drar fordel av denne observasjonen for å redusere beregningskraften nødvendig for utførelse av flerkopikorrelasjonsoperasjonen. På grunn av dette omfatter fremgangsmåten i henhold til oppfinnelsen et trinn som omfatter en operasjon for sammenkjeding av spektra. Hensikten med denne operasjonen er å erstatte det mottatte signalet og hver dopplerkopi med et syntetisk signal hvilket spektrum er innenfor et bånd tilsvarende det tidligere defi-nerte båndet Beff. For hver kopi, omfatter det syntetiske spektret alle frekvenskomponentene med nivå ulik null til det originale signalet. Illustrasjonen i figur 4 bidrar til å tydeliggjøre, på en enkel måte, funksjonen til trinnet for sammenkjeding av spektrene til dopplerkopiene.
Figur 4 viser et forstørret begrenset utsnitt av det virkelige spektret til signalet som vises i figur 2. I denne figuren, vises kun tre spektrallinjer 41, som kan være hvilke som helst av tre linjer tatt fra linjene som danner spektret til et signal som, for eksempel, signalet sendt fra sonaren eller én av dets dopplerkopier. Figur 4 viser også et begrenset snitt av hele gruppen med frekvenskanaler 42 definert av frekvensanalyseoperasjonen gjennomført på signalet. Denne frekvensanalysen kan, for eksempel, gjennomføres av en FFT-type operasjon.
Sammenkjedingstrinnet gjennomføres samtidig for hver dopplerkopi, hvor hver har et forskjellig linjespektrum. Det kan deles opp i et bestemt antall operasjoner. En første operasjon omfatter bestemmelse av frekvenskanaler 42 som sammenfaller med spektrallinjene 41, med kjent posisjon, som danner spektret til kopien. Spektralbåndet som okkuperes av hver av linjene dekkes dermed av en gruppe 43 med nabofrekvenskanaler hvor antallet avhenger særlig av oppløsningen til frekvensanalysen som gjennomføres initielt. Denne gruppen med dopplerkanaler 43 omfatter det minimum antall kanaler som er nødvendig og deres spektralbredde er, i praksis, hovedsakelig lik den til den aktuelle linjen 41.
Generelt sett, som illustrert i figur 4, er spektralbåndet som er okkupert av en linje 41 i spektret til det analyserte signalet ikke nødvendigvis sentrert på midtfrekvensen i en frekvenskanal. Videre er spektralbredden til en linje større enn bredden til en dopplerkanal. Dette er grunnen til at det er nødvendig, for å dekke båndet som er okkupert av en spektrallinje, å betrakte flere nabofrekvenskanaler. Minimum antall kanaler som er nødvendig bestemmes av bredden til den aktuelle linjen og av størrelsen på dopplerkanalene.
Gruppene med nabo-dopplerkanaler 43 som korresponderer med hver av linjene bestemmes dermed, og selve sammenkjedingsoperasjonen omfatter eliminering, fra spektralrepresentasjonen til signalet, av komponentene som ligger i de mellom-liggende frekvenskanalene som ikke tilhører noen av de dannede gruppene, og om-danning av en spektralrepresentasjon av signalet ved å sidestille de gjenstående komponentene. Denne representasjonen kan betraktes å være den til det effektive spektret til den aktuelle kopien.
Som vist i figur 5, dannes denne spektralrepresentasjonen av en samling av linjer med bredde p/T, atskilt fra hverandre av et frekvensintervall b/T hovedsakelig lik p/T. Dette intervallet er, i virkeligheten, definert av antall frekvenskanaler bevart rundt hver spektrallinje. Hele spektret som dermed dannes har en spektralbredde hovedsakelig lik pB/N.
Sammenkjedingstrinnet slik som beskrevet gjør det dermed fordelaktig mulig å redusere spektret til signalet eller signalene som skal behandles. Fremgangsmåten i henhold til oppfinnelsen anvender sammenkjedingen av spektra slik som illustrert i figur 6 som viser et flytdiagram for trinnene omfattet av fremgangsmåten.
I henhold til oppfinnelsen, gjennomgår det mottatte signalet en spektraldekomposi-sjon 32, for eksempel av type FFT. Spektraldekomposisjonen av det mottatte signalet er fortrinnsvis identisk med den som ble brukt for å danne spektrene til kopiene av det mottatte signalet. Videre er de sammenkjedede spektrene 63 til frekvensko-piene av det sendte signalet tilgjengelig som lagrede data, hvor hvert sammenkjedet spektrum omfatter en gruppe med gitte frekvenskanaler. Disse kanalene identi-fiseres, for eksempel, av et tall som representerer deres sekvens i kanalsenene 42 som dekker båndet B til signalet.
Spektret 33 til det mottatte signalet 32 går deretter gjennom et trinn for å danne sammenkjedede spektra 62, hvor hvert sammenkjedet spektrum 62 kombinerer frekvenskomponentene som ligger i de samme frekvenskanalene som de som er tilbakeholdt for å danne det sammenkjedede spektret til dopplerkopien 63 som den må korreleres med.
Fremgangsmåten omfatter deretter et trinn for interspektralberegning 64, gjen-nomført på en konvensjonell måte, som omfatter, som før, beregning av, på en kjent måte, interspektret som er resultatet av hvert produkt av et syntetisk spektrum til det mottatte signalet og det tilhørende konjugerte spektret til den lagrede dopplerkopien 35 av det sendte signalet. Dermed tar fremgangsmåten i henhold til oppfinnelsen fordelaktig i bruk hvert produkt av spektra, ikke på de originale spektrene, men på de sammenkjedede spektrene. Da de sammenkjedede spektrene omfatter et antall N2spektralkomponenter, mindre enn antallet Ni fra de originale spektrene, medfører flerkopikorrelasjonsoperasjonen en beregningsbelastning som er mindre enn tilfellet ved en konvensjonell fremgangsmåte.
Det oppnådde interspektret 65 gjennomgår så et trinn for frekvens-tid-transformasjon 66, for eksempel av type invers FFT, for å gjenopprette tidssignale-ne som representerer korrelasjonsproduktene. Denne operasjonen gjennomføres samtidig på alle interspektraene beregnet fra de sammenkjedede spektrene til det mottatte signalet og fra dopplerkopiene. Denne inverstransformasjonen ligner, i prinsippet, på de som er implementert i kjente fremgangsmåter av tidligere teknikk. Men, det bør bemerkes at, når det gjelder transformasjon ved bruk av invers FFT, gjennomføres operasjonen på interspektra som omfatter færre komponenter, hvilket reduserer størrelsen på FFT-en som skal gjennomføres og dermed antall nødvendige beregninger. Disse FFT-operasjonene utføres på antall N2frekvens-punkter i et interspektrum eller, i praksis, på et antall 2P punkter som tilsvarer annen potens som er det nærmeste dette antallet. Antall tilgjengelige punkter gjen-nomføres, i det siste tilfellet, ved addisjon av punkter som er lik null. Dermed, som for trinnet for interspektralkorrelasjon, resulterer effekten av de sammenkjedede resultatene i en vesentlig reduksjon av antall operasjoner som må gjennomføres, hvor reduksjonens størrelsesorden er gitt av forholdet N/p. Spektralanalysen kan gjennomføres mer generelt med en DFT (Diskret Fourier Transform) som ikke nød-vendigvis anvendes på et antall punkter lik annen potens.
Fremgangsmåten i henhold til oppfinnelsen, som beskrevet ovenfor, medfører derfor at antall aritmetiske operasjoner som er nødvendige for å gjennomføre doppleranalysen av det mottatte signalet reduseres betraktelig, så lenge det sendte signalet har et linjespektrum og en veiet distribusjon som gjør at spektralbredden til hver linje kan begrenses. Reduksjonen anvendes fordelaktig i trinnene for korrelasjon og frekvens-tid-transformasjon som er de avgjørende trinnene i fremgangsmåten når det gjelder beregningsbelastning.
Denne fremgangsmåten har fordelen av å være spesielt tilpasset signalene som kommer fra utsendinger i form av et utbrudd med N etterfølgende pulser over tid, hvor hver puls frekvensmoduleres over et bånd B med en midtfrekvens f0.

Claims (5)

1. Fremgangsmåte for dopplerprosessering av ekkoer av sendte signaler, hvor disse sendte signalene har et linjespektrum, karakterisert vedat den gjennomfører en flerkopikorrelasjon av det mottatte signalet, og denne flerkopikorrelasjonen omfatter minst føl-gende trinn: - et trinn (32) for tid-frekvens-transformasjonen av det mottatte signalet, - et trinn (61) for dannelse av syntetiske spektra ved hjelp av sammenkjeding, med utgangspunkt i spektret til det mottatte signalet, - et trinn (64) for beregningen av interspektralproduktet, med andre ord produktet av hvert genererte spektrum med det sammenkjedede konjugerte spektret til den korresponderende dopplerkopien, og - et trinn (66) for frekvens-tid-transformasjonen av interspektralproduktene.
2. Fremgangsmåte ifølge krav 1, karakterisert vedat, for hver dopplerkopi, fremskaffes det sammenkjedede spektret til nevnte kopi ved å velge, for hver spektrallinje som danner signalet, frekvenskanalen som ligger nærmest posisjonen til linjen, ved også å velge nabofrekvenskanalene, ved å eliminere de andre frekvenskanalene og ved å sammenkjede frekvensområdene som korresponderer med de valgte frekvenskanalene.
3. Fremgangsmåte ifølge enten krav 1 eller krav 2,karakterisert vedat hvert syntetisk spektrum fremskaffes ved å velge og ved å sammenkjede frekvensområdene som korresponderer med de valgte frekvenskanalene for å danne det sammenkjedede spektret til den korresponderende dopplerkopien.
4. Fremgangsmåte ifølge hvilket som helst av kravene 1 til 3,karakterisert vedat frekvens-tid-transformasjonene til interspektralproduktene gjennomføres ved hjelp av en invers FFT utført på et antall N2punkter lik antall frekvenskanaler som danner dopplerkopiene.
5. Fremgangsmåten i krav 4, karakterisert vedat den inverse FFT-en gjennomføres på en gruppe M punkter dannet av de N2punktene som korresponderer med frekvenskanalene til det sammenkjedede spektret og hvor M - N2null-punkter legges til slik at M er lik kvadratet umiddelbart over N2.
NO20053337A 2004-07-09 2005-07-08 Fremgangsmåte for prosessering av ekkoer fra utsendte signaler. NO336331B1 (no)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
FR0407706A FR2872919B1 (fr) 2004-07-09 2004-07-09 Traitement coherent rapide pour codes a spectre de raies periodiques

Publications (3)

Publication Number Publication Date
NO20053337D0 NO20053337D0 (no) 2005-07-08
NO20053337L NO20053337L (no) 2006-01-10
NO336331B1 true NO336331B1 (no) 2015-08-03

Family

ID=34940279

Family Applications (1)

Application Number Title Priority Date Filing Date
NO20053337A NO336331B1 (no) 2004-07-09 2005-07-08 Fremgangsmåte for prosessering av ekkoer fra utsendte signaler.

Country Status (10)

Country Link
US (1) US7286443B2 (no)
EP (1) EP1615049B1 (no)
JP (1) JP4899031B2 (no)
KR (1) KR101175087B1 (no)
AT (1) ATE371875T1 (no)
DE (1) DE602005002197T2 (no)
ES (1) ES2292064T3 (no)
FR (1) FR2872919B1 (no)
NO (1) NO336331B1 (no)
SG (1) SG119299A1 (no)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2872919B1 (fr) * 2004-07-09 2006-09-29 Thales Sa Traitement coherent rapide pour codes a spectre de raies periodiques
FR2900504B1 (fr) * 2006-04-26 2009-11-20 Thales Sa Procede pour optimiser l'alimentation d'une antenne lineaire d'emission remorquee pour emettre en mode omnidirectionnel.
JP4890195B2 (ja) * 2006-10-24 2012-03-07 日本電信電話株式会社 ディジタル信号分波装置及びディジタル信号合波装置
DE102008040248A1 (de) * 2008-07-08 2010-01-14 Robert Bosch Gmbh Verfahren und Vorrichtung zum Bestimmen einer Geschwindigkeit eines Objekts
FR2951561B1 (fr) * 2009-10-20 2011-12-09 Areva T & D Sas Procede de detection de la position d'un front d'onde correspondant a un evenement dans un signal recu par un detecteur
JP5047314B2 (ja) * 2010-01-15 2012-10-10 富士フイルム株式会社 有機電界発光素子
DE102010044742A1 (de) * 2010-09-08 2012-03-08 Atlas Elektronik Gmbh Verfahren und Vorrichrung zur Bestimmung einer aus dem Doppler-Effekt resultierenden Doppler-Frequenzverschiebung
PT106909B (pt) * 2013-04-26 2017-06-01 Cuf - Quim Ind S A Processo e equipamento para remoção de impurezas de compostos aromáticos nitrados
DE102017107634A1 (de) * 2017-04-10 2018-10-11 Atlas Elektronik Gmbh Verarbeitungseinheit für eine sonaranlage zum verarbeiten von hydrophonsignalen sowie sonaranlage und verfahren
KR102057271B1 (ko) * 2018-07-24 2019-12-18 국방과학연구소 음향 신호 확장 장치 및 방법
CN113534115B (zh) * 2021-05-31 2023-09-29 中国船舶重工集团公司第七一五研究所 一种主被动联合处理的声纳目标高精度线谱提取方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2653564B1 (fr) 1989-10-20 1992-01-24 Thomson Csf Procede de formation de voies pour sonar.
US5245587A (en) * 1990-12-14 1993-09-14 Hutson William H Multi-dimensional signal processing and display
US5490516A (en) * 1990-12-14 1996-02-13 Hutson; William H. Method and system to enhance medical signals for real-time analysis and high-resolution display
FR2675587B1 (fr) 1991-04-16 1993-07-16 Thomson Csf Procede de compensation du mouvement de l'antenne pour un sonar.
US5557637A (en) * 1994-09-24 1996-09-17 Glynn; Thomas W. Convolutional ambiguity multiple access (CAMA) transmission system
FR2754907B1 (fr) 1996-10-22 1999-01-08 Thomson Marconi Sonar Sas Procede de mesure du decalage doppler dans un systeme de detection utilisant des codes ambigus
FR2778748B1 (fr) * 1998-05-15 2004-08-27 Thomson Marconi Sonar Sas Procede de detection d'objets mobiles par sonar actif
JP3199240B2 (ja) * 1999-03-16 2001-08-13 防衛庁技術研究本部長 ドップラー補正を行う航走雑音用音響測位装置及び方法
JP3367462B2 (ja) * 1999-04-20 2003-01-14 日本電気株式会社 アクティブソーナー及びその目標検出方法
FR2872919B1 (fr) * 2004-07-09 2006-09-29 Thales Sa Traitement coherent rapide pour codes a spectre de raies periodiques

Also Published As

Publication number Publication date
EP1615049A1 (fr) 2006-01-11
SG119299A1 (en) 2006-02-28
FR2872919A1 (fr) 2006-01-13
KR20060050017A (ko) 2006-05-19
JP4899031B2 (ja) 2012-03-21
ATE371875T1 (de) 2007-09-15
FR2872919B1 (fr) 2006-09-29
EP1615049B1 (fr) 2007-08-29
KR101175087B1 (ko) 2012-08-21
DE602005002197D1 (de) 2007-10-11
ES2292064T3 (es) 2008-03-01
US20070008818A1 (en) 2007-01-11
JP2006023306A (ja) 2006-01-26
NO20053337D0 (no) 2005-07-08
NO20053337L (no) 2006-01-10
DE602005002197T2 (de) 2008-05-21
US7286443B2 (en) 2007-10-23

Similar Documents

Publication Publication Date Title
NO336331B1 (no) Fremgangsmåte for prosessering av ekkoer fra utsendte signaler.
US5440311A (en) Complementary-sequence pulse radar with matched filtering and Doppler tolerant sidelobe suppression preceding Doppler filtering
Callow Signal processing for synthetic aperture sonar image enhancement
JP4365909B2 (ja) パルス反転ドップラー超音波診断画像処理方法及び装置
JP4091067B2 (ja) ドップラー受信ビーム生成器システムのための方法とシステム
US5570691A (en) Method and apparatus for real-time, concurrent adaptive focusing in an ultrasound beamformer imaging system
US6087981A (en) Method for pulse compression with a stepped frequency waveform
CN101452073B (zh) 基于多发多收频分雷达的宽带信号合成方法
US6139501A (en) Coincident tissue and motion ultrasonic diagnostic imaging
JP2001133544A (ja) レーダ装置及びそのコヒーレント積分方法
US6072419A (en) Method for the processing of the reception signal of a deramp type synthetic aperture radar
US6184820B1 (en) Coherent pulse radar system
EP0749585A1 (en) Improvements in digital time-delay acoustic imaging
US11971470B2 (en) Beamforming hardware accelerator for radar systems
GB2355529A (en) Signal processing and scanning sonar system
US6856577B1 (en) Method and apparatus for forming a reception beam
CN104854473A (zh) 使用多普勒聚焦的亚奈奎斯特雷达处理
FR3050104A1 (no)
FR3079059A1 (fr) Filtrage adaptatif de fouillis dans une imagerie par ultrasons basée sur une force de rayonnement acoustique
CN114609595A (zh) 一种频分正交mimo雷达信号处理方法
JP4008830B2 (ja) レーダ信号処理装置
US20230408640A1 (en) System and method for radar range walk compensation in radar systems
KR20150058682A (ko) 표적 속도에 따른 도플러 효과를 보상하는 고속 lfm 표적 검출 방법 및 장치
CN117741582B (zh) 一种基于多维域编码的阵列雷达抗主瓣干扰方法及系统
RU2782249C1 (ru) Устройство цифровой обработки сигналов в импульсно-доплеровской рлс с компенсацией миграции целей по дальности за один период излучения и приема пачки радиоимпульсов

Legal Events

Date Code Title Description
MM1K Lapsed by not paying the annual fees