注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

AlexYoung

做好大家都必须要做的事,发展自己感兴趣的事,结束圣人所不齿的事...

 
 
 

日志

 
 

多普勒展宽截面的求取@课程练习  

2013-06-26 00:03:16|  分类: 课程练习 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

初始化,对已知量赋值.............................................................................................. 1

多普勒函数值............................................................................................................ 1

多普勒展宽截面........................................................................................................ 2

 

% 多普勒函数的数值积分及多普勒展宽截面的求取

% 用于计算截面中的多普勒展开截面

 

初始化,对已知量赋值

clear                                  %清除变量

k=1.38E-23/(1.6*10^-19);               %玻尔兹曼常数, 原本单位为J/K,此处需要转换为eV/K

u=1.66E-27;                            %原子质量单位, kg

A=238;                                 %U238相对原子质量

CapitalGamma=0.02752;                   %6.67eV处的U238原子的总能级宽度

En=6.67;

E0=6.67;

T = 293;

%Ea=k*T;  %最可几能量,在麦克斯韦分布中的一个最可几能量值

%Ec=En+Ea;

%y=2/CapitalGamma*(Ec-E0);

delta=sqrt(4*k*T*En/A);  %多普勒宽度

Theta=CapitalGamma/delta %多普勒函数的第一个参数

x=2/CapitalGamma*(En-E0) %多普勒函数的第二个参数

 

 

Theta =

 

    0.5170

 

 

x =

 

     0

 

 

 

多普勒函数值

psi=quadl(@(s)beijihanshu(s,Theta,x),-100,100)

 

 

psi =

 

    0.3501

 

 

 

多普勒展宽截面

DopplerSigma=20000*0.026/(0.026+0.00152)*sqrt(E0/En)*psi;

vpa(DopplerSigma,5)

 

 

ans =

 

6615.4

 

 

 

 

  评论这张
 
阅读(129)| 评论(0)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2018