dcoefficientvst.m

来自「用于计算脉冲激光照射下物体的温度分布」· M 代码 · 共 24 行

%由扩散系数与不同温度下浓度的关系得到扩散系数与温度关系
load D_Dcoeff.mat
load N_Dcoeff.mat
load T_Dcoeff.mat
mark=['+' '*' 'o' 's' '+' '*' 'o' '-'];
tt=[300:10:1000]+273.15;
T=T++273.15;
text=['','','','',''];
for i=1:4,
   
 p=polyfit(1./T,log10([D(:,i)]'),1);
 pp(i,:)=p;
 plot(T,D(:,i),mark(i),tt,10.^(polyval(pp(i,:),1./tt,'b')));
 hold on
% text(i,50)=num2str(N(i))
end
title('The diffusion coefficient D versus T in different Zn concentration '),Xlabel('T/K'),Ylabel('D/cm^2/s')
hold off
text1=[num2str(N(1)) '(1/cm^3)'];
text2=[num2str(N(2)) '(1/cm^3)'];
text3=[num2str(N(3)) '(1/cm^3)'];
text4=[num2str(N(4)) '(1/cm^3)'];
legend(text1,'',text2,'',text3,'',text4,'');