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

Click For Summary
Deconvolution in signal analysis is essential for quantifying blood flow in the brain, particularly when dealing with concentration functions of contrast solutions. The user seeks to find h(t) from the convolution of two known functions, f(t) and g(t). The discussion suggests using the Fourier Transform method, where the Fourier Transform of f(t) is divided by that of g(t) to obtain H(ω), followed by an inverse Fourier Transform. A common issue arises when the results show a peak followed by a valley, which can occur due to differences in the parameters of the functions or the presence of unit step functions. Understanding these nuances is crucial for accurate deconvolution results.
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 F(\omega) by G(\omega) to get H(\omega) 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/)
 
Question: A clock's minute hand has length 4 and its hour hand has length 3. What is the distance between the tips at the moment when it is increasing most rapidly?(Putnam Exam Question) Answer: Making assumption that both the hands moves at constant angular velocities, the answer is ## \sqrt{7} .## But don't you think this assumption is somewhat doubtful and wrong?

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
949
  • · 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