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

AlexYoung

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

 
 
 

日志

 
 

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

2009-08-12 11:51:58|  分类: 课程练习 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

当然如果使用内联函数的话,就可以不用建立两个M文件了,现在可以单建立一个运行M文件即可:

clc;clear;
%对应组的衰变常数依次为: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

%建立内联函数,表征反应性方程
inhourequ=inline('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');
%对内联函数进行求解
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'

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

历史上的今天

评论

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

页脚

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