电力系统算法

zibuyu_2017

function y = curve(x,y,n)
L=size(x);
5 S" [6 b+ _( T1 k3 r0 Fm=L(2);
7 G/ c1 x' {5 J: Y( vfor i=1:m;
. b0 {6 Z2 Z% l9 [5 }3 t    D(i,1)=1;
. U1 t/ F" m3 w& q    for j=2:n+1;
: H/ G3 B1 H! C" h        D(i,j)=D(i,j-1)*x(i);
# b7 ?2 M# C# x8 q& Q( ?& w    end
end
% l! H; m3 X! i/ b. y; SA=D'*D;
) o& y; g* |5 g% `( Z- Jb=D'*y';
S=A\b;
for i=1:n+1
    R(i)=S(n+2-i);
end
y=poly2str(R,'x');
6 K' r7 F' S' ]& w9 kYY=D*S;
sum=0;
max=0;
for i=1:m
7 V0 ]9 `3 }. D' X    D(i)=abs(YY(i)-y(i));
    if max<D(i)
' f& A7 P7 \) _        max=D(i);
+ p5 z, @5 N, [5 q5 w    end
1 `) L3 V  ~1 {6 p- {6 M/ _6 o8 y( \    sum=sum+D(i)^2;
: Z8 K8 @/ v/ ]3 v4 z/ `* eend
sum=sqrt(sum);
7 U# B! u5 ^: g! S* cdisp(y);
# S& M5 b8 N' z4 f1 `; I$ H% W+ sdisp('均方误差为');
/ c% K: z- A- D, u9 ~2 ]! @disp(sum);
disp('最大误差为');
' }" ~8 I4 ?+ J/ @4 a: D& udisp(max);
end

8 d% z8 I  Y2 ]3 d7 O8 ?
& I" i2 T' D: C

zibuyu_2017

强烈支持楼主ing……