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

AlexYoung

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

 
 
 

日志

 
 

阶跃扰动时点堆模型动态方程的解(3)  

2009-08-12 13:45:35|  分类: 课程练习 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

%点堆模型动态方程在反应性阶跃引入的解
tic;clc;clear;
l=0.0001;%l--the neutron lifetime(中子寿命)
lambda=[0.0127   0.0317   0.1155   0.3106   1.3966   3.8760];

%λ--the decay constant(对应组的衰变常数)
beta=[0.000247 0.001358 0.001222 0.002645 0.000832 0.000169];

%β--the delayed neutron fraction(对应组的缓发中子份额)
delta=0.001;

%δ--the initial reactivity insertion(初始状态有效增值因数阶跃扰动量)
%X=l*omega/(1+l*omega)+omega/(1+l*omega).*(beta./(omega+lambda));
%rho=sum(X);%ρ(t)--随时间变化的反应性
inhourequ=inline('sum(l*omega/(1+l*omega)+omega/(1+l*omega).*(beta./(omega+lambda)))-delta','omega');
%================================================================
%ρ=lω/(1+lω)+ω/(1+lω)Σ[β/(ω+λ)]--the reactivity equation
%================================================================
omega(1)=fzero('inhourequ',0);
omega(2)=fzero('inhourequ',[-lambda(2)+eps -lambda(1)-eps]) ;
%omega(2)=fzero('inhourequ',[-0.0316 -0.0128]);
omega(3)=fzero('inhourequ',[-lambda(3)+eps -lambda(2)-eps]) ;
%omega(3)=fzero('inhourequ',[-0.1154 -0.0318]);
omega(4)=fzero('inhourequ',[-lambda(4)+eps -lambda(3)-eps]) ;
%omega(4)=fzero('inhourequ',[-0.3105 -0.1156]);
omega(5)=fzero('inhourequ',[-lambda(5)+eps -lambda(4)-eps]) ;
%omega(5)=fzero('inhourequ',[-1.3965 -0.3107]);
omega(6)=fzero('inhourequ',[-lambda(6)+eps -lambda(5)-eps]) ;
%omega(6)=fzero('inhourequ',[-3.8759 -1.3967]);
omega(7)=fzero('inhourequ',[-1e3 -lambda(6)-eps]) ;
%omega(7)=fzero('inhourequ',[-60 -3.8761]);
omega=omega'
toc;
%omega =

%    0.0188
%   -0.0141
%   -0.0622
%   -0.1867
%   -1.2097
%   -3.7650
%  -55.5793

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

历史上的今天

评论

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

页脚

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