最小二乘法 MATLAB

有了前面对最小二乘法估计学习过，了解了原理（参考：一元线性回归模型：最小二乘估计），于是用MATLAB看了一下效果，后期再对RANSAC学习一下。

[plain]  view plain  copy

 print ?

1. close all  
2. clear,clc  
3.   
4. x=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18];  
5. y=[3.75 20.95 19.95 39.3 28.45 36.85 54.40 71.55 68.00 73.85 78.65 73.75 88.25 92.3 85.9 100.2 109.2 114.9];  
6.   
7. len = length(x)  
8.   
9. ave_x = sum(x(1:len))/len  
10. ave_y = sum(y(1:len))/len  
11. A = 0  
12. B = 0  
13. hold on;  
14. for i = 1:length(x)  
15.     plot(x(i),y(i),'rs','MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize',6)  
16.     %text(x(i),y(i),['(',num2str(x(i),'%5.2f'),',',num2str(y(i),'%5.2f'),')'],'EdgeColor','red','BackgroundColor',[.7 .9 .7],'VerticalAlignment','bottom');  
17.     A = A +(x(i)-ave_x)*(y(i)-ave_y)  
18.     B = B +(x(i)-ave_x)*(x(i)-ave_x)  
19. end  
20. b = A/B  
21. a = ave_y-b*ave_x  
22. xx = 0:1:20  
23. yy = a+b*xx  
24. plot(xx,yy,'r-')  

结果如下：

来源： http://blog.youkuaiyun.com/chentravelling/article/details/50817774