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

AlexYoung

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

 
 
 

日志

 
 

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

2009-08-11 20:54:53|  分类: 课程练习 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

经过一番化简后,得到反应性方程:...

对于给定条件下的解过程如下,借助MATLAB数学工具:

一、首先,建立一个M文件,取名为inhourequ.m,用于描述反应性方程在给定条件下的形式

function y=inhourequ(x)
y=0.0001*x/(1+0.0001*x)+1/(1+0.0001*x)*(0.000247*x/(x+0.0127)...
    +0.001385*x/(x+0.0317)+0.001222*x/(x+0.1155)...
    +0.002645*x/(x+0.3106)+0.000832*x/(x+1.3966)+0.000169*x/(x+3.8760))-0.001;
%对应组的衰变常数依次为:0.0127     0.0317     0.1155     0.3106     1.3966     3.8760
%对应组的缓发中子份额为:0.000247 0.001358 0.001222 0.002645 0.000832 0.000169
%假设的中子寿命: l=0.0001s
%初始状态有效增值因数阶跃扰动量:delta(k)=0.001

二、接下来,在MATLAB命令窗口里输入以下命令,或者是单独建立一个M文件

clc;clear;
w(1)=fzero('inhourequ',0);
w(2)=fzero('inhourequ',[-0.0316 -0.0128]);
w(3)=fzero('inhourequ',[-0.1154 -0.0318]);
w(4)=fzero('inhourequ',[-0.3105 -0.1156]);
w(5)=fzero('inhourequ',[-1.3965 -0.3107]);
w(6)=fzero('inhourequ',[-3.8759 -1.3967]);
w(7)=fzero('inhourequ',[-60 -3.8761]);
w=w'

三、所得的解便是在以上绿色标识的条件下的7个解,之后就可以进行后续相关的讨论了...

w =

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

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

历史上的今天

评论

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

页脚

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