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

AlexYoung

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

 
 
 

日志

 
 

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

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

  下载LOFTER 我的照片书  |

Contents

  • 初始化,对已知量赋值
  • 多普勒函数值
  • 多普勒展宽截面
% 多普勒函数的数值积分及多普勒展宽截面的求取
% 用于计算截面中的多普勒展开截面

初始化,对已知量赋值

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
 


Published with MATLAB? 7.13

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

历史上的今天

评论

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

页脚

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