蒙特卡罗积分法在求解数值积分中的应用

对于一个实际结构,由有限元法离散化处理后,动力学方程可写为:

对以上方程式(10),考虑某一模态的振动,并略去下标i可写为:



将上式对时间求导,并将t,t+2△(其中2△为时间步长)带入t0,最后通过抛物线法则计算式中的卷积则可得到




应用上述递推公式,以前一时刻来求后一时刻的结果。计算不重复。当aij求出后,以后在时域中的步进求解只是一些简单的数组相乘。计算速度很快。

%可计算无阻尼和阻尼强迫震动

%输入可为数组或函数荷载

disp('扩大四倍后位移随时间的变化');

title('加载时间增加四倍时的位移位移时程曲线');

disp('荷载作用时间内位移变化');

disp('荷载作用时间内位移最大值');

disp('荷载作用时间内速度变化');

disp('荷载作用时间内速度最大值');

disp('荷载作用时间内加速度变化');

disp('荷载作用时间内加速度最大值');

title('荷载作用时间内位移时程曲线');

title('荷载作用时间内速度时程曲线');

title('荷载作用时间内加速度时程曲线');

disp('扩大四倍后位移随时间的变化');

disp('扩大四倍后位移最大值');

disp('扩大四倍后速度随时间的变化');

disp('扩大四倍后速度最大值');

disp('扩大四倍后加速度随时间的变化');

disp('扩大四倍后加速度最大值');

title('加载时间增加四倍时的位移时程曲线');

title('加载时间增加四倍时的速度时程曲线');

title('加载时间增加四倍时的加速度时程曲线');

}

Bessel函数源于微分方程,含Bessel函数的积分一直是物理、化学等诸多工程应用领域中的一个热点问题。在大多数情况下,这类积分无法得到解析解,需进行数值求积;而在振荡因子很大时Bessel函数呈高振荡性,传统的数值积分方法往往效率低,甚至失去效用。本文应用近年来快速发展的高振荡函数数值积分方法理论,结合Luke、Watson等人对Bessel函数的研究成果,充分利用Maple以及MATLAB的计算功能,给出两类Bessel变换的Filon-type积分方法,并将此算法以可视化界面形式予以实现,可供需要解决含Bessel函数数值积分的科学工作者使用参考。 全文由四章组成。 本文第一章综述了Bessel函数的来历背景,以及高振荡函数数值积分特别是含有Bessel函数的数值积分的经典方法和最新发展,并指出它们的成功与不足。 第二章给出了两类含驻点Bessel函数积分的基于特殊函数的Filon-type方法及误差估计;重点研究了Moments的求解。 第三章将本文提出来的含有Bessel函数的数值积分方法以MATLAB7.0.1为平台予以实现,为用户提供一个计算这类积分的简洁的、可视化的积分工具。 第四章,运用论文提出的方法计算几个积分,作出误差图。得到的误差以及与Maple计算时间的对比展示了本文方法的高效性。

【学位授予单位】:中南大学
【学位授予年份】:2007

支持CAJ、PDF文件格式


纪刚,张纬康,卢晓平;[J];海军工程大学学报;2004年06期
翁爱华,王雪秋;[J];西北地震学报;2003年03期
中国硕士学位论文全文数据库
孙胜先;褚玉明;;[J];安徽工学院学报;1990年04期
刘芳;[J];安徽师范大学学报(自然科学版);1999年03期
刘芳;[J];安徽教育学院学报(自然科学版);1999年02期
郑明;吴晓龙;赵宝生;;[J];辽宁科技大学学报;2010年03期
吕金城;原新生;;[J];安阳师范学院学报;2006年02期
王伟,郭尚平,简水生;[J];北方交通大学学报;1996年01期
刘式适;[J];北京大学学报(自然科学版);1990年03期
付遵涛,刘式适,付彩霞;[J];北京大学学报(自然科学版);1998年05期
杨争光,苏东林,李梅,吕善伟;[J];北京航空航天大学学报;2004年07期
中国博士学位论文全文数据库
徐凌洋;[D];中国农业科学院;2011年
程瑶;[D];辽宁工程技术大学;2011年
刘希强;[D];中国工程物理研究院北京研究生部;2002年
中国硕士学位论文全文数据库
史秀宝;[D];辽宁工程技术大学;2009年
王建勇;[D];辽宁工程技术大学;2009年
赵勇;[D];西安电子科技大学;2011年
王蓉;[D];西安电子科技大学;2011年
郭汉伟,何建国,尹家贤,刘培国,刘克成;[J];国防科技大学学报;2000年04期
李吉松;朴化荣;;[J];物探化探计算技术;1993年02期
邵先喜,尹传存;[J];数学研究与评论;2001年03期
王秀友;[J];淮北煤炭师范学院学报(自然科学版);2005年02期
冯兆生;[J];安徽建筑工业学院学报;1994年02期
杨文熊,顾尔祚,朱敏;[J];应用数学和力学;1990年08期
中国重要会议论文全文数据库
张前安;吴逢铁;郑维涛;;[A];中国光学学会2011年学术大会摘要集[C];2011年
卞磊;唐少强;;[A];北京力学会第18届学术年会论文集[C];2012年
中国博士学位论文全文数据库
张军勇;[D];中国工程物理研究院;2011年
中国硕士学位论文全文数据库
王海丽;[D];南京航空航天大学;2012年
}

我要回帖

更多关于 蒙特卡罗积分 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信