SE514875C2 - Förfarande och anordning för konstruktion av digitala filter - Google Patents
Förfarande och anordning för konstruktion av digitala filterInfo
- Publication number
- SE514875C2 SE514875C2 SE9903160A SE9903160A SE514875C2 SE 514875 C2 SE514875 C2 SE 514875C2 SE 9903160 A SE9903160 A SE 9903160A SE 9903160 A SE9903160 A SE 9903160A SE 514875 C2 SE514875 C2 SE 514875C2
- Authority
- SE
- Sweden
- Prior art keywords
- power density
- high frequency
- spectral
- spectral power
- frequency resolution
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims description 20
- 230000003595 spectral effect Effects 0.000 claims abstract description 84
- 238000012935 Averaging Methods 0.000 claims abstract description 10
- 230000001629 suppression Effects 0.000 claims abstract description 9
- 238000012546 transfer Methods 0.000 claims abstract description 9
- 230000006870 function Effects 0.000 claims description 16
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims 2
- 230000000694 effects Effects 0.000 claims 2
- 238000005530 etching Methods 0.000 claims 1
- 238000003384 imaging method Methods 0.000 claims 1
- 238000004611 spectroscopical analysis Methods 0.000 claims 1
- 238000001228 spectrum Methods 0.000 abstract description 14
- 230000011218 segmentation Effects 0.000 abstract description 3
- 238000013461 design Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 3
- 101100129499 Arabidopsis thaliana MAX2 gene Proteins 0.000 description 2
- 230000004075 alteration Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 102100022002 CD59 glycoprotein Human genes 0.000 description 1
- 101000897400 Homo sapiens CD59 glycoprotein Proteins 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001373 regressive effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Computational Linguistics (AREA)
- Quality & Reliability (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- Mathematical Physics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Noise Elimination (AREA)
- Picture Signal Circuits (AREA)
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
Fx) C 5 30 534: 875 2 SUMMERING Ett syftemål för föreliggande uppfinning är reducering eller eliminering av dessa nackdelar hos den kända tekniken.
Detta syftemäl uppnås i enlighet med de bifogade patentkraven.
Kort uttryckt anpassar föreliggande uppfinning dynamiskt segmentens längder och positioner till den aktuella formen av effektspektrum för talsig- nalen. Topparna och dalarna för spektrum bestäms, och metoden säkerstäl- ler att topparna ej delas mellan olika segment när segmenten fördelas över frekvensdomänen. Företrädesvis täcks varje topp av ett segment som är centrerat på toppen. Segmentlängden styrs företrädesvis av det mänskliga hörselsystemets frekvenskarakteristika.
Denna metod har fördelen att den reducerar komplexiteten i filterberäkning» en utan att offra noggrannhet vid de betydelsefulla spektrumtopparna.
Vidare reducerar metoden också spektrumvariansen från ram till ram, vilket förbättrar talkvaliteten.
KORT BESKRIVNING AV RITNINGARNA Uppfinningen samt ytterligare syftemäl och fördelar som uppnås med denna förstås bäst genom hänvisning till nedanstående beskrivning och de bifogade ritningarna, i vilka: Fig. 1 är ett diagram som illustrerar den spektrala effekttäthetsskatt- ningen för en brusig talsignal; Fig. 2 är ett diagram som illustrerar den spektrala effekttäthetsskatt- ningen i fig. 1 efter segmentering och medelvärdesbildning i enlighet med den kända tekniken; Fig. 3 är ett diagram liknande diagrammet i fig. 1; lO 15 20 25 5,13! 875* Fig. 4 är ett diagram som illustrerar den spektrala effekttäthetsskatt- ningen i ñg. 1 efter segmentering och medelvärdesbildning i enlighet med föreliggande uppfinning; Fig. 5 är ett flödesschema som illustrerar en exemplifierande utföringsf form av förfarandet i enlighet med föreliggande uppfinning; och Fig. 6 är ett blockschema som illustrerar en exemplífierande utförings- form av en filterkonstruktions- och ñltreríngsanordning i enlighet med föreliggande uppfinning.
DETALJERAD BESKRIVNING I vissa tillämpningar bestäms filtret i frekvensdomänen. Vid exempelvis telefontillämpningar används ofta brusundertryckning baserad på spektral subtraktion (se [2, 3}). I detta fall bestäms filtret som funktion H60) aV fre- kvens: ß H(w) = l-óißàgilï CDXMI) där O, ß, Ö är lwflslïafltcf och Cl>v(a>) och ClDX (co) är skattníngar av den spek- trala effekttätheten för rent brus resp. brusigt tal. Detta uttryck erhålles ur modellen: x[n] = y[n] + v[n] där v[n] är brussignalen, x[n] är den brusiga talsignalen och y[n] är den önskade signalen. En skattning pä den önskade signalen y[n] erhålls genom applicering av filtret som representeras av H(o)) på den brusiga signalen x[n].
(U UI 30 533 ”875 4 Ett annat exempel på en tillämpning där filtret bestäms i frekvensdomånen är en frekvensselektiv icke-linjär processor för ekosläckning. I detta fall definie- ras filtret av funktionen: H = f(x,<ï>e) där ClDXUu) representerar en skattníng på den spektrala effekttätheten för en signal x[n] innehållande residualeko och (lâeßu) representerar en skattning av den spektrala effekttätheten för residualekosignalen e[n]. Ett exempel på en lämplig funktion f utgörs av ovanstående funktion för brusundertryck» A ning, varvid CDx(c0) nu representerar den spektrala effekttäthetsskattningen av den residualeko innehållande signalen och Cl>v(a)) ersätts med den spektrala effekttäthetsskattningen e(cu) aV fesidualekot- Ffltfeï Hlffl) baseras på modellen: Xlfll = ylfll + fin] där y[n] är den önskade signalen. En skattning av den önskade signalen y[n] erhålls genom applicering av filtret som representeras av HW) på den fesídlla» leko innehållande signalen x[n]. l ovanstående exempel beskrivs filtret av en reell-värd frekvenskontinuerlig överföringsfunktíon Hgo). Denna funktion samplas för erhållande av en frekvensdiskret överföringsfunktion H[k]. Detta steg utförs i typfallet när skattningarna baseras på parametriska skattningsmetoder. Det är dock även möjligt att erhålla den frekvensdiskreta överföringsfunktionen H[k] direkt, t.ex. genom användande av periodogrambaserade skattningsmetoder. En fördel med parametriska skattningsmetoder är att skattningarna i typfallet har lägre varians från ram till ram än skattningar från periodogrambaserade metoder. lO 15 20 25 s14 a1s Följande beskrivning kommer att begränsas till brusundertryckning, men det inses att samma principer även kan användas i andra tillämpningar, såsom ekosläckning.
De frekvensdiskreta spektrala effektskattningarna Cl>x[k], QÅDÅk] är initialt kända med hög frekvensupplösning, i typfallet 128 eller 256 frekvenstappar.
Fig. 1 är ett diagram som illustrerar den spektrala effekttäthetsskattningen d^>x[k] för en brusig talsignal, i detta fall med 256 frekvenstappar. Detta spektrum erhålls ur en parametrisk skattningsmetod (i typfallet en autore~ gressiv modell, t.ex. av ordning lO). Om filteröverföringsfunktionen H[k] skall bestämmas med samma upplösning kommer detta att lägga en tung beräk- ningsbörda på den hårdvara/ mjukvara som implementerar den ovan beskrivna brusundertryckningsalgoritmen. I enlighet med [1] delas därför frekvensområdet upp i segment med konstant längd och bildas ett medel- värde av Cl>x[k] inom varje segment, såsom illustreras i fig. 2. Detta medel- värde Éxpegmenz] används istället för filbx[k] vid beräkning av H[k]. En på liknande sätt segmenterad och medelvärdesbildad skattning švßegmeny] används istället för v[k]. På detta sätt kan ett enda värde på H[k], beteck- nat 17 [segment] och definierat av T a ß (DV [segmenflï fI-[segment] = 1 - 5 Å CD x [segment] användas för ett helt segment och ej endast för ett värde k. En nackdel med denna metod är att de betydelsefulla topparna på grund av segmentens fixa längd och läge kan delas mellan flera segment, såsom illustreras vid MAX2 och MAXS i fig. 2. Detta leder till en dålig upplösning av dessa toppar.
Eftersom topparna vidare kan skifta sin position från talram till talram är de ibland uppdelade och ibland inte uppdelade, vilket leder till mycket irrite- rande ”musikbrus” (music noise).
U! 20 ša4 avs Såsom illustreras av den exemplifierande algoritmen som presenteras nedan löser föreliggande uppfinning detta problem genom att dynamiskt anpassa längden och positionen av segmenten till formen av den aktuella skattningen Cl>x[k]. Kort uttryckt börjar algoritmen med att bestämma lokala maxima (toppositioner) och lokala minima (dalpositioner) för ÖXUQ. Sedan centreras ett segment på varje maximum och fördelas segment mellan topparna för att täcka dalarna. Vid lokalisering av maxima och minima är en parametrisk spektrumskattning baserad på autoregression särskilt attraktiv, eftersom ett sådant spektrum garanterat har maximalt M/2 toppar, där M är modellord» ningen.
Företrädesvis anpassas längden av de enskilda segmenten till egenskaperna för det mänskliga hörselsystemet, som har studerats i [4]. På basis av denna information kan följande relation erhållas mellan segmentcentrum fc och segmentlängd (för ett frekvensområde på 256 tappar och en samp1ingsfre~ kvens på 8000 Hz, vilket ger en frekvensupplösning på 31.25 Hz/tapp): '93 Hz (3 tappar), o s f; s 906 Hz 155 Hz (5 tappar), 906 < fc S 1417 Hz 218 Hz (7 tappar), l4l7 281 Hz (9 tappar), 1812 < fc S 2250Hz 343 Hz (1 1 tappar), 2250 < fc S 2593 Hz 406 Hz (13 tappar), 2593 < fc S 2937 Hz 468 Hz (15 tappar), 2937 < fc S 3250 Hz x53l Hz (17 tappar), 3250 < fc S 4000 Hz segment(fc) = En omvandling till den disktreta frekvensdomänen k ger: 10 15 20 15:14 sin 7 B tappar, O S kc S 29 5 tappar, 29 < kc S 45 7 tappar, 45 < kc S 58 9 tappar, 58 < kc S 72 k -_- Segmenq Û 111 tappar, 72 < k, s 83 13 tappar, 83 < kc S 94 15 tappar, 94 < kc S 104 \l7 tappar, 104 < kc S 128 Genom användning av denna relation kan följande algoritm, vilken även illustreras i ñg. 5, användas för bestämning av en segmenterad och medel- värdesbildad ñlteröverföringsfunktion H[k] med dynamiskt bestämda segmentlängder och positioner: Sl: S2: S3: S4: S5: S6: S7: S8: S9: S10: Hämta nästa signalblock av x[n] Bestäm ClJXUc] för signalblocket Bestäm lokala maxima och minima av Cl)x[k] Sätt kmax till k-värdet för första maximum Sätt kmm till k-värdet för första minimum Sätt kc=kmax Bestäm medelvärdet av Cl>x[k] i segment(k<;) centrerat på ke Bestäm medelvärdet av Cl), [k] i samma segment Bestäm Ü [segment] genom användning av medelvärdena ClDÄk] and ätt/f] Om kc-segmenflkc) /2>kmm, utförs steg S9-S10 Sätt kc=kc-segment(kc) (det gamla värdet kc används för segment(kc)) Bestäm medelvärdet av fl>x[k] i segment(kc) centrerat på ke Bestäm medelvärdet av Cl>v[k] i samma segment Bestäm Ü [segment] genom användning av d^>x[k] och Cl>v[k] Gå till S8 U! [\) ro 30 5:11: is 15 8 S11: Sätt kc=kmax S 12: Sätt kinm till k-värdet för nästa minimum S13: Om kc+segment(kc)/2 S14: Sätt kc=kC-segment(kc) (det gamla värdet på ke används för seg- ment(kc)) S15: Bestäm medelvärdet av ClDÄk] i segment(kc) centrerat på k: Bestäm medelvärdet av ClDv[k] i samma segment Bestäm Üßegment] genom användning ClDÄk] och Cl>v[k] Gå till S13 S16: Om kmax är sista maximum, gå till S1 S17: Sätt kmax till k-värdet för nästa maximum och gå till S6 Genom användning av denna algoritm på spektrum i fig. 1 produceras det segmenterade och medelvårdesbíldade spektrumet i fig. 4. Såsom en illust- ration är lokala maxima i fig. 4 placerade vid MAXlz k = 20 MAX2: k = 41 MAXS: k = 73 och lokala mínima är placerade vid: MINI: k= O MlN2: k= 31 MIN3: k= 61 MIN4: k=l28 Applícering av ovanstående algoritm på exempelvis maximum nr. 2 vid k=41, ger ett segment på 5 tappar centrerat på k=41, två segment på 5 tappar till vänster om maximivärdet och två segment på 7 tappar till höger om maximi- värdet. En jämförelse mellan fig. 2 och ñg. 1 visar att segment som täcker en topp alltid är centrerade på toppen. Vidare noteras att lägre frekvenser 10 20 25 30 - 514: 8.75 resulterar i kortare segment. Såsom noterats ovan kommer algoritmen att garantera att segmenten alltid är centrerade på topparna och att segment- bredden anpassas till topparnas läge när dessa byter position från ram till fam .
Fig. 6 är ett blockschema som illustrerar en exemplífierande utföringsform av en filterkonstruktionsanordning i enlighet med föreliggande uppfinning, i detta fall avsedd för brusundertryckning genom spektral subtraktion. En ström av brusiga talsampel x[n] leds till en buffert 10, som samlar ihop ett block eller en ram av sampel. En spektrumestimator 12 bestämmer AR- parametrarna för detta block och använder dessa parametrar för bestämning av skattningen Cl>x[k] av den spektrala effekttätheten för det aktuella blocket av den brusiga talsignalen x[n]. I typfallet innehåller denna skattning 128 eller 256 sampel. En max-min-detektor 14 söker igenom skattningen med avseende på lokala maxima och minima. Positionerna för lokala maxima och minima leds till en segmentfördelare 18, som fördelar segmenten i enlighet med förfarandet som beskrivits och nu hänvisning till fig. 5. Segmentpositio- nerna och -längderna leds till en medelvärdesbildare 20. Medelvärdesbilda- ren 2O mottager samplen i skattningen ClDXUc] och bildar medelvärdet i varje specificerat segment. Under ett block utan tal insamlas ett block av bak- grundsbrussígnalen v[n] i en buffert 22. En spektrumestimator 24 bestäm- mer AR-parametrarna för detta block och använder dessa parametrar för bestämning av skattningen <í>v[k] av den spektrala effekttätheten för detta block av bakgrundsgrussignalen v[n]. Denna skattning har samma antal sampel som skattningen fll>x[k]. Segmentpositionerna och -längderna från segmentfördelaren 18 leds även till en annan medelvärdesbildare 26.
Medelvärdesbildaren 26 mottagare samplen i skattningen <í>v[k] och bildar medelvärdet i varje specificerat segment., Båda medelvärdesbildama 20, 26 leder medelvärdena i varje segment till en filterberäknare 28, som bestäm- mer ett värde av filteröverföringsfunktionen för varje segment. Detta produ- cerar ett segmenterat filter 17 [k] representerat av blocket 30. Detta filter leds (fl (_ 'J Såå 875 lO till en ingång av en multiplicerare 32. Signalblocket i bufferten lO leds även till ett F FT-block (Fast Fourier Transform) 34, som transformerar blocket till frekvensdomänen. Längden av det transformerade blocket är densamma som längden av det segmenterade filtret. Den transformerade signalen leds till en annan ingång av multipliceraren 32, där den multipliceras med det seg- menterade filtret. Slutligen transformeras den filtrerade signalen tillbaka till tidsdomänen i ett inverst FFT-block 36.
I typfallet implementeras de olika blocken i fig. 6 genom en eller flera mikroprocessorer eller mikro/ signalprocessor kombinationer. De kan dock även implementeras genom en eller flera ASICs (application specific integra- ted circuits).
En liknande struktur kan användas för icke-linjär filtrering vid ekosläck- ning. I detta fall representerar x[n} den residualeko innehållande signalen och ersätts v[n] med en skattning av residualekot e[n]. En annan skillnad är att i detta fall härrör skattningarna ClDXUr] och åklk] från samma talrarn (vid brusundertryfckning genom spektral subtraktion betraktas brusspektrum vara stationärt och sker skattningen under talpauser).
Fackrnannen inser att olika modifieringar och förändringar kan göras vid föreliggande uppfinning utan avvikelse från dess ram, som definieras av de bifogade patentkraven. lO [ll [3] [4] su envis 1 1 REFERENCES U.S. Patent Nr 5,839,101 (A. Vähitalo et al).
J.S. Lim och A.V. Oppenheim, “Enhancement and bandwídth com- pression of noisy speech”, Proc. of the IEEE, Vol. 67, Nr. 12, 1979, sid. 1586-1604.
S.F. Boll, “Suppression of acoustic noise in speech using spectral subtraction”, IEEE Trans. on Acoustics, Speech and Signal Processing, Vol. ASSP-27, Nr. 2, 1979, sid.1l3-120.
U. Zölser, “Digital audio signal processing”, John Wiley ön Sons, Chichester, U.K., 1997, sid. 252-253.
Claims (16)
1. Förfarande för konstruktion av digitala filter för brusundertryckning genom spektral subtraktion, innefattande stegen bestämning av en spektral effekttäthetsskattning med hög frekvens- upplösning av brusigt tal ur ett brusígt talsignalblock, bildande av en styckevis konstant spektral effekttäthetsskattning av det brusiga talet genom medelvärdesbildning av effekttätheter i frekvenstappseg- ment av den spektrala effekttäthetsskattningen med hög frekvensupplösníng av det brusiga talet, bestämning av en spektral effekttäthetsskattning med hög frekvens- upplösning av bakgrundsbrus ur ett bakgrundsbrussignalblock, bildande av en styckevis konstant spektral effekttäthetsskattning av bakgrundsbrus genom medelvärdesbildning av effekttätheter i frekvenstapp- segment av den spektrala effekttäthetsskattningen med hög frekvensupplös- ning av bakgrundsbrus, bestämning av en styckevis konstant digital filteröverföringsfunktion genom användning av spektral subtraktion baserad på nämnda styckevis konstanta spektrala effekttäthetsskattning av brusigt tal och nämnda stycke- vis konstanta spektrala effekttäthetsskattning av bakgrundsbrus, känne- tecknat av anpassning av längden av individuella segment till formen av den spek- trala effekttäthetsskattningen med hög frekvensupplösning av brusigt tal; och användning av samma segment för både den spektrala effekttäthets~ skattningen med hög frekvens av brusigt tal och den spektrala effekttäthets- skattningen med hög frekvensupplösning av bakgrundsbrus.
2. Förfarande enligt krav 1, kännetecknat av centrering av segment på lokala maxima i den spektrala effekttäthetsskattningen med hög frekvens- upplösning av brusigt tal.
3. Förfarande enligt krav l eller 2, kännetecknat av ökning av segment- längden för höga frekvenser i enlighet med det mänskliga hörselsystemet, 10 15 20 25 30 t VSÉ-išii it87f5 13
4. Förfarande för konstruktion av digitala ñlter för icke-linjär ekosläckning, innefattande stegen bestämning av en spektral taleffekttäthetsskattning med hög frekvens- upplösning ur ett residualeko innehållande talsignalblock, bildande av en styckevis konstant spektral taleffekttäthetsskattning ge- nom medelvärdesbildning av effekttätheter inom frekvenstappsegment av nämnda spektrala taleffekttäthetsskattning med hög frekvensupplösning, bestämning av en spektral resídualekoeffekttäthetsskattning med hög frekvensupplösning ur ett ekosignalblock, bildande av en styckevis konstant spektral residualekoeffekttäthetss- kattning genom medelvärdesbildning av effekttätheter inom frekvenstappseg- ment av nämnda spektrala residualekoeffekttäthetsskattning med hög frekvensupplösning, bestämning av en styckevis konstant digital filteröverföringsfunktion genom användning av nämnda styckevis konstanta spektrala taleffekttäthets- skattning och nämnda styckevis konstanta spektrala residualekoeffekt- täthetsskattning, kännetecknat av anpassning av längden av individuella segment till formen av nämnda spektrala taleffekttäthetsskattning med hög frekvensupplösning; och användning av samma segment för både nämnda spektrala taleffekttät- hetsskattning med hög frekvensupplösning och nämnda spektrala residual- ekoeffektäthetsskattning med hög frekvensupplösning.
5. Förfarande enligt krav 4, kännetecknat av centrering av segment pä lokala maxima av nämnda spektrala taleffekttäthetsskattning med hög frekvensupplösning.
6. Förfarande enligt krav 4 eller 5, kännetecknat av ökning av segment- längden för höga frekvenser i enlighet med det mänskliga hörselsystemet. 30 514 375 14
7. Förfarande för konstruktion av digitala filter, innefattande stegen bestämning av en spektral effekttäthetsskattning med hög frekvens~ upplåsning av ett insignalbloek, bildande av en styckevis konstant spektral effekttäthetsskattning genom medelvärdesbildning av effekttätheter inom frekvenstappsegment av nämnda spektrala effekttäthetsskattning med hög frekvensupplösning, bestämning av en styckevis konstant digital filteröverföringsfunktion genom användning av nämnda styckevis konstanta spektrala effekttäthets- skattning, kännetecknat av anpassning av längden av individuella segment till formen av nämnda spektrala effekttäthetsskattning med hög frekvensupplösning.
8. F örfarande enligt krav 7, kännetecknat av centrering av segment på lokala maxima av nämnda spektrala effekttäthetsskattning med hög frekvens- upplösning.
9. Anordning för konstruktion av digitala filter för brusundertryckning genom spektral subtraktion, innefattande organ för bestämning en spektral effekttäthetsskattning med hög frekvensupplösning av brusigt tal ur ett brusigt talsignalblock, organ för bildande av en styckevis konstant spektral effekttäthetsskattv ning av brusigt tal genom medelvärdesbildning av effekttätheter inom fre~ kvenstappsegment av nämnda spektrala effekttäthetsskattning med hög frekvensupplösning av brusigt tal, organ för bestämning av en spektral effekttäthetsskattning med hög frekvensupplösning av bakgrundsbrus ur ett bakgrundsbrussignalbloek, organ för bildande av en styckevis konstant spektral effekttäthetsskatt- ning av bakgrundsbrus genom medelvärdesbildning av effekttätheter inom frekvenstappsegment av nämnda spektrala effekttäthetsskattning med hög frekvensupplösning av bakgrundsbrus, organ för bestämning av en styckevis konstant digital ñlteröverförings- funktion genom användning av spektral subtraktion baserad på nämnda styckevis konstanta spektrala effekttäthetsskattning av brusigt tal och lO 15 20 25 '30 5114 8.75 15 nämnda styckevis konstanta spektrala effekttäthetsskattning av bakgrunds- brus, kännetecknad av organ (14, 16, 18) för anpassning av längden av individuella segment till formen av nämnda spektrala effekttäthetsskattning med hög frekvensupplös- ning av brusigt tal; och organ (18) för användning av samma segment för både nämnda spek- trala effekttäthetsskattning med hög frekvensupplösning av brusigt tal och nämnda spektrala effekttäthetsskattning med hög frekvensupplösning av bakgrundsbrus.
10. Anordning enligt krav 9, kännetecknad av organ (18) för centrering av segment på lokala maxima av nämnda spektrala effekttäthetsskattning med hög frekvensupplösning av brusigt tal.
11. Anordning enligt krav 9 or 10, kännetecknad av organ (18) för ökning av segmentlängden för höga frekvenser i enlighet med det mänskliga hörselsy- SÉCIIICÉ.
12. Anordning för konstruktion av digitala filter för icke-linjär ekosläckning, innefattande organ för bestämning av en spektral effekttäthetsskattning med hög frekvensupplösning av tal från ett residualeko innehållande talsignalblock, organ för bildande av en styekevis konstant spektral taleffekttäthets- skattning genom medelvärdesbildning av effekttätheter inom frekvenstapp- segment av nämnda spektrala effekttäthetsskattning med hög frekvensupp- lösning av tal, organ för bestämning av en spektral effekttäthetsskattning med hög frekvensupplösning av residualeko ur ett ekosignalblock, organ för bildande av en styckevis konstant spektral effekttäthetsskatt- ning av residualeko genom medelvärdesbildning av effekttätheter inom frekvenstappsegment av nämnda spektrala effekttäthetsskattning med hög frekvensupplösning av residualeko, i J (fl l\) UT 534! 875 16 organ för bestämning av en styckevis konstant digital filteröverförings» funktion genom användning av nämnda styckevis konstanta spektrala taleffekttäthetsskattning och nämnda styckevis konstanta spektrala effekttät- hetsskattning av residualeko, kännetecknad av organ (14, 16, 18) för anpassning av längden av individuella segment till formen av nämnda spektrala effekttäthetsskattning med hög frekvensupplös- ning av tal; och organ (18) för användning av samma segment för både nämnda spek- trala effekttäthetsskattning med hög frekvensupplösning av tal och nämnda spektrala effekttäthetsskattning med hög frekvensupplösning av residualeko.
13. Anordning enligt krav 12, kännetecknad av organ (18) för centrering av segment på lokala maxima av nämnda spektrala effekttäthetsskattning med hög frekvensupplösning av tal.
14. Anordning enligt krav 12 or 13, kännetecknad av organ (20) för ökning av segmentlängden för höga frekvenser i enlighet med det mänskliga hörselsy- SÉCITLCI.
15. Anordning för konstruktion av digitala filter, innefattande organ för bestämning av en spektral effekttäthetsskattning med hög frekvensupplösning ur ett ínsignalblock, organ för bildande av en styckevis konstant spektral effekttäthetsskatt- ning genom medelvärdesbildning av effekttätheter inom frekvenstappsegment av nämnda spektrala effekttäthetsskattning med hög frekvensupplösning, organ för bestämning av en styckevis konstant digital filteröverförings- funktion genom användning av nämnda styckevis konstanta spektrala effekttäthetsskattning, kännetecknad av organ (14, 16, 18) för anpassning av längden av individuella segment till formen av nämnda spektrala effekttäthetsskattning med hög frekvensupplös- ning. v51! :v87 5 1 7
16. Anordning enligt krav 15, kännetecknad av organ (18) för centrering av segment på lokala maxima av nämnda spektrala effekttäthetsskattning med hög frekvensupplösning.
Priority Applications (8)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
SE9903160A SE514875C2 (sv) | 1999-09-07 | 1999-09-07 | Förfarande och anordning för konstruktion av digitala filter |
AT00957195T ATE445260T1 (de) | 1999-09-07 | 2000-08-23 | Digitales filter und entwicklungsmethode zur geräuschunterdrückung mittels spektralersubtraktion |
CN00812582.1A CN1197243C (zh) | 1999-09-07 | 2000-08-23 | 一种利用频谱减法抑制噪声的数字滤波器设计方法和装置 |
DE60043115T DE60043115D1 (de) | 1999-09-07 | 2000-08-23 | Digitales filter und entwicklungsmethode zur geräuschunterdrückung mittels spektralersubtraktion |
AU68848/00A AU6884800A (en) | 1999-09-07 | 2000-08-23 | Digital filter design method and apparatus for noise suppression by spectral substraction |
PCT/SE2000/001609 WO2001018961A1 (en) | 1999-09-07 | 2000-08-23 | Digital filter design method and apparatus for noise suppression by spectral substraction |
EP00957195A EP1221197B1 (en) | 1999-09-07 | 2000-08-23 | Digital filter design method and apparatus for noise suppression by spectral substraction |
US09/656,295 US6564184B1 (en) | 1999-09-07 | 2000-09-06 | Digital filter design method and apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
SE9903160A SE514875C2 (sv) | 1999-09-07 | 1999-09-07 | Förfarande och anordning för konstruktion av digitala filter |
Publications (3)
Publication Number | Publication Date |
---|---|
SE9903160D0 SE9903160D0 (sv) | 1999-09-07 |
SE9903160L SE9903160L (sv) | 2001-03-08 |
SE514875C2 true SE514875C2 (sv) | 2001-05-07 |
Family
ID=20416885
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
SE9903160A SE514875C2 (sv) | 1999-09-07 | 1999-09-07 | Förfarande och anordning för konstruktion av digitala filter |
Country Status (8)
Country | Link |
---|---|
US (1) | US6564184B1 (sv) |
EP (1) | EP1221197B1 (sv) |
CN (1) | CN1197243C (sv) |
AT (1) | ATE445260T1 (sv) |
AU (1) | AU6884800A (sv) |
DE (1) | DE60043115D1 (sv) |
SE (1) | SE514875C2 (sv) |
WO (1) | WO2001018961A1 (sv) |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6859531B1 (en) * | 2000-09-15 | 2005-02-22 | Intel Corporation | Residual echo estimation for echo cancellation |
EP1278185A3 (en) * | 2001-07-13 | 2005-02-09 | Alcatel | Method for improving noise reduction in speech transmission |
US7187719B2 (en) * | 2003-06-18 | 2007-03-06 | Mindspeed Technologies, Inc. | Method and system for data rate optimization in a digital communication system |
CN100336307C (zh) * | 2005-04-28 | 2007-09-05 | 北京航空航天大学 | 接收机射频系统电路内部噪声的分配方法 |
CN100370519C (zh) * | 2005-10-17 | 2008-02-20 | 西安交通大学 | 一种增强电子喉语音的方法及其系统 |
US20070237341A1 (en) * | 2006-04-05 | 2007-10-11 | Creative Technology Ltd | Frequency domain noise attenuation utilizing two transducers |
US8761387B2 (en) | 2006-05-04 | 2014-06-24 | Mindspeed Technologies, Inc. | Analog transmit crosstalk canceller |
US8300849B2 (en) * | 2007-11-06 | 2012-10-30 | Microsoft Corporation | Perceptually weighted digital audio level compression |
CN101904097B (zh) * | 2007-12-20 | 2015-05-13 | 艾利森电话股份有限公司 | 噪声抑制方法和设备 |
KR20090110244A (ko) * | 2008-04-17 | 2009-10-21 | 삼성전자주식회사 | 오디오 시맨틱 정보를 이용한 오디오 신호의 부호화/복호화 방법 및 그 장치 |
CN101615538B (zh) * | 2008-06-27 | 2013-01-30 | 上海亿盟电气自动化技术有限公司 | 一种脱扣器滤波方法 |
CN101320566B (zh) * | 2008-06-30 | 2010-10-20 | 中国人民解放军第四军医大学 | 基于多带谱减法的非空气传导语音增强方法 |
CN101778183B (zh) * | 2009-01-13 | 2013-06-05 | 华为终端有限公司 | 一种残留回声抑制方法及设备 |
US8229126B2 (en) * | 2009-03-13 | 2012-07-24 | Harris Corporation | Noise error amplitude reduction |
CN102667928B (zh) * | 2009-12-25 | 2013-06-12 | 三菱电机株式会社 | 噪声消除装置以及噪声消除方法 |
US8600312B2 (en) * | 2010-01-25 | 2013-12-03 | Qualcomm Incorporated | Method and apparatus for spectral sensing |
CN102158197A (zh) * | 2010-02-11 | 2011-08-17 | 北京国能子金电气技术有限公司 | 变频率信号的自适应低通数字滤波器和方法 |
US9648421B2 (en) | 2011-12-14 | 2017-05-09 | Harris Corporation | Systems and methods for matching gain levels of transducers |
CN103021420B (zh) * | 2012-12-04 | 2015-02-25 | 中国科学院自动化研究所 | 一种基于相位调整和幅值补偿的多子带谱减法的语音增强方法 |
CN103632677B (zh) * | 2013-11-27 | 2016-09-28 | 腾讯科技(成都)有限公司 | 带噪语音信号处理方法、装置及服务器 |
CN104143337B (zh) * | 2014-01-08 | 2015-12-09 | 腾讯科技(深圳)有限公司 | 一种提高音频信号音质的方法和装置 |
EP3125429A1 (en) * | 2015-07-28 | 2017-02-01 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Partitioned block frequency domain adaptive filter device comprising adaptation modules and correction modules |
CN105812995B (zh) * | 2016-03-09 | 2018-11-09 | 西安航空学院 | 一种噪音滤除方法和装置 |
CN105877706A (zh) * | 2016-03-31 | 2016-08-24 | 济南大学 | 一种基于改进谱减法的心音增强方法 |
CN105827218B (zh) * | 2016-05-17 | 2018-04-06 | 中国电子科技集团公司第四十一研究所 | 一种对fft数据实施的噪声带限滤波方法 |
US10481831B2 (en) * | 2017-10-02 | 2019-11-19 | Nuance Communications, Inc. | System and method for combined non-linear and late echo suppression |
EP3667420A1 (en) * | 2018-12-13 | 2020-06-17 | ASML Netherlands B.V. | Method of measuring a parameter of a lithographic process, metrology apparatus |
CN111239839B (zh) * | 2020-02-10 | 2021-04-09 | 山东大学 | 一种针对频率域电磁法的频谱密度计算方法及装置 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4658426A (en) * | 1985-10-10 | 1987-04-14 | Harold Antin | Adaptive noise suppressor |
IT1184023B (it) * | 1985-12-17 | 1987-10-22 | Cselt Centro Studi Lab Telecom | Procedimento e dispositivo per la codifica e decodifica del segnale vocale mediante analisi a sottobande e quantizzazione vettorariale con allocazione dinamica dei bit di codifica |
CN1062963C (zh) * | 1990-04-12 | 2001-03-07 | 多尔拜实验特许公司 | 用于产生高质量声音信号的解码器和编码器 |
DK173156B1 (da) | 1994-10-05 | 2000-02-14 | Damixa As | Etgrebsblandingsarmatur |
SE505156C2 (sv) * | 1995-01-30 | 1997-07-07 | Ericsson Telefon Ab L M | Förfarande för bullerundertryckning genom spektral subtraktion |
FI100840B (sv) | 1995-12-12 | 1998-02-27 | Nokia Mobile Phones Ltd | Brusdämpare och förfarande för dämpning av bakgrundsbrus i brusigt tal samt en mobilteleapparat |
US5729694A (en) * | 1996-02-06 | 1998-03-17 | The Regents Of The University Of California | Speech coding, reconstruction and recognition using acoustics and electromagnetic waves |
FR2748342B1 (fr) * | 1996-05-06 | 1998-07-17 | France Telecom | Procede et dispositif de filtrage par egalisation d'un signal de parole, mettant en oeuvre un modele statistique de ce signal |
-
1999
- 1999-09-07 SE SE9903160A patent/SE514875C2/sv not_active IP Right Cessation
-
2000
- 2000-08-23 AU AU68848/00A patent/AU6884800A/en not_active Abandoned
- 2000-08-23 AT AT00957195T patent/ATE445260T1/de not_active IP Right Cessation
- 2000-08-23 DE DE60043115T patent/DE60043115D1/de not_active Expired - Lifetime
- 2000-08-23 CN CN00812582.1A patent/CN1197243C/zh not_active Expired - Fee Related
- 2000-08-23 WO PCT/SE2000/001609 patent/WO2001018961A1/en active Application Filing
- 2000-08-23 EP EP00957195A patent/EP1221197B1/en not_active Expired - Lifetime
- 2000-09-06 US US09/656,295 patent/US6564184B1/en not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
WO2001018961A1 (en) | 2001-03-15 |
ATE445260T1 (de) | 2009-10-15 |
SE9903160D0 (sv) | 1999-09-07 |
EP1221197A1 (en) | 2002-07-10 |
EP1221197B1 (en) | 2009-10-07 |
DE60043115D1 (de) | 2009-11-19 |
SE9903160L (sv) | 2001-03-08 |
CN1373930A (zh) | 2002-10-09 |
CN1197243C (zh) | 2005-04-13 |
AU6884800A (en) | 2001-04-10 |
US6564184B1 (en) | 2003-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
SE514875C2 (sv) | Förfarande och anordning för konstruktion av digitala filter | |
CN110634497B (zh) | 降噪方法、装置、终端设备及存储介质 | |
KR101120679B1 (ko) | 이득-제한된 잡음 억제 | |
JP2003534570A (ja) | 適応ビームフォーマーにおいてノイズを抑制する方法 | |
US6343268B1 (en) | Estimator of independent sources from degenerate mixtures | |
JP2002510930A (ja) | 多重非相関化法を用いた未知の混在ソースの分離 | |
SE505156C2 (sv) | Förfarande för bullerundertryckning genom spektral subtraktion | |
US20060072766A1 (en) | Reverberation removal | |
JP2005518118A (ja) | 周波数解析のためのフィルタセット | |
JP3588030B2 (ja) | 音声区間判定装置及び音声区間判定方法 | |
JP4816711B2 (ja) | 通話音声処理装置および通話音声処理方法 | |
Dharanipragada et al. | MVDR based feature extraction for robust speech recognition | |
JP2014126856A (ja) | 雑音除去装置及びその制御方法 | |
FanChiang et al. | Low complexity formant estimation adaptive feedback cancellation for hearing aids using pitch based processing | |
JP4434813B2 (ja) | 雑音スペクトル推定方法、雑音抑圧方法および雑音抑圧装置 | |
JP4965891B2 (ja) | 信号処理装置およびその方法 | |
JP2000330597A (ja) | 雑音抑圧装置 | |
JPH04245300A (ja) | 雑音除去装置 | |
KR19990028308A (ko) | 전력 스펙트럼밀도 추정방법 및 장치 | |
Agcaer et al. | Optimization of amplitude modulation features for low-resource acoustic scene classification | |
Yamashita et al. | Spectral subtraction iterated with weighting factors | |
JP2002023790A (ja) | 音声特徴量抽出装置 | |
Chinaev et al. | Fast and accurate audio resampling for acoustic sensor networks by polyphase-Farrow filters with FFT realization | |
JP2003044077A (ja) | 音声特徴量抽出方法と装置及びプログラム | |
KR101096091B1 (ko) | 음성 분리 장치 및 이를 이용한 단일 채널 음성 분리 방법 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
NUG | Patent has lapsed |