电力系统算法

zibuyu_2017

function y = curve(x,y,n)
L=size(x);
% m/ m) T$ K7 m' k3 S1 w4 q# Qm=L(2);
for i=1:m;
    D(i,1)=1;
    for j=2:n+1;
& G* h' V1 {: w/ a! N" R        D(i,j)=D(i,j-1)*x(i);
. L; f8 j& ]5 d& I. E    end
end
A=D'*D;
$ A% R3 a/ Q& I5 K/ }  Db=D'*y';
S=A\b;
. ^; i: U0 t5 I  q, `$ ^0 ~for i=1:n+1
    R(i)=S(n+2-i);
: t: A! t5 Y9 M4 o) [' f7 _) pend
7 G7 q6 c1 J% r$ J: V8 y- zy=poly2str(R,'x');
' @; L' s- u& B2 h2 y2 K- q: @YY=D*S;
sum=0;
5 s, c& A) W0 Jmax=0;
for i=1:m
8 ^9 f1 t) E3 D+ I, O5 x9 L    D(i)=abs(YY(i)-y(i));
8 H+ Z! C& e; n  d    if max<D(i)
' A( @7 L/ x. m" l& h        max=D(i);
6 ~: G( ?0 v- N" W6 J& Q    end
    sum=sum+D(i)^2;
9 e9 U& t/ Z- K* z' ?) Tend
9 [; y2 Y$ k% Csum=sqrt(sum);
: X  E6 _3 ]# A, |' h: Z- V6 Rdisp(y);
disp('均方误差为');
0 ~/ O* I7 m3 D, ~1 Edisp(sum);
disp('最大误差为');
disp(max);
' i% @8 q# t7 ^! T: E1 Qend
) U* [4 r5 w8 E4 O' ~9 y- @) A
% r7 l% m6 d$ g4 I0 }1 g" M8 \
& p) D/ c+ k8 o" E: J+ k

zibuyu_2017

强烈支持楼主ing……