2013电工杯数学建模A题论文（国家二等奖）(8)

代null的前后2N个数据

for i=1:M

if (s(i)==-1000)&(i-N>0) msum=0; for j=-N:N if j~=0 msum=msum+s(i+j);

end

s(i)=msum/2/N; end elseif (s(i)==-1000)&(i-N<=0)

msum=0; for j=1:N msum=msum+s(i+j);

end

s(i)=msum/N; elseif (s(i)==-1000)&(i+N>M)

msum=0; for j=-N:1 msum=msum+s(i+j);

end

s(i)=msum/N; end end q=q+s; end plot(q);

NL=12;Nt=length(q)/NL;k=1;%NL12*5=60s一分钟

for i=(1:Nt)*NL s1(k)=q(i); k=k+1; end

T=length(s1);N=2; % N滑动平均个数 2N+1个

ysum=0;k=1; for i=N+1:N:T-N for j=-N:N

ysum=ysum+s1(i+j);

end

yk(k)=ysum/(2*N+1); ym(k)=s1(i)-yk(k); ysum=0;k=k+1; end figure plot(yk) figure plot(ym)

n=length(ym) yav=sum(ym)/n

ys=sqrt(sum((ym-yav).^2)/(n-1)) ymaxmin=max(ym)-min(ym) %%灰色预测模型（滚动预测未来7天数据）

clc,clear Da=42;

x0=[3610 3853 3705 4097 2963 2230 1781 1601 1692 1627 1470 1450 1366 1060 1335 1064]';

yz=zeros(length(x0),Da);yep=yz;ydelt=yz;yrho=zeros(length(x0)-1,Da);

for i=1:Da

n=length(x0);

lamda=x0(1:n-1)./x0(2:n); range=min(max(lamda')); x1=cumsum(x0);

B=[-0.5*(x1(1:n-1)+x1(2:n)),ones(n-1,1)];

Y=x0(2:n); u=B\\Y;

x=dsolve('Dx+a*x=b','x(0)=x0');

x=subs(x,{'a','b','x0'},{u(1),u(2),x0(1)});

yuce1=subs(x,'t',[0:n-1]); y=vpa(x,6);

yuce=[x0(1),diff(yuce1)]; epsilon=x0'-yuce;

delta=abs(epsilon./x0');

rho=1-(1-0.5*u(1))/(1+0.5*u(1))*lamda';

yz(:,i)=yuce'; yep(:,i)=epsilon'; ydelt(:,i)=delta'; yrho(:,i)=rho'; x0=yuce'; end

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库2013电工杯数学建模A题论文（国家二等奖）(8)在线全文阅读。