【Matlab】实现趋势面分析

%%读取所在路径下数据文件中sheet1的C2：M54区域数据并存储到双精度矩阵num，
%将降水p、经纬度x、y以及x、y的变式读取到num
num = xlsread('/T_data.xls',1,'C2:M54');

P=num(:,1);%%提取降水数据
x=num(:,3);
y=num(:,4);
X_1=num(:,2:4);%%提取经纬度数据x、y,命名为矩阵X_1,用于一次趋势面模型
X_2=num(:,2:7);
X_3=num(:,2:11);
%%数据预处理处理结束

%%方法1求解回归系数矩阵
%%求解一次趋势面模型回归系数矩阵A_1，利用正规方程组矩阵形式反推
% X_1ts=X_1.';%%矩阵转置
% A_1=(X_1ts*X_1)^(-1)*X_1ts*P;%%求解回归系数矩阵的公式
% %