搜索
人大经济论坛 标签 模型 相关日志

tag 标签: 模型经管大学堂:名校名师名课

相关日志

分享 金融计量学—统计量
accumulation 2015-4-25 13:38
AIC最小原则是判定模型好坏标准之一,犹如R2(R平方)一样。 AIC和SC(舒瓦茨信息)常常一并作为判断模型拟合程度的标准之一,特别是在滞后阶数的选择上。比如说,一个VAR(向量自回归模型),经济理论往往无法确定滞后阶数,这时往往采用AIC或者SC最小原则,即观察不同的阶数的VAR模型,哪个模型的AIC或者SC值最小就选用哪个模型进行分析。 AIC、SC都会在模型参数中给出。除了R2、AIC、 SC之外,常用的判断标准还有Lg(极大似然法则)等。这些法则主要用在同一模型不同滞后阶数选择的判断上。 AIC和BIC是同一个指标,一般用于选择模型,也就是模型的比较优劣 他们的不同之处在于 AIC=-2 ln( L ) + 2 k 中文名字:赤池信息量 akaike information criterion BIC= -2 ln( L ) + ln(n)*k 中文名字:贝叶斯信息量 bayesian information criterion HQ= -2 ln( L ) + ln(ln(n))*k hannan-quinn criterion 构造这些统计量所遵循的统计思想是一致的,就是在考虑拟合残差的同时,依自变量个数施加“惩罚”。 但,倘若因此就说它们是同一个指标,恐怕还是有些不妥,毕竟“惩罚”的力度还是不尽相同的。 此外,这些信息量的用途不仅限于选择模型,还能用于选取合适的变换等等。而在那些时候,这些信息量又是另一个模样,也就是说它们有许多变体。因此,它们也被称为AIC准则、BIC准则等等。它们中的每一个体现的都是一系列的标准,而非单独的一个简单式子。
个人分类: 金融工程|0 个评论
分享 裂变模型—IMSL程序库
accumulation 2015-4-17 00:33
裂变模型—IMSL程序库
DBCONF函数 DBCONF函数
个人分类: 裂变模型|0 个评论
分享 裂变模型—DBCONF函数
accumulation 2015-4-17 00:29
1. 编译器本身的问题(FORTRAN编译器是否存在计算的上限); 2. 调用其它 IMSL 函数(IMSL函数中计算约束下最小值问题的函数共有8个); 3. 有约束和无约束(如果用无约束的最值函数,是否会产生物理问题); 4. 自己编(算法在FORTRAN说明书中已有); 5. LINUX 是否能调用 IMSL ; 6. C++ —— FORTRAN(通过修改已有的C++、Fortran函数构造出FORTRAN的约束下最值问题函数) ;
个人分类: 裂变模型|0 个评论
分享 美林投资时钟与经济周期模型
accumulation 2015-4-14 14:01
1. 在经济周期的四个阶段中,股票与债券的不同表现;增长下降(经济衰退)—债券(贴现率)优于股票;增长上升(经济复苏与增长)—股票优于债券; 2. 核心结论:增长与通胀变动方向共可组成四种组合,增长变动与股票价值变动同向,与债券价值变动方向无关;通胀变动与股票价值变动反向,也与债券价值变动反向;(美林的投资时钟 - 统计结果) 3. 公理化推演的结果—在增长向下通胀向下阶段,股票受益于通胀的下行但同时受制于增长,债券确定性的受益于通胀下降,因而此阶段股票的表现要差于债券,此阶段对应于投资时钟的衰退阶段; 4. 在增长向上通胀向下的阶段,股票同时受益于增长与通胀的变化,债券仅受益于通胀的下降,因而此阶段股票的表现优于债券,此阶段对应于投资时钟的复苏阶段; 5. 当增长向上通胀向上的阶段,股票受益于增长但同时受制于通胀,债券仅受通胀的影响。因而此阶段股票的表现要优于债券,此阶段对应于投资时钟的过热阶段; 6. 当增长向下通胀向上的阶段,股票同时受制于增长与通胀,债券仅受通胀的影响,因而此阶段股票的表现仍要差于债券,此阶段对应于投资时钟的滞涨阶段;
个人分类: 金融工程|0 个评论
分享 裂变模型
accumulation 2015-4-10 14:21
裂变模型
势能曲面—搜索最优裂变路径 裂变路径 Brosa模型—多模式无规颈断裂模型—质量分布 裂变后现象 激发能 丰中子核 结合能的宏观模型 壳修正的HF方法 对称裂变与非对称裂变 轻锕系元素裂变的三种模式
个人分类: 裂变模型|0 个评论
分享 计算物理基础
accumulation 2015-4-8 01:36
1.建立模型—计算方法—语言编程—结果分析; 2.物理学、数值分析、计算机程序设计; 3.数值计算方法; 4.s=solve('cos(2*x)+sin(x)=1') 5.syms a b c x f=a*x^2+b*x+c solve(f) 6.解常微分方程:dsolve('Dx=-a*x') 7.画正弦曲线:x=0:0.1:6; A= ; B=sin(A); plot(x,B(1,:),x,B(2,:)) 8.矩阵的输入:A= 求逆:inv(A); 9.AX=b: A= ; b= ; X=A \ b 10.sum det prod eig max min median sort mean std trace rank cumsum cumprob kron cross
个人分类: 金融学|0 个评论
分享 液滴模型—STATA
accumulation 2015-4-7 02:14
. gen delta=1 if mod(A,2)==0 mod(A-Z,2)==0 (255 missing values generated) . replace delta=0 if mod(A,2)==1 (162 real changes made) . replace delta=-1 if mod(A,2)==0 mod(A-Z,2)==1 (93 real changes made) . gen dripmodel=15.835*A-18.33*A^(2/3)-0.714*Z^2/A^(1/3)-92.8*(A/2-Z)^2/A+11.2*delta/A^(1/2)
个人分类: 原子核物理|0 个评论
分享 CAPM模型
accumulation 2015-4-6 18:54
3.1 决定资产期望回报的是边际风险 马可维兹均值方差分析的一个重要理念是“有得必有失”。 换句话说,风险与收益是对等的。由于投资者喜欢更高的期望回报率,但厌恶更高的回报率波动率。因此,一个很直观的想法是, 波动率(或是波动方差)越高的资产,期望回报率越高(当前价格越低)。 但CAPM 实际上证伪了这一观点。 决定资产期望回报率的不是资产回报率的波动率,而是资产回报率与市场组合波动的相关性。 为什么会有这样违反直觉的结论? 原因在于决定价格的不是总效用,而是边际效用。 也就是说, 决定资产期望回报率的是资产的边际风险,而不是总风险。 水对人非常重要,但为什么价格那么低?原因是水的数量非常庞大,从而将水带给人的边际效用压得非常低。所以水的价格也就很低。而从前面的证明中我们看到, 一项资产添加进市场组合后,带给投资者的边际效用取决于该资产的期望回报率和它与市场组合的协方差。 所以, 在资产定价的时候,决定性的因素是资产与市场组合的协方差,而非资产本身的回报率方差。一项资产可能本身回报率的波动非常大,但它和市场组合的相关性很小。这样它的期望回报率就会比较低。相应的,当前的资产价格就会比较高。 现在我们可以回答上节课一开始提出的一个随堂思考问题。一个汽车制造商和一个新药开发商, 如果预期利润一样,谁的股价更高?在没有学过CAPM 之前,我们的回答可能是汽车制造商的股价更高。 因为相比新药开发商,汽车制造商利润的确定性更大,风险更小。由于人们厌恶风险,所以会给汽车制造商更高的股票估价。 但在学习了CAPM 之后,我们知道正确的答案应该是新药开发商的股价更高。因为药品开发带来的利润波动率与整个经济,也就是市场组合的相关性不大,会远远小于汽车制造商利润波动性与整个经济的相关性。 因此,新药开发商会有更低的期望回报。相应的,它当前的股价估值就应该高于汽车制造商。所以, 在重庆啤酒乙肝疫苗开发失败之前,投资者将重庆啤酒的股价推高到令人乍舌的程度是有道理的。
个人分类: 金融工程|0 个评论
分享 CAPM模型
accumulation 2015-4-6 18:42
1. 均衡时的市场组合 通过上节课的分析我们知道, 所有理性投资者都应该以市场组合M 的形式持有风险资产。 两个问题自然产生 :(1)在均衡的时候,组合M 是什么样的?(2)在均衡的时候,投资者这样的行为会产生什么样的资产价格状态? 问题1 的答案相对简单。 市场组合M 包含了市场上所有的风险资产。其中各类资产的权重占比就是市场上各类风险资产价值的比率。换言之,市场组合就是整个市场。 有人可能会认为这个结果相当惊人。因为根据资产的期望收益与波动状况,完全有可能我们用均值方差分析方法计算出来的组合权重与市场中各类资产的组合比重不一致。甚至,我们有可能计算出来某些资产的组合权重为0 甚至为负数。但这样的思路还停留在个人选择的层面上。 最终的经济状态是由经济均衡所决定的。如果某些资产的收益波动状况不太好,让投资者对其购买量(决定于投资者的投资组合权重)低于其供给量,那就意味着这类资产存在着供过于求的问题,它的价格一定会下跌。而价格的下跌会提升其预期回报率。价格会下跌至正好平衡这类资产的供给与需求水平。 所以,在均衡时,各类资产的价格(回报)会处在恰当的水平,保证每样资产在整个资产市场中的比重就是理性投资者组合优化所会选择的组合权重。这样一来,投资者会选择的市场组合就变成了整个市场。 这个逻辑简单说起来,就是如果市场组合与整个市场不同,则某些资产的供需一定不平衡,从而引发价格调整,最终使得市场组合一定会等同于整个市场。 这就回答了前面随堂思考题3 中提出的问题。思路走到这里,就自然地把我们引到了前面提出的问题2。我们可以合理地猜测, 在均衡的时候,不同资产的价格(回报率)之间一定存在着某种特定的关系。否则风险资产的市场就无法出清。 CAPM 表明,这一猜测是正确的。以下,我们来严格证明这一点。
个人分类: 金融工程|0 个评论
分享 金融计量模型
accumulation 2015-4-4 19:05
金融计量模型的比较:如果这些模型可以解释数据反映的现象,就比较AIC或BIC指标。 (1)自回归模型(Autoregressive Model):根据预测变量的历史数据,通过自回归进行预测,即y(t)=β0+β1*y(t-1)+β2*y(t-2)+u(t),模型中含有n个历史数据,则该自回归模型称为n阶自回归模型;相应的参数β1、β2等称为冲击倾向(冲击乘数);自回归模型的根据是时间序列的相关性,即记忆性;自回归模型的预测能力很强; (2)分布滞后模型(Distributed Lag Model):用与被解释变量y相关的解释变量x的时间序列数据对y进行预测,即y(t)=β0+β1*x(t)+β2*x(t-1)+u(t),例如,y(t)代表时期t的经济增长率,x(t)代表时期t的投资量,分布滞后模型相当于将静态模型中的解释变量替代为解释变量的历史时间序列;分布滞后模型具有明确的经济学、金融学理论解释; (3)自回归分布滞后模型(Autoregressive Distributed Lag Model):将以上两种模型结合起来,即构造模型y(t)=β0+β1*x(t)+β2*x(t-1)+β3*y(t-1)+β4*y(t-2)+u(t),自回归分布滞后模型综合了自回归模型和分布滞后模型的优点; (4)GARCH模型(General Autoregressive Conditional Heteroskedasticity Model):用长期平均方差、n-1期及以前的方差估计值与n-1期及以前的残差平方去估计n期的方差的模型,模型的解释变量中含有p个方差估计量,q个残差平方,则该模型称为GARCH(p,q);GARCH模型常用于估计波动率。
个人分类: 金融工程|0 个评论
分享 线性回归模型的建模假设(LINE)
yatouha 2015-4-2 12:25
• 线性 LINEARITY : 因变量均数 m 与 X 间呈直线关系 • 独立 INDEPENDENCE 每一观察值之间彼此独立(在给定自变量 x 的条件下,残差的条件期望值为零,本假设又称零均值假设);该条件实际上隐含着:各随机误差项互不相关;随机误差项与相应的自变量 x 不相关。 • 给定 X 时, Y 正态分布 ( Normal ) : 对于任何给定的 X, Y 服从正态分布, 均数为 m Y|X , 标准差为 s Y|X , 该条件实际隐含着误差项服从正态分布(随机误差,即残差服从均值为零,方差为 s 2 的正态分布)。 • 等方差 (Equal variance) 或标准差相等: 对于任何 X 值,随机变量 Y 的标准 差 s Y|X 相等;隐含着不论 x 取何值, y 都具有相同的方差(对于所有的自变量 x ,残差 的条件方差为 s 2 ,且 s 为常数)。
个人分类: 统计|33 次阅读|0 个评论
分享 裂变物理模型问题
accumulation 2015-3-27 14:19
1.用球坐标描写的形变参量 2.用柱坐标描写的形变参量(Funny Hills) 3.小形变描述—轴对称椭球 4.卵形线 5.非轴对称形状 6.形变参量β 7.液滴模型计算的裂变位垒高度公式 8.有限力程模型 9.位能曲面的壳修正与HF方法 10.位能曲面的微观计算
个人分类: 裂变模型|0 个评论
分享 有限力程模型
accumulation 2015-3-27 11:52
1.原子核宏观模型与微观修正 2.小液滴模型:考虑库仑能和压缩能引起的核密度的不均匀性对能量公式修正的液滴模型 3.核的能量:核子间相互作用的位能和核子的动能 4.表面能的两种具体体现: (1)量子力学原理测不准关系的边界约束引起表面核子动能的增加 (2)外界并无核子,减少了核力的作用,导致位能的增加 5.位能在表面的减少是与表面的形状有关的 6.鞍点形状接近断点,中间有一个很细的颈部,这些地方表面能较小,用同一表面张力计算,就高估了表面能,因而提高了裂变位垒 7.表面能的力程修正 8.库仑能弥散层的修正 9.β稳定线上原子核裂变位垒的计算 10.汤川阱+指数阱的核势能形式
个人分类: 裂变模型|0 个评论
分享 液滴模型的位能曲面
accumulation 2015-3-27 11:42
1.位能曲面关于形变参量的函数 2.可裂变参量 3.核能量随形变的变化与宏观模型的选择 4.Myers-Swiatecki质量公式 5.描述核形变的球坐标模型—球谐函数展开 6.表面能在核为球形时最小,随形变而增大;而库仑能则与此相反,球形核库仑能最大,随形变而减少;正是这两种相互对抗的因素的竞争,导致裂变位垒和鞍点的形成 7.四极形变、八极形变与十六极形变 8.原子核能否连续演化到中子星? 9.极值:峰、谷与鞍点—Businaro-gallone点 10.裂变位垒与核对裂变不稳定 11.液滴模型计算的裂变位垒高度 12.液滴模型计算裂变位垒的三个缺点
个人分类: 裂变模型|0 个评论
分享 裂变势能曲面的计算
accumulation 2015-3-26 18:41
本工作从裂变势能曲面和Brosa 模型两方面入手,研究与核能应用相关的锕系元素的核数据理论计算问题。 用宏观-微观模型的方法计算锕系元素的裂变势能曲面。核的表面用五维广义Lawrence 形状描述。采用Strutinsky 的壳修正方法计算形变核的势能。宏观模型能采用了两组不同的公式,即Myers-Swiatecki 公式和LSD 公式。独立粒子势 采用了变形核Woods-Saxon 势和折叠Yukawa 势两种。 设计并实现了两种以分水岭算法为基础的搜索算法。一种是模拟降水算法。将原本用于处理二维图像的分水岭算法移植到了五维曲面上。另一种是模拟泛洪算法。通过引入“岸线”概念,将算法每次迭代所需访问的区域局限在那些状态 有变化的点周围,大大提高了计算速度。 独立编写了势能曲面计算和搜索算法的程序。计算了U 和Pu 元素共23 个同位素核的400 万格点的势能曲面,并从中搜索得到最优裂变路径和裂变势垒。用Brosa 模型(多模式无规颈断裂模型)计算了低能中子诱发U 和Pu 元素裂变的碎片质量分布,动能分布和累计产额。计算裂变碎片的质量分布和动能分布时,将裂变系统的断前形状的几个主要参量与裂变后现象的物理量和分布相联系,从而确定裂变核断前形状。计算累计产额时,使用了大型核反应程序Talys计算了各个反应道的截面。
个人分类: 裂变模型|0 个评论
分享 原子核液滴模型参数数值含义
accumulation 2015-3-26 10:28
1. 液滴模型的 5 个参数大小相差很大,其数值大小的含义是什么? 液滴模型的 5 个参数: (1) 体积能参数 aV=15.835MeV :在不扣除增加一个核子带来的其它效应时,增加一个核子将使体积能为原子核结合能贡献 15.835MeV ;根据实验上 B/A 为常数及原子核的不可压缩性,推导出体积能,体积能与体积成正比,即与 A 成正比,体积能是原子核结合能的主导项; (2) 表面能参数 -aS=-18.33MeV :液滴模型中,原子核存在表面,表面上的核子没有四周核子的包围,因此表面核子的结合能将会弱一些;因此表面能前面的参数为负值,即从体积能中减去表面核子所高估的结合能;表面能的作用与体积能的作用相反,表面能的作用是使得原子核有尽可能小的表面积;表面能参数的数值含义是:计算原子核结合能时,核子数 A 增加导致的每单位 A^ ( 2/3 )的增加,需要扣除 aS=18.33MeV 的表面核子所高估的结合能; (3) 库仑能参数 -aC=-0.714MeV :原子核中存在质子,质子之间的库仑斥力是导致原子核不稳定的因素,库仑能的作用与体技能的作用相反,库仑能降低了原子核的稳定性,即减少了原子核的结合能;库仑能参数的含义为每单位 Z^2/(A^(1/3)) 的增加,将会降低 0.714MeV 的原子核结合能; (4) 对称能参数 -aa=-92.80MeV :原子核中的质子、中子有数目相等的趋势,如果质子、中子数目不相等,在不考虑多余中子抵消库仑斥力的情况下,原子核的稳定性将降低;因此,对称能项反映了质子数与中子数偏离对原子核稳定性的削弱;对称能是一项量子效应,根据泡利不相容原理,在中子、质子对称相处的情形下,能填充的单粒子能级更低;参数 -aa 的含义为,每单位( A/2-Z ) ^2/A 所带来的偏差,将使得原子核结合能降低 92.80MeV ; (5) 对能参数 aP=11.2MeV :在原子核的组成中,有中子、质子各自成对相处的趋势,同类核子成对相处时结合能增大,原子核更稳定;不成对时,结合能就会小些,原子核稳定性降低;对能是一种量子效应,全同费米子的波函数是交换反对称的,当成对的同类核子相距很近,它们的空间波函数对称而自旋波函数反对称时,这样的成对核子有较强的吸引作用,总自旋为零,附加的对能使结合能增大;以奇 A 核的对能为基准,偶偶核会多具有 11.2MeV*A^(-1/2) 的结合能而更稳定,而奇奇核会减少 -11.2MeV*A^(-1/2) 的结合能而更不稳定;
个人分类: 原子核物理|0 个评论
分享 原子核液滴模型对称能形式与费米气体模型
accumulation 2015-3-26 10:27
1. 液滴模型的对称能项采用( A/2-Z)**2/A, 而不是 |A/2-Z| 等其他形式,你猜测原因何在? 对称能采用( A/2-Z)**2/A 形式可由费米气体模型推出: 根据原子核费米气体模型: 中子和质子的最大动量: 费米能量与费米动量之间的非相对论关系: 每个核子的平均动能: 原子核的总平均动能: 取 Z-N= δ,令 Z+N=A 为定值,则 假定δ /A 为小量,则利用二项式展开式 可以将原子核的总动能展开为: 展开式第一项即可理解为体积能,第二项即可理解为对称能; 因此, Z=N 时原子核的能量最小,原子核最稳定; 因为费米气体模型只考虑泡利原理,而没有考虑核子之间的剩余相互作用,因此费米气体模型的理论值只有实验值的一半; 原子核液滴模型的推广公式中,也存在用其它形式表示对称能的情形,如: Myers Swiateeki 公式: Danos Gillet 公式(将对称能与对能结合起来处理): Zeng , Cheng Yang 公式: 修正的液滴模型公式在结合能计算、β稳定线位置、可裂变度随 A 变化、库仑能查变化规律、同位旋多重态质量方程等方面都有改进;
个人分类: 原子核物理|0 个评论
分享 裂变物理模型—目录
accumulation 2015-3-24 01:39
第一章 引言……………………………………………………………………..1 1.1 研究工作的意义……………………………….…...…………………...2 1.2 裂变反应机制 ……………………………………..…...………………...2 1.3 裂变后现象 ………………………………………………….…………...3 1.4 结合能的宏观模型 ……………………………………………………...5 1.5 结合能的宏观-微观修正方法 …………………………………..….…...6 1.6 裂变势能曲面 ……………………………………………………….…...7 1.7 裂变理论模型概述……………………………………….……………...8 1.7.1 统计模型…………………………………………………………...8 1.7.2 输运理论…………………………………………………………...8 1.7.3 多模式模型………………………………………………………...9 1.7.4 无规颈断裂模型………………………………………………….10 1.7.5 多模式无规颈断裂模型 ……………………………..…………...10 1.8 裂变产额的系统学……………………………………………………...11 1.9 TALYS 程序 …………………………………………………………….12 1.10 本工作的研究内容……………………………………………………12 第二章 原子核形变能的计算 …………………………………………….…...15 2.1 核形状描述………………………………………………………..…....15 2.2 宏观模型……………………………………………………………….17 2.2.1 Myers-Swiatecki 液滴模型能公式 …………………………………...17 2.2.2 LSD 公式……………………………………………………….....18 2.3 独立粒子势……………………………………………………….........19 2.3.1 核势的自旋无关部分V1………………………………………...20 2.3.2 自旋轨道耦合项………………………………………………....21 2.3.3 库伦势……………………………….…………………………...21 2.4 计算单粒子能级…………………………………………………….....22 2.5 Hamilton 矩阵元…………………………………………..….………...24 2.6 结合能的宏观-微观模型方法 ……………………………………........25 2.6.1 壳修正………………………………………………………........25 2.6.2 对修正…………………………………………………………....28 2.7 数值算法……………………………………………………….............29 第三章 势能曲面搜索算法研究 ……………………………………………....31 3.1 搜索算法的设计思路……………………………………………….....31 3.2 模拟降水算法………………………………………………………......33 3.3 模拟泛洪算法……………………………………………………….....36 3.4 算法的其它问题讨论………………………………………………….39 第四章 势能曲面计算和搜索结果 …………………………………………....41 4.1 核在裂变路径上的形态演化 ……………………………………….....41 4.2 不同模型的计算结果比较………………………………………….....43 4.3 多模式裂变路径 …………………………………………….………....45 4.4 锕系元素的裂变位垒高度计算结果 ……………………………….....47 第五章 U 同位素裂变产额的Brosa 模型研究………………………….…...51 5.1 Brosa 模型概述…………………………………………………….…...51 5.2 低能中子诱发235U 裂变的碎片质量分布和动能分布 ……………...55 5.3 低能中子诱发238U 裂变的累积产额与能量关系 …………………...62 第六章 总结和展望………………………………………….………………...67 6.1 总结………………………………………………………………..…...67 6.2 结论和创新……………………………………..………..……..……...67 6.3 展望………………………………………………..…………………...68 参考文献…………………………………………………………..….………...71 附录 Hamilton 矩阵元推导过程…………………………………….………...75 A.1 椭球谐振子函数…………………………………………...…..……...75 A.2 动能算符矩阵元………………………………………….…………...75 A.3 和 的矩阵元…………………………………………….…………...77 A.4 自旋-轨道耦合能量算符的矩阵元…………………………………...78
个人分类: 裂变模型|0 个评论
分享 裂变物理模型—摘要
accumulation 2015-3-24 01:31
第一种是 模拟降水算法 。通过 将二维图像处理领域的分水岭算法 的基本思想推广到高维势能曲面中, 设计并实现了用于高维势能曲面搜索的模拟降水算法 。在本工作之前,还未见模拟降水的算法应用于势能曲面搜索的先例。模拟降水算法可以给出完整的 最优裂变路径 , 基态位置、鞍点和断点都在这条路径上 。另一种是改进的 模拟泛洪算法 。Mller 小组在工作中采用了简单的模拟泛洪算法,他们的算法的缺点是非常耗费时间,每次迭代都需要遍历所有的网格点。 本工作首次引入“岸线”概念,使得每次迭代所需访问的区域限定在那些状态有变化的点周围,通过优化算法使程序的执行速度提高了一个数量级。本工作独立编写了 高维势能曲面计算程序和搜索算法的程序 。这两个程序都提供了不同的理论模型和算法供选择,势能曲面计算程序中可以选择不同宏观模型和独立粒子势;搜索算法程序可以指定采用模拟降水算法或模拟泛洪算法。程序的有效性的验证是通过对各种输入和输出进行分析检验来进行的。 计算了U和Pu 元素一系列同位素的裂变位垒 ,并与其他人的工作了比较,提出了进一步改进程序的方向。 用 Brosa 模型(多模式无规颈断裂模型) 计算了 低能中子诱发U 和Pu 元素裂变的碎片质量分布,动能分布和累计产额 。计算裂变碎片的质量分布和动能分布时,将 裂变系统的断前形状的几个主要参量 与 裂变后现象 的物理量和分布相联系,从而确定裂变核断前形状。计算累计产额时, 使用了大型核反应程序Talys计算了各个反应道的截面 。 关键词:裂变位垒,势能曲面,裂变产额,Brosa 模型
个人分类: 裂变模型|0 个评论
分享 裂变物理模型—摘要
accumulation 2015-3-24 01:25
独立编写了 宏观-微观模型 的 高维势能曲面计算程序 。核的表面用 广义Lawrence 形状 描述。采用Strutinsky 的 宏观势能(液滴模型能)加微观壳修正的方法计算形变核的势能 。计算液滴模型能时采用了两组不同的公式,即 Myers-Swiatecki 公式和LSD 公式 。 独立粒子势也采用了变形核Woods-Saxon 势和折叠Yukawa 势两种 。比较不同理论模型的计算结果,确定了LSD 公式加折叠汤川势的模型组合能得出最合理的结果。这是第一次将LSD 公式和折叠汤川势与Lawrence 形状 结合 来计算裂变截面。 计算U 和Pu 同位素的势能曲面时,使用了 多达400 万格点的大型网格 ,保证了计算结果包含了裂变路径上所有的关键细节的信息。为了 在高维曲面上搜索基态、鞍点、断点、裂变路径等关键位置 ,借鉴数学形态学和数字图像处理中的基本理论,设计并实现了 两种以分水岭算法为基础的新型搜索算法 。
个人分类: 裂变模型|0 个评论
GMT+8, 2026-2-13 05:08