电力系统算法

zibuyu_2017

function y = curve(x,y,n)
. B  |7 Y: A7 T! r7 R# EL=size(x);
  t7 _, z% \6 A- e2 v3 `3 C1 Tm=L(2);
for i=1:m;
/ Z) q. J, _( T( e& F) L    D(i,1)=1;
5 M/ J" r' V$ W0 |4 L" k+ I- v    for j=2:n+1;
        D(i,j)=D(i,j-1)*x(i);
    end
1 v7 p8 ~/ Q) x* `5 v9 \7 Iend
, f; \9 q; y9 j3 T' Q% m" a# H3 n2 }A=D'*D;
b=D'*y';
S=A\b;
for i=1:n+1
! @" y0 `2 P" b    R(i)=S(n+2-i);
end
y=poly2str(R,'x');
YY=D*S;
sum=0;
' `7 _- w% K4 L" z" ymax=0;
* O! t3 A( @, Q& M9 o# u  h% ]for i=1:m
    D(i)=abs(YY(i)-y(i));
    if max<D(i)
/ f9 l  S$ w& V! g. A        max=D(i);
    end
; P7 a- Z- o9 t8 a& q% @0 e    sum=sum+D(i)^2;
end
sum=sqrt(sum);
disp(y);
  ~( ?/ v+ n4 D; \* j: Cdisp('均方误差为');
disp(sum);
1 {% c% T# e: {( ]: J& R! P8 Idisp('最大误差为');
disp(max);
! o, o9 c( P( \$ @9 y5 Q' [end


1 R4 s; J/ ]* ~

zibuyu_2017

强烈支持楼主ing……