load('convex_subpixel_ptmean_03_SNR_30_sigma_1_Nb_2_Train')
load('最终vca提取mat')
duanyuan=[E_t,E_minus]
wavelength=wavelength/1000;
sum=0
for i=1:4
sum=sum+acos(dot(duanyuan(1:211,i),end_member(1:211,i))/(norm(duanyuan(1:211,i))*norm(end_member(1:211,i))))*180/pi
end
plot_sig( 'VCA',end_member,duanyuan,wavelength,101,1)