离散试验数据点的正交多项式最小二乘拟合

function a=ZJZXEC（xym）
if（length（x） == length（y））
    n = length（x）;  
else
    disp（‘x和y的维数不相等！‘）;
    return;
end                  %维数检查

syms v;
d = zeros（1m+1）;
q = zeros（1m+1）;
alpha = zeros（1m+1）;
for k=0:m
   px（k+1）=power（vk）;
end                      %x的幂多项式
B2 = [1];
d（1） = n;
for l=1:n
    q（1） = q（1） + y（l）;
    alpha（1） = alpha（1） + x（l）;
end
q（1） = q（1）/d（1）;
alpha（1） = alpha（1）/d（1）;
a（1） = q（1）;
B1 = [-alpha（1） 1];
for l=1:n
    d（2） = d（2） + （x（l）-alpha（1））^2;
    q（2） = q（2） + y（l）*（x（l）-alpha（1））;
    alpha（2） = alpha（2） + x（l）*（x（l）-alpha（1））^2;
end
q（2） = q（2）/d（2）;
alpha（2） = alpha（2）/d（2）;
a（1） = a