benzhengmo.m

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

load ('1999090105.txt');X1999090105=X1999090105';X1999090105=X1999090105(:); 
load ('1999090610.txt');X1999090610=X1999090610';X1999090610=X1999090610(:); 
load ('1999091115.txt');X1999091115=X1999091115';X1999091115=X1999091115(:); 
load ('1999091620.txt');X1999091620=X1999091620';X1999091620=X1999091620(:); 
load ('1999092125.txt');X1999092125=X1999092125';X1999092125=X1999092125(:); 
load ('1999092630.txt');X1999092630=X1999092630';X1999092630=X1999092630(:); 

x1=zeros(30,24);
for m=1:5
    for n=1:24
        x1(m,n)=sum(X1999090105((m-1)*1440+(n-1)*60+1:(m-1)*1440+n*60))/60;
    end
end
for m=6:10
    for n=1:24
        x1(m,n)=sum(X1999090610((m-6)*1440+(n-1)*60+1:(m-6)*1440+n*60))/60;
    end
end
for m=11:15
    for n=1:24
        x1(m,n)=sum(X1999091115((m-11)*1440+(n-1)*60+1:(m-11)*1440+n*60))/60;
    end
end
for m=16:20
    for n=1:24
        x1(m,n)=sum(X1999091620((m-16)*1440+(n-1)*60+1:(m-16)*1440+n*60))/60;
    end
end
for m=21:25
    for n=1:24
        x1(m,n)=sum(X1999092125((m-21)*1440+(n-1)*60+1:(m-21)*1440+n*60))/60;
    end
end
for m=26:30
    for n=1:24
        x1(m,n)=sum(X1999092630((m-26)*1440+(n-1)*60+1:(m-26)*1440+n*60))/60;
    end
end
% plot(x1(1,:))
mean_x1= mean(x1')';
x1=x1-repmat(mean_x1,1,size(x1,2));  % 0 均值
% 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(521);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:5   % for j=1:5
% figure(2)
subplot(5,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