View Single Post
Jul30-05, 02:04 PM   #5
 
Quote by tramoman
I am trying to model the behaviour of a reactor vessel when some
particles inside the container strike the inner wall using FINITE
ELEMENT MODELLING software.

I reproduced several random particles hitting the wall at different
times and position across the surface of the inside wall.

Now I would like to get rid of the TRANSFER FUNCTION of the container
from the random particles data using the deconvolution to extract only
the particle informations.

I apply an inpulse inside the vessel (the middle) and I acquire the
output on the outer surface.

THe PROBLEM: How can I deconvolve the spectra in a right way?
In fact the TRANSFER FUNCTION of the vessel changes according with the
position of the impact.

The tranfer function of the vessel will be completely different when
considering the position A or position B.


| |
posB> | |
| |
posA> [itex]| |>[/itex] output
| |
| |






Can somebody help me please?
What kind of DSP tool should I use?
it is really important for me to solve this problem.

Thanks a lot
M
Do you know the transfer function?If yes, I think dividing your data by an averaged form of it should do.Else,if you can assume that the collisions are completely random(white spectrum) and the impulse response of the vessel 'minimum phase' ---find the power spectrum from autocorrelation of data(so that you have the amplitude spectrum).Phase spectrum can then be uniquely found(for 'minimum phase' wavelets)--for the formula see any DSP book or Jon F. Claerbout's book on geophysical signal processing.