%*****************************
%
%畫出相位分布的偽彩色圖
%
%*****************************
III=1:I-1;
JJJ=1:J-1;
[IX,JY]=meshgrid(JJJ,III); % 劃分網格
pcolor(IX,JY,EI(III,JJJ));% 偽彩色顯示
% caxis([0,10000])
colorbar('horiz') %顏色對應表
% hold on
% contour(IX,JY,EI,30,'k')
% hold off
% axis square   %顯示正方形網格區域
% axis equal
axis ij      
view(-90,90)  % 矯正視角
shading interp% 消除網格
% colormap(gray)