matlab練習程序（三種方法解最小二乘）

解最小二乘的方法有很多，這里給出常見的三種方法實現。

一是一般方法，之前博客一般都用這種方法。

二是svd分解法，之前有用過svd(見這里，這里，這里和這里)解其他問題，但是沒用來解過最小二乘。

三是qr分解法，這個好像沒用過。

這里主要總結記錄一下實現方法。

代碼如下：

clear all;
close all;
clc;

a=2;b=2;c=-3;d=1;e=2;f=30;   %系數         
[x,y]=meshgrid(0:0.1:30);
z=a*x.^2+b*y.^2+c*x.*y+d*x+e*y +f;      %原始模型     
mesh(x,y,z)
hold on;

X=x(:);
Y=y(:);
Z=z(:);
X=X(1:100:end);
Y=Y(1:100:end);
Z=Z(1:100:end);
plot3(X,Y,Z,'ro')

XX=[X.^2 Y.^2 X.*Y X Y ones(length(X),1)];
YY=Z;

%一般方法
C1=inv(XX'*XX)*XX'*YY;                                

%svd分解法
Y=zeros(6,1);
[U,S,V]=svd(XX);
B=U'*YY;
for i=1:6
    Y(i)=B(i)/S(i,i);
end
C2=V*Y;

%qr分解法
[Q,R] = qr(XX);
C3 = pinv(R)*(Q'*YY); 

C1'
C2'
C3'

三種方法結果是一致的：

模型如下：