用surf画三维图
figure
min_v=min(min(BSR_all));
max_v=max(max(BSR_all));
% step=(max_v-min_v)/1000;
surf(ball,Nall,BSR_all)
xlabel('Fractional bandwidth')
ylabel('Number of antennas')
zlabel('BSR')
xlim([ball(1) ball(end)])
ylim([Nall(1) Nall(end)])
grid on
box on
title('BSR')
colormap Jet
hc=colorbar;
set(hc,'ylim',[min_v max_v])