search for: a_j

...{ double sum = 0; for( int k = 0; k < 2*N; k++ ) { sum += in[k]*cos(M_PI*(2*j+1+N)*(2*k+1)/8/N); } out[j] = sum; } } > It just a normal shifted mdct, like the one used in mpeg layer 3. > it's equations are: > (time domain: a_j, freq domain: b_k; w_j is the window) > forward: > b_k = sum(j=0..4N-1) a_j * w_j * cos( (2k+1)(2j+1+-N) / 8N ) // forward > a_j = w_j * sum(k=0..2N-1) a_j * cos( (2k+1)(2j+1+-N) / 8N ) // backward (missing PI in the cos, and a_j should be b_k in summation) > where +- means: I'm n...

...ing else from a text terminal with gcc and vi right now). Hope you like it, Segher -- Suppose we want to calculate a MDCT (encode side): (x_ is input (time) samples, w_ is window, b_ is output (freq) samples): b_k = \sum_{j=0}^{4N-1} x_j w_j \cos{(2k+1)(2j+1-2N) \over 8N} \pi After setting a_j := w_{-j-1+N} x{-j-1+N} + w_{j+N} x{j+N}, j < N a_j := w_{j-1-N} x{j-1-N} - w_{j+3N} x{j+3N}, otherwise (or something similar; the cos'es should agree) we get b_k = \sum_{j=0}^{2N-1} a_j \cos{(2k+1)(2j+1) \over 8N} \pi This is a DCT-IV. For decode side, we first compute DCT-IV, and...

Hi, Can someone tell me which MDCT and invMDCT equation uses? I implemented the invMDCT one given in eusipco.corrected.ps file (handed out by Monty way back) and it produces different time domain samples. I tried both the FFT method and the slow way directly from the equation and couldn't reproduce the results from the original code. This leads me to believe that the forward MDCT used in

I should'nt have sent the last mail so fast. Same problem with \eqn{u_j = a_j + b x + c x^2, \quad j=1, \ldots, r-1} {u[j] = a[j] + b*x + c*x^2 j = 1,\dots,r-1} I thought the problem in the first case could have to do with the use use of \mbox{} (with the braces) within the arguments of \eqn, but here there are none braces in the arguments of \eqn{}{}. Another c...

suppose something like probability(passing test) is driven by 1. fixed effects -- sex 2. district effects - district funding 3. school effects - neighborhood income, racial composition, % two parent families, ... 4. class effects - teacher quality measurement, 5. individual random effects - IQ. how would such a model be setup in lmer? I can't find much discussion on the

...t the case of a single group. The model is: Y_i= a +bX_i + error where I indexes the different values of X and Y in this group . The a priori hypothesis of the slope is: b=K. This is easily tested with a t-test (b-K=0). Now imagine that there are j groups. For each group j the model is: Y_ij= a_j + b_jX_ij + error. Both the intercepts (a) and the slopes (b) are allowed to vary between groups. The a priori (null) hypothesis of interest involved the between-group values of the slopes and is: b_j=Kj where Kj is specified a priori for each group j based on theoretical considerations but whos...

Hi all, I have a problem using rdirichlet{gtools}. For Dir( a1, a2, ..., a_n), its mode can be found at $( a_i -1)/ ( \sum_{i}a_i - n)$; The means are $a_i / (\sum_{i} a_i ) $; I tried to study the above properties using rdirichlet from gtools. The code are: ############## library(gtools) alpha = c(1,3,9) #totoal=13 mean.expect = c(1/13, 3/13, 9/13) mode.expect = c(0, 2/10, 8/10) #

...+7(M-BUI > # M+6$M=RUE+68M9RUQ+34M-BTW+6DM82TS=0HC($!H5TER/3MX=GA1;&]Z>7!C > # M9G%B?FI'/58J4W!;,G)R4CA6<6%@7#U*,FMF67E8,VI'8G(]/CM[=7IX='AX > # MW<L,14(2SWS1$J0=[8?[[?=T-T!2LK,S>W,S5W4;[TXLD4"CT:]/-^"JC)-> > # M$F?5E]ZP_WJ^^^0W^-$'@Y>'A_J)UOKET<':;_<ST/KHZ&"P]_+@Y>&!UGOX > # M=/1$'S[Y'7XJ6P:%UD_^27^J#?U'L;WMT4/[OV<*_XC^#UG_P^'1P3?]_Q_K > # M_SRX-;,X,;]ZC;W!X.C@X*OZ/]C;/R#]#UX.#_?V!["%O8/!X.43/?BF_]_\ > # MYYGV:M:]7O-YEA>Z,0U]QO_2@[B,C57/]'%A@M)$(ST)2OW6A'IXI(>#...