%%%%%%%%%%%%%%%%%%%%%%%%%%%GPS高程異常曲面擬合程序%%%%%%%%%
%?k=input（‘已知點個數(shù)k=‘）
%?c=input（‘檢核（擬合）點數(shù)c=‘）
clear?all;
time=clock;??????????????????????????????????%?計時開始
load?study3.txt;
load?work3.txt;
k=10;
c=30;
%%%%%%%%%%%%%%%%%%%%利用10個已知點求擬合公式的系數(shù)
Xk=study3（:1:3）;
Xc=work3（:1:2）;
%?Xx=zeros（k?3）
for?i=1:k
????Xx（i1）=1;
????Xx（i2）=Xk（i1）;
????Xx（i3）=Xk（i2）;
????Xx（i4）=Xk（i1）*Xk（i1）;
????Xx（i5）=Xk（i1）*Xk（i2）;
????Xx（i6）=Xk（i2）*Xk（i2）;???
end
????A=Xx\Xk（1:k3）?;
????
????A
?%%%%%%%%%%%%%%%%%%%用擬合公式求工作數(shù)據(jù)的高程異常值
for?j=1:c
????z（j1）=A（11）*1+A（21）*Xc（j1）+A（31）*Xc（j2）+A（41）*Xc（j1）*Xc（j1）+A（51）*Xc（j1）*Xc（j2）+A（61）*Xc（j2）*Xc（j2）;
end
%?%z
tt=work3（:3）;
%?t