NO831719L - Fremgangsmaate og anordning for syntetisering av et kontinuerlig estimatsignal ut fra bruddstykker av et gaussisk signal fremkommet ved ultralyd-dopplermaaling paa en fluidumstroem - Google Patents

Fremgangsmaate og anordning for syntetisering av et kontinuerlig estimatsignal ut fra bruddstykker av et gaussisk signal fremkommet ved ultralyd-dopplermaaling paa en fluidumstroem

Info

Publication number
NO831719L
NO831719L NO831719A NO831719A NO831719L NO 831719 L NO831719 L NO 831719L NO 831719 A NO831719 A NO 831719A NO 831719 A NO831719 A NO 831719A NO 831719 L NO831719 L NO 831719L
Authority
NO
Norway
Prior art keywords
signal
filter
fragment
frequency
parameter
Prior art date
Application number
NO831719A
Other languages
English (en)
Inventor
Bjoern A J Angelsen
Kjell Kristoffersen
Original Assignee
Vingmed As
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 Vingmed As filed Critical Vingmed As
Priority to NO831719A priority Critical patent/NO831719L/no
Priority to GB08411946A priority patent/GB2142753B/en
Priority to IT20891/84A priority patent/IT1173985B/it
Priority to JP59094403A priority patent/JPH0616782B2/ja
Priority to DE19843417568 priority patent/DE3417568A1/de
Priority to FR848407409A priority patent/FR2549957B1/fr
Publication of NO831719L publication Critical patent/NO831719L/no
Priority to US07/213,474 priority patent/US4934373A/en

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/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S15/586Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/24Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
    • G01P5/241Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by using reflection of acoustical waves, i.e. Doppler-effect

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Multimedia (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Measuring Volume Flow (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Description

Denne oppfinnelse angår en fremgangsmåte og en anordning
for syntetisering av et kontinuerlig estimat-signal ut fra bruddstykker av et Gaussisk signal fremkommet ved ultralyd-dopplermåling på en fluidumstrøm. Oppfinnelsen er utviklet primært med sikte på anvendelse ved ultralyd-blodstrømmåling i levende biologiske strukturer. Således kan fremgangsmåten og anordningen inngå i et komplett system for undersøkelse av blodsystemer i levende biologiske strukturer, som beskrevet i samtidig norsk patentsøknad nr..1 denne forbindelse skal det også nevnes at foreliggende oppfinnelse er basert på hoved-prinsipper for syntetisering som angitt i norsk patentsøknad nr. 82.1245.
Som det vil fremgå av den følgende beskrivelse og slik det tildels fremgår av forannevnte to norske patentsøknader, er det bl.a. i forbindelse med blodstrømmålinger av stor betydning å kunne syntetisere et tilnærmet Gaussisk signal med tidsvariabel autokorrelasjonsfunksjon til bruk ved intermitterende ultralyd-dopplermåling på strømmende fluida, herunder blod som strømmer i levende vevstrukturer, for dannelse av et kontinuerlig estimat av et målt dopplersignal som opprinnelig i prinsippet fremkommer som et kontinuerlig signal.
Nærmere angivelser av fremgangsmåten og anordningen samt
de nye og særegne trekk ifølge oppfinnelsen, er opptatt i patent-kravene.
Oppfinnelsen skal i det følgende forklares nærmere under henvisning til tegningene, hvor:
Figur 1 viser eksempler på bruddstykker av et dopplersignal
slik det kan forekomme på utgangen av et høypassfilter i en dopplerprosessor for blodstrømmåling,
Figur 2 viser i prinsippet grunnstrukturen for en syntetisator basert på et eksiteringssignal i form av bredbåndet eller hvit støy, Figur 3 er et blokkskjema for en syntetisator med et transversalt filter, for syntese av et stasjonært signal ut fra bruddstykker som f.eks. vist på Figur 1, Figur 4 illustrerer multiplikasjon av et gitt signal med en veiefunksjon for dannelse av filterkoeffisienter til bruk i filteret på Figur 3, Figur 5 er et blokkskjema for en syntetisator i henhold til
anordningen ifølge oppfinnelsen, for et ikke-stasjonært komplekst signal,
Figur 6 viser en oppsplitting av en kompleks multiplikasjon
som inngår i blokkskjemaet på Figur 5, til reelle ope-rasjoner ,
Figur 7 viser eksempler på veiefunksjoner for blanding av signaler fra to filtre som vist på Figur 5,
og
Figur 8 viser et eksempel på veiefunksjon ved bruk av bare
ett filtet i blokkskemaet på Figur 5.
Figur 9 viser et blokkskjema som er et alternativ til ut- førelsen på Figur 5, og Figur 10 viser et generalisert blokkskjema som innbefatter utførelsen på Fig. 5 såvel som utførelsen på Fig. 9.
Den her omhandlede syntetisator er generelt innrettet til ikke-
syntese av et stasjonært Gaussisk signal ut fra bruddstykker av et annet foreliggende Gaussisk signal slik at de stokastiske egenskaper av det syntetiserte signal gitt ved autokorrela-sjon-funksjonen, approksimerer de stokastiske egenskaper til det foreliggende signal. Metoden kan som nevnt med fordel an-vendes ved blodstrømmåling basert på dopplereffekten for til-bakespredt ultralyd fra blod. Dopplersignalet vil i dette tilfelle være et Gaussisk signal og ved regelmessige avbrudd i dopplermålingen for korte intervaller, vil det være behov for å erstatte det direkte målte dopplersignal med et estimert signal, enten hele tiden eller deler av tiden, slik som forklart i norsk patentsøknad nr. 82.1245.
Det syntetiserte signal som fremkommer ved den i det føl-gende beskrevne metode, vil kunne benyttes som estimat på basis av bruddstykker av det direkte målte signal.
Det forutsettes at det foreligger bruddstykker på lengde T av dopplersignalet fra den dybden det gjelder, fra utgangen av høypassfilteret, med regelmessige intervaller slik som vist i Fig. 1. Disse bruddstykker benyttes til å beregne koeffisienter i et filter som påtrykkes bredbåndet (tilnærmet hvit) støy eller en annen hensiktsmessig eksitasjon, f.eks. et pulstog. Utgangen av dette filteret er da et tilnærmet Gaussisk signal, og dette benyttes som syntetisk signal. Grunnstrukturen i syn-tetisatoren er illustrert i Fig. 2.
De stokastiske egenskaper til et ikke-stasjonært Gaussisk signal er beskrevet av signalets autokorrelasjonsfunksjon R(t^, t^)• Dersom signalet er stasjonært, blir denne en funksjon av Det kan ^a defineres et effektspektrum for signalet som Fourier-transformasjonen av autokorrelasjons-funksjonen.
For ikke-stasjonære signaler kan man beregne et korttids-spektrum over så kort tid at signalet er tilnærmet stasjonært. Korttidsspektra for forskjellige utfall av samme prosess-ensemble vil være noe ulike på grunn
av stokastisk usikkerhet i spektralestimeringen. Det kan fore-tas ensemblemidling over alle korttidsspektrene og dette vil gi hastighetsfordelingen i observasjonsvolumet foldet med det spek-tralvindu som benyttes og blodets transittidsvindu gjennom ob-servas jonsvolumet .
Her skal først beskrives hvordan man kan syntetisere et stasjonært Gaussisk signal med tilnærmet samme spektrum som et gitt stasjonært Gaussisk signal, ut fra et bruddstykke av det gitte signal. Blokkskjema for en syntetisator som gjør dette, er vist i Fig. 3 for et reelt signal. N sampler av signalet veies med en veiefunksjon w^(n) slik at man får koeffisientene
Dette er illustrert i Fig. 4. Veiefunksjonen er av sam-
me type som benyttes ved spektralestimering for å redusere sidelobe-nivået og kan f.eks. være et Hamming- eller Hanning-vindu. Den benyttes her for å redusere sidelobe-nivået i spektret til det syntetiserte signal.
Koeffisientene a^(n) benyttes i et transversalt filter
som påtrykkes bredbåndet (tilnærmet hvit) støy v(n), som vist i
Fig. 3.. x(n) er det syntetiserte signal. z indikerer lagring og forsinkelse av signalet på ett skritt i n. v(n) kan være bredbåndet (tilnærmet hvit) Gaussisk støy, men bredbåndet binær støy eller annen eksitasjonskilde kan også benyttes hvis filteret har mange koeffisienter som beskrevet ovenfor. På grunn av sentral-grense-teoremet vil x(n) da bli tilnærmet Gaussisk. Dette har den fordel at man i filteret bare får multiplikasjon med - 1, og forsinkelsene kan utføres av datavipper eller digi-tale skiftregistre. Effektspekteret til det syntetiserte signal er
hvor F{ } indikerer Fourier-transf ormas jonen, W^(ta) er Fourier-transf ormas jonen til w ,G er effektspekteret til x og<*>in-
dikerér foldning i frekvensplanet. Koeffisientene som er gitt ovenfor, kan også endres gjennom visse typer lineær transfor-masjoner (all-pass-operasjon) uten at effektspekteret i det syntetiserte signal endres. Dette kan ha den fordel at fil-terets impulsrespons gjøres symmetrisk slik at en reduserer antallet multiplikasjoner til det halve, men det har også den ulempe at det må utføres en lineær transformasjon på koeffisientene først.
Syntese av komplekse Gaussiske signaler kan foregå på samme måte. (n) består da av en reell og en imaginær del. Koeffisientene a^(n) vil tilsvarende få en reell og en imagi-når del og for reell støy vil man få et komplekst syntetisert
Oi
signal x(n). Dersom det benyttes kompleks støy, vil korrela-sjonsegenskapene til x(n) ligne mer på dem til et dopplersignal.
To eksempler på syntetisatorer for et ikke-stasjonært, komplekst Gaussisk signal skal beskrives i det følgende under henvisning til Fig. 5,henholdsvis Fig. 9. Det forutsettes at bruddstykker av det opprinnelige komplekse Gaussiske signal x^(n) foreligger som vist i Fig. 1. De beskrevne syntetisatorer er prinsipielt av samme type som vist på Fig. 3 for stasjonære signaler, men filterkoeffisientene gjøres tidsvariable, slik det fremgår av det følgende. Tidsvariasjonen i koeffisientene beregnes på grunnlag av flere bruddstykker av signalet x^(n). For beregningen av tidsvariasjonen i koeffisientene utnyttes det forhold at tidsvariasjonen i fluidum-strømmens hastighet er båndbegrenset. Båndbredden av dopplersignalet endrer seg også forholdsvis langsomt, mens senter-frekvensen og maksimal frekvens undergår hurtigere endringer. Med andre ord endres frekvensspekterets posisjon hurtigere.
De to utførelseseksempler som er'vist, har som felles hovedtrekk at man først syntetiserer et signal med tilnærmet riktig båndbreddevariasjon og derefter flytter dette til det aktuelle område i frekvensplanet, ved multiplikasjon med et egnet signal.
Ifølge Fig. 5 beregnes for hvert bruddstykke x^ (n) en karakteristisk spektral parameter uk ved hjelp av en be-reghingsinnretning 161. Samtidig lagres bruddstykket i en lagringsinnretning 162. Parameteren uk kan f.eks. være maksi mal, midlere eller roten av midlere kvadratisk vinkelfrekvens. Eftersom båndbredden av signalet varierer langsomt, vil for-skjellen mellom de nevnte frekvensparametre variere langsomt. Alle disse kan derfor benyttes med tilfredsstillende resultat, unntatt når det i dopplersignalet finnes rester av signal fra vev som beveger seg. Isåfall er den maksimale vinkelfrekvens å foretrekke fordi denne påvirkes lite av signaler fra vev.
Det lagrede bruddstykket, x\ (n) , multipliseres så i en første multiplikasjbnsinnretning 163 med e win som genereres i blokk 164. Ved dette flyttes spektret av xi(n) ned rundt null (i det følgende kalt basisbåndet) og man,.får fjernet endringen i korttidsspektret gitt av uk . Den komplekse multiplikasjon kan utføres som i Fig. 6. Resultatet multipliseres så med en vindusfunksjon w^(n) i blokk 165, på samme måte som i Fig. 4. Dette danner filterkoeffisienter for syntese av et signal som ligger i basisbåndet, på samme måte som vist i Fig. 3. Blokken 165 utgjør en annen multiplikasjonsinnretning hvis funksjon eventuelt kan kombineres med den multiplikasjon som skjer i den første multiplikasjonsinnretning 163, idet innretningen 165 vil kunne ansees inkorporert i blokk 164,
til en mer generell genereringsinnretning.
I Fig. 5 er det illustrert to filtre 166 og 167 for syntese av signalet i basisbåndet. Koeffisientene lastes inn i de respektive filtre for annet hvert bruddstykke av det opprinnelige signal. Utgangen av filtrene veies med veiefunksjoner ved hjelp av ytterligere multiplikasjonsinnretninger 168 og 169 som illustrert på Fig. 7. Dette gjør at bidraget efter summasjonen fra de to filtrene 166 og 167 blir null når endringen av koeffisienter i vedkommende filter skjer.
På grunn av den gradvise endringen i vindusfunksjonen vil en
få en gradvis overgang fra at det syntetiserte signal er gitt av koeffisientene fra ett bruddstykke til det er gitt av koeffisientene for neste bruddstykke. I Fig. 7 er det vist vinduer med lineært økende flanker, men andre vinduer kan benyttes, f.eks. Hamming-vindu o.l. Vinduene må være null når inn-lesning av koeffisienter i vedkommende filter foregår. Denne tiden kan dog gjøres meget kort (10 - 50 ys).
På basis av de karakteristiske vinkelfrekvenser co. for
i
flere bruddstykker av det opprinnelige signal, estimeres en kon-
tmuerlig variabel vinkelfrekvens co(n) i en estimeringsenhet 171 på Fig. 5. Signalet i basisbåndet efter blokk 170 multipliseres så i blokk 173 med e<iw>(<n>)<n>generert i blokk 172. Ved dette flyttes spekteret fra basisbåndet til det aktuelle området.
En enkel metode for estimering av to(n) er ved lineær interpolasjon mellom uk og Den kan også dannes ved filtrering av en rekke verdier av uk , siden endringen i karakteristisk vinkelfrekvens er båndbegrenset.
Ved dannelse av co(n) ved lineær interpolasjon mellom co.
og + må man forsinke innlesningen av koeffisientene i filtrene 166 og 167 med ett trinn slik at foreligger når man laster inn koeffisientene fra x (n). Dette gir en forsinkelse mellom syntetisert signal og opprinnelig signal, og anvendelsen avgjør hvorvidt dette kan tolereres.
En forenklet estimator kan også benytte bare ett filter istedenfor de to filtre på Fig. 5. Ett forslag til veiefunk-
sjon for dette tilfellet er vist i Fig. 8.
Den alternative utførelse på Fig. 9 adskiller seg fra
den på Fig. 5 ved at den første multiplikasjon av signalet med exp (-iuKn) er utelatt. Til gjengjeld må utgangen fra filtrene 266 og 267 multipliseres med to forskjellige signaler som representerer de instantane avvik mellom den i løpende tid eStimerte spektrale parameter ^<t>Oo(n) og de tilsvarende spektrale parametre u^jog <^2j+l^or ^ver av ^ ilterkoef f isientene.
Nærmere bestemt viser Fig. 9 en lagringsinnretning 262 for bruddstykkesignalene. Videre er det anordnet en beregningsinnretning 261 og en estimeringsenhet 271, idet disse blokker 261, 262 og 271 svarer til de respektive blokker 161, 162 og 171 på Fig. 5. Efter lagringsinnretningen 262 følger en multi-plikas jonsinnretning 263 som mottar utgangssignalet fra innretningen 262 og som har en annen inngang for påtrykning av en veie- eller vindusfunksjon wf(n). Det fremgår således at bruddstykkesignalene x^her brukes direkte til å danne filterkoeffisientene for den følgende filteranordning bestående av filtrene 266 og 267. På utgangen av filtrene er det anordnet ytterligere multiplikasjonsinnretninger henholdsvis 273 a og 273 b som også mottar respektive signaler
exp i(u)(n)-^2j^n°9exPi(w(n)- w2j+l^n ^ra i1111176^11!11?6116 272 a henholdsvis 272 b, som mottar frekvensparameteren u)(n) fra estimeringsenheten 271. Dermed blir filter-utgangssignalene flyttet tilbake til frekvensområdet for det opprinnelige signal. Derefter blir signalene veiet i multi-plikas jonsinnretninger 268 og 269 som også befinner seg på utgangen av filtrene. Disse multiplikasjoner kunne imidlertid også skje umiddelbart efter filtrene, dvs. foran innret-ningene 273 a og 273 b. Endelig følger i blokk 270 en summering av signalene fra de to filtre. Rekkefølgen av disse funksjoner kan tildels ombyttes slik som allerede antydet,
og Fig. 5 viser en rekkefølge hvor summeringen i blokk 170 går forut for transponeringen i blokk 173. I praksis blir utførelsen på Fig. 9 ansett for å være mest fordelaktig da den gir en enklere oppbygning av anordningen. Imidlertid er det klart at det ønskede resultat blir oppnådd såvel med ut-førelsen på Fig. 5 som med utførelsen på Fig. 9 og det er følgelig et spørsmål om hva som er praktisk hensiktsmessig, avhengig bl.a. av tilgjengelig teknologi og komponenter, hvilket alternativ som skal velges.
Fig. 10 viser et blokkskjema som er generelt for de to alternative utførelser på Fig. 5 og 9. Hvis man på Fig. 10 setter Ato^ = 0, fremkommer utførelsen på Fig. 9, og dersom man setter Auk = uk, fremkommer utførelsen på Fig. 5. I dette, tilfelle vil blandesignalet efter filtrene være det samme for begge filtre slik at multiplikasjonen kan utføres efter summeringen. I praksis vil man enten velge Aco^= 0 eller Acb. = u). fordi dette resulterer i bare to multiplikasjoner,
ilcJ
Endelig kan man også tenke seg benyttet flere enn to filtre for dermed å ta bedre vare på variasjonene i signalets båndbredde. Isåfall benyttes alternativet med Auk = uk slik at det bare trenges to komplekse multiplikasjoner.

Claims (17)

1. Fremgangsmåte for syntetisering av et kontinuerlig estimatsignal av et Gaussisk signal fremkommet ved ultralyd-dopplermåling på en fluidumstrøm, ut fra foreliggende bruddstykker av det Gaussiske signal, karakterisert ved at bruddstykkesignalene legges til grunn for bestem-melse av tidsvariable koeffisienter i en styrt filteranordning som påtrykkes et bredbåndet inngangssignal, slik at det fremkommer et syntetisert signal med tilnærmet riktige bånd-breddevariasjoner som det Gaussiske signal, og at det syntetiserte signal ved multiplikasjon med et annet signal avledet fra bruddstykkesignalene, flyttes til det aktuelle frekvens-område, slik at det syntetiserte kontinuerlige estimatsignal blir dannet med stokastiske egenskaper som tilnærmet svarer til de stokastiske egenskaper ved det foreliggende signal.
2. Fremgangsmåte ifølge krav 1, karakterisert ved a) at bruddstykke-signalene først lagres b) at det samtidig beregnes en karakteristisk frekvensparameter (uk ) for hvert bruddstykke, c) at det lagrede bruddstykkesignal benyttes som grunnlag for å danne koeffisientene i den styrte filteranordning, som omfatter ett eller flere filtre og som med det bredbåndete inngangssignal gir et syntetisert.filterutgangssig-nal som er et tilnærmet Gaussisk signal, d) at hvert filter-utgangssignal multipliseres med en veie-eller vindusfunksjon for utjevning av overgangene mellom på-følgende bruddstykker, e) at utgangssignaler fra eventuelt to eller flere filtre summeres efter multiplikasjonen, f) at det parallelt med signalbehandlingen i trinnene a)-c) ovenfor på basis av frekvensparameteren (uk ) for hvert bruddstykke av det opprinnelige Gaussiske signal, dannes en parameter (u)(n)) for hvert tidspunkt, som er et estimat av frekvensparameteren til det opprinnelige signal mellom bruddstykkene og derved representerer de forholdsvis hurtige variasjoner i beliggenheten av det opprinnelige signals frekvensspektrum, og g) at den sistnevnte estimerte frekvensparameter direkte eller indirekte kombineres med filter-utgangssignalet før eller efter multiplikasjonen og eventuell summering, slik at det syntetiserte, kontinuerlige estimatsignal blir dannet.
3. Fremgangsmåte ifølge krav 2, karakterisert ved at den karakteristiske frekvensparameter kom- . bineres med det lagrede bruddstykkesignal for transponering av dette til et basisbånd ved eller omkring frekvens-nullpunktet, og at kombineringen av den nevnte estimerte f rekvensparameter ('co (n) ) med f ilter-utgangssignalet. skjer vindus- eftervmultiplikasjonen og eventuell summering, for tilbaketransponering til frekvensområdet for det opprinnelige signal (Fig. 5)
4. Fremgangsmåte ifølge krav 2, karakterisert ved at bruddstykkesignalene (x^ ) brukes direkte til å danne filterkoeffisientene og at det tilsvarende filter-utgangssignal blandes med et signal exp ii(cu(n) - co.)n} dannet av den estimerte f rekvensparameter (co(n)).
5. Fremgangsmåte ifølge et av de foregående krav for anvendelse ved blodstrømmåling på levende biologiske strukturer, karakterisert ved at den karakteristiske f rekvensparameter (co^ ) er den maksimale frekvens i bruddstykke-signalene.
6. Fremgangsmåte ifølge et av de foregående krav for anvendelse ved blodstrømmåling på levende biologiske strukturer, karakterisert ved at den karakteristiske frekvensparameter (uk) er den midlere frekvens i bruddstykke-signalene.
7. Fremgangsmåte ifølge et av de foregående krav for anvendelse ved blodstrømmåling på levende biologiske strukturer, karakterisert ved at den karakteristiske f rekvensparameter (co^ ) er kvadratroten av den midlere kvadra-tiske frekvens i bruddstykke-signalene.
8. Fremgangsmåte ifølge et av de foregående krav, karakterisert ved at basisbånd-signalet for påfølgende bruddstykker vekselvis benyttes til å danne koeffisienter til et første og et annet styrt filter i en anordning av to parallelle filtre hvis utgangssignaler fortrinnsvis hver for seg multipliseres med en veie- eller vindusfunksjon for utjevning av overgangen mellom påfø lgende bruddstykker.
9. Fremgangsmåte ifølge krav 8, karakterisert ved at veie- eller vindusfunksjonen efter de(t) styrte filter(e) har monotont økende resp. avtagende forløp, fortrinnsvis med en verdi lik null når filteranordningens filterkoeffisienter endres.
10. Fremgangsmåte ifølge et av de foregående krav, karakterisert ved at bruddstykkesignalene eller signaler avledet av disse multipliseres med en veiefunksjon og benyttes sdm koeffisienter i en filteranordning bestående av ett eller flere parallelle transversale filtre.
11. Fremgangsmåte ifølge krav 10, karakterisert ved at koeffisientene underkastes en lineær allpasstransformasjon. .
12. Fremgangsmåte ifølge et av de foregående krav, karakterisert ved at tidsvariasjonen i filterkoeffisientene beregnes på grunnlag av to eller flere bruddstykker.
13. Fremgangsmåte ifølge et av de foregående krav, karakterisert ved at estimeringen av den kontinuerlig variable frekvensparameter skjer ved interpolasjon mellom den karakteristiske frekvensparameter for to eller flere påfølgende bruddstykker.
14. Anordning for utførelse av fremgangsmåten ifølge krav 3, karakterisert ved at den omfatter: en lagringsinnretning (162) for bruddstykkesignaler, en beregningsinnretning (161) for beregning og lagring av en karakteristisk spektral parameter for hvert bruddstykke, en innretning (164) til å generere et signal som repre senterer den spektrale parameter, en første multiplikasjonsinnretning (163) som har inn- ganger forbundet med utgangen av lagringsinnretningen (162) resp. genereringsinnretningen (164), en annen multiplikasjonsinnretning (165) hvis ene inngang er forbundet med utgangen av den første multiplikasjonsinnretning, og hvis annen inngang er innrettet til å påtrykkes en (veie-) vindusfunksjon, - en filteranordning (166, 167) som mottar koeffisienter avgitt fra utgangen av den annen multiplikasjonsinnretning (165) og minst ett kontinuerlig eksitasjonssi <q> nal fra en eksitasjonskilde, for frembringelse av et syntetisert, kontinuerlig utgangssignal, en tredje multiplikasjonsinnretning (168, 169) på utgangen av filteranordningen (166, 167) for påtrykning av en veie- eller vindusfunksjon, en estimeringsenhet (171) som på basis av frekvenspara meteren (ok ) for hvert bruddstykke av det opprinnelige Gaussiske signal danner en parameter (uj(n)) for hvert tidspunkt, som er et estimat av frekvensparameteren til det opprinnelige signal mellom bruddstykkene og derved representerer de forholdsvis hurtige variasjoner i beliggenheten av det opprinnelige signals frekvensspektrum, og en enhet (172) som danner den komplekse eksponensial- funksjon av ovennevnte f rekvensparameter (co(n)) multiplisert med tiden og en multiplikasjonsenhet (173) for tilbaketransponering av det utjevnede filter-utgangssignal slik at dette legges til frekvensområdet for det opprinnelige signal og derved danner det syntetiserte, kontinuerlige estimatsignal.
15. Anordning for utfø relse av fremgangsmåten ifølge krav 4, karakterisert ved at den omfatter: - en lagringsinnretning (262) for bruddstykkesignaler, en beregningsinnretning (261) for beregning og lagring av en karakteristisk spektral parameter for hvert bruddstykke, en inngangs-multiplikasjonsinnretning (263) hvis ene inn gang er forbundet med utgangen av lagringsinnretningen (262) og hvis annen inngang er innrettet til å påtrykkes en veie- eller vindusfunksjon, en filteranordning (266, 267) som mottar koeffisienter avgitt fra inngangs-multiplikasjonsinnretningen (263) og minst ett kontinuerlig eksitasjonssiqnal fra en eksitasjonskilde, for frembringelse av et syntetisert kontinuerlig utgangssignal, en utgangs-multiplikasjonsinnretning (268, 269) efter ut gangen av filteranordningen (266, 267) for påtrykning av en veie- eller vindusfunksjon, en estimeringsenhet (271) som-på basis av frekvens parameteren (co.) for hvert bruddstykke av det opprinnelige Gaussiske signal danner en parameter (co(n)) for hvert tidspunkt, som er et estimat av frekvensparameteren til det opprinnelige signal mellom bruddstykkene og derved representerer de forholdsvis hurtige variasjoner i beliggenheten av det opprinnelige signals frekvensspektrum, - en innretning (272 a, b) som er forbundet med utgangen av estimeringsenheten (271) og danner et signal exp{i(co(n) - co^ )n} av den estimerte parameter (co(n)), og en ytterligere multiplikasjonsinnretning (273 a, b) for an eller efter utgangs-multiplikasjonsinnretningen for blanding av det sistnevnte signal med filterutgangs- signalet, slik at dette flyttes til frekvensområdet for det opprinnelige signal og derved danner det syntetiserte, kontinuerlige estimatsignal.
16. Anordning ifølge krav 14 eller 15, karakterisert ved at filteranordningen består av to parallelle filtre (166, 167, 266, 267) og at bruddstykkesignalene eller signaler avledet av disse benyttes til vekselvis å danne filterkoeffisienter for det ene og det annet filter.
17. Anordning ifølge krav 14, 15 eller 16, karakterisert ved at filteranordningen består av ett eller flere transversale filtre hvis koeffisienter bestemmes av bruddstykkesignalene eller signaler avledet av disse, multiplisert med en vindusfunksjon.
NO831719A 1983-05-13 1983-05-13 Fremgangsmaate og anordning for syntetisering av et kontinuerlig estimatsignal ut fra bruddstykker av et gaussisk signal fremkommet ved ultralyd-dopplermaaling paa en fluidumstroem NO831719L (no)

Priority Applications (7)

Application Number Priority Date Filing Date Title
NO831719A NO831719L (no) 1983-05-13 1983-05-13 Fremgangsmaate og anordning for syntetisering av et kontinuerlig estimatsignal ut fra bruddstykker av et gaussisk signal fremkommet ved ultralyd-dopplermaaling paa en fluidumstroem
GB08411946A GB2142753B (en) 1983-05-13 1984-05-10 Synthesizing continuous estimate signal
IT20891/84A IT1173985B (it) 1983-05-13 1984-05-11 Procedimento ed apparecchiatura per sintetizzare un segnale stimato continuo da segmenti di un segnale gaussiano ottenuti mediante misurazione doppler ultrasonica su un flusso di fluido
JP59094403A JPH0616782B2 (ja) 1983-05-13 1984-05-11 連続した推定信号を実時間で合成する方法及び装置
DE19843417568 DE3417568A1 (de) 1983-05-13 1984-05-11 Verfahren und vorrichtung zur synthetisierung eines kontinuierlichen schaetzsignals aus segmenten eines bei ultraschall-doppler-messungen aus einer fluessigkeitsstroemung erhaltenen gaussschen signals
FR848407409A FR2549957B1 (fr) 1983-05-13 1984-05-14 Procede et appareil pour synthetiser un signal continu estime a partir de segments d'un signal gaussien resultant de la mesure doppler par ultra-sons de l'ecoulement d'un fluide
US07/213,474 US4934373A (en) 1983-05-13 1988-06-30 Method and apparatus for synthesizing a continuous estimate signal from segments of a Gaussian signal provided by ultra sonic doppler measurement on a fluid flow

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
NO831719A NO831719L (no) 1983-05-13 1983-05-13 Fremgangsmaate og anordning for syntetisering av et kontinuerlig estimatsignal ut fra bruddstykker av et gaussisk signal fremkommet ved ultralyd-dopplermaaling paa en fluidumstroem

Publications (1)

Publication Number Publication Date
NO831719L true NO831719L (no) 1984-11-14

Family

ID=19887096

Family Applications (1)

Application Number Title Priority Date Filing Date
NO831719A NO831719L (no) 1983-05-13 1983-05-13 Fremgangsmaate og anordning for syntetisering av et kontinuerlig estimatsignal ut fra bruddstykker av et gaussisk signal fremkommet ved ultralyd-dopplermaaling paa en fluidumstroem

Country Status (7)

Country Link
US (1) US4934373A (no)
JP (1) JPH0616782B2 (no)
DE (1) DE3417568A1 (no)
FR (1) FR2549957B1 (no)
GB (1) GB2142753B (no)
IT (1) IT1173985B (no)
NO (1) NO831719L (no)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5016641A (en) * 1989-11-13 1991-05-21 Advanced Technology Laboratories, Inc. Spectral interpolation of ultrasound Doppler signal
JP2678124B2 (ja) * 1993-01-08 1997-11-17 アロカ株式会社 超音波ドプラ診断装置
JPH06342022A (ja) * 1993-03-17 1994-12-13 Hewlett Packard Co <Hp> ゲートスペクトル分析方法
US5433206A (en) * 1994-06-29 1995-07-18 Elscint, Ltd. System and a method for simultaneous, real time ultrasound imaging of biological tissue and measuring of blood flow velocity
US5476097A (en) * 1994-10-13 1995-12-19 Advanced Technology Laboratories, Inc. Simultaneous ultrasonic imaging and Doppler display system
JP3657706B2 (ja) * 1996-09-11 2005-06-08 株式会社日立メディコ 超音波ドプラ診断装置
US5808195A (en) * 1997-05-21 1998-09-15 Ads Environmental Services Arrangement for determining liquid velocity versus depth utilizing historical data
KR100248065B1 (ko) * 1997-12-30 2000-03-15 윤종용 디지털 신호 처리장치 및 그 방법
CN101161205B (zh) * 2006-10-13 2011-07-06 深圳迈瑞生物医疗电子股份有限公司 多普勒血流声音信号抗混叠的方法及装置
US10330508B2 (en) * 2015-12-09 2019-06-25 Texas Instruments Incorporated Ultrasonic flowmeter using windowing of received signals

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4111055A (en) * 1977-10-11 1978-09-05 Del Mar Avionics Pulse echo ultrasonic system
US4182173A (en) * 1978-08-23 1980-01-08 General Electric Company Duplex ultrasonic imaging system with repetitive excitation of common transducer in doppler modality
FR2447041A1 (fr) * 1979-01-19 1980-08-14 Inst Nat Sante Rech Med Perfectionnements aux velocimetres doppler a bruit pseudo-aleatoires
US4324258A (en) * 1980-06-24 1982-04-13 Werner Huebscher Ultrasonic doppler flowmeters
US4407293A (en) * 1981-04-24 1983-10-04 Diasonics, Inc. Ultrasound imaging apparatus for providing simultaneous B-scan and Doppler data
NO150015C (no) * 1981-11-13 1984-08-08 Vingmed As Fremgangsmaate ved blodstroemhastighetsmaaling med ultralyd, kombinert med ekko-amplitudeavbildning, for undersoekelse av levende biologiske strukturer
US4441368A (en) * 1982-04-19 1984-04-10 General Electric Company Method and means for determining ultrasonic wave attenuation in tissue
NO831718L (no) * 1983-05-13 1984-11-14 Vingmed As Fremgangsmaate og apparat ved blodstroem-hastighetsmaaling med ultralyd for dannelse av todimensjonal avbildning av blodets hastighet

Also Published As

Publication number Publication date
FR2549957A1 (fr) 1985-02-01
FR2549957B1 (fr) 1989-12-08
GB2142753A (en) 1985-01-23
JPS6035221A (ja) 1985-02-23
IT8420891A0 (it) 1984-05-11
GB2142753B (en) 1986-10-29
IT8420891A1 (it) 1985-11-11
DE3417568A1 (de) 1984-11-29
US4934373A (en) 1990-06-19
GB8411946D0 (en) 1984-06-13
IT1173985B (it) 1987-06-24
DE3417568C2 (no) 1993-02-11
JPH0616782B2 (ja) 1994-03-09

Similar Documents

Publication Publication Date Title
US4848354A (en) Method and apparatus for investigating a circulatory system in living biological structures
EP2570782B1 (en) Doppler measurement instrument and doppler measurement processing method
US20060195279A1 (en) Method and system for the detection and/or removal of sinusoidal interference signals in a noise signal
NO831719L (no) Fremgangsmaate og anordning for syntetisering av et kontinuerlig estimatsignal ut fra bruddstykker av et gaussisk signal fremkommet ved ultralyd-dopplermaaling paa en fluidumstroem
EP1884196A1 (en) Velocity measuring method and velocity measuring device using the same
JPH065192B2 (ja) 振動制御装置
CN108324319A (zh) 用于无失真多波束超声接收波束形成的***和方法
US20080156106A1 (en) Method and apparatus for dealiasing Doppler signals
JPH11142425A (ja) 流速測定装置および超音波診断装置
AU2018301988B2 (en) Method and device for quantifying medium viscoelasticity
EP3127487A1 (en) Ultrasonic imaging system
JPS60153849A (ja) 媒体の構造を特徴付けする方法及びこの方法を実施する装置
Hein et al. A real-time ultrasound time-domain correlation blood flowmeter. I. Theory and design
JP4369427B2 (ja) ドプラ速度検出装置及びそれを用いた超音波診断装置
US5537344A (en) High-speed processing apparatus and method, signal analyzing system, and measurement apparatus and method
JPH06254095A (ja) 複素mtiフィルタ
JP4698003B2 (ja) 超音波診断装置
US5574674A (en) Fourier transform processing for digital filters or other spectral resolution devices
JP5612984B2 (ja) 超音波診断装置
JP2019000158A (ja) 超音波診断装置
JP5559655B2 (ja) 超音波診断装置
EP1853171B1 (en) Method and device for determining the motion vector of tissues in a biological medium
JP2011217898A (ja) 超音波診断装置
JP2003250802A (ja) 超音波診断装置
JP2012147894A (ja) 超音波診断装置