Fitting every fifth datapoint one after another

Question

J.P. 2016-11-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/311193-fitting-every-fifth-datapoint-one-after-another

评论： J.P. 2016-11-8

Hello!

I want to do a linear regression for a bunch of data, using every fifth or third (maybe tenth) Datapoint. The data will be vary in size and therefor I need to do a function that varies as well.

I did the poly1-fit to the first 13 data already by typing it for each fit. Can you think of any way to solve it via function maybe a loop to solve this seccessively for any size of x- or y-datapoints.

I will be happy for any kind of response that helps me to solve this problem.

Thanks a lot!

The result should look like this for more datapoints:

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

dbmn 2016-11-8

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/311193-fitting-every-fifth-datapoint-one-after-another#answer_242456

在 MATLAB Online 中打开

I have a solution that avoids the loops altogehter, but might not be as readable

n = 5; % Number of datapoints
x1 = reshape(x, n, []); % please make sure, that your number will be divisable into n Elements
y1 = reshape(y, n, []);
a = (x1-x1(1,:))\(y1-y1(1,:)); % we use mldivide here and correct for the origin
a = a(1,:);
b = y1(1,:) - a.*x1(1,:); % here we shifted back

b is the offset and a the slope of your curves

3 个评论
显示 1更早的评论隐藏 1更早的评论

J.P. 2016-11-8

Hi! Thanks for your help. But your solution gives this Error:

Error using - Matrix dimensions must agree.

Error in polyfit_Inet (line 6) a = (x1-x1(1,:))\(y1-y1(1,:)); % we use mldivide here and correct for the origin

J.P. 2016-11-8

在 MATLAB Online 中打开

Okay, I changed the parts a little bit and it works!!!

n = 6; % Number of datapoints
x1 = reshape(x, n, []); % please make sure, that your number will be divisable into n Elements
y1 = reshape(y, n, []);
a = (x1(1)-x1(1,n))\(y1(1)-y1(1,n)); % we use mldivide here and correct for the origin
a = a(1,:);
b = y1(1,:) - a.*x1(1,:); % here we shifted back
figure()
hold on
hdl1 = plot(x1, y1, 'b.');
hdl2 = plot(x1, a.*x1(2,:)+b, 'r-');