记得数学还是自己喜欢的科目,可熟人发了这么一个求解,要求代码实现。结果把自己整的不会了。
参看他人文档“”Y=AX+B推导求解
以及ppt:中科院的线性代数ppt
1 大家熟悉的拟合法
2矩阵法
学会了吧??
3 例题
原图来源于 第一个引文的WellP.c's blog,若侵权请告知删除。
function [ Xpara ]=FitPlane( input_X , input_Y , input_Z)
% FileName : FitPlane
% CreatDate : 2018/7/10
% Describe : Least square fitting plane equation
% OutputPara : Xpara
% Note : Plane equation: z=a*x+b*y+c
% Author : wellp
X = input_X';
Y = input_Y';
I = ones(size(input_X'));
A = [X, Y, I];
b = input_Z';
Xpara = (A' * A) ^ -1 * A' * b;
end