电力系统算法

zibuyu_2017

function y = curve(x,y,n)
% g/ ?( ^6 E% c' h/ `) p, d; uL=size(x);
  W8 r8 h& Z" b! n8 s! U0 um=L(2);
for i=1:m;
2 P# W1 n/ `' r. }; D/ P) c- p, F    D(i,1)=1;
    for j=2:n+1;
        D(i,j)=D(i,j-1)*x(i);
    end
end
2 z. a, v" E; t' J: RA=D'*D;
b=D'*y';
9 K0 q* s4 e4 d+ OS=A\b;
for i=1:n+1
9 I% K7 t# G& C. k% F, }    R(i)=S(n+2-i);
  x: N& N; \  a( T) tend
y=poly2str(R,'x');
YY=D*S;
" M  P% w- ~7 w& J, Msum=0;
) ]- t, K. K0 Nmax=0;
for i=1:m
3 X' m' i/ O% I    D(i)=abs(YY(i)-y(i));
    if max<D(i)
        max=D(i);
    end
    sum=sum+D(i)^2;
7 K2 D; c5 _0 }$ U8 ^# p" y% b3 Aend
sum=sqrt(sum);
: \- F- r7 w0 Bdisp(y);
disp('均方误差为');
disp(sum);
disp('最大误差为');
( |. \0 r7 `! X8 a* k/ Z' w# Ndisp(max);
( n% Q, t( |/ X) N) O! Z4 ^% `2 Pend
" ^  E8 _1 C& L1 V8 q

zibuyu_2017

强烈支持楼主ing……