统计物理双原子分子定容热容图像（平动热容、转动热容、振动热容）

考虑经典极限条件下的双原子分子定容热容（经典统计求解）

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');

最终出图的结果如下：（最后将图像转化为对数坐标即可）