Glauber模型数值积分方法的matlab实现

JuST2333333

已于 2022-11-19 11:37:54 修改

阅读量441

点赞数 1

文章标签： matlab

于 2022-11-19 11:37:01 首次发布

本文链接：https://blog.csdn.net/qq_50452777/article/details/127934632

版权

本文介绍了一篇关于粒子物理课程实践作业的内容，重点是使用matlab实现Glauber模型的数值积分方法。作业涉及对Ncoll和Npart公式的二重积分，该积分过程参考了特定的数值积分技术，如Simpson公式，并假设inelastic nucleon-nucleon cross section为42。实验结果显示了积分的计算效果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

这是粒子物理课程的实践作业

其实就是对着Ncoll和Npart的公式数值积分，参考这篇的2.3部分

Glauber Modeling in High-Energy Nuclear Collisions | Annual Review of Nuclear and Particle Science (annualreviews.org)

这里的二重积分稍微有点麻烦，参考的这篇

二重积分的Simpson公式及其误差估计.pdf (book118.com)

clc;clear;close all
r=linspace(0,20,100);
rhoA1=1./(1+exp((r-6.38)./0.535));
rhoB1=1./(1+exp((r-6.38)./0.535));%画图观察过，这俩在15fm左右就降为零了
A=197;
B=197;

%%%开始归一化
funcA=@(r) 1./(1+exp((r-6.38)./0.535)).*4.*pi.*r.^2;
funcB=@(r) 1./(1+exp((r-6.38)./0.535)).*4.*p