用先验获取知识的代码

for i=1:1000
ZeroV(i)=acos(dot(E_t,X(1:211,indices(i)))/(norm(E_t)*norm(X(1:211,indices(i)))))*180/pi
end
[n,im]=min(ZeroV);
cols=im; %cols 为最小值所在列数
pix_col=indices(cols)