matlab绘制玻色-爱因斯坦凝聚的内能曲线_matlab bose-einstein distribution-CSDN博客

本文链接：https://blog.csdn.net/m0_46808930/article/details/129686095

前言

玻色–爱因斯坦凝态（Bose-Einstein condensate）又称玻色–爱因斯坦凝聚态，是玻色子原子在冷却到接近绝对零度所呈现出的一种气态的、超流性的物质状态（物态）。这种大量具有玻色统计性质的粒子，如同原子“凝聚”到同一状态，称为玻色–爱因斯坦凝聚（Bose-Einstein condensation，BEC）。

执行文件

clc,clear;
% BEC
syms z n T;
y1 = symsum(1/(n^(5/2)),n,1,1000);
y2 = symsum(1/(n^(3/2)),n,1,1000);

U1 = (3/2)*T^(5/2)*(y1/y2);   
U_class = 3/2*T;

xx1 = 0:0.1:1;
xx2 = 2:0.2:3;
xx = [xx1,xx2]; 
y11 = [];
y22 = [];
j = 1;
p = 1;
for i=0:0.1:1
    y11(j) = (3/2)*i^(5/2)*(y1/y2);
    j = j + 1;
end

for k = 2:0.2:3
   y22(p) = (3/2)*k;
   p = p + 1;
end
yy = [y11,y22];

xi = 1:0.2:2;
yi = interp1(xx,yy,xi,'spline'); 

% 绘图
plot(xi,yi,'r','LineWidth',1.5),hold on;
fplot(U1,[0,1],'r','LineWidth',1.5);hold on;
fplot(U_class,[0,3],'b','LineWidth',1.5);hold on;
legend('BE','BE','Classical particle','Location','northwest')
title('Internal energy of Bose Einstein condensation');
xlabel('(T/Tc)/K');
ylabel('(U)/NK');