基于EDA数据包络分析的综合评价代码.zip

%function y=zhangyong(x0) clear clc format long g x0=[2804.5 2924 2991 3301 3508]; n=size(x0,2); for i=1:n x1(i)=sum(x0(1:i)); end YN=x0(2:n)'; B=[-0.5.*(x1(1:n-1)+x1(2:n)).' ones(n-1,1)]; cap_a=inv(B'*B)*B'*YN; a=cap_a(1);b=cap_a(2); cap_x1(1)=x0(1); for k=1:n-1 cap_x1(k+1)=(x0(1)-b/a)*exp(-1*a*k)+b/a; end cap_x0(1)=cap_x1(1); for k=1:n-1 cap_x0(k+1)=cap_x1(k+1)-cap_x1(k); end a b y=cap_x0 c1=x0(1)-b/a c2=b/a s=30; error=cap_x0-x0; error_mean=mean(error); error_cov=sqrt(cov(error)/n); origin_cov=sqrt(cov(x0)/n); C=error_cov/origin_cov P=sum((error-error_mean)<0.6745*origin_cov)/n model_x0(1)=x0(1); model_x1(1)=x0(1); for t=1:s model_x1(t+1)=c1*exp(-1*a*t)+c2; model_x0(t+1)=model_x1(t+1)-model_x1(t); end model_x0