solve_two_temperature.m

激光与物质的相互作用 双温方程算法 matlab源程序代码

function dy =Solve_Two_Temperature(t,y)
%t 代表时间变化量
%y 代表所求的电子及晶格温度
dy = zeros(2,1);    % a column vector
g=4.2e17;%耦合系数
Ce=194;%电子热容
Ci=3.6e6;%晶格热容
R=0.431;%反射率
delta=24e-9;%趋肤深度
tao=70e-15;%脉宽
J=4e-9/(pi*20e-12);%能量密度
A=-g/Ce;%系数
 
S=0.94*(1-R)/(tao*delta)*J*exp(-4*log(2)*(t/tao-2).^2);%体激光源项
%dt = zeros(2,1); 
% dt(1)=0.0;
% dt(2)=700e-15;

% S(1)=0.94*(1-R)/(tao*delta)*J*(exp(-4*log(2)*(t/tao-2).^2)+exp(-4*log(2)*((t-dt(1))/tao-2).^2));%双脉冲体激光源项
% S(2)=0.94*(1-R)/(tao*delta)*J*(exp(-4*log(2)*(t/tao-2).^2)+exp(-4*log(2)*((t-dt(2))/tao-2).^2));
%dy(1) =A*1/y(1)*(y(1)-y(2))+S(1);%方程1 ；y1 电子温度；y2 晶格温度
dy(1) =A*1/y(1)*(y(1)-y(2))+S;
dy(2) =g*(y(1)-y(2))/Ci;%方程2