estimation of autocorrelation matrix or eigenvalue spread.

Started by mk in comp.dsp15 years ago 12 replies

Hi everyone, I am experimenting with some equalization schemes and to judge the performance, I need to estimate the eigenvalue spread of the...

Hi everyone, I am experimenting with some equalization schemes and to judge the performance, I need to estimate the eigenvalue spread of the input to the equalizer. Currently I am working with a single input stream of many thousands of samples, and I am assuming that my data is more or less uniformly distributed. To get the eigenvalues, I generate R = Sum( x[k]*x[k]T) where x[k] is a vector...


Pitch recognition of a musical note on a smart phone

Started by mzincali in comp.dsp12 years ago 6 replies

Hi. With limited resources such as slower CPUs, code size and RAM, how best to detect the pitch of a musical note, similar to what a tuner...

Hi. With limited resources such as slower CPUs, code size and RAM, how best to detect the pitch of a musical note, similar to what a tuner would do? Should I use: - Kiss FFT - FFTW - autocorrelation - zero crossing analysis - octave-spaced filters other? Thank you in advance. -mz


Levinson algo. wtih complex coef. required

Started by herbert arthur in comp.dsp18 years ago 7 replies

Hi, I am desperatly looking for a generalized Levinson algorithm that can cope with complex autocorrelation values and giving complex filter...

Hi, I am desperatly looking for a generalized Levinson algorithm that can cope with complex autocorrelation values and giving complex filter coefficients... I would like to implement it in C, and i have only found the classical LPC analysis but it can not deal with complex coefficients... Anyone can help ? Thanks in advance


error caused by Wien-Hammerstein system (Bussgang?)

Started by Peter Mairhofer in comp.dsp5 years ago

Hi, I have the following setup: x ---> [ G(z) ]-----> [ (.)^3 ]-----> [ H(z) ]-----> y where x is a Gaussian input with E(x)=0 and...

Hi, I have the following setup: x ---> [ G(z) ]-----> [ (.)^3 ]-----> [ H(z) ]-----> y where x is a Gaussian input with E(x)=0 and autocorrelation Rxx. For simplicity, I can set Rxx=sigma^2*delta(tau) and G(z)=1, i.e. a pure Hammerstein system driven with white Gaussian noise (I would like to generalize afterwards if possible): x ---> [ (.)^3 ]-----> [ H(z) ]-----> y Now I would lik


Noise estimation

Started by Doug in comp.dsp11 years ago 10 replies

Hello, Is there a method to estimate the statistical properties of additive noise from a received signal observation. Assume we have y(k)...

Hello, Is there a method to estimate the statistical properties of additive noise from a received signal observation. Assume we have y(k) = x(k) + n(k) and we observe the received signal y(k). I would like to know if there is a way to determine the autocorrelation of the noise process n(k), that is Rnn. This is not a speech problem. x(k) is a digital PAM signal that is corrupted ...


ARMA model statistical properties

Started by Thomas Arildsen in comp.dsp14 years ago 8 replies

I have an ARMA process with known parameters. Are there general expressions for variance and possibly autocorrelation of the output for such...

I have an ARMA process with known parameters. Are there general expressions for variance and possibly autocorrelation of the output for such processes? I can find the above for an order-(1,1) process in for example Shanmugan & Breipohl's "Random Signals" (Wiley, 1988), but is it possible for higher orders? If it helps, I would be happy to confine the processess in question to order (p,p...


Detection of frequency with minimum latency

Started by devidasbhonde in comp.dsp15 years ago 16 replies

Hi , Rune said in thread of "Basics of Spectral Density and Autocorrelation" that it is possible to detect a single sine wave in a data...

Hi , Rune said in thread of "Basics of Spectral Density and Autocorrelation" that it is possible to detect a single sine wave in a data stream. How to do it ? I have a low frequency time series data (7Hz to 50Hz) I need to get its sprectum with minimum latency. and with 1Hz resolution. What should I do ? Regards and thanks, Devidas


Simple Question on Random Variables and Random Processes

Started by Randy Yates in comp.dsp16 years ago 16 replies

Let X(t) be a stationary white noise process with autocorrelation function R_{xx}(\tau) = N0 * \delta(\tau). Let V be a random...

Let X(t) be a stationary white noise process with autocorrelation function R_{xx}(\tau) = N0 * \delta(\tau). Let V be a random variable determined from X(t) by V = \int_{-\infty}^{+\infty} X(t) dt. How do you compute statistics on V, like E[V^2]? Yes, this is a homework problem. -- % Randy Yates % "Watching all the days go by... %% Fuquay-Varina, NC ...


LPC using cumulants

Started by shamailsaeed in comp.dsp13 years ago 2 replies

Hi All, I want to implement LPC using cumulants (instead of LPC using autocorrelation ) algorithms in C. Where can I find that algorithms. some...

Hi All, I want to implement LPC using cumulants (instead of LPC using autocorrelation ) algorithms in C. Where can I find that algorithms. some formulas or c code.


LMS maximum step size

Started by cfy30 in comp.dsp11 years ago 13 replies

Hi all, w(n+1) = w(n) + mu*e(n)*x(n) Can anyone tell me how to determine the maximum mu? I read some materials that mu should be less than...

Hi all, w(n+1) = w(n) + mu*e(n)*x(n) Can anyone tell me how to determine the maximum mu? I read some materials that mu should be less than 2/lambda_max, where lambda_max is the maximum eigen value of the autocorrelation matrix of the input. Take the 50Hz noise cancellation filter in Widrow's paper as an example, "Adaptive Noise Cancelling : Principles and Applications". How to determine t...


GSM AMR optimization

Started by Yaroslav in comp.dsp16 years ago 2 replies

Hi, I am optimizing GSM AMR reference code. The optimization needs not to be bit-exact, in fact, some quality degradation is permitted. I...

Hi, I am optimizing GSM AMR reference code. The optimization needs not to be bit-exact, in fact, some quality degradation is permitted. I have done all the obvious things like primitives inlining, autocorrelation approximation, etc. Are there any other mathematical shortcuts to make the codec faster? Right now, the codec is 25% slower that I need it to be. TIA, ~Yaroslav


Scaling factor in inverse fourier transform

Started by Titi Anggono in comp.dsp14 years ago 3 replies

Hi all, I tried to calculate the autocorrelation from FFT using C. As far as I know the input for the inverse FT is the power spectrum density...

Hi all, I tried to calculate the autocorrelation from FFT using C. As far as I know the input for the inverse FT is the power spectrum density of the data X(t) PSD = (Xreal(f)^2 + Ximag(f)^2) / N, N is the length of the real number After calculating the inverse FT, I scaled the result with 1/N but the result is not 1 (less than 1). I wonder if I made mistake with the scaling factor or ...


algorithms for location versus for frequency?

Started by Ron N. in comp.dsp16 years ago 26 replies

I suppose all characterization of signals inside some buffer of samples is a form of statistical signal processing. When asking about the...

I suppose all characterization of signals inside some buffer of samples is a form of statistical signal processing. When asking about the frequency of some signal in comp.dsp, all sorts of algorithms are mentioned, forms of autocorrelation (or squared difference, etc.), or fft/dft usage (windowed, zero padded, phase differentiated, etc.). Also, there seem to be lots of textbooks on these ...


Plot of Autocorrelation Error using LMS in Matlab

Started by Anonymous in comp.dsp15 years ago

Hi !! I hope somebody can help !! N = 8000; % Number of Samples n = 0:N-1; sigmax = 0.01 L = 12 f1 = 101.7 f2 = 142.4 f3 =...

Hi !! I hope somebody can help !! N = 8000; % Number of Samples n = 0:N-1; sigmax = 0.01 L = 12 f1 = 101.7 f2 = 142.4 f3 = 231.5 fs = 1000 K = 20 % Number of realizations cumensq = zeros(K,N) en1 = zeros(K,N) % length of filter coefficients wn = zeros(1,L); a = 0.0001; muhat = 1.00 wmat = zeros(N+1,L); % 8001 filters wmat(1,:)=wn; for kk = 1: K % ensemble process xn =...


Quadriphase barker sequences

Started by Anonymous in comp.dsp16 years ago 2 replies

Hello, I am looking for a quadriphase Barker sequence (a sequence whose elements are from the set {1,j,-1,-j} and which has a...

Hello, I am looking for a quadriphase Barker sequence (a sequence whose elements are from the set {1,j,-1,-j} and which has a periodic autocorrelation function whose maximum sidelobe magnitude


should first coef of spiking decon Wiener filter be positive?

Started by Anonymous in comp.dsp5 years ago

Spiking decon has been used in geophysics for over 50 years. An FIR Wiener filter is computed by Levinson recursion. The inputs are the...

Spiking decon has been used in geophysics for over 50 years. An FIR Wiener filter is computed by Levinson recursion. The inputs are the autocorrelation of the trace, and an impulse (1,0,0,0,...) as desired signal. It is based on certain assumptions, such as the wavelet being minimum phase and stationary (which are only approximations in reality). Now I have been looking at an open-source pac...


Different Definitions for Inverse Fourier Transform

Started by Randy Yates in comp.dsp15 years ago 16 replies

It's not completely standard, but I've usually seen a factor of 1/(2*\pi) in front of the integral for the inverse Fourier transform. I believe...

It's not completely standard, but I've usually seen a factor of 1/(2*\pi) in front of the integral for the inverse Fourier transform. I believe the theore says that this factor must be somewhere in the round-trip journey from forward transform to revers. It could be 1/sqrt(2*\pi) in front of both, or whatever. [papoulis] has the present in his definition of the autocorrelation function as ...


Sine Wave autocorrelation, interpolation of phase

Started by eduardoG26 in comp.dsp4 years ago 11 replies

Hallo. For the purpose of measuring complex impedances i need to compare the phase of two copies of a sinewave over a number of periods. The...

Hallo. For the purpose of measuring complex impedances i need to compare the phase of two copies of a sinewave over a number of periods. The sinewave is generated in the same ┬ÁC. Frequency is known and stable. The original and the shifted signals are sampled by a double synchr. ADC. The phase shift is the base for calculation of the complex Z of a load. I have implemented an (auto-)correlatio...


each PMEPR has a distinct waveform?

Started by hype...@gmail.com in comp.dsp14 years ago

hi all, while i am reading abt Peak Power reduction from Hanzo book pg 336. he described that a AA ( Aperiodic Autocorrelation) and its...

hi all, while i am reading abt Peak Power reduction from Hanzo book pg 336. he described that a AA ( Aperiodic Autocorrelation) and its transfroms ( actually circular shifted version of original waveform ) completely characterize PMEPR. so if one AA and its transforms have one PMEPR and one distinct waveform then can i conclude that the relation bet PMEPR and Distinct Waveform is one...


IIR Wiener Filter in Matlab

Started by SlvrDreams in comp.dsp15 years ago 1 reply

I am attempting to code an IIR wiener filter for system identification in Matlab, but I am currently having some difficulty with the theory......

I am attempting to code an IIR wiener filter for system identification in Matlab, but I am currently having some difficulty with the theory... Can someone please clarify? To find the spectral factorization, my understanding is that I take the roots of the autocorrelation sequence, then select those whose magnitude fall within the unit circle? Does anyone have any example IIR Wiener Matlab sc...