高程异常代码

%%%%%%%%%%%%%%%%%%%%%%%%%%%GPS高程异常曲面拟合程序%%%%%%%%%
% k=input（‘已知点个数k=‘）
% c=input（‘检核（拟合）点数c=‘）
clear all;
time=clock;                                  % 计时开始
load study3.txt;
load work3.txt;
k=10;
c=30;
%%%%%%%%%%%%%%%%%%%%利用10个已知点求拟合公式的系数
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
 %%%%%%%%%%%%%%%%%%%用拟合公式求工作数据的高程异常值
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