noc_dst.m

平时练习时的一些仿真程序

StartTime=[90 04 09 00 00];
EndTime=[90 04 11 23 00];
HDataY=AnyChooseY(StartTime,EndTime);
% t=datenum(EndTime(1),EndTime(2),EndTime(3))-datenum(StartTime(1),StartTime(2),StartTime(3))+1;
col=length(HDataY)/60;
x1=zeros(1,col);
% for m=1:t
    for n=1:col
         x1(n)=sum(HDataY((n-1)*60+1:n*60))/60;
    end
% end
% mean_x1=mean(x1')';
% x1=x1-repmat(mean_x1,1,size(x1,2));  % 0 均值
x1=x1-x1(1);
% plot(x1(1,:))

v=x1'*x1;
[F D]=eig(v);
D1=diag(D);
[D2 I]=sort(D1);
[m n]=size(F);
I=I([n:-1:1]);
D2=D2(I);                %把特征值按从大到小排列
for j=1:n
    F2(:,j) = F(:,I(j));
end

figure
subplot(121);plot(F2(:,1))
% subplot(523);plot(F2(:,2))
% subplot(525);plot(F2(:,3))
% subplot(527);plot(F2(:,4))
% subplot(529);plot(F2(:,5))

A=x1*F2;
w=D2/sum(D2);
f1=A(:,1)*F2(:,1)';f2=A(:,2)*F2(:,2)';f3=A(:,3)*F2(:,3)';

for j=1:1   % for j=1:5
% figure(2)
subplot(1,2,2*j)
hold on
plot(x1(j,:))
plot(f1(j,:),'r')
plot(f2(j,:),'g')
plot(f3(j,:),'y')
plot(f1(j,:)+f2(j,:)+f3(j,:),'k')
% hold off
end

DstData=AnyDst(StartTime,EndTime);
DstData=DstData-DstData(1);
% DstData=AnyDst([99 3 4 4 0],[99 3 5 22 0]);
plot(DstData,'c')