电力系统算法

zibuyu_2017

function y = curve(x,y,n)
$ r2 L* U- d* D2 T) G- m+ w; EL=size(x);
" F7 j: d9 D( \$ ?6 E9 O0 ~8 e6 Z" Wm=L(2);
2 p" m* t% X# Y3 c$ d2 ^5 [for i=1:m;
    D(i,1)=1;
    for j=2:n+1;
3 [- ~/ R# x+ a0 d4 y+ R        D(i,j)=D(i,j-1)*x(i);
    end
8 w& q. _, u0 m- c- K3 [& l0 bend
A=D'*D;
0 m) C' o6 n  I2 {b=D'*y';
S=A\b;
) j9 \3 y! [+ w5 S5 cfor i=1:n+1
, _, m7 f6 {1 Z/ X" `) \    R(i)=S(n+2-i);
0 O4 x& }1 u8 C! w) pend
y=poly2str(R,'x');
YY=D*S;
9 C+ J+ n2 t9 y: ~# ^5 Isum=0;
( U& t1 S( s% j% wmax=0;
* V+ L) d! N/ {, B( d1 a$ L# u  H; jfor i=1:m
9 B3 K7 A4 u/ @( o1 r    D(i)=abs(YY(i)-y(i));
: P" A( {  m/ c5 a1 V' [6 \" [    if max<D(i)
        max=D(i);
    end
    sum=sum+D(i)^2;
9 t2 c$ m" [( ~8 a5 e9 s! ~* oend
- G9 Y' X* }$ o% u: F3 ]sum=sqrt(sum);
disp(y);
! L. ^, k2 ^9 {. T% w* w* _) Odisp('均方误差为');
disp(sum);
disp('最大误差为');
disp(max);
end
5 H0 J. y7 `3 A7 R
; e: Z& ~7 Y  E) ]0 k$ J* d4 k: k$ N
2 r. ?9 K+ ~8 u; `4 g5 G

zibuyu_2017

强烈支持楼主ing……