Deconvolution of a polynomial and exponential function

I got the same problem with following link:

https://www.mathworks.com/matlabcentral/answers/602500-deconvolution-of-a-polynomial-and-exponential-function

I got trouble using Bayesian deconvolution

here I already got wt(z) and a'(z) with same one dimension size array

I have no idea how to got iterative.

If need my data, I will give my data. thanks a lot ><

[larry liu]

this is the original equation, I need to solve R(z) by deconvolution

[larry liu]

here is my w(z) and a' and z axis file

[高飞 支]

I have the same trouble have you solved it?can you give me some advices?

Answer 1

[Yash]

Hi Larry,

I understand that you are facing issues while using the "deconv" function in MATLAB. "deconv(y,h)" deconvolves a vector "h" out of a vector "y" using polynomial long division, and returns the quotient "x" and remainder "r" such that "y = conv(x,h) + r".

Given that you already have a deconvolution function in your case, you may not need to use the "deconv" function. Instead, you can make an initial guess for "R(z)" and then obtain the approximation iteratively using a loop as follows:

n = length(df2);
w = Wzinstead';
a = df2';

R = randi(1000,n,1)/1000;     % Initial guess, random values

nIter = 10;
for i=1:1:nIter
    R1=R*(corr(w, a./(conv(w,R,'same'))));
    R=R1;
end

Hope this helps you address the issue.