电力系统算法

zibuyu_2017

function y = curve(x,y,n)
L=size(x);
0 O* |- Z& x+ s6 _4 Km=L(2);
for i=1:m;
) q. g* C8 c  G. Z$ e, O    D(i,1)=1;
    for j=2:n+1;
        D(i,j)=D(i,j-1)*x(i);
    end
end
A=D'*D;
b=D'*y';
S=A\b;
for i=1:n+1
    R(i)=S(n+2-i);
0 J, D0 |* V$ n$ |3 X6 T4 Mend
- N- J! ^, H6 `8 `y=poly2str(R,'x');
. Y3 d; O! }& ?% iYY=D*S;
8 k& k% g( ?* t1 F, osum=0;
' @% _. b/ Z; w7 w& q% Vmax=0;
2 q4 w  }) A7 C! ]0 n0 A0 nfor i=1:m
    D(i)=abs(YY(i)-y(i));
4 t' o1 t8 {( ?* u, q1 i, y    if max<D(i)
0 H4 b& l& ]4 |9 H        max=D(i);
8 P0 N0 t7 T. n: M0 r    end
3 J& H3 U" ?( _4 y6 C" Y. x    sum=sum+D(i)^2;
! _) k( ?! t& T4 t  |0 Rend
. l# G  F$ n* Q8 tsum=sqrt(sum);
disp(y);
disp('均方误差为');
disp(sum);
) m$ K- j9 Q! r6 G  w: w/ Tdisp('最大误差为');
" N. P3 n5 m6 y1 V" h  E+ I, wdisp(max);
* I6 r. x  x' F# \% h, }end
/ q; G% E" g& ]/ j4 |1 H

9 u$ n. R$ O% \

zibuyu_2017

强烈支持楼主ing……