% matlab用最小二乘法解下列矛盾方程組x1-x2=1;-x1+x2=2;2x1-2x2=3;-3x1+x2=4
% 對於方程組Ax=b,A為n×m矩陣,如果A列滿秩,且n>m。
% 則方程組沒有精確解,此時稱方程組為超定方程組。
% 線性超定方程組經常遇到的問題是資料的曲線擬合。
% 對於超定方程,在MATLAB中,利用左除命令(x=A\b)來尋求它的最小二乘解。
clc; clear all;
A = [1 -1
-1 1
2 2
-3 1];
B = [1; 2; 3; 4];
A\B
結果
ans =
-0.3654
1.5192
對於本題:
a=[
476.1905 238.0952 1428.571
514.2857 102.8571 1028.571
1016.949 112.9944 1016.949
1108.972 138.6215 970.3504
827.5862 275.8621 1379.31
1759.398 270.6767 947.3684 ];
A=a(:,1:2);
b=a(:,3);
A\b
%即可求得M,N的值
0.0166
5.2825
% matlab用最小二乘法解下列矛盾方程組x1-x2=1;-x1+x2=2;2x1-2x2=3;-3x1+x2=4
% 對於方程組Ax=b,A為n×m矩陣,如果A列滿秩,且n>m。
% 則方程組沒有精確解,此時稱方程組為超定方程組。
% 線性超定方程組經常遇到的問題是資料的曲線擬合。
% 對於超定方程,在MATLAB中,利用左除命令(x=A\b)來尋求它的最小二乘解。
clc; clear all;
A = [1 -1
-1 1
2 2
-3 1];
B = [1; 2; 3; 4];
A\B
結果
ans =
-0.3654
1.5192
對於本題:
a=[
476.1905 238.0952 1428.571
514.2857 102.8571 1028.571
1016.949 112.9944 1016.949
1108.972 138.6215 970.3504
827.5862 275.8621 1379.31
1759.398 270.6767 947.3684 ];
A=a(:,1:2);
b=a(:,3);
A\b
%即可求得M,N的值
ans =
0.0166
5.2825