matlab代码实现热度图（以数学建模2011A题为例）

x=xlsread('11adata.xls','sheet1','B4:B322'); y=xlsread('11adata.xls','sheet1','C4:C322'); z=xlsread('11adata.xls','sheet1','D4:D322');%高度 area=xlsread('11adata.xls','sheet1','E4:E322'); temp=xlsread('11adata.xls','sheet2'); metal=temp(:,2:end); uij=min(x(:,:)):(max(x(:,:))-min(x(:,:)))/100:max(x(:,:)); vij=min(y(:,:)):(max(y(:,:))-min(y(:,:)))/100:max(y(:,:)); [xi,yj]=meshgrid(uij,vij); figure(1); z2=griddata(x,y,metal(:,6),xi,yj,'cubic'); h=pcolor(xi,yj,z2); hold on; for i=1:length(area(:,1)) if(area(i)==1) h1=plot3(x(i),y(i),z(i),'k*'); elseif(area(i)==2) h2=plot3(x(i),y(i),z(i),'k+'); elseif(area(i)==3) h3=plot3(x(i),y(i),z(i),'kx'); elseif(area(i)==4) h4=plot3(x(i),y(i),z(i),'ko'); elseif(area(i)==5) h5=plot3(x(i),y(i),z(i),'k.'); end end for i=1:length(area(:,1)) if(area(i)==1) h1=plot3(x(i),y(i),metal(i),'r*'); elseif(area(i)==2) h2=plot3(x(i),y(i),metal(i),'y+'); elseif(area(i)==3) h3=plot3(x(i),y(i),metal(i),'gx'); elseif(area(i)==4) h4=plot3(x(i),y(i),metal(i),'mo'); elseif(area(i)==5) h5=plot3(x(i),y(i),metal(i),'k.'); end end %contour(x,y,z,6); hold off; title('金属元素浓度分布图'); set(h,'edgecolor','none','facecolor','interp'); colorbar; legend([h1,h2,h3,h4,h5],'生活区','工业区','山区','主干道路区','公园绿地区',2);