squar_least.m

来自「常微分方程数值解法Runge_Kutta解法」· M 代码 · 共 32 行

function s=squar_least(x,y,n)
if nargin<3 n=1;end
phi2=zeros(n+1,n+1);
h=size(x);





for i=1:n+1
    for j=1:n+1
        phi2(i,j)=0;
        for k=1:h(2);    
        phi2(i,j)= phi2(i,j)+(x(k)^(j-1)*x(k)^(i-1));
        end
    end
end



phif=zeros(n+1,1);
for i=1:n+1
    phif(i)=0;
    for k=1:h(2)
    phif(i)= phif(i)+x(k)^(i-1)*y(k);
    end
end


s=phi2\phif;