Robustness Techniques for Speech Recognition

Size: px
Start display at page:

Download "Robustness Techniques for Speech Recognition"

Transcription

1 Robustness Techniques for peech Recognition Berin Chen, 2005 References: 1. X. Huang et a. poken Language Processing (2001). Chapter J. C. Junqua and J. P. Haton. Robustness in Automatic peech Recognition (1996), Chapters 5, T. F. Quatieri, Discrete-Time peech igna Processing (2002), Chapter 13

2 Introduction Cassification of peech Variabiity in Five Categories Pronunciation Variation Intra-speaker variabiity Linguistic variabiity Inter-speaker variabiity peaker-independency peaker-adaptation peaker-dependency Robustness Enhancement Variabiity caused by the environment Variabiity caused by the contet Contet-Dependent Acoustic Modeing peech - Berin Chen 2

3 Introduction (cont.) The Diagram for peech Recognition Acoustic Processing Linguistic Processing peech signa Feature Etraction Likeihood computation Linguistic etwork Decoding Recognition resuts Acoustic mode mode Language mode mode Leicon Importance of the robustness in speech recognition peech recognition systems have to operate in situations with uncontroabe acoustic environments The recognition performance is often degraded due to the mismatch in the training and testing conditions Varying environmenta noises, different speaker characteristics (se, age, diaects), different speaking modes (styistic, Lombard effect), etc. peech - Berin Chen 3

4 Introduction (cont.) If a speech recognition system s accuracy does not degrade very much under mismatch conditions, the system is caed robust AR performance is rather uniform for Rs greater than 25dB, but there is a very steep degradation as the noise eve increases 25dB Es Es 2.5 = 10 og 10 => = 10 E E 316 igna energy measured in time domain, e.g.: E s = 1 T 1 T * [] n s [] n s n = 0 Various noises eist in varying rea-word environments periodic, impusive, or wide/narrow band peech - Berin Chen 4

5 Introduction (cont.) Therefore, severa possibe robustness approaches have been deveoped to enhance the speech signa, its spectrum, and the acoustic modes as we Environmenta compensation processing (feature-based) Acoustic mode adaptation (mode-based) Inherenty robust acoustic features (both mode- and featurebased) Discriminative acoustic features peech - Berin Chen 5

6 peech - Berin Chen 6 The oise Types [ ] [ ] [ ] [ ] ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) { } ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) () ( ) ( ) ( ) ( ) () spectrum :, or power spectrum :, or cos 2 2 Re * = + = = + + = + = + = P P P P P H H H H H X H X m n m h m s m nn hh ss H X θ h[m] s[m] n[m] [m] A mode of the environment.

7 Additive oises Additive noises can be stationary or non-stationary tationary noises uch as computer fan, air conditioning, car noise: the power spectra density does not change over time (the above noises are aso narrow-band noises) on-stationary noises Machine gun, door sams, keyboard cicks, radio/tv, and other speakers voices (babbe noise, wide band nose, most difficut): the statistica properties change over time peech - Berin Chen 7

8 Additive oises (cont.) peech - Berin Chen 8

9 Convoutiona oises Convoutiona noises are mainy resuted from channe distortion (sometimes caed channe noises ) and are stationary for most cases Reverberation, the frequency response of microphone, transmission ines, etc. peech - Berin Chen 9

10 oise Characteristics White oise The power spectrum is fat nn =,a condition equivaent to different sampes being uncorreated, R nn [ m] = qδ [ m] White noise has a zero mean, but can have different distributions We are often interested in the white Gaussian noise, as it resembes better the noise that tends to occur in practice Coored oise ( ) q The spectrum is not fat (ike the noise captured by a microphone) Pink noise A particuar type of coored nose that has a ow-pass nature, as it has more energy at the ow frequencies and ros off at high frequency E.g., the noise generated by a computer fan, an air conditioner, or an automobie peech - Berin Chen 10

11 Musica oise oise Characteristics (cont.) Musica noise is short sinusoids (tones) randomy distributed over time and frequency That occur due to, e.g., the drawback of origina spectra subtraction technique and statistica inaccuracy in estimating noise magnitude spectrum Lombard effect A phenomenon by which a speaker increases his voca effect in the presence of background noise (the additive noise) When a arge amount of noise is present, the speaker tends to shout, which entais not ony a high ampitude, but aso often higher pitch, sighty different formants, and a different cooring (shape) of the spectrum The vowe portion of the words wi be overemphasized by the speakers peech - Berin Chen 11

12 Robustness Approaches

13 Three Basic Categories of Approaches peech Enhancement Techniques Eiminate or reduce the noisy effect on the speech signas, thus better accuracy with the originay trained modes (Restore the cean speech signas or compensate for distortions) The feature part is modified whie the mode part remains unchanged Mode-based oise Compensation Techniques Adjust (changing) the recognition mode parameters (means and variances) for better matching the testing noisy conditions The mode part is modified whie the feature part remains unchanged Inherenty Robust Parameters for peech Find robust representation of speech signas ess infuenced by additive or channe noise Both of the feature and mode parts are changed peech - Berin Chen 13

14 Assumptions & Evauations Genera Assumptions for the oise The noise is uncorreated with the speech signa The noise characteristics are fied during the speech utterance or vary very sowy (the noise is said to be stationary) The estimates of the noise characteristics can be obtained during non-speech activity The noise is supposed to be additive or convoutiona Performance Evauations Inteigibiity, quaity (subjective assessment) Distortion between cean and recovered speech (objective assessment) peech recognition accuracy peech - Berin Chen 14

15 pectra ubtraction (). F. Bo, 1979 A peech Enhancement Technique Estimate the magnitude (or the power) of cean speech by epicity subtracting the noise magnitude (or the power) spectrum from the noisy magnitude (or power) spectrum Basic Assumption of pectra ubtraction The cean speech s [ m] is corrupted by additive noise n[ m] Different frequencies are uncorreated from each other s [ m] and n[ m] are statisticay independent, so that the power spectrum of the noisy speech [ m] can be epressed as: P ( ) = P ( ) + P ( ) X To eiminate the additive noise: P ( ) = P ( ) P ( ) X We can obtain an estimate of P ( ) using the average period of M frames that known to be just noise: 1 M ( ) 1 Pˆ = P ( ),i i= 0 M frames peech - Berin Chen 15

16 pectra ubtraction (cont.) Probems of pectra ubtraction s [ m] and n[ m] are not statisticay independent such that the cross term in power spectrum can not be eiminated ( ) is possiby ess than zero Pˆ Introduce musica noise when ( ) P ( ) eed a robust endpoint (speech/noise/sience) detector P X peech - Berin Chen 16

17 pectra ubtraction (cont.) Modification: oninear pectra ubtraction () Pˆ P PX ( ) P ( ), if PX ( ) P ( ) ( ) = P ( ), otherwise ( ) and P ( ): smoothed noisy and noise spectrum X or Pˆ P φ P ( ) ( ) ( ) ( ) ( ) X φ, if PX > φ + β P ( ) = β P ( ), otherwise ( ) and P ( ) X : smoothed noisy and noise ( ): a non - inear function according to R spectrum peech - Berin Chen 17

18 pectra ubtraction (cont.) pectra ubtraction can be viewed as a fitering operation Power pectrum Pˆ ( ) = PX ( ) P ( ) P = ( ) ( ) P 1 X P ( ) = P X ( ) 1 + X 1 R ( ) 1 = P X ( ) ( R P ( ) P ( ) ( ) + P ( ) = P P ( ) ( ) (supposed that : instantane ous R ) P X ( ) P ( ) + P ( ) ) The time varying suppressio n fiter is given approimat ey by : H ( ) = R ( ) 1/ 2 pectrum Domain peech - Berin Chen 18

19 Wiener Fitering A peech Enhancement Technique From the tatistica Point of View The process [ m] is the sum of the random process s[ m] and the additive noise process n[ m] [ m] = s[ m] + n[ m] Find a inear estimate ŝ [ m] in terms of the process [ m] : Or to find a inear fiter h [ m] such that the sequence ŝ[ m] = [ m] h[ m] minimizes the epected vaue of ( ŝ[ m] s[ m] ) 2 [ m ] ŝ [ m ] oisy peech [ m ] s ˆ = [ m ] h [ m ] = A inear fiter h[n] = h [] [ m ] Cean peech peech - Berin Chen 19

20 Wiener Fitering (cont.) Minimize the epectation of the squared error (MME estimate) Minimize F k F h [ m] h[ ] [ m ] [] k s[ m] [ m k] = h[ ] [ m ] [ m k] k k = k = k = R s ss s s = 0 = E s [ m] [ m k] = h[ ] [ m ] [ m k] [ m] ( s[ m k] + n[ m k] ) = h[] [ m ] [ m k] s[ m] s[ m k] + s[ m] n[ m k] = h[ ] R[ k ] k = = [] k = h[] k R[] k Rs [ n] and R [ n] ( ) = H ( ) ( ) sequences of s n = = = k= = 2 k= and Take Fourier transform Take summation for k [] [ n] s [ m] n[ m] and are statisticay independent! : are respective y the autocorre ation peech - Berin Chen 20

21 Wiener Fitering (cont.) Minimize the epectation of the squared error (MME estimate) Q ss ( ) = H ( ) ( ) H ( ) = ss ( ) ( ) = ss ss ( ) ( ) + ( ) nn, is caed the noncausa Wiener fiter (where ( ) = ( ) + ( ) ss nn ) peech - Berin Chen 21

22 Wiener Fitering (cont.) The time varying Wiener Fiter aso can be epressed in a simiar form as the spectra subtraction H ( ) = ss = 1 + ( ) ss ( ) + ( ) P P nn P -1 ( ) 1 = 1 ( ) + R( ) = P ( ) ( ) + P ( ) -1, ( R ( ) = P P ( ) ( ) : instantane ous R ) vs. Wiener Fiter: 1. Wiener fiter has stronger attenuation at ow R region 2. Wiener fiter does not invoke an absoute threshoding 10 og P P ( ) ( ) peech - Berin Chen 22

23 Wiener Fitering (cont.) Wiener Fitering can be reaized ony if we know the power spectra of both the noise and the signa A chicken-and-egg probem Approach - I : Ephraim(1992) proposed the use of an HMM where, if we know the current frame fas under, we can use its mean spectrum as ( ) or P ( ) In practice, we do not know what state each frame fas into either Weight the fiters for each state by a posterior probabiity that frame fas into each state ss peech - Berin Chen 23

24 Approach - II : Wiener Fitering (cont.) The background/noise is stationary and its power spectrum can be estimated by averaging spectra over a known background region For the non-stationary speech signa, its time-varying power spectrum can be estimated using the past Wiener fiter (of previous frame) Pˆ H ( t, ) = PX ( t, ) H( t 1, ) Pˆ ( ) ( t, ) t, = Pˆ ( t, ) + P ( ), ( t : frameinde, H( ) : Wiener fiter) ~ P ( t, ) = P ( t, ) H( t, ) X The initia estimate of the speech spectrum can be derived from spectra subtraction ometimes introduce musica noise peech - Berin Chen 24

25 peech - Berin Chen 25 Wiener Fitering (cont.) Approach - III : ow down the rapid frame-to-frame movement of the object speech power spectrum estimate by appy tempora smoothing ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ),,,, ˆ, ˆ, in, ˆ repace to, use Then, ˆ 1 1, ~, α α P t P t P t H P t P t P t H t P t P t P t P t P + = + = + = ) ) ) )

26 Wiener Fitering (cont.) Cean peech oisy peech Enhanced oise peech Using Approach III α = 0.85 Other more compicate Wiener fiters peech - Berin Chen 26

27 The Effectives of Active oise peech - Berin Chen 27

28 Cepstra Mean ormaization (CM) A peech Enhancement Technique and sometimes caed Cepstra Mean ubtraction (CM) CM is a powerfu and simpe technique designed to hande convoutiona (Time-invariant inear fitering) distortions [ n] = s[ n] h[ n] ( ) ( ) H ( ) X = X = og H = og + og H = + H ( + H ) = C CH 1 T 1 = ( C t ) + CH = C CH CX = C + Time Domain pectra Domain Cepstra Domain 1 T 1 C = C t and CX + t = 0 t = 0 T T if the training and testing speech materias were recored from two different channes () 1 = C( + H(1) ) = C + CH(1), Testing : CX ( 2) = C( + H( 2 ) ) = C CH( 2 ) Training : CX + Log Power pectra Domain CX CX ( 1) CX ( 1) = C C ( 2 ) CX ( 2 ) = C C The spectra characteristics of the microphone and room acoustics thus can be removed! Can be eiminated if the assumption of zero-mean speech contribution! peech - Berin Chen 28

29 Cepstra Mean ormaization (cont.) ome Findings Interesting, CM has been found effective even the testing and training utterances are within the same microphone and environment Variations for the distance between the mouth and the microphone for different utterances and speakers Be carefu that the duration/period used to estimate the mean of noisy speech Why? Probematic when the acoustic feature vectors are amost identica within the seected time period peech - Berin Chen 29

30 Cepstra Mean ormaization (cont.) Performance For teephone recordings, where each ca has different frequency response, the use of CM has been shown to provide as much as 30 % reative decrease in error rate When a system is trained on one microphone and tested on another, CM can provide significant robustness peech - Berin Chen 30

31 Cepstra Mean ormaization (cont.) CM has been shown to improve the robustness not ony to varying channes but aso to the noise White noise added at different Rs ystem trained with speech with the same R (matched Condition) Cepstra deta and deta-deta features are computed prior to the CM operation so that they are unaffected. peech - Berin Chen 31

32 Cepstra Mean ormaization (cont.) From the other perspective We can interpret CM as the operation of subtracting a ow-pass tempora fiter d[ n], where a the T coefficients are identica and equa to 1, which is a high-pass tempora fiter T Aeviate the effect of conventiona noise introduced in the channe Tempora (Moduation) Frequency peech - Berin Chen 32

33 Cepstra Mean ormaization (cont.) Rea-time Cepstra ormaization CM requires the compete utterance to compute the cepstra mean; thus, it cannot be used in a rea-time system, and an approimation needs to be used Based on the above perspective, we can impement other types of high-pass fiters CX t = α CX t + α ( 1 ) CX t 1, ( CX t :cepstra mean) peech - Berin Chen 33

34 RATA Tempora Fiter Hyneck Hermansky, 1991 A peech Enhancement Technique RATA (Reative pectra) Assumption The inguistic message is coded into movements of the voca tract (i.e., the change of spectra characteristics) The rate of change of non-inguistic components in speech often ies outside the typica rate of change of the voca tact shape E.g. fi or sow time-varying inear communication channes A great sensitivity of human hearing to moduation frequencies around 4Hz than to ower or higher moduation frequencies Effect RATA uppresses the spectra components that change more sowy or quicky than the typica rate of change of speech peech - Berin Chen 34

35 RATA Tempora Fiter (cont.) The IIR transfer function 1 C ( ) ~ ( ) z z H z = = 0.1z C ( z) 1 MFCC stream 3 z 2z z c [ t] c ~ [ t] H(z) H(z) 4 ew MFCC stream Frame inde H(z) An other version z H ( z) = z z RATA has a peak at about 4Hz (moduation frequency) 2 z 4 c ~ [ t] = 0.98 c ~ [ t 1] c[ t] c[ t 1] 0.1 c[ t 3] c[ t 4] moduation frequency 100 Hz peech - Berin Chen 35

36 Retraining on Corrupted peech A Mode-based oise Compensation Technique Matched-Conditions Training Take a noise waveform from the new environment, add it to a the utterance in the training database, and retrain the system If the noise characteristics are known ahead of time, this method aow as to adapt the mode to the new environment with reativey sma amount of data from the new environment, yet use a arge amount of training data peech - Berin Chen 36

37 Retraining on Corrupted peech (cont.) Muti-stye Training Create a number of artificia acoustica environments by corrupting the cean training database with noise sampes of varying eves (30dB, 20dB, etc.) and types (white, babbe, etc.), as we as varying the channes A those waveforms (copies of training database) from mutipe acoustica environments can be used in training peech - Berin Chen 37

38 Mode Adaptation A Mode-based oise Compensation Technique The standard adaptation methods for speaker adaptation can be used for adapting speech recognizers to noisy environments MAP (Maimum a Posteriori) can offer resuts simiar to those of matched conditions, but it requires a significant amount of adaptation data MLLR (Maimum Likeihood Regression) can achieve reasonabe performance with about a minute of speech for minor mismatch. For severe mismatches, MLLR aso requires a arger amount of adaptation data peech - Berin Chen 38

39 igna Decomposition Using HMMs A Mode-based oise Compensation Technique Recognize concurrent signas (speech and noise) simutaneousy Parae HMMs are used to mode the concurrent signas and the composite signa is modeed as a function of their combined outputs Three-dimensiona Viterbi earch oise HMM (especiay for non-stationary noise) Computationay Epensive for both Training and Decoding! Cean speech HMM peech - Berin Chen 39

40 Parae Mode Combination (PMC) A Mode-based oise Compensation Technique By using the cean-speech modes and a noise mode, we can approimate the distributions obtained by training a HMM with corrupted speech peech - Berin Chen 40

41 Parae Mode Combination (cont.) The steps of tandard Parae Mode Combination (Log- orma Approimation) Cepstra domain μ Σ c c Σ μ Cean speech HMM s oisy speech HMM s μˆ c Σˆ c = C 1 μ c = C Σ ( C 1 c 1 ) T c μˆ = C μˆ c Σˆ = C Σˆ C Constraint: the estimate of variance is positive Log-spectra domain μ Σ μˆ Σˆ Σ ij μ = ep μ i = μ μ i ( + Σ 2) In inear spectra domain, the distribution is ognorma T j i ii [ ep( Σ ) 1] Because speech and noise are independent and additive in the inear spectra domain ij 1 Σˆ ii ( ˆ μi ) og( + 1) ˆ μi = og 2 ˆ μi 2 Σˆ ij = og ( ) Σˆ + 1 ˆ μ ˆ μ i ij j Σ Linear spectra domain μ Σ μ ˆ = gμ + μ~ ˆ 2 = g μˆ Σˆ Σ ~ + Σ oise HMM s μ~ ~ Σ Log-norma approimation (Assume the new distribution is ognorma) peech - Berin Chen 41

42 Parae Mode Combination (cont.) Modification-I: Perform the mode combination in the Log- pectra Domain (the simpest approimation) Log-Add Approimation: (without compensation of variances) ( ( ep μ ) ep( μ ) ˆ μ = og + ~ The variances are assumed to be sma A simpified version of Log-orma approimation Reduction in computationa oad Modification-II: Perform the mode combination in the Linear pectra Domain (Data-Driven PMC, DPMC, or Iterative PMC) Use the speech modes to generate noisy sampes (corrupted speech observations) and then compute a maimum ikeihood of these noisy sampes This method is ess computationay epensive than standard PMC with comparabe performance peech - Berin Chen 42

43 Parae Mode Combination (cont.) Modification-II: Perform the mode combination in the Linear pectra Domain (Data-Driven PMC, DPMC) Cean peech HMM oise HMM oisy peech HMM Cepstra domain Appy Monte Caro simuation to draw random cepstra vectors (for eampe, at east 100 for each distribution) Generating sampes Linear spectra domain Domain transform peech - Berin Chen 43

44 Parae Mode Combination (cont.) Data-Driven PMC peech - Berin Chen 44

45 Vector Tayor eries (VT) P. J. Moreno,1995 A Mode-based oise Compensation Technique VT Approach imiar to PMC, the noisy-speech-ike modes is generated by combining of cean speech HMM s and the noise HMM Unike PMC, the VT approach combines the parameters of cean speech HMM s and the noise HMM ineary in the ogspectra domain P = P P + P Power spectrum X X ( ) ( ) H ( ) ( ) = og( P ( ) PH ( ) + P ( ) ) P ( ) ( ) ( ) = og P PH 1 + P ( ) PH ( ) ogp ( ) ogp ( ) ogph ( ) = ogp ( ) + ogp ( ) + og1+ e = = + H + H + f H Log Power spectrum ( ) ( ) H + e on-inear function ( ) ( ) ( ) H, H,, wheref, H, = og1+ e + og1 Is a vector function peech - Berin Chen 45

46 peech - Berin Chen 46 Vector Tayor eries (cont.) The Tayor series provides a poynomia representation of a function in terms of the function and its derivatives at a point Appication often arises when noninear functions are empoyed and we desire to obtain a inear approimation The function is represented as an offset and a inear term ( ) ( ) ( )( ) ( )( ) ( ) ( )( ) ( ) n n n o f n f f f f R R f ! : =

47 Vector Tayor eries (cont.) Appy Tayor eries Approimation f df ( ) ( ) ( 0, H 0, 0 ),,,, ( H f 0 H d ) df ( 0, H 0, 0 ) ( ) df ( 0, H 0, 0 ) + + ( H H ) +... dh VT-0: use ony the 0th-order terms of Tayor eries VT-1: use ony the 0th- and 1th-order terms of Tayor eries f 0, H 0, 0 is the vector function evauated at a particuar vector point If VT-0 is used E u Σ [ X ] = E[ + H + f (,H, )] = u [ ( )] s + uh + E f,h, u + u + E[ f ( u,u,u )] s h s h n u + u + f ( u,u,u ) ( X s Σ s h + Σ ( ) h (if s h n and H To get the cean speech statistics is aso Gaussian) u Σ are independent) 0 u s Σ + g + s d If the channe fiter is inear - time invariant, 0-th order VT we can regard it as a bias (constant), g, in the og power spectrum domain f ( u,g,u ) ( X is aso Gaussian) s n 0 peech - Berin Chen 47

48 Vector Tayor eries (cont.) peech - Berin Chen 48

49 Retraining on Compensated Features A Mode-based oise Compensation Technique that aso Uses enhanced Features (processed by, CM, etc.) Combine speech enhancement and mode compensation peech - Berin Chen 49

50 Principa Component Anaysis Principa Component Anaysis (PCA) : Widey appied for the data anaysis and dimensionaity reduction in order to derive the most epressive feature Criterion: for a zero mean r.v. R, find k (k ) orthonorma vectors {e 1, e 2,, e k } so that (1) var(e 1 T )=ma 1 (2) var(e it )=ma i subject to e i e i-1 e 1 1 i k {e 1, e 2,, e k } are in fact the eigenvectors of the covariance matri (Σ ) for corresponding to the argest k eigenvaues Fina r.v y R k : the inear transform (projection) of the origina r.v., y=a T A=[e 1 e 2 e k ] data Principa ais peech - Berin Chen 50

51 Principa Component Anaysis (cont.) peech - Berin Chen 51

52 Principa Component Anaysis (cont.) Properties of PCA The components of y are mutuay uncorreated E{y i y j }=E{(e it ) (e jt ) T }=E{(e it ) ( T e j )}=e it E{ T } e j =e it Σ e j = λ j e it e j =0, if i j the covariance of y is diagona The error power (mean-squared error) between the origina vector and the projected is minimum =(e 1T )e 1 + (e 2T )e 2 + +(e kt )e k + +(e T )e =(e 1T )e 1 + (e 2T )e 2 + +(e kt )e k (ote : R ) error r.v : - = (e k+1t )e k+1 + (e k+2t )e k+2 + +(e T )e E((- ) T (- ))=E((e k+1t ) e k+1t e k+1 (e k+1t ))+ +E((e T ) e T e (e T )) =var(e k+1t )+ var(e k+2t )+ var(e T ) = λ k+1 + λ k+2 + +λ minimum peech - Berin Chen 52

53 PCA Appied in Inherenty Robust Features Appication 1 : the inear transform of the origina features (in the spatia domain) Origina feature stream t Frame inde z t = A T t A T A T A T A T The coumns of A are the first k eigenvectors of Σ transformed feature stream z t Frame inde peech - Berin Chen 53

54 peech - Berin Chen 54 PCA Appied in Inherenty Robust Features (cont.) Appication 2 : PCA-derived tempora fiter (in the tempora domain) The effect of the tempora fiter is equivaent to the weighted sum of sequence of a specific MFCC coefficient with ength L sid aong the frame inde quefrency Frame inde B 1 (z) B 2 (z) B n (z) L z k (1) z k (2) z k (3) Origina feature stream t The impuse response of B k (z) is one of the eigenvectors of the covariance for z k ( ) ( ) ( ) ( ) () ( ) () ( ) ( ) n m y m y m y m y K k K n k n n n K k K k K k K k ), ( ), (,2) (,1) ( ), ( ), (,2) (,1) ( ) (3, ) (3, (3,2) (3,1) ) (2, ) (2, (2,2) (2,1) ) (1, ) (1, (1,2) (1,1) 2 1 L L M M M M L M M L M M M M M M z k (n)=[ y k (n) y k (n+1) y k (n+2) y k (n+l-1)] T ( ) + = + = L n k n L μ zk z ( ) ( ) ( ) ( ) + = + = Σ L n T k k k k k n n L z z z z z μ μ ( ) ( ) ( ) n k n k T k z e 1, ˆ = The eement in the new feature vector From Dr. Jei-wei Hung

55 PCA Appied in Inherenty Robust Features (cont.) The frequency responses of the 15 PCA-derived tempora fiters From Dr. Jei-wei Hung peech - Berin Chen 55

56 PCA Appied in Inherenty Robust Features (cont.) Appication 2 : PCA-derived tempora fiter Mismatched condition Fiter ength L=10 Matched condition From Dr. Jei-wei Hung peech - Berin Chen 56

57 PCA Appied in Inherenty Robust Features (cont.) Appication 3 : PCA-derived fiter bank Power spectrum obtained by DFT h 1 h 2 h 3 h k is one of the eigenvectors of the covariance for k From Dr. Jei-wei Hung peech - Berin Chen 57

58 PCA Appied in Inherenty Robust Features (cont.) Appication 3 : PCA-derived fiter bank From Dr. Jei-wei Hung peech - Berin Chen 58

59 Linear Discriminative Anaysis Linear Discriminative Anaysis (LDA) Widey appied for the pattern cassification In order to derive the most discriminative feature Criterion : assume w j, μ j and Σ j are the weight, mean and covariance of cass j, j=1. Two matrices are defined as: Between - cass covariance : Within - cass covariance : Find W=[w 1 w 2 w k ] such that T W bw Wˆ = arg ma W T W W The coumns w j of W are the eigenvectors of w -1 B having the argest eigenvaues w w b = = j = 1 j = 1 w w j j Σ ( μ μ)( μ μ) j j j T peech - Berin Chen 59

60 Linear Discriminative Anaysis (cont.) The frequency responses of the 15 LDA-derived tempora fiters From Dr. Jei-wei Hung peech - Berin Chen 60

61 Minimum Cassification Error Minimum Cassification Error (MCE): Genera Objective : find an optima feature presentation or an optima recognition mode to minimize the epected error of cassification The recognizer is often operated under the foowing decision rue : C(X)=C i if g i (X,Λ)=ma j g j (X,Λ) Λ={λ (i) } i=1 M (M modes, casses), X : observations, g i (X,Λ): cass conditioned ikeihood function, for eampe, g i (X,Λ)=P(X λ (i) ) Traditiona Training Criterion : find λ (i) such that P(X λ (i) ) is maimum (Maimum Likeihood) if X C i This criterion does not aways ead to minimum cassification error, since it doesn't consider the mutua reationship between different casses For eampe, it s possibe that P(X λ (i) ) is maimum but X C i peech - Berin Chen 61

62 Minimum Cassification Error (cont.) Threshod τ P ( LR ( k ) KW k C k ) P LR ( k ) k ( ) KW k C k Type I error Type II error LR ( k ) Eampe showing histograms of the ikeihood ratio when keyword KW and KW k C k k C k LR ( k ) Type I error: Fase Rejection Type II error: Fase Aarm/Fase Acceptance peech - Berin Chen 62

63 Minimum Cassification Error (cont.) Minimum Cassification Error (MCE) (Cont.): One form of the cass miscassification measure : d d d i i i () i ( X ) g X, λ 1 M 1 ( ) ( ( () i + og ep g X, λ ) α ) = ( X ) 0 impies a miscassification (error = 1) ( X ) < 0 impies a correct cassification (error = 0) A continuous oss function is defined as foows : i ( X, Λ) = ( d ( X )) i X C where the sigmoid function ( d ) j i Cassifier performance measure : L i 1 = 1 + ep ( γ d + θ ) ( Λ) = E [ L( X, Λ) ] = ( X, Λ) δ( X C ) X X M i= 1 i i 1 α X C i peech - Berin Chen 63

64 Minimum Cassification Error (cont.) Using MCE in mode training : Find Λ such that Λ ˆ = arg min L Λ ( Λ) = arg min E [ L( X, Λ) ] Λ the above objective function in genera cannot be minimized directy but the oca minimum can be achieved using gradient decent agorithm ( ) L Λ wt+ = wt ε, w : an arbitrary parameter of w Using MCE in robust feature representation ˆ ( f ) f = arg mine L( f ( X ), Λ ) f ote: X 1 Λ f X [ ] : a transform of the originafeature X whie feature presentation is changed, the mode is aso changed accordingy peech - Berin Chen 64

Robustness Techniques for Speech Recognition

Robustness Techniques for Speech Recognition Robustness Techniques for peech Recognition Berin Chen, 2004 References: 1. X. Huang et a. poken Language Processing (2001). Chapter 10 2. J. C. Junqua and J. P. Haton. Robustness in Automatic peech Recognition

More information

Robustness Techniques for Speech Recognition

Robustness Techniques for Speech Recognition Robustness Techniques for peech Recognition Berin Chen Department of Computer cience & Information Engineering ationa Taiwan orma University References: 1. X. Huang et a. poken Language rocessing (2001).

More information

ASummaryofGaussianProcesses Coryn A.L. Bailer-Jones

ASummaryofGaussianProcesses Coryn A.L. Bailer-Jones ASummaryofGaussianProcesses Coryn A.L. Baier-Jones Cavendish Laboratory University of Cambridge caj@mrao.cam.ac.uk Introduction A genera prediction probem can be posed as foows. We consider that the variabe

More information

FRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA)

FRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA) 1 FRST 531 -- Mutivariate Statistics Mutivariate Discriminant Anaysis (MDA) Purpose: 1. To predict which group (Y) an observation beongs to based on the characteristics of p predictor (X) variabes, using

More information

A Brief Introduction to Markov Chains and Hidden Markov Models

A Brief Introduction to Markov Chains and Hidden Markov Models A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,

More information

A Robust Voice Activity Detection based on Noise Eigenspace Projection

A Robust Voice Activity Detection based on Noise Eigenspace Projection A Robust Voice Activity Detection based on Noise Eigenspace Projection Dongwen Ying 1, Yu Shi 2, Frank Soong 2, Jianwu Dang 1, and Xugang Lu 1 1 Japan Advanced Institute of Science and Technoogy, Nomi

More information

Bayesian Learning. You hear a which which could equally be Thanks or Tanks, which would you go with?

Bayesian Learning. You hear a which which could equally be Thanks or Tanks, which would you go with? Bayesian Learning A powerfu and growing approach in machine earning We use it in our own decision making a the time You hear a which which coud equay be Thanks or Tanks, which woud you go with? Combine

More information

Stochastic Variational Inference with Gradient Linearization

Stochastic Variational Inference with Gradient Linearization Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,

More information

Melodic contour estimation with B-spline models using a MDL criterion

Melodic contour estimation with B-spline models using a MDL criterion Meodic contour estimation with B-spine modes using a MDL criterion Damien Loive, Ney Barbot, Oivier Boeffard IRISA / University of Rennes 1 - ENSSAT 6 rue de Kerampont, B.P. 80518, F-305 Lannion Cedex

More information

Radar/ESM Tracking of Constant Velocity Target : Comparison of Batch (MLE) and EKF Performance

Radar/ESM Tracking of Constant Velocity Target : Comparison of Batch (MLE) and EKF Performance adar/ racing of Constant Veocity arget : Comparison of Batch (LE) and EKF Performance I. Leibowicz homson-csf Deteis/IISA La cef de Saint-Pierre 1 Bd Jean ouin 7885 Eancourt Cede France Isabee.Leibowicz

More information

An Algorithm for Pruning Redundant Modules in Min-Max Modular Network

An Algorithm for Pruning Redundant Modules in Min-Max Modular Network An Agorithm for Pruning Redundant Modues in Min-Max Moduar Network Hui-Cheng Lian and Bao-Liang Lu Department of Computer Science and Engineering, Shanghai Jiao Tong University 1954 Hua Shan Rd., Shanghai

More information

Target Location Estimation in Wireless Sensor Networks Using Binary Data

Target Location Estimation in Wireless Sensor Networks Using Binary Data Target Location stimation in Wireess Sensor Networks Using Binary Data Ruixin Niu and Pramod K. Varshney Department of ectrica ngineering and Computer Science Link Ha Syracuse University Syracuse, NY 344

More information

II. PROBLEM. A. Description. For the space of audio signals

II. PROBLEM. A. Description. For the space of audio signals CS229 - Fina Report Speech Recording based Language Recognition (Natura Language) Leopod Cambier - cambier; Matan Leibovich - matane; Cindy Orozco Bohorquez - orozcocc ABSTRACT We construct a rea time

More information

Bayesian Unscented Kalman Filter for State Estimation of Nonlinear and Non-Gaussian Systems

Bayesian Unscented Kalman Filter for State Estimation of Nonlinear and Non-Gaussian Systems Bayesian Unscented Kaman Fiter for State Estimation of Noninear and Non-aussian Systems Zhong Liu, Shing-Chow Chan, Ho-Chun Wu and iafei Wu Department of Eectrica and Eectronic Engineering, he University

More information

SVM: Terminology 1(6) SVM: Terminology 2(6)

SVM: Terminology 1(6) SVM: Terminology 2(6) Andrew Kusiak Inteigent Systems Laboratory 39 Seamans Center he University of Iowa Iowa City, IA 54-57 SVM he maxima margin cassifier is simiar to the perceptron: It aso assumes that the data points are

More information

Explicit overall risk minimization transductive bound

Explicit overall risk minimization transductive bound 1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

AST 418/518 Instrumentation and Statistics

AST 418/518 Instrumentation and Statistics AST 418/518 Instrumentation and Statistics Cass Website: http://ircamera.as.arizona.edu/astr_518 Cass Texts: Practica Statistics for Astronomers, J.V. Wa, and C.R. Jenkins, Second Edition. Measuring the

More information

Supervised i-vector Modeling - Theory and Applications

Supervised i-vector Modeling - Theory and Applications Supervised i-vector Modeing - Theory and Appications Shreyas Ramoji, Sriram Ganapathy Learning and Extraction of Acoustic Patterns LEAP) Lab, Eectrica Engineering, Indian Institute of Science, Bengauru,

More information

Turbo Codes. Coding and Communication Laboratory. Dept. of Electrical Engineering, National Chung Hsing University

Turbo Codes. Coding and Communication Laboratory. Dept. of Electrical Engineering, National Chung Hsing University Turbo Codes Coding and Communication Laboratory Dept. of Eectrica Engineering, Nationa Chung Hsing University Turbo codes 1 Chapter 12: Turbo Codes 1. Introduction 2. Turbo code encoder 3. Design of intereaver

More information

Improved Min-Sum Decoding of LDPC Codes Using 2-Dimensional Normalization

Improved Min-Sum Decoding of LDPC Codes Using 2-Dimensional Normalization Improved Min-Sum Decoding of LDPC Codes sing -Dimensiona Normaization Juntan Zhang and Marc Fossorier Department of Eectrica Engineering niversity of Hawaii at Manoa Honouu, HI 968 Emai: juntan, marc@spectra.eng.hawaii.edu

More information

(This is a sample cover image for this issue. The actual cover is not yet available at this time.)

(This is a sample cover image for this issue. The actual cover is not yet available at this time.) (This is a sampe cover image for this issue The actua cover is not yet avaiabe at this time) This artice appeared in a journa pubished by Esevier The attached copy is furnished to the author for interna

More information

Statistical Learning Theory: A Primer

Statistical Learning Theory: A Primer Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO

More information

A Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c)

A Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c) A Simpe Efficient Agorithm of 3-D Singe-Source Locaization with Uniform Cross Array Bing Xue a * Guangyou Fang b Yicai Ji c Key Laboratory of Eectromagnetic Radiation Sensing Technoogy, Institute of Eectronics,

More information

Problem set 6 The Perron Frobenius theorem.

Problem set 6 The Perron Frobenius theorem. Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator

More information

Proceedings of Meetings on Acoustics

Proceedings of Meetings on Acoustics Proceedings of Meetings on Acoustics Voume 9, 23 http://acousticasociety.org/ ICA 23 Montrea Montrea, Canada 2-7 June 23 Architectura Acoustics Session 4pAAa: Room Acoustics Computer Simuation II 4pAAa9.

More information

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.

More information

Nonlinear Analysis of Spatial Trusses

Nonlinear Analysis of Spatial Trusses Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes

More information

Alberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain

Alberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain CORRECTIONS TO CLASSICAL PROCEDURES FOR ESTIMATING THURSTONE S CASE V MODEL FOR RANKING DATA Aberto Maydeu Oivares Instituto de Empresa Marketing Dept. C/Maria de Moina -5 28006 Madrid Spain Aberto.Maydeu@ie.edu

More information

Automobile Prices in Market Equilibrium. Berry, Pakes and Levinsohn

Automobile Prices in Market Equilibrium. Berry, Pakes and Levinsohn Automobie Prices in Market Equiibrium Berry, Pakes and Levinsohn Empirica Anaysis of demand and suppy in a differentiated products market: equiibrium in the U.S. automobie market. Oigopoistic Differentiated

More information

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy

Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation

More information

Reliability: Theory & Applications No.3, September 2006

Reliability: Theory & Applications No.3, September 2006 REDUNDANCY AND RENEWAL OF SERVERS IN OPENED QUEUING NETWORKS G. Sh. Tsitsiashvii M.A. Osipova Vadivosto, Russia 1 An opened queuing networ with a redundancy and a renewa of servers is considered. To cacuate

More information

Translation Microscopy (TRAM) for super-resolution imaging.

Translation Microscopy (TRAM) for super-resolution imaging. ransation Microscopy (RAM) for super-resoution imaging. Zhen Qiu* 1,,3, Rhodri S Wison* 1,3, Yuewei Liu 1,3, Aison Dun 1,3, Rebecca S Saeeb 1,3, Dongsheng Liu 4, Coin Ricman 1,3, Rory R Duncan 1,3,5, Weiping

More information

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA

T.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network

More information

ESTIMATION OF SAMPLING TIME MISALIGNMENTS IN IFDMA UPLINK

ESTIMATION OF SAMPLING TIME MISALIGNMENTS IN IFDMA UPLINK ESTIMATION OF SAMPLING TIME MISALIGNMENTS IN IFDMA UPLINK Aexander Arkhipov, Michae Schne German Aerospace Center DLR) Institute of Communications and Navigation Oberpfaffenhofen, 8224 Wessing, Germany

More information

Convolutional Networks 2: Training, deep convolutional networks

Convolutional Networks 2: Training, deep convolutional networks Convoutiona Networks 2: Training, deep convoutiona networks Hakan Bien Machine Learning Practica MLP Lecture 8 30 October / 6 November 2018 MLP Lecture 8 / 30 October / 6 November 2018 Convoutiona Networks

More information

STA 216 Project: Spline Approach to Discrete Survival Analysis

STA 216 Project: Spline Approach to Discrete Survival Analysis : Spine Approach to Discrete Surviva Anaysis November 4, 005 1 Introduction Athough continuous surviva anaysis differs much from the discrete surviva anaysis, there is certain ink between the two modeing

More information

Nearest Neighbor Discriminant Analysis for Robust Speaker Recognition

Nearest Neighbor Discriminant Analysis for Robust Speaker Recognition INTERSPEECH 2014 Nearest Neighbor Discriminant Anaysis for Robust Speaker Recognition Seyed Omid Sadjadi, Jason W. Peecanos, Weizhong Zhu IBM Research, Watson Group {sadjadi,jwpeeca,zhuwe}@us.ibm.com Abstract

More information

Research of Data Fusion Method of Multi-Sensor Based on Correlation Coefficient of Confidence Distance

Research of Data Fusion Method of Multi-Sensor Based on Correlation Coefficient of Confidence Distance Send Orders for Reprints to reprints@benthamscience.ae 340 The Open Cybernetics & Systemics Journa, 015, 9, 340-344 Open Access Research of Data Fusion Method of Muti-Sensor Based on Correation Coefficient

More information

Discrete Techniques. Chapter Introduction

Discrete Techniques. Chapter Introduction Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, we as various

More information

Gokhan M. Guvensen, Member, IEEE, and Ender Ayanoglu, Fellow, IEEE. Abstract

Gokhan M. Guvensen, Member, IEEE, and Ender Ayanoglu, Fellow, IEEE. Abstract A Generaized Framework on Beamformer esign 1 and CSI Acquisition for Singe-Carrier Massive MIMO Systems in Miimeter Wave Channes Gokhan M. Guvensen, Member, IEEE, and Ender Ayanogu, Feow, IEEE arxiv:1607.01436v1

More information

SVM-based Supervised and Unsupervised Classification Schemes

SVM-based Supervised and Unsupervised Classification Schemes SVM-based Supervised and Unsupervised Cassification Schemes LUMINITA STATE University of Pitesti Facuty of Mathematics and Computer Science 1 Targu din Vae St., Pitesti 110040 ROMANIA state@cicknet.ro

More information

DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE

DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE Yury Iyushin and Anton Mokeev Saint-Petersburg Mining University, Vasiievsky Isand, 1 st ine, Saint-Petersburg,

More information

The EM Algorithm applied to determining new limit points of Mahler measures

The EM Algorithm applied to determining new limit points of Mahler measures Contro and Cybernetics vo. 39 (2010) No. 4 The EM Agorithm appied to determining new imit points of Maher measures by Souad E Otmani, Georges Rhin and Jean-Marc Sac-Épée Université Pau Veraine-Metz, LMAM,

More information

An Information Geometrical View of Stationary Subspace Analysis

An Information Geometrical View of Stationary Subspace Analysis An Information Geometrica View of Stationary Subspace Anaysis Motoaki Kawanabe, Wojciech Samek, Pau von Bünau, and Frank C. Meinecke Fraunhofer Institute FIRST, Kekuéstr. 7, 12489 Berin, Germany Berin

More information

Scalable Spectrum Allocation for Large Networks Based on Sparse Optimization

Scalable Spectrum Allocation for Large Networks Based on Sparse Optimization Scaabe Spectrum ocation for Large Networks ased on Sparse Optimization innan Zhuang Modem R&D Lab Samsung Semiconductor, Inc. San Diego, C Dongning Guo, Ermin Wei, and Michae L. Honig Department of Eectrica

More information

Approach to Identifying Raindrop Vibration Signal Detected by Optical Fiber

Approach to Identifying Raindrop Vibration Signal Detected by Optical Fiber Sensors & Transducers, o. 6, Issue, December 3, pp. 85-9 Sensors & Transducers 3 by IFSA http://www.sensorsporta.com Approach to Identifying Raindrop ibration Signa Detected by Optica Fiber ongquan QU,

More information

Methods for Ordinary Differential Equations. Jacob White

Methods for Ordinary Differential Equations. Jacob White Introduction to Simuation - Lecture 12 for Ordinary Differentia Equations Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Micha Rewienski, and Karen Veroy Outine Initia Vaue probem exampes Signa

More information

Fast Blind Recognition of Channel Codes

Fast Blind Recognition of Channel Codes Fast Bind Recognition of Channe Codes Reza Moosavi and Erik G. Larsson Linköping University Post Print N.B.: When citing this work, cite the origina artice. 213 IEEE. Persona use of this materia is permitted.

More information

Power Control and Transmission Scheduling for Network Utility Maximization in Wireless Networks

Power Control and Transmission Scheduling for Network Utility Maximization in Wireless Networks ower Contro and Transmission Scheduing for Network Utiity Maximization in Wireess Networks Min Cao, Vivek Raghunathan, Stephen Hany, Vinod Sharma and. R. Kumar Abstract We consider a joint power contro

More information

Asymptotic Gains of Generalized Selection Combining

Asymptotic Gains of Generalized Selection Combining Asymptotic Gains of Generaized Seection Combining Yao a Z Wang S Pasupathy Dept of ECpE Dept of ECE Dept of ECE Iowa State Univ Iowa State Univ Univ of Toronto, Emai: mayao@iastateedu Emai: zhengdao@iastateedu

More information

Active Learning & Experimental Design

Active Learning & Experimental Design Active Learning & Experimenta Design Danie Ting Heaviy modified, of course, by Lye Ungar Origina Sides by Barbara Engehardt and Aex Shyr Lye Ungar, University of Pennsyvania Motivation u Data coection

More information

Wavelet Methods for Time Series Analysis. Wavelet-Based Signal Estimation: I. Part VIII: Wavelet-Based Signal Extraction and Denoising

Wavelet Methods for Time Series Analysis. Wavelet-Based Signal Estimation: I. Part VIII: Wavelet-Based Signal Extraction and Denoising Waveet Methods for Time Series Anaysis Part VIII: Waveet-Based Signa Extraction and Denoising overview of key ideas behind waveet-based approach description of four basic modes for signa estimation discussion

More information

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia

More information

Discrete Techniques. Chapter Introduction

Discrete Techniques. Chapter Introduction Chapter 3 Discrete Techniques 3. Introduction In the previous two chapters we introduced Fourier transforms of continuous functions of the periodic and non-periodic (finite energy) type, as we as various

More information

Classification of Multi-Frequency Polarimetric SAR Images Based on Multi-Linear Subspace Learning of Tensor Objects

Classification of Multi-Frequency Polarimetric SAR Images Based on Multi-Linear Subspace Learning of Tensor Objects Remote Sens. 2015, 7, 9253-9268; doi:10.3390/rs70709253 OPEN ACCESS remote sensing ISSN 2072-4292 www.mdpi.com/journa/remotesensing Artice Cassification of Muti-Frequency Poarimetric SAR Images Based on

More information

Efficiently Generating Random Bits from Finite State Markov Chains

Efficiently Generating Random Bits from Finite State Markov Chains 1 Efficienty Generating Random Bits from Finite State Markov Chains Hongchao Zhou and Jehoshua Bruck, Feow, IEEE Abstract The probem of random number generation from an uncorreated random source (of unknown

More information

Optimality of Inference in Hierarchical Coding for Distributed Object-Based Representations

Optimality of Inference in Hierarchical Coding for Distributed Object-Based Representations Optimaity of Inference in Hierarchica Coding for Distributed Object-Based Representations Simon Brodeur, Jean Rouat NECOTIS, Département génie éectrique et génie informatique, Université de Sherbrooke,

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

A GENERAL METHOD FOR EVALUATING OUTAGE PROBABILITIES USING PADÉ APPROXIMATIONS

A GENERAL METHOD FOR EVALUATING OUTAGE PROBABILITIES USING PADÉ APPROXIMATIONS A GENERAL METHOD FOR EVALUATING OUTAGE PROBABILITIES USING PADÉ APPROXIMATIONS Jack W. Stokes, Microsoft Corporation One Microsoft Way, Redmond, WA 9852, jstokes@microsoft.com James A. Ritcey, University

More information

LINEAR DETECTORS FOR MULTI-USER MIMO SYSTEMS WITH CORRELATED SPATIAL DIVERSITY

LINEAR DETECTORS FOR MULTI-USER MIMO SYSTEMS WITH CORRELATED SPATIAL DIVERSITY LINEAR DETECTORS FOR MULTI-USER MIMO SYSTEMS WITH CORRELATED SPATIAL DIVERSITY Laura Cottateucci, Raf R. Müer, and Mérouane Debbah Ist. of Teecommunications Research Dep. of Eectronics and Teecommunications

More information

BDD-Based Analysis of Gapped q-gram Filters

BDD-Based Analysis of Gapped q-gram Filters BDD-Based Anaysis of Gapped q-gram Fiters Marc Fontaine, Stefan Burkhardt 2 and Juha Kärkkäinen 2 Max-Panck-Institut für Informatik Stuhsatzenhausweg 85, 6623 Saarbrücken, Germany e-mai: stburk@mpi-sb.mpg.de

More information

MC-CDMA CDMA Systems. Introduction. Ivan Cosovic. Stefan Kaiser. IEEE Communication Theory Workshop 2005 Park City, USA, June 15, 2005

MC-CDMA CDMA Systems. Introduction. Ivan Cosovic. Stefan Kaiser. IEEE Communication Theory Workshop 2005 Park City, USA, June 15, 2005 On the Adaptivity in Down- and Upink MC- Systems Ivan Cosovic German Aerospace Center (DLR) Institute of Comm. and Navigation Oberpfaffenhofen, Germany Stefan Kaiser DoCoMo Euro-Labs Wireess Soution Laboratory

More information

Title: Crack Diagnostics via Fourier Transform: Real and Imaginary Components vs. Power Spectral Density

Title: Crack Diagnostics via Fourier Transform: Real and Imaginary Components vs. Power Spectral Density Tite: rack Diagnostics via ourier Transform: ea and maginary omponents vs. Power pectra Density Authors: Leonid M. Geman Victor Giurgiutiu van V. Petrunin 3 rd nternationa Workshop on tructura Heath Monitoring

More information

Appendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model

Appendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model Appendix of the Paper The Roe of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Mode Caio Ameida cameida@fgv.br José Vicente jose.vaentim@bcb.gov.br June 008 1 Introduction In this

More information

Improving the Accuracy of Boolean Tomography by Exploiting Path Congestion Degrees

Improving the Accuracy of Boolean Tomography by Exploiting Path Congestion Degrees Improving the Accuracy of Booean Tomography by Expoiting Path Congestion Degrees Zhiyong Zhang, Gaoei Fei, Fucai Yu, Guangmin Hu Schoo of Communication and Information Engineering, University of Eectronic

More information

Do Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix

Do Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix VOL. NO. DO SCHOOLS MATTER FOR HIGH MATH ACHIEVEMENT? 43 Do Schoos Matter for High Math Achievement? Evidence from the American Mathematics Competitions Genn Eison and Ashey Swanson Onine Appendix Appendix

More information

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After

More information

Kernel pea and De-Noising in Feature Spaces

Kernel pea and De-Noising in Feature Spaces Kerne pea and De-Noising in Feature Spaces Sebastian Mika, Bernhard Schokopf, Aex Smoa Kaus-Robert Muer, Matthias Schoz, Gunnar Riitsch GMD FIRST, Rudower Chaussee 5, 12489 Berin, Germany {mika, bs, smoa,

More information

A proposed nonparametric mixture density estimation using B-spline functions

A proposed nonparametric mixture density estimation using B-spline functions A proposed nonparametric mixture density estimation using B-spine functions Atizez Hadrich a,b, Mourad Zribi a, Afif Masmoudi b a Laboratoire d Informatique Signa et Image de a Côte d Opae (LISIC-EA 4491),

More information

Cryptanalysis of PKP: A New Approach

Cryptanalysis of PKP: A New Approach Cryptanaysis of PKP: A New Approach Éiane Jaumes and Antoine Joux DCSSI 18, rue du Dr. Zamenhoff F-92131 Issy-es-Mx Cedex France eiane.jaumes@wanadoo.fr Antoine.Joux@ens.fr Abstract. Quite recenty, in

More information

Unconditional security of differential phase shift quantum key distribution

Unconditional security of differential phase shift quantum key distribution Unconditiona security of differentia phase shift quantum key distribution Kai Wen, Yoshihisa Yamamoto Ginzton Lab and Dept of Eectrica Engineering Stanford University Basic idea of DPS-QKD Protoco. Aice

More information

Strain Energy in Linear Elastic Solids

Strain Energy in Linear Elastic Solids Strain Energ in Linear Eastic Soids CEE L. Uncertaint, Design, and Optimiation Department of Civi and Environmenta Engineering Duke Universit Henri P. Gavin Spring, 5 Consider a force, F i, appied gradua

More information

C. Fourier Sine Series Overview

C. Fourier Sine Series Overview 12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a

More information

A Generalized Framework on Beamformer Design and CSI Acquisition for Single-Carrier Massive MIMO Systems in Millimeter-Wave Channels

A Generalized Framework on Beamformer Design and CSI Acquisition for Single-Carrier Massive MIMO Systems in Millimeter-Wave Channels 1 A Generaized Framework on Beamformer esign and CSI Acquisition for Singe-Carrier Massive MIMO Systems in Miimeter-Wave Channes Gokhan M. Guvensen, Member, IEEE and Ender Ayanogu, Feow, IEEE Abstract

More information

EXEMPLAR-BASED VOICE CONVERSION USING NON-NEGATIVE SPECTROGRAM DECONVOLUTION

EXEMPLAR-BASED VOICE CONVERSION USING NON-NEGATIVE SPECTROGRAM DECONVOLUTION 8th ISCA Speech Synthesis Worshop August 31 September 2, 2013 Barceona, Spain EXEMPLAR-BASED VOICE CONVERSION USING NON-NEGATIVE SPECTROGRAM DECONVOLUTION Zhizheng Wu 1,2, Tuomas Virtanen 3, Tomi Kinnunen

More information

A. Distribution of the test statistic

A. Distribution of the test statistic A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch

More information

6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7

6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7 6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17 Soution 7 Probem 1: Generating Random Variabes Each part of this probem requires impementation in MATLAB. For the

More information

Traffic data collection

Traffic data collection Chapter 32 Traffic data coection 32.1 Overview Unike many other discipines of the engineering, the situations that are interesting to a traffic engineer cannot be reproduced in a aboratory. Even if road

More information

Soft Clustering on Graphs

Soft Clustering on Graphs Soft Custering on Graphs Kai Yu 1, Shipeng Yu 2, Voker Tresp 1 1 Siemens AG, Corporate Technoogy 2 Institute for Computer Science, University of Munich kai.yu@siemens.com, voker.tresp@siemens.com spyu@dbs.informatik.uni-muenchen.de

More information

General Tensor Discriminant Analysis and Gabor Features for Gait Recognition

General Tensor Discriminant Analysis and Gabor Features for Gait Recognition > PAI-057-005< Genera ensor Discriminant Anaysis and Gabor Features for Gait Recognition Dacheng ao, Xueong Li, Xindong Wu 2, and Stephen J. aybank. Schoo of Computer Science and Information Systems, Birkbeck

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

Discriminant Analysis: A Unified Approach

Discriminant Analysis: A Unified Approach Discriminant Anaysis: A Unified Approach Peng Zhang & Jing Peng Tuane University Eectrica Engineering & Computer Science Department New Oreans, LA 708 {zhangp,jp}@eecs.tuane.edu Norbert Riede Tuane University

More information

Structural Control of Probabilistic Boolean Networks and Its Application to Design of Real-Time Pricing Systems

Structural Control of Probabilistic Boolean Networks and Its Application to Design of Real-Time Pricing Systems Preprints of the 9th Word Congress The Internationa Federation of Automatic Contro Structura Contro of Probabiistic Booean Networks and Its Appication to Design of Rea-Time Pricing Systems Koichi Kobayashi

More information

Primal and dual active-set methods for convex quadratic programming

Primal and dual active-set methods for convex quadratic programming Math. Program., Ser. A 216) 159:469 58 DOI 1.17/s117-15-966-2 FULL LENGTH PAPER Prima and dua active-set methods for convex quadratic programming Anders Forsgren 1 Phiip E. Gi 2 Eizabeth Wong 2 Received:

More information

Lecture Note 3: Stationary Iterative Methods

Lecture Note 3: Stationary Iterative Methods MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or

More information

Partial permutation decoding for MacDonald codes

Partial permutation decoding for MacDonald codes Partia permutation decoding for MacDonad codes J.D. Key Department of Mathematics and Appied Mathematics University of the Western Cape 7535 Bevie, South Africa P. Seneviratne Department of Mathematics

More information

A Novel Learning Method for Elman Neural Network Using Local Search

A Novel Learning Method for Elman Neural Network Using Local Search Neura Information Processing Letters and Reviews Vo. 11, No. 8, August 2007 LETTER A Nove Learning Method for Eman Neura Networ Using Loca Search Facuty of Engineering, Toyama University, Gofuu 3190 Toyama

More information

The influence of temperature of photovoltaic modules on performance of solar power plant

The influence of temperature of photovoltaic modules on performance of solar power plant IOSR Journa of Engineering (IOSRJEN) ISSN (e): 2250-3021, ISSN (p): 2278-8719 Vo. 05, Issue 04 (Apri. 2015), V1 PP 09-15 www.iosrjen.org The infuence of temperature of photovotaic modues on performance

More information

From Margins to Probabilities in Multiclass Learning Problems

From Margins to Probabilities in Multiclass Learning Problems From Margins to Probabiities in Muticass Learning Probems Andrea Passerini and Massimiiano Ponti 2 and Paoo Frasconi 3 Abstract. We study the probem of muticass cassification within the framework of error

More information

Support Vector Machine and Its Application to Regression and Classification

Support Vector Machine and Its Application to Regression and Classification BearWorks Institutiona Repository MSU Graduate Theses Spring 2017 Support Vector Machine and Its Appication to Regression and Cassification Xiaotong Hu As with any inteectua project, the content and views

More information

Hierarchical Fault Response Modeling of Analog/RF Circuits. Gurusubrahmaniyan Subrahmaniyan Radhakrishnan

Hierarchical Fault Response Modeling of Analog/RF Circuits. Gurusubrahmaniyan Subrahmaniyan Radhakrishnan Hierarchica Faut Response Modeing of Anaog/RF Circuits by Gurusubrahmaniyan Subrahmaniyan Radhakrishnan A Thesis Presented in Partia Fufiment of the Requirements for the Degree Master of Science Approved

More information

Source and Relay Matrices Optimization for Multiuser Multi-Hop MIMO Relay Systems

Source and Relay Matrices Optimization for Multiuser Multi-Hop MIMO Relay Systems Source and Reay Matrices Optimization for Mutiuser Muti-Hop MIMO Reay Systems Yue Rong Department of Eectrica and Computer Engineering, Curtin University, Bentey, WA 6102, Austraia Abstract In this paper,

More information

Proceedings of Meetings on Acoustics

Proceedings of Meetings on Acoustics Proceedings of Meetings on Acoustics Voume 19, 213 http://acousticasociety.org/ ICA 213 Montrea Montrea, Canada 2-7 June 213 Underwater Acoustics Session 1pUW: Seabed Scattering: Measurements and Mechanisms

More information

Statistics for Applications. Chapter 7: Regression 1/43

Statistics for Applications. Chapter 7: Regression 1/43 Statistics for Appications Chapter 7: Regression 1/43 Heuristics of the inear regression (1) Consider a coud of i.i.d. random points (X i,y i ),i =1,...,n : 2/43 Heuristics of the inear regression (2)

More information

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z Chapter W3 Mechanica Systems II Introduction This companion website chapter anayzes the foowing topics in connection to the printed-book Chapter 3: Lumped-parameter inertia fractions of basic compiant

More information

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017 In-pane shear stiffness of bare stee deck through she finite eement modes G. Bian, B.W. Schafer June 7 COLD-FORMED STEEL RESEARCH CONSORTIUM REPORT SERIES CFSRC R-7- SDII Stee Diaphragm Innovation Initiative

More information

arxiv: v1 [cs.db] 1 Aug 2012

arxiv: v1 [cs.db] 1 Aug 2012 Functiona Mechanism: Regression Anaysis under Differentia Privacy arxiv:208.029v [cs.db] Aug 202 Jun Zhang Zhenjie Zhang 2 Xiaokui Xiao Yin Yang 2 Marianne Winsett 2,3 ABSTRACT Schoo of Computer Engineering

More information

hole h vs. e configurations: l l for N > 2 l + 1 J = H as example of localization, delocalization, tunneling ikx k

hole h vs. e configurations: l l for N > 2 l + 1 J = H as example of localization, delocalization, tunneling ikx k Infinite 1-D Lattice CTDL, pages 1156-1168 37-1 LAST TIME: ( ) ( ) + N + 1 N hoe h vs. e configurations: for N > + 1 e rij unchanged ζ( NLS) ζ( NLS) [ ζn unchanged ] Hund s 3rd Rue (Lowest L - S term of

More information

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS TONY ALLEN, EMILY GEBHARDT, AND ADAM KLUBALL 3 ADVISOR: DR. TIFFANY KOLBA 4 Abstract. The phenomenon of noise-induced stabiization occurs

More information