matlab绘制经典气体、玻色气体和费米气体的化学势随温度的变化曲线
前言
经典气体、玻色气体和费米气体的化学势随温度的变化曲线是热力学与统计物理中老生常谈的话题了,具有典型的代表性和对比性,下面简要介绍如何使用matlab绘制相关曲线,以及其算法演示。对于使用不同的画图软件和程序基本算法思路不变,但是一些程序命令需要读者根据软件的情况自行修改!
运行执行文件:
clc,clear;% Classical particle
k = 1.380649*10^(-23);
h = 1.05457266*10^(-34);
syms T
u_classcial = T*log(T^(-3/2)*(2*pi)^(3/2));
fplot(u_classcial,[0,22],'r','LineWidth',1);hold on;% Fermion
syms e2 u2
%ft = (e2.^(1/2))./(exp((e2-u2)./T)+1);
%N2 = integral(@(e2)ft,0,inf);
%fz = @(u2)N2-13.96;
t = 0.1:0.1:22;
i = 1;
u_fermi = [];
for T = 0.1:0.1:22u_fermi(i) = fzero(@(u2)integral(@(e2) (e2.^(1/2))./(exp((e2-u2)./T)+1),0,inf)-13.96,1);i = i+1;
end
plot(t,u_fermi,'b','LineWidth',1);hold on;% Boson
fz0=g(3/2,1);
Tc=fz0^(-2/3)*2*pi;
z=0.15:0.05:1;
t=[];
for i=1:length(z)
t(i)=g(3/2,z(i));
end
t=t.^(-2/3)*2*pi;
u_bose = log(z).*t;% Plot
plot(t,u_bose,'c','LineWidth',1);hold on;
plot(linspace(0,Tc,10),linspace(0,0,10),'c','LineWidth',1)
%line([Tc,Tc], [7,-41],'LineStyle','--','color', 'g','LineWidth',0.7);hold on;
%line([4,22], [0,0],'LineStyle','-.', 'color', 'm','LineWidth',0.1);
title('Change of chemical potential energy of classical, Bose and Fermi gases with temperature');
legend('Classic','Fermi','Bose','Location','best')
xlabel('Temperature (T/K)');
ylabel('Chemical potential (μ/ J*mol^(-1) )');
调用到的function文件:(玻色爱因斯坦函数)
function g=g(v,z)
syms x
g=double(symsum(z^x/x^v,x,1,inf));
end
绘制完成后的图像演示:
感谢阅读,有帮助点赞收藏支持呀~
有问题可以私信~
matlab绘制经典气体、玻色气体和费米气体的化学势随温度的变化曲线相关推荐
- MATLAB数学经典建模之风扇特性:流量 随 压比函数值变化的图形 (2 维图形)
clc % 清屏 clear all; % 删除workplace变量 close all; % 关掉显示图形窗口 format short warning off H = 11; %飞行高度 Ma ...
- nyquist图怎么画matlab,用MATLAB绘制Nyquist图
<用MATLAB绘制Nyquist图>由会员分享,可在线阅读,更多相关<用MATLAB绘制Nyquist图(9页珍藏版)>请在人人文库网上搜索. 1.用MATLAB绘制Nyqu ...
- matlab绘制圆,且求解两个圆的交点坐标
前言 这里简单记录下对于matlab绘制圆,以及求解两圆交点的坐标的实现(包括C++). 一 绘制圆 circle.m文件: function [] = circle(x, y, r, color) ...
- matlab系统的根轨迹,实验五 利用MATLAB绘制系统根轨迹
<实验五 利用MATLAB绘制系统根轨迹>由会员分享,可在线阅读,更多相关<实验五 利用MATLAB绘制系统根轨迹(6页珍藏版)>请在人人文库网上搜索. 1.实验五 利用MAT ...
- Matlab 绘制三维立体图(以地质异常体为例)
Matlab 绘制三维立体图(以地质异常体为例) 参考文章: (1)Matlab 绘制三维立体图(以地质异常体为例) (2)https://www.cnblogs.com/yangwenbo214/p ...
- matlab绘制频散曲线,Matlab绘制频散曲线程序代码.docx
Matlab绘制频散曲线程序代码.docx 下载提示(请认真阅读)1.请仔细阅读文档,确保文档完整性,对于不预览.不比对内容而直接下载带来的问题本站不予受理. 2.下载的文档,不会出现我们的网址水印. ...
- matlab绘制X,Y二维散点图并标出序号
matlab绘制X,Y二维散点图并标出序号 标签: matlab 2016-11-22 12:20 193人阅读 评论(0) 收藏 举报 分类: matlab(1) 版权声明:本文为博主原创文章, ...
- 如何使用MATLAB绘制不同类型的二维图形
如何使用MATLAB绘制不同类型的二维图形 听语音 | 浏览:360 | 更新:2014-11-30 18:10 | 标签:matlab 1 2 3 4 5 分步阅读 如何用MATLAB绘制各种二维图 ...
- 用matlab绘制抛物线y的x平方,我刚刚学习MATLAB,想画一下(x^2+y^2-1)^3=x^2*y^3这个曲线的图像,该怎么画呢,谢谢大家了。...
syms x y zezsurf(x,y,1-x-2*y)如果希望去除网格可以再加上shading interp如果想让坐标轴显示比例相等,可以再加上axis equalwww.mh456.com防采 ...
最新文章
- Collection集合List、Set
- 彩色图像--色彩空间 YIQ 、YUV 、YCbCr 、YC1C2 和I1I2I3
- 利用Sql Server2005发送邮件
- dropdownList级联刷新gridView
- 批处理学习之Bat命令——获取当前盘符、当前目录、上级目录
- spring-第八篇之容器中的bean的生命周期
- 正确理解MySQL中的where和having的区别
- oracle jde优势介绍,JDE 的扫盲知识介绍。。。
- 网上订餐php论文,php032网上订餐系统
- win7 安装英文语言包
- Sprint 敏捷开发
- 从外包、互联网到国有企业,再到研究生拟录取,三年时间让我悟出了自己人生的意义
- matlab模拟塞曼图谱,塞曼效应理论与实验分析.doc
- pc端的微信dat文件解码转换成普通图片
- MySQL数据表中的数据单表查询
- win10Edge浏览器或第三方浏览器网页字体模糊解决方法
- linux can命令详解,Linux CAN编程详解
- python中换行符用法_python换行符是什么?
- CSDN的markdown编辑器使用(含表格背景色、文本颜色字体字号、锚点跳转、表情符号等)
- oracle 19c INS-20802