nihe.m

来自「自己编写的多项式拟合程序」· M 代码 · 共 46 行

load  b02.dat
b=b02.';
b0=b-mean(b);
N=length(b0);

%n阶多项式拟合
n=5;
for i=1:2*n
    k(i)=0;   
    for j=1:N
        k(i)=k(i)+j^i/N;
    end
end
for i=1:n+1
    ky(i)=0;
     for j=1:N
       ky(i)=ky(i)+j^(i-1)*b(j)/N;
    end
end

for i=1:n+1
    for j=1:n+1
        if (i==1)&&(j==1) 
            kk(i,j)=1;
        else
        kk(i,j)=k((i-1)+j-1);
        end
    end
end
a=kk\ky';    %采用直接法解方程组

%拟和值
nihezhi=0;
for i=1:n+1
    nihezhi=nihezhi+a(i)*(1:N).^(i-1);
end 

plot(1:length(b),b,1:length(nihezhi),nihezhi);

%figure(2);
%subplot(2,1,1);
%plot(b);
%subplot(2,1,2);
%plot(nihezhi);