最小二乘拟合matlab程序

%线性最小二乘逼近的求法
x=input（‘输入x观测值：‘）;
y=input（‘输入y观测值：‘）;
n=input（‘输入阶次n：‘）;
m=length（x）;
b=0;
A=0;
if n>m
    disp（‘阶次应小于观测值个数‘）;
else
for i=1:m
    for j=0:n
        A（ij+1）=x（i）^j;
    end
end
G=A‘*A;
for j=0:n
    d=0;
    for i=1:m
        d=d+y（i）*x（i）^j;
    end
    b（j+1）=d;
end