Analyzing a Pulse: Averaging Spectra & Data Window Criteria

mediocre
Messages
14
Reaction score
0
Is this method good for analyzing a pulse?
Why is the averaging of the spectra good for noise cleaning?
What's the criteria for defining how much windows in a set of data one should take?

Anyone who clarifies this has a big thanks from me.
 
Mathematics news on Phys.org
Welch's method divides the data string into short segments, applies a window to each prior to the Fourier transform, averages the transforms, and takes the squared magnitude to produce a power spectrum. It reduces the variance of the DFT spectral estimate (just as Bartlett's method does) and at the same time smooths the estimate in the way windowing always does, at the expense of degraded resolution (also in the way windowing always does). Another advantage over the regular Bartlett periodogram is that the power spectrum is always positive. Most spectrum analyzers use this method.

Answering your other questions requires a deeper understanding of spectral analysis. There are a dozen well-established methods of spectral analysis. DFT-based ones are the simplest and most common but not necessarily best. Model-based methods (AR, MA, ARMA, etc.) are better if you have a model and the time to develop the specific analysis. Non-linear methods (Capon method or minimum variance distortionless response, maximum likelihood, MUSIC, etc.) have far better resolution but are complicated and require a high state of knowledge to use wisely. Wavelets and so-called fractional Fourier transform methods are good for signals with components of varying finite durations (speech, say). Maximum entropy and Bayesian methods are at the top in terms of accuracy, resolution, freedom from artifacts, and also in terms of difficulty and required knowledge. They are seldom used.

The wording of your questions suggests that you are not knowledgeable in these methods. Therefore, go ahead and use a simple FFT-based method. You have plenty of company--probably 90+% of engineers know of and use nothing else.
 
Well,i didn't expect this much information :) ... but of course,your observation is very much correct.
I had some assignments for class and i just wanted to know,on a whim,why these methods are suited for my task.
MATLAB help points to the Introduction to Spectral analysis from Stoica and Moises,ill guess ill start there and then perhaps another time be more specific...But thanks anyway!
 
Sure thing.

I can answer your last question in a general way, anyway. If your signal is repetitive and low in jitter (example: output of a radar transmitter), then averaging a lot of data is helpful. Overlap the windows so that all segments receive equal waiting on average. If the signal is variable (maybe receiving UHF radar pulses through the ionosphere during high solar activity), then you'll get a lot of smearing. This might be good if you want to characterize the average propagation channel. If you want to determine the spectrum of an individual pulse, on the other hand, the smearing is bad. Instead use a short data record.

For a given data sequence, the trade between many short segments and fewer longer ones is the trade between lower variance in the spectral estimate, and finer spectral resolution.
 
I'll be honest and say that i am not following everything you've said but at least you gave some pointers.

And to be direct my problem was determining the amplitude response for a horizontal seimograph.
In the picture on the left is the one obtained by forcing it with a current(via induction).While the right one is given by analyzing the spectral density of a pulse and then multiplying it with phase velocity squared to get the amplitude response.

My problem is the difference that gets bigger with higher frequencies.

Anyway,thanks for the effort,its much appreciated
 

Attachments

  • praktikum.png
    praktikum.png
    12.9 KB · Views: 589
Insights auto threads is broken atm, so I'm manually creating these for new Insight articles. In Dirac’s Principles of Quantum Mechanics published in 1930 he introduced a “convenient notation” he referred to as a “delta function” which he treated as a continuum analog to the discrete Kronecker delta. The Kronecker delta is simply the indexed components of the identity operator in matrix algebra Source: https://www.physicsforums.com/insights/what-exactly-is-diracs-delta-function/ by...
Fermat's Last Theorem has long been one of the most famous mathematical problems, and is now one of the most famous theorems. It simply states that the equation $$ a^n+b^n=c^n $$ has no solutions with positive integers if ##n>2.## It was named after Pierre de Fermat (1607-1665). The problem itself stems from the book Arithmetica by Diophantus of Alexandria. It gained popularity because Fermat noted in his copy "Cubum autem in duos cubos, aut quadratoquadratum in duos quadratoquadratos, et...
I'm interested to know whether the equation $$1 = 2 - \frac{1}{2 - \frac{1}{2 - \cdots}}$$ is true or not. It can be shown easily that if the continued fraction converges, it cannot converge to anything else than 1. It seems that if the continued fraction converges, the convergence is very slow. The apparent slowness of the convergence makes it difficult to estimate the presence of true convergence numerically. At the moment I don't know whether this converges or not.
Back
Top