电力系统算法

zibuyu_2017

function y = curve(x,y,n)
L=size(x);
m=L(2);
for i=1:m;
  P: q; L' T  _3 d    D(i,1)=1;
    for j=2:n+1;
- E& v7 ?2 I# \. R% k        D(i,j)=D(i,j-1)*x(i);
1 X3 T7 j! V4 d    end
end
$ F+ H# }( r9 W8 zA=D'*D;
b=D'*y';
S=A\b;
8 X6 X; l9 q  K/ s2 Xfor i=1:n+1
8 b3 f3 K) |) E    R(i)=S(n+2-i);
. R+ L( |- C6 G" R! S0 nend
5 ?, e: E& n2 N4 Uy=poly2str(R,'x');
YY=D*S;
& c( @& ]* x( D4 P8 D: msum=0;
max=0;
for i=1:m
0 \" R# g6 R2 |6 Q3 C$ |0 O    D(i)=abs(YY(i)-y(i));
$ v# b6 y9 B' y2 c! y; |- T7 F    if max<D(i)
$ z* d- r5 q+ R: _) ?' Y        max=D(i);
% w& p' n6 |( l: [% R) O. Y: ~    end
+ |5 P# `+ o" G0 T2 S1 [    sum=sum+D(i)^2;
" ~2 j; B5 L& ]2 {% _9 G: e2 Aend
sum=sqrt(sum);
disp(y);
disp('均方误差为');
disp(sum);
* P  J% a$ Z# F( X* r& ^. O: ?disp('最大误差为');
" s  o# u6 F: s6 O$ O  B# n, ]disp(max);
. ]# c  A0 n6 _& a' p  I) Vend
' w1 F" W1 Y( j9 [

zibuyu_2017

强烈支持楼主ing……