电力系统算法

zibuyu_2017

function y = curve(x,y,n)
0 L& N5 O9 C; tL=size(x);
: l; J3 e$ A  C% Q1 C; i. Qm=L(2);
for i=1:m;
7 I5 L# c1 G9 t6 g0 {+ C    D(i,1)=1;
    for j=2:n+1;
        D(i,j)=D(i,j-1)*x(i);
    end
6 \: s' {# d' jend
% K# ^; \( P" [) k2 d) E$ c5 `A=D'*D;
b=D'*y';
1 G9 U9 J6 x7 _S=A\b;
; n! g) G: ]5 K/ Y0 ?for i=1:n+1
    R(i)=S(n+2-i);
end
y=poly2str(R,'x');
1 w8 N; G- @2 l( sYY=D*S;
/ k. |  E  i0 e. M, ]7 `- c+ x! F9 [sum=0;
( }9 b5 k# p/ k0 J- _max=0;
* H+ b7 {) I  U6 }8 Gfor i=1:m
6 W, T8 [; o# U* O% C" A+ a    D(i)=abs(YY(i)-y(i));
    if max<D(i)
7 ^) M+ B. j: y$ Y8 s( L8 _        max=D(i);
0 o& e" a/ r: @; T/ [  ~. a/ r! X    end
+ r( Q  I6 u6 Q5 y    sum=sum+D(i)^2;
+ T6 b" r) h8 T0 [8 l, [, ]8 |end
' ^: F: O5 {, msum=sqrt(sum);
& @+ r$ y/ c& B4 J' V0 bdisp(y);
disp('均方误差为');
+ k) q4 Z0 J/ ^, j& [. Z, \8 ndisp(sum);
, w: |1 y( i$ H% s7 R/ w/ Gdisp('最大误差为');
disp(max);
' c9 P# @* w0 ^end
/ U' m/ q& Q9 e! x% V
7 `3 e) W. {# \- h4 \% w3 ]
8 C0 v" ?6 e+ w( A

zibuyu_2017

强烈支持楼主ing……