加权最小二乘法源码

加权最小二乘法源码（使用最原始的的定义公式编码）

function p=zxechf（xywn）
%w是每点的权重，n是最高次幂 p是输出的系数
figure（1）
plot（xy‘+‘）; %做出原始散点图     
legend（‘数据点（xiyi）‘）
xlabel（‘x‘）; 
ylabel（‘y‘）; 
title（‘数据点（xiyi）的散点图‘）
hold on


A=zeros（（n+1）（n+1））;
for i=1:（n+1）%行  循环n+1次
    for j=1:（n+1）%列  循环n+1次
        A（ij）=sum（w.*x.^（（i-1）+（j-1）））;%给矩阵A赋值
    end
end


B=zeros（（n+1）1）;
for i=1:（n+1）%循环n+1次
    
B（i1）=sum（w.*x.^（i-1）.*y）;%给矩阵B赋值
end

p=A\B;              %系数的列矩阵  是由低次向高次顺序 排列

%poly2str %X是表示一个多项式的向量 %Y多项式的字符串表示
C=zeros（1（n+1））;
C=p‘;  
X=fliplr（C）; %对调顺序，由高次到低次排列
            %

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1356  2015-07-17 21:43  jia_quan_zui_xiao_er_cheng_fa.m