电力系统算法

zibuyu_2017

function y = curve(x,y,n)
L=size(x);
m=L(2);
for i=1:m;
    D(i,1)=1;
. o2 }, A4 {9 S2 z- y9 S- W) }    for j=2:n+1;
        D(i,j)=D(i,j-1)*x(i);
8 L9 W/ j+ r; r- i  n    end
; [, j% X* }( ^) P" Eend
+ v, i  l: q8 ?8 YA=D'*D;
5 A; S: R) W" a9 G/ @, o) X9 s/ R* zb=D'*y';
S=A\b;
* m3 o# \1 u" R. ~6 N9 hfor i=1:n+1
    R(i)=S(n+2-i);
end
% g3 R( k, d: Y$ B/ Q( py=poly2str(R,'x');
YY=D*S;
2 P! _# R+ g+ J* x: Y; v: p- Dsum=0;
$ l4 m, f( Y( b5 rmax=0;
for i=1:m
* m2 a- I$ V; V. x  B    D(i)=abs(YY(i)-y(i));
    if max<D(i)
) {" S7 q. g  K) p6 J        max=D(i);
8 T* Y, ?- o# o( o' s9 m+ U    end
    sum=sum+D(i)^2;
end
sum=sqrt(sum);
+ c) u/ ^) u4 a7 W3 q9 P) rdisp(y);
disp('均方误差为');
. e& @+ C3 r4 G  l% r! }3 Z+ }$ U' ydisp(sum);
; n, i$ {5 b8 L# A; }disp('最大误差为');
: o9 D" w; r+ L5 Bdisp(max);
5 F9 \* V. E0 n! o& ~end
; N6 P# w) T9 Y/ u* @" s# l, I
4 G- ~$ J) X. j  q

zibuyu_2017

强烈支持楼主ing……