Why Does Deconvolution Result in a Peak Followed by a Valley?

Click For Summary
SUMMARY

The discussion focuses on the challenges of performing deconvolution in signal analysis, specifically in the context of quantifying blood flow in the brain. The user is attempting to deconvolve two functions, f(t) and g(t), which describe the concentration of a contrast solution, using the discrete Fourier transformation. The resulting function h(t) often exhibits a peak followed by a valley, raising questions about the validity of the deconvolution process. Key insights suggest that discrepancies in parameters c and m may affect the deconvolution outcome.

PREREQUISITES
  • Understanding of deconvolution techniques in signal processing
  • Familiarity with Fourier Transform and its applications
  • Knowledge of mathematical modeling of functions in signal analysis
  • Experience with blood flow quantification methods in medical imaging
NEXT STEPS
  • Study the principles of deconvolution in signal processing
  • Learn about the implications of parameter discrepancies in deconvolution
  • Explore advanced Fourier Transform techniques for signal analysis
  • Investigate the use of unit step functions in convolution and deconvolution
USEFUL FOR

Researchers in medical imaging, signal processing engineers, and professionals involved in quantifying physiological signals will benefit from this discussion.

PaulPaul
Messages
1
Reaction score
0
Sorry, but I don't know where the topic about signal analysis should go..?
----------------------------------------------------------------------
I'm doing research about quantifying blood flow in the brain. Basically I need to know how to perform deconvolution (I think.)

I have two functions that describe the concentration of contrast solution:

f(t) = a(t-to)^b exp[ -(t-to)/c ]
g(t) = d(t-t1)^k exp[ -(t-t1)/m ]

where a,b,c,d,k,m,to, and t1 are known parameters (found by fitting the function to measurements.)

now say f = h * g (* - convolution)
How can I find h(t)? What type of result can I expect?

I have tried using the discrete Fourier transformation. But I often find h to be a peak followed by a small valley. I don't understand why this would be the deconvolution.

Please help. I'm stuck.
(http://s153.photobucket.com/albums/s235/s1020099/)
 
Physics news on Phys.org
PaulPaul said:
Sorry, but I don't know where the topic about signal analysis should go..?

try the USENET newsgroup comp.dsp .

----------------------------------------------------------------------
I'm doing research about quantifying blood flow in the brain. Basically I need to know how to perform deconvolution (I think.)

I have two functions that describe the concentration of contrast solution:

f(t) = a(t-to)^b exp[ -(t-to)/c ]
g(t) = d(t-t1)^k exp[ -(t-t1)/m ]

where a,b,c,d,k,m,to, and t1 are known parameters (found by fitting the function to measurements.)

now say f = h * g (* - convolution)
How can I find h(t)? What type of result can I expect?


if c is not the same number as m, i don't think you can do it. well, it might depend a little on what if f() and g() have unit step functions applied to them.

but the basic idea is to compute the Fourier Transform of f(t) and g(t), divide [itex]F(\omega)[/itex] by [itex]G(\omega)[/itex] to get [itex]H(\omega)[/itex] and then inverse Fourier transform that result.

I have tried using the discrete Fourier transformation. But I often find h to be a peak followed by a small valley. I don't understand why this would be the deconvolution.

Please help. I'm stuck.
(http://s153.photobucket.com/albums/s235/s1020099/)
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 6 ·
Replies
6
Views
5K
  • · Replies 23 ·
Replies
23
Views
3K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K
Replies
8
Views
2K
Replies
1
Views
2K