考虑经典极限条件下的双原子分子定容热容(经典统计求解)
clc,clear;
% 双原子分子热熔
% 平动热熔
Cvt = 3/2;
% 转动热熔
syms l t;
Zr = symsum((2*l+1)*exp(-30*l*(l+1)/t),l,0,1000);
temp = diff(log(Zr));
Ur = t^2*temp;
Cvr = diff(Ur,t,1);
%振动热熔
syms n t;
Zv = symsum(exp(-(3000*(n+1/2)/t)),n,0,1000);
temp2 = diff(log(Zv));
Uv = t^2*temp2;
Cvv = diff(Uv,t,1);
% 绘图
Cv_fin = Cvt + Cvr + Cvv;
fplot(Cv_fin,[0,5000],'r','LineWidth',2);hold on;
fplot([2.5,2.5],[0,5000],'b--','LineWidth',1.5);hold on;
fplot([3.5,3.5],[0,5000],'b--','LineWidth',1.5);hold on;
line([30,30], [1.5,3.6], 'color', 'g','LineWidth',1);hold on;
line([3000,3000], [1.5,3.6], 'color', 'g','LineWidth',1);
title('双原子分子定容热容');
xlabel('温度/K');
ylabel('定容热容/NK');
最终出图的结果如下:(最后将图像转化为对数坐标即可)