§¹i häc quèc gia Hµ néi Tr−êng ®¹i häc c«ng nghÖ NguyÔn Träng Kh¸nh Xö lý tù thÝch nghi theo b¨ng con vμ øng dông LuËn v¨n th¹c sÜ Hµ Néi - 2006 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com §¹i häc quèc gia Hµ Néi tr−êng ®¹i häc c«ng nghÖ NguyÔn Träng Kh¸nh Xö lý tù thÝch nghi theo b¨ng con vμ øng dông Ngµnh: C«ng nghÖ ®iÖn tö – viÔn th«ng Chuyªn ngµnh: Kü thuËt v« tuyÕn ®iÖn tö vµ th«ng tin liªn l¹c M· sè: 2.00 LuËn v¨n th¹c sÜ Ng−êi h−íng dÉn khoa häc: GS. TSKH Huúnh H÷u TuÖ Hµ Néi - 2006 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Môc lôc B¶ng ký hiÖu c¸c ch÷ viÕt t¾t. 2 Ch−¬ng 1: Tãm t¾t lý thuyÕt xö lý tù thÝch nghi. M¹ch läc tuyÕn tÝnh tèi −u.
CÊu tróc cña c¸c m¹ch läc thÝch nghi. M¹ch läc FIR Wiener. Xö lý tÝn hiÖu sè ®a tèc ®é. Bé gi¶m tèc ®é mÉu (Decimation).
Bé t¨ng tèc ®é mÉu (Interpolation). Thay ®æi tèc ®é mÉu b»ng hÖ sè ph©n sè h÷u tØ. C¸c yªu cÇu vÒ tÝnh to¸n. C¸c thuËt to¸n tù thÝch nghi vµ øng dông.
Ph−¬ng ph¸p gi¶m b−íc nhanh nhÊt. ThuËt to¸n toµn ph−¬ng trung b×nh tèi thiÓu LMS. ThuËt to¸n b×nh ph−¬ng tèi thiÓu ®Ö quy RLS. 18 Ch¦¬ng 2: CÊu tróc dμn läc vμ c¸c phÐp to¸n.
Kh¸i niÖm dµn läc sè. Quan hÖ gi÷a c¸c tÝn hiÖu trong dµn läc QMF. Dµn läc QMF kh«ng chång phæ. 29 Ch−¬ng 3: Xö lý tù thÝch nghi theo b¨ng con vμ ¸p dông vμo bμi to¸n läc nhiÔu ©m vμ lo¹i tiÕng vang.
Mét sè ch−¬ng tr×nh m« pháng c¬ së .Ph©n tÝch d¹ng tÝn hiÖu nhiÔu ©m häc. T¸c dông cña viÖc t¨ng gi¶m mÉu. Läc thÝch nghi nhiÔu ©m häc. Läc thÝch nghi nhiÔu ©m häc trªn toµn d¶i th«ng.
Läc thÝch nghi nhiÔu ©m häc dïng b¨ng con. Bé khö tiÕng vang tù thÝch nghi. Bé khö tiÕng vang tù thÝch nghi toµn d¶i. Bé khö tiÕng vang tù thÝch nghi theo b¨ng con.
61 Tμi liÖu tham kh¶o. 68 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com B¶ng ký hiÖu c¸c ch÷ viÕt t¾t SNR Signal to noise ratio AWGN Additive white gaussian noise FIR Finite Impulse Response IIR Infinite Impulse Response BER Bit error ratio LMS Least Mean Square RLS Recursive Least Square BPSK Binary pulse shift keying QPSK Quadrature pulse shift keying QAM Quadrature amplitude modulation QMF Quadrature Mirror Filter Bank OFDM Orthogonal frequency-division multiplexing UMF Uniformly modulated filter bank DFT Discrete Fourier Transform 1 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Më ®Çu Chóng ta ®Òu biÕt r»ng lo¹i bá t¹p ©m vµ khö tiÕng vang lµ c«ng viÖc rÊt cÇn thiÕt trong viÖc ®¶m b¶o chÊt l−îng trong lÜnh vùc xö lý ©m thanh [1]. Th«ng th−êng th× chÊt l−îng ©m thanh sÏ ®−îc tèt h¬n nhê bé läc tiÕng vang thÝch nghi (AEC) [2]. CÊu tróc cña bé läc thÝch nghi th«ng dông nhÊt th−êng lµ cÊu tróc LMS hoÆc RLS.
TÊt nhiªn lµ ®èi víi nh÷ng nhµ xö lý tÝn hiÖu sè chuyªn nghiÖp th× cã rÊt nhiÒu ph−¬ng ph¸p thÝch nghi kh¸c, ®−îc sö dông cho c¸c øng dông phï hîp. Giíi h¹n cña cña luËn v¨n nµy chØ lµ 2 ph−¬ng ph¸p LMS vµ RLS. VÊn ®Ò ®Æt ra ®ã lµ ®é phøc t¹p cña viÖc xö lý trªn toµn d¶i th«ng. B»ng c¸ch chia d÷ liÖu ©m thanh thµnh c¸c b¨ng con nhá phï hîp, th× ®é phøc t¹p cña hÖ thèng ®−îc gi¶m theo hÖ sè lµ sè b¨ng con.
NÕu chóng ta ¸p dông ph−¬ng ph¸p ph©n tÝch phæ rêi r¹c DFT hoÆc c¸c dµn läc lÊy mÉu ë tÇn sè giíi h¹n th× t¹i c¸c b¨ng con sÏ xuÊt hiÖn hiÖn t−îng chång phæ vµ ¶nh h−ëng tíi sù héi tô cña c¸c bé läc thÝch nghi. VÊn ®Ò nµy cã thÓ xö lÝ tèt nÕu sö dông tÇn sè lÊy mÉu cao h¬n mµ ë ®©y nghÜa lµ lo¹i bá bít c¸c mÉu d÷ liÖu (sè mÉu bá ®i gi÷a 2 lÇn lÊy mÉu ph¶i nhá h¬n sè b¨ng con). Vµ ®ã chÝnh lµ träng t©m cña luËn v¨n nµy: xö lý tù thÝch nghi theo b¨ng con vµ ¸p dông trong mét bµi to¸n cô thÓ xö lý lo¹i bá tiÕng vang. LuËn v¨n nµy tr×nh bµy tõng b−íc qu¸ tr×nh võa ®−îc ®Ò cËp ®Õn, cÊu tróc bao gåm c¸c ch−¬ng sau: Ch−¬ng 1: Tãm t¾t lý thuyÕt xö lý tù thÝch nghi Ch−¬ng 2: Tr×nh bµy cÊu tróc dµn läc vµ c¸c phÐp to¸n Ch−¬ng 3: Xö lý tù thÝch nghi b¨ng con ¸p dông trong bµi to¸n lo¹i bá tiÕng vang.
PhÇn kÕt luËn tæng kÕt l¹i nh÷ng kÕt qu¶ ®¹t ®−îc trong luËn v¨n. 2 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Ch−¬ng 1: Tãm t¾t lý thuyÕt xö lý tù thÝch nghi 1. M¹ch läc tuyÕn tÝnh tèi −u ThuËt ng÷ “mạch läc” dïng ®Ó chØ c¸c hÖ thèng cã kh¶ n¨ng kh«i phôc l¹i d¹ng cña c¸c thµnh phÇn tÇn sè cña tÝn hiÖu lèi vµo ®Ó t¹o ra tÝn hiÖu lèi ra tháa m·n c¸c yªu cÇu mong muèn. M¹ch läc cã thÓ tuyÕn tÝnh hoÆc phi tuyÕn.
Ta chØ xÐt c¸c m¹ch läc tuyÕn tÝnh, tøc c¸c m¹ch läc cã tÝnh chÊt tu©n theo nguyªn lý chång chÊt. Gi¶n ®å khèi cña bµi to¸n läc thèng kª Trong ®ã : x[n] : tÝn hiÖu vµo d [n] : tÝn hiÖu mong muèn e[n] : tÝn hiÖu sai sè TÝn hiÖu lèi vµo lµ mét d·y thêi gian rêi r¹c x[n] , m¹ch läc ®−îc ®Æc tr−ng bëi ®¸p øng xung h[n] , cßn tÝn hiÖu lèi ra ë thêi ®iÓm n lµ y[n]. Lèi ra nµy ®−îc sö dông ®Ó x¸c ®Þnh mét ®¸p øng mong muèn d [n]. Trong tr−êng hîp m¹ch läc thÝch nghi, c¸c hÖ sè cña m¹ch läc ph¶i ®−îc chän lùa sao cho d·y tÝn hiÖu mong muèn cã d¹ng phï hîp nhÊt víi tÝn hiÖu lèi vµo.
§iÒu nµy cã thÓ thùc hiÖn nÕu d·y tÝn hiÖu sai sè e[n] héi tô vÒ 0 nhanh nhÊt. §Ó lµm ®−îc ®iÒu nµy ta ph¶i tèi −u hãa mét hµm sai sè ®−îc x¸c ®Þnh theo ph−¬ng ph¸p thèng kª hoÆc theo ph−¬ng ph¸p quyÕt ®Þnh. ¾ §èi víi ph−¬ng ph¸p thèng kª, hµm sai sè ®−îc sö dông lµ gi¸ trÞ toµn ph−¬ng trung b×nh cña tÝn hiÖu sai sè e[n]. NÕu tÝn hiÖu vµo vµ tÝn hiÖu mong muèn lµ nh÷ng tÝn hiÖu dõng th× viÖc cùc tiÓu hãa sai sè toµn ph−¬ng trung b×nh ®−a ®Õn m¹ch läc Wiener – m¹ch läc tèi −u theo toµn ph−¬ng 3 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com trung b×nh.
HÇu hÕt c¸c thuËt to¸n thÝch nghi lµ ¸p dông cho c¸c lo¹i m¹ch läc Wiener. §Ó thiÕt kÕ m¹ch läc Wiener cÇn ph¶i biÕt tr−íc c¸c tÝnh chÊt thèng kª cña c¸c tÝn hiÖu c¬ së. C¸c tÝn hiÖu nµy ®−îc gi¶ thiÕt lµ ergodic, nghÜa lµ tÝn hiÖu dõng vµ trung b×nh theo thêi gian b»ng trung b×nh thèng kª. ¾ §èi víi ph−¬ng ph¸p quyÕt ®Þnh, c¸ch chän hµm sai sè lµ mét tæng träng sè cña tÝn hiÖu sai sè toµn ph−¬ng.
ViÖc cùc tiÓu hãa hµm nµy dÉn ®Õn mét m¹ch läc tèi −u ®èi víi d·y d÷ liÖu ®· cho 1. CÊu tróc cña c¸c m¹ch läc thÝch nghi CÊu tróc th−êng ®−îc sö dông trong m¹ch läc thÝch nghi [3, 4] lµ cÊu tróc ngang: x[n] x[n − 1] x[n − N + 1] • • Z −1 • Z −1 • Z −1 w0 [n] w1 [n] wN −1 [n] ⊗ ⊗ ⊗ ⊗ ⊕ y[n ] − ThuËt to¸n thÝch nghi ⊕ e[n] + d [n] H×nh 1. CÊu tróc cña mét m¹ch läc ngang thÝch nghi Trong ®ã x[n] : tÝn hiÖu lèi vµo, d [n] : tÝn hiÖu mong muèn, y[n] : tÝn hiÖu lèi ra, e[n] : tÝn hiÖu sai sè Ph−¬ng tr×nh sai ph©n m« t¶ quan hÖ gi÷a lèi ra vµ lèi vµo: N −1 y[n] = ∑ wk [n]x[n − k ] (1.1) k =0 4 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com ë ®©y N lµ chiÒu dµi cña m¹ch läc cßn wk [n] lµ c¸c hÖ sè cña m¹ch läc, c¸c hÖ sè nµy cã thÓ thay ®æi theo thêi gian vµ ®−îc ®iÒu khiÓn b»ng thuËt to¸n thÝch nghi. M¹ch läc FIR Wiener Ta chØ xÐt tr−êng hîp c¸c hÖ sè gi¸ trÞ thùc.
BiÓu diÔn lèi vµo vµ c¸c hÖ sè m¹ch läc d−íi d¹ng vÐc t¬ ta cã: w = [w0 w1 K wN −1 ] T (1.3) Ta cã lèi ra cña m¹ch läc: N −1 y[n] = ∑ wk [n]x[n − k ] = wT x[n] (1.5) §èi víi m¹ch läc Wiener, hµm hiÖu n¨ng ®−îc chän lµ sai sè toµn ph−¬ng trung b×nh: [ ξ = E e[n] 2 ] (1.7) Ta ®Þnh nghÜa vÐc t¬ t−¬ng quan chÐo bËc Nx1: p = E [x[n]d [n]] = [ p0 p N −1 ] T p1 p2 K (1.8) Vµ ma trËn tù t−¬ng quan: 5 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.9) ⎢ ⎥ ⎢ M M M K M ⎥ ⎢rN −1,0 rN −1,1 rN −1, 2 K rN −1, N −1 ⎥⎦ ⎣ Chó ý r»ng: E [d [n]xT [n]] = pT vµ wT p = pT w ta cã: ξ = E [d 2 [n]] − 2wT p + wT Rw (1.10) §Ó thu ®−îc d·y hÖ sè øng víi hµm phÝ tæn ξ cã gi¸ trÞ cùc tiÓu ta cÇn ph¶i gi¶i hÖ ph−¬ng tr×nh ®−îc t¹o thµnh tõ ®¹o hµm bËc nhÊt cña ξ ®èi víi mçi hÖ sè wi b»ng kh«ng, tøc lµ: ∂ξ = 0 víi i = 0, 1, 2, K N − 1 ∂wi ViÕt d−íi d¹ng ma trËn: Δξ = 0 ë ®©y ∇ lµ to¸n tö vi ph©n ®−îc x¸c ®Þnh nh− mét vÐc t¬ cét: ⎡ ∂ ⎤ ⎢ ∂w ⎥ ⎢ 0 ⎥ ⎢ ∂ ⎥ ∇ = ⎢ ∂w1 ⎥ ⎢ M ⎥ ⎢ ∂ ⎥ ⎢ ⎥ ⎣⎢ ∂wN −1 ⎦⎥ §Ó t×m c¸c ®¹o hµm riªng cña ξ ®èi víi c¸c hÖ sè wi cña m¹ch läc, tr−íc hÕt ta biÓu diÔn hµm phÝ tæn ξ thµnh d¹ng sau: ξ = E [d 2 [n]] − 2∑ pk wk + ∑ ∑ wk wm rkm N −1 N −1 N −1 (1.11) k =0 k =0 m =o N −1 N −1 N −1 N −1 N −1 N −1 Ta cã: ∑ ∑ wk wm rkm = ∑ ∑ wk wm rkm + wi ∑ wk rki + wi ∑ wm rkm + wi2 rki k =0 m=o k =0 m=0 k =0 m=0 k ≠i m≠i k ≠i m≠i ∂ξ N −1 Tõ ®ã: = −2 pi + ∑ wk (rki + rik ) víi i=0, 1, …,N-1 (1.12) ∂wi k =0 6 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Ta thÊy: rki = E [x[n − k ]x[n − i ]] = Φ xx [i − k ] (1.13) ë ®©y: Φ xx [i − k ] lµ hµm tù t−¬ng quan cña x[n] t−¬ng tù: rik == Φ xx [k − i ] Do tÝnh chÊt ®èi xøng cña hµm tù t−¬ng quan Φ xx [k ] = Φ xx [− k ] ta thu ®−îc: rki = rik ∂ξ N −1 Tõ ®ã: = −2 pi + 2∑ rik wk víi i = 0, N − 1 (1.