# wiener deconvolution

Started by November 3, 2009
```Hi every one

I am working on Wiener deconvolution. my aim is to deconvolve the
refelectivity function. I am using threshold method in inverse filtering
for the deconvolution.

Hf = H.*(abs(H)>0)+1/gamma*(abs(H)==0);
iHf = 1./Hf;
iHf = iHf.*(abs(H)*gamma>1)+gamma*abs(Hf).*iHf.*(abs(Hf)*gamma<=1);

W = iHf.*Yyy./(Yyy+ sigma); %Inverse and denoise filters in cascade

the code is working when gamma = 1 and signal =0;

But when I change the parameters gamma there is no deconvolved output.

can anybody help me in estmating the parameter gamma.
or please mention refference to find out the gamma;

```
```On 3 Nov, 19:15, "ramjanagam" <ramamoorthy...@yahoo.com> wrote:
> Hi every one
>
> I am working on Wiener deconvolution. my aim is to deconvolve the
> refelectivity function. I am using threshold method in inverse filtering
> for the deconvolution.
>
> Hf = H.*(abs(H)>0)+1/gamma*(abs(H)==0);
> iHf = 1./Hf;
> iHf = iHf.*(abs(H)*gamma>1)+gamma*abs(Hf).*iHf.*(abs(Hf)*gamma<=1);
>
> W = iHf.*Yyy./(Yyy+ sigma); %Inverse and denoise filters in cascade
>
> the code is working when gamma = 1 and signal =0;
>
> But when I change the parameters gamma there is no deconvolved output.

If the algorithm breaks down when you change one parameter, it
is likely an implementation blunder somewhere. First of all make
sure your impleentation is correct with respect to the algorithm
you have chosen. Only after you have verified this, does it make