Linear best fit model
2 views (last 30 days)
Show older comments
Hello Community,
I have a query about best fit model.
From a large data set, I could generate a linear best fit line, which has a specific slope.
Now, now I want to compare how deviated this best fit is from another line having a different slope.
How can I force a specific slope value to draw another fit to this dataset?
Thanks in advance,
Ani
5 Comments
Answers (1)
William Rose
on 14 Jul 2022
Here is an example in which the fixed slope is 1.3 and the best fit slope is approximately 1.0.
x=0:1:100;
N=length(x); %number of data points
y=x+5*randn(1,N); %noisy data
% Find best-fit straight line
P = polyfit(x,y,1);
y1=P(1)*x+P(2);
rmse1=sqrt(sum((y-y1).^2)/N);
% Fit straight line with fixed slope
% The straight line should pass through the mean of the data
fixedslope=1.3;
xm=mean(x); ym=mean(y);
intercept=ym-fixedslope*xm;
y2=fixedslope*x+intercept;
rmse2=sqrt(sum((y-y2).^2)/N);
%% Display results on console
fprintf('Fit 1 (best fit) : y1=%.3fx%+.3f. RMS Error=%.2f\n',P(1),P(2),rmse1)
fprintf('Fit 2 (fixed slope): y2=%.3fx%+.3f. RMS Error=%.2f\n',fixedslope,intercept,rmse2)
%% Plot results
figure; plot(x,y,'rx',x,y1,'-g',x,y2,'-b')
xlabel('X'); ylabel('Y'); grid on; legend('data','fit 1','fit 2')
0 Comments
See Also
Categories
Find more on Interpolation in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!