代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