nplot.m

用于计算脉冲激光照射下物体的温度分布

%Nplot.m
format short;
t1=num2str(t(l)*1e0);
t22=t(fix(l/5.5));
t2=num2str(t22*1e0);

z0=num2str(z(1));
z1=num2str(z(2));
z2=num2str(z(10));
z3=num2str(z(10));

r1=num2str(r(1));
r2=num2str(r(3));
r3=num2str(r(2));


subplot(2,2,1),
plot(z([1:end]),[Ntemp([L1+1:L1+L1])]'),xlabel('Depth z/(um)'),ylabel('Zn Concentration N(um^-3)'),
title('the Zn Concentration along z axis'),
%legend('The profile of zn concentration in the time t=500us')
hold on
j=1;
plot(z([1:end]),[N([L1*(j-1)+1:L1*(j-1)+L1],fix(l/5.5))]',':r'),xlabel('Depth z/(um)'),ylabel('Zn Concentration N(um^-3)'),
title('the Zn Concentration along z axis'),
text1=['t=' t1 's'];
text2=['t=' t2 's'];
legend(text1,text2),
hold off
subplot(2,2,2),
axis([0 zd 1 10^8]),
semilogy(z([1:end]),[Ntemp([L1+1:L1+L1])]'),xlabel('Depth z/(um)'),ylabel('Zn Concentration N(um^-3) in logarithm'),
title('the Zn Concentration along z axis'),
%legend('The profile of zn concentration in the time t=500us')
hold on
j=2;
axis([0 zd 1 10^8]),
semilogy(z([1:end]),[N([L1*(j-1)+1:L1*(j-1)+L1],fix(l/5.5))]',':r'),xlabel('Depth z/(um)'),ylabel('Zn Concentration N(um^-3) in logarithm'),
title('the Zn Concentration along z axis'),
text1=['t=' t1 's'];
text2=['t=' t2 's'];
legend(text1,text2),
hold off
subplot(2,2,3),
plot(r,[N([2:L1:L1*(M1-1)+2],l)]'),xlabel('r/(um)'),ylabel('Zn Concentration N(um^-3)'),
title('the Zn Concentration along z axis'),
%legend('The profile of zn concentration in the time t=500us')
hold on
plot(r,[N([2:L1:L1*(M1-1)+2],fix(l/5.5))]',':r'),xlabel('r/(um)'),ylabel('Zn Concentration N(um^-3)'),
title('the Zn Concentration along z axis'),
i=10;
plot(r,[N([i:L1:L1*(M1-1)+i],fix(l))]','-.b'),xlabel('r/(um)'),ylabel('Zn Concentration N(um^-3)'),
plot(r,[N([i:L1:L1*(M1-1)+i],fix(l/5.5))]','+b'),xlabel('r/(um)'),ylabel('Zn Concentration N(um^-3)'),
text1=['t=',t1, 's,z=',z1, 'um'];
text2=['t=' t2 's,z=' z1 'um'];
text3=['t=' t1 's,z=' z2 'um'];
text4=['t=' t2 's,z=' z2 'um'];
legend(text1,text2 ,text3,text4)
hold off
subplot(2,2,4),
plot(t([1:tn]),N(2,[1:end])),xlabel('diffusion t/(s)'),ylabel('Zn Concentration N(um^-3)'),
title('the Zn Concentration versus t'),
hold on
%j=3;
%plot(tt([1:tn]),N(L1*(j-1)+30,[1:end]),'r:'),xlabel('diffusion t/(s)'),ylabel('Zn Concentration N(um^-3)'),
j=2;
plot(t([1:tn]),N(L1*(j-1)+1,[1:end]),'g:'),xlabel('diffusion t/(s)'),ylabel('Zn Concentration N(um^-3)'),
j=2;
plot(t([1:tn]),N(L1*(j-1)+10,[1:end]),'r-.'),xlabel('diffusion t/(s)'),ylabel('Zn Concentration N(um^-3)'),
text1=['r=' r1 'um,z=' z1 'um'];
%text2=['r=' r2 'um,z=' z2 'um'];
text3=['r=' r3 'um,z=' z0 'um'];
text4=['r=' r3 'um,z=' z3 'um'];
legend(text1,text3,text4);
hold off