Adaptive Filter | ADPCM | ADSP | ADSP-2181 | Aliasing | AMR | Anti-Aliasing | ARMA | Autocorrelation | AutoCovariance | Beamforming | Bessel | Blackfin | Butterworth | C6713 | CCS | Chebyshev | CIC Filter | Circular Convolution | Code Composer Studio | Comb Filter | Compression | Convolution | Cross Correlation | DCT | Decimation | Deconvolution | Demodulation | DM642 | DSP Boards | DSP/BIOS | DTMF | Echo Cancellation | Equalization | Equalizer | ETSI | EZLITE (Ez-kit Lite) | FFT | FFTW | FIR Filter | Fixed Point | FSK | G.711 | G.723 | G.729 | Gaussian Noise | Goertzel | GPIO | Hilbert Transform | IFFT | IIR Filter | Interpolation | Invariance | JTAG | Kalman | Laplace Transform | Levinson | LPC | McBSP | MIPS | Modulation | MPEG | Multirate | Notch Filter | Nyquist | OFDM | Oversampling | Pink Noise | Pitch | PLL | Polyphase | QAM | QDMA | Quantization | Quantizer | Radar | Random Noise | Reed Solomon | Remez | Resampling | RTDX | Sampling | Sharc | TI C6711 | Undersampling | Viterbi | Wavelets | White Noise | Wiener Filter | Windowing | XDS510PP | Z Transform

Hi All,

I am working with injecting contrast/dye and measuring the resulting
enhancement (resulting from the injection) in a patient.  

The injection is a square input (e.g., 20 ml @ 4 ml/sec for 5 sec).  Now,
I measure the resulting enhancement (OUTPUT) at a certain location.  Now,
this output is related to injection via the transfer function of the
patient.  

so, OUTPUT = TRANSFER FUNCTION of patient (convolved with) INPUT

Now, if we have to find the enhancement in the patient for an input of 4
ml @ 4ml/sec for 1 sec in time domain using the above data in time domain.
 In other words, how can I find the impluse response of the patient in the
time domain using deconvolution.  How do we do it in Matlab?  Is there any
Matlab function that can do this.  

Does this make sense?

Thankyou very much for taking your time!!

Best Regards,
NU




		
This message was sent using the Comp.DSP web interface on
www.DSPRelated.com

______________________________
New DSP Code Snippets Section now Live.   Learn more about the reward program for contributors here.

NU wrote:
> Now, if we have to find the enhancement in the patient for an input of 4
> ml @ 4ml/sec for 1 sec in time domain using the above data in time domain.
>  In other words, how can I find the impluse response of the patient in the
> time domain using deconvolution.
To deconvolve atleast one of the two operands should be known.
>How do we do it in Matlab?  Is there any
> Matlab function that can do this.
Yes. There is a Matlab function for deconvolution. But again you must
know one of the operands i.e. either transfer function or the input.
The function is 'deconv'.

______________________________
New DSP Code Snippets Section now Live.   Learn more about the reward program for contributors here.

> OUTPUT = TRANSFER FUNCTION of patient (convolved with) INPUT

This is only true if the patient is a linear system, which is unlikely.


(This is a good example of the opposite of anthropomorphism, which was
a thread here recently... the suggested name was mechanicism.)

______________________________
New DSP Code Snippets Section now Live.   Learn more about the reward program for contributors here.