电力系统算法

zibuyu_2017

function y = curve(x,y,n)
L=size(x);
m=L(2);
8 |% \9 r9 }. n! t$ u, }; Ifor i=1:m;
    D(i,1)=1;
    for j=2:n+1;
        D(i,j)=D(i,j-1)*x(i);
: P' F3 w6 T: G: C    end
& L1 }4 p* O  t# H5 r" w* _: bend
& S$ S/ D+ M- `0 oA=D'*D;
! y8 r$ o/ K, q& Ab=D'*y';
, \( _4 l9 F2 Y0 o( LS=A\b;
# i' m5 v% W6 y+ p1 m) pfor i=1:n+1
    R(i)=S(n+2-i);
end
* d. V! }0 \5 R4 A2 jy=poly2str(R,'x');
+ Y) X0 u1 x7 r" B% FYY=D*S;
sum=0;
( D- f- f( j5 s8 n# s2 {$ emax=0;
for i=1:m
4 Z: f8 S4 c% x; w9 P    D(i)=abs(YY(i)-y(i));
3 w! L/ c( P) m    if max<D(i)
" @" ^) Z" s2 N' H) w1 O        max=D(i);
    end
    sum=sum+D(i)^2;
end
. E& i' e, |# ?sum=sqrt(sum);
disp(y);
disp('均方误差为');
3 |9 x( b1 A8 A2 v7 gdisp(sum);
disp('最大误差为');
3 a% @, |+ U+ kdisp(max);
' U% R8 Y7 G# W$ H4 Qend
- n, x3 R& }3 J% P4 M0 t  n/ S. E
* T& J- q. \' Q0 `! u- ^* F
0 t0 W! V& g0 w( A

zibuyu_2017

强烈支持楼主ing……