svd_sdosvd.m

来自「svd分析的例子」· M 代码 · 共 36 行

 %奇异值分解子程序 sdosvd.m 
sum=0; 
 
[eigenassay,singularvalue,eigengene]=svd(normdata,0); %svd  的紧凑格式 
eigengene=eigengene';%A=usv' 
 
for i=1:n(2) 
  singularv(i)=singularvalue(i,i); 
  sum=sum+singularv(i); 
end 
singularv=singularv/sum; 
clear sum 
 
i=1:n(2); 
subplot(2,2,1) 
bar(i,singularv) 
xlabel('component') 
ylabel('relative variance') 
title('(a)') 
clear i 
 
time=0:10:times; 
subplot(2,2,2) 
plot(time,eigengene(1,:),'o-') 
xlabel('time(min)') 
title('(b) first eigengene') 
 
subplot(2,2,3) 
plot(time,eigengene(2,:),'o-') 
xlabel('time(min)') 
title('(c) second eigengene') 
 
subplot(2,2,4) 
plot(time,eigengene(3,:),'o-') 
xlabel('time(min)') 
title('(d) third eigengene')