电力系统算法

zibuyu_2017 发表于 2017-8-10 19:15:53

function y = curve(x,y,n)
L=size(x);
m=L(2);
for i=1:m;
D(i,1)=1;
for j=2:n+1;
D(i,j)=D(i,j-1)*x(i);
end
end
A=D'*D;
b=D'*y';
S=A\b;
for i=1:n+1
R(i)=S(n+2-i);
end
y=poly2str(R,'x');
YY=D*S;
sum=0;
max=0;
for i=1:m
D(i)=abs(YY(i)-y(i));
if max<D(i)
max=D(i);
end
sum=sum+D(i)^2;
end
sum=sqrt(sum);
disp(y);
disp('均方误差为');
disp(sum);
disp('最大误差为');
disp(max);
end

zibuyu_2017 发表于 2017-8-10 19:17:47