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

AlexYoung

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

 
 
 

日志

 
 

热中子裂变时的裂变中子能谱  

2009-05-18 13:43:26|  分类: 课程练习 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

tic
clear;clf;clc;
E=linspace(0,10,101);
f=0.453*exp(-1.036*E).*sinh(sqrt(2.29*E));
plot(E,f);grid on;
xlabel('能量E/MeV');
ylabel('{\chi}(E)/MeV^{-1}');
title('\fontsize{18}\fontname{隶书}^{235}U核热中子裂变时的裂变中子能谱');
legend('裂变中子能谱分布函数\chi(E)=0.453e^{-1.036E}sinh[(2.29E)^{0.5}]');
f1=inline('0.453*exp(-1.036*x).*sinh(sqrt(2.29*x))');%建立内联函数;
I=quadl(f1,0,10)
%求解数值积分,验证此概率密度函数的正确性,结果为I=1.0001符合概率密度函数的完整性;

%但是在数值积分时上式中的积分限不能够为无穷大,若为inf时显示的结果为I=NaN;

f2=inline('x.*(0.453*exp(-1.036*x).*sinh(sqrt(2.29*x)))');
M=quadl(f2,0,10)%结果为M=1.9717;


%下面以符号函数的积分验证该概率密度函数的正确性以及裂变中子的平均能量

syms e
q=0.453*exp(-1.036*e).*sinh(sqrt(2.29*e));
I2=numeric(int(q,0,inf))%I2=1.0012
M2=numeric(int(e*q,0,inf))%M2=1.9836
%通过比较还是符号函数积分更准确;
toc

热中子裂变时的裂变中子能谱 - 青楼薄名 - yangwankui1573的博客

 

  评论这张
 
阅读(328)| 评论(1)
推荐 转载

历史上的今天

评论

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

页脚

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