三维荧光 衍射峰处理 画图

function [V,bl]=ygjf(eem_cor,delt) %% %清除瑞利散射与拉曼散射后的效果图 [m,n]=size(eem_cor); EX=eem_cor(1,2:n); EM=eem_cor(2:m,1); FL=eem_cor(2:m,2:n); FL=reshape(FL',m-1,n-1,1); contourf(EM,EX,reshape(FL,n-1,m-1)); colormap cool A=[250 540 200 490]; axis(A); xlabel('Em(nm)') ylabel('Ex(nm)') line([A(1) A(2)],[250 250],'LineStyle','-','color','k') line([330 330],[A(3) 250],'LineStyle','-','color','k') line([380 380],[A(3) A(4)],'LineStyle','-','color','k') text((A(1)+330)/2,(250+A(3))/2,'Ⅰ','color','w','FontSize',10,'fontangle','italic','fontweight','bold'); text((330+380)/2,(250+A(3))/2,'Ⅱ','color','w','FontSize',10,'fontangle','italic','fontweight','bold'); text((380+A(2))/2,(250+A(3))/2,'Ⅲ','color','w','FontSize',10,'fontangle','italic','fontweight','bold'); text((380+A(1))/2,(250+A(4))/2,'Ⅳ','color','w','FontSize',10,'fontangle','italic','fontweight','bold'); text((380+A(2))/2,(250+A(4))/2,'Ⅴ','color','w','FontSize',10,'fontangle','italic','fontweight','bold'); colorbar set(gcf, 'color', 'w'); figure contour(EM,EX,reshape(FL,n-1,m-1),20); axis tight colormap(jet); colorbar shading interp axis(A); xlabel('Em. (nm)') ylabel('Ex. (nm)') figure pcolor(EM,EX,reshape(FL,n-1,m-1)); line([A(1) A(2)],[250 250],'LineStyle','-','color','k') line([330 330],[A(3) 250],'LineStyle','-','color','k') line([380 380],[A(3) A(4)],'LineStyle','-','color','k') text((A(1)+330)/2,(250+A(3))/2,'Ⅰ','color','w','FontSize',10,'fontangle','italic','fontweight','bold'); text((330+380)/2,(250+A(3))/2,'Ⅱ','color','w','FontSize',10,'fontangle','italic','fontweight','bold'); text((380+A(2))/2,(250+A(3))/2,'Ⅲ','color','w','FontSize',10,'fontangle','italic','fontweight','bold'); text((380+A(1))/2,(250+A(4))/2,'Ⅳ','color','w','FontSize',10,'fontangle','italic','fontweight','bold'); text((380+A(2))/2,(250+A(4))/2,'Ⅴ','color','w','FontSize',10,'fontangle','italic','fontweight','bold'); axis tight colormap(jet); colorbar shading interp axis(A); xlabel('Em. (nm)') ylabel('Ex. (nm)') figure; surfc(EM,EX,reshape(FL,n-1,m-1)) shading interp axis tight B=[A 0 600]; axis(B); xlabel('Em(nm)') ylabel('Ex(nm)') view(0,90) colorbar set(gcf, 'color', 'w'); figure; view(-100,30) surfc(EM,EX,reshape(FL,n-1,m-1)) shading interp axis tight xlabel('Em. (nm)') ylabel('Ex. (nm)') zlabel('fl') colorbar set(gcf, 'color', 'w'); figure; view(-100,30) surfc(EM,EX,reshape(FL,n-1,m-1)) shading interp axis tight camlight left lighting phong xlabel('Em. (nm)') ylabel('Ex. (nm)') zlabel('fl') colorbar set(gcf, 'color', 'w'); figure; contour(EM,EX,reshape(FL,n-1,m-1)) axis tight xlabel('Em. (nm)') ylabel('Ex. (nm)') zlabel('fl') colorbar axis(A); line([A(1) A(2)],[250 250],'LineStyle','-','color','k') line([330 330],[A(3) 250],'LineStyle','-','color','k') line([380 380],[A(3) A(4)],'LineStyle','-','color','k') text((A(1)+330)/2,(250+A(3))/2,'I','color','k','FontSize',10,'fontangle','italic','fontweight','bold'); text((330+380)/2,(250+A(3))/2,'II','color','k','FontSize',10,'fontangle','italic','fontweight','bold'); text((380+A(2))/2,(250+A(3))/2,'III','color','k','FontSize',10,'fontangle','italic','fontweight','bold'); text((380+A(1))/2,(250+A(4))/2,'IV','color','k','FontSize',10,'fontangle','italic','fontweight','bold'); text((380+A(2))/2,(250+A(4))/2,'V','color','k','FontSize',10,'fontangle','italic','fontweight','bold'); h1=annotation('textarrow',[0.27 .3],[.2 .25]); set(h1,'string','PR','color','k','FontSize',10,'fontangle','italic','fontweight','bold'); h2=annotation('textarrow',[0.29 .32],[.5 .55]); set(h2,'string','FL','color','k','FontSize',10,'fontangle','italic','fontweight','bold'); set(gcf, 'color', 'w'); %% for i=1:1:m-1 if EM(i)==A(1) m1=i+1; elseif EM(i)==A(2) m2=i+1; end end for j=1:1:n-1 if EX(j)==A(3) n1=j+1; elseif EX(j)==A(4) n2=j+1; end end new_eem=eem_cor(m1:m2,n1:n2); mm1=(330-A(1))/delt; nn1=(250-A(3))/10; mm2=(380-330)/delt+mm1; nn2=(250-A(3))/10; mm3=(A(2)-380)/delt+mm2; nn3=(250-A(3))/10; mm4=(380-A(1))/delt; nn4=(A(4)-250)/10; mm5=(A(2)-380)/delt+mm4; nn5=(A(4)-250)/10; V1=0; V2=0; V3=0; V4=0; V5=0; for i=1:1:mm1 for j=1:1:nn1 V1=V1+(new_eem(i,j)+new_eem(i,j+1)+new_eem(i+1,j)+new_eem(i+1,j+1))/4*delt*10; end end for i=mm1+1:1:mm2 for j=1:1:nn2 V2=V2+(new_eem(i,j)+new_eem(i,j+1)+new_eem(i+1,j)+new_eem(i+1,j+1))/4*delt*10; end end for i=mm2+1:1:mm3 for j=1:1:nn3 V3=V3+(new_eem(i,j)+new_eem(i,j+1)+new_eem(i+1,j)+new_eem(i+1,j+1))/4*delt*10; end end for i=1:1:mm4 for j=nn3+1:1:nn4 V4=V4+(new_eem(i,j)+new_eem(i,j+1)+new_eem(i+1,j)+new_eem(i+1,j+1))/4*delt*10; end end for i=mm4+1:1:mm5 for j=nn3+1:1:nn5 V5=V5+(new_eem(i,j)+new_eem(i,j+1)+new_eem(i+1,j)+new_eem(i+1,j+1))/4*delt*10; end end bl=zeros(5,1); V=[V1 V2 V3 V4 V5]'; for i=1:1:5 bl(i)=V(i)/sum(V); end