Fitting a convolution

11 次查看（过去 30 天）
Niles Martinsen
Niles Martinsen 2012-6-7
Hi
I have the following piece of code:
------------------------------------------------------------------------------------------
dataX = -4:1:4;
dataY = [0 -1 -10 -40 -55 10 40 10 1];
figure(1)
plot(dataX, dataY, '*');
x = -2:0.01:2;
plot(-4:0.01:4, 1*conv(exp(-x).*heaviside(x),sin(x)), dataX, dataY, '*')
mdl = @(a, x)(a(1)*conv(exp(-x).*heaviside(x),sin(x)));
par=[1];
[fitted_par, r, J, cov, mse] = nlinfit(dataX, dataY, mdl, par);
------------------------------------------------------------------------------------------
What I am trying to do is to fit the function (a convolution) to the data set. However, I get an error due to vector size mismatch. I'm not sure what is going on here: The fitting routine should not care about how many data points I have?
Regards, Niels.

请先登录，再回答此问题。

回答（2 个）

Andrei Bobrov
Andrei Bobrov 2012-6-8
Try
mdl = @(a, x)a*conv(exp(-x).*heaviside(x),sin(x),'same');
[fitted_par, r, J, cov, mse] = nlinfit(dataX, dataY, mdl, 1);
  1 个评论
Niles Martinsen
Niles Martinsen 2012-6-8
Hi
Thanks for replying. It runs now, but the fit gives me the amplitude 48.17, which is a horrible guess. Is something wrong with the program, or is it just a difficult function to fit?
Best regards,
Niels.

请先登录，再进行评论。

Niles Martinsen
Niles Martinsen 2012-6-8
Andrei's suggestion works, but when I increase the resolution of x, then the convolution changes as well. Is this just a property on conv, or have I made an error?

类别

在 Help Center 和 File Exchange 中查找有关 Correlation and Convolution 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by