Linear best fit model
2 次查看(过去 30 天)
显示 更早的评论
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 个评论
回答(1 个)
William Rose
2022-7-14
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 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Linear and Nonlinear Regression 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!