residue.m

多项式二次曲线拟合程序！用于实验设计！另有画图函数并且可以进行方差分析1

load x.mat
load y.mat
data_X=x;
data_Y=y;
[m,n]=size(x);
for i=1:m
    Y_estimated(1,i)=A+data_X(i,:)*B+data_X(i,:)*C*data_X(i,:)';
end

y_mean=mean(data_Y);
m=9
n=17
%计算回归
U=sum((Y_estimated-y_mean).^2);
a1=U/m
%计算残差
Q=sum((data_Y-Y_estimated').^2);
b1=Q/(n-m-1)
%计算总方差
total=U+Q;
total_freedom=size(data_Y,1)-1
%计算方差比
F=a1/b1
%计算复相关系数
R=sqrt(U/total)