《计算物理习题》由会员分享，可在线阅读，更多相关《计算物理习题（6页珍藏版）》请在金锄头文库上搜索。
1、第二章第二章 习题习题（1）采用线性同余法（参见公式(2.2.3)）产生伪随机数。取，和5a 1c 16m .记录下产生出的前 20 数，它产生数列的周期是多少？01x （2）取，和，用线性同余法产生出三维数组137a 187c 256m 01x 和二维数组，然后分别绘出其三维和二维分布图形。12,nnn 1,nn （3）用“投针法”计算出圆周率的数值，画出程序流程框图，并编写程序。（4）已知电子在物质中的作用截面，试写出电子在物质电子对光电总compton层中相互作用的抽样程序框图和程序。（5）编写一个程序按照产生随机数序列，并绘图表明其分布满足分布ln1 i密度函数。 其它, 00, 0,
2、)(xexfx（6）轻子的平均寿命为，试写出 N 个轻子在实验室系中以速度运s13104 . 3v动的飞行距离的抽样程序框图和程序。（7）写出各向同性分布的角度 抽样程序（） 。,dddsin（8）如分布密度函数为， （其中，为整数） ，试写出抽样nyxxenyxf ),(nyx, 0, 1程序框图和程序。 （9）证明 Breit-Wigner 分布22 0)(1)(xxxf可以通过抽样得到 。0()iixxcot （10） 归一化黑体辐射频谱为)()1(15)(44Tkhxdxexdxxfx 其中证明如下抽样步骤得到的抽样分布满足上面的分布，求出它的抽样效率。抽样步骤：让 L 等于满足下面不
3、等式的整数 的最小值，l4 14 1190ljj 然后置，其中为 0，1 区间均匀分布的伪随机数。 )ln(1 5432Lxi（11） 对正则高斯分布抽样:.dxxdxxp 222exp21)(（12） Gamma 函数的一般形式为0,)!1()(1xdxexnadxxfaxnn抽样证明其抽样方法可以为.).ln(1 121nna（13）分布的一般形式为20,)2/(21)(2/12/ 2/xdxexndxxfxn n抽样证明其抽样方法可以为, 其中为标准正态分布的 n 个独立抽样值. niix12nxxx,.,21（14） 选择偏倚分布密度函数，用蒙特卡洛重要抽样法求积分xexg)(.02/
4、3dxexx（15） 编写一个程序，采用 Metropolis 随机游走的方法产生按高斯分布， 222/exp)(xAxf21的随机点。抽样中常数的值需要知道吗？试决定接受点与试探步数之比，到达A平衡分布的时间与最大试探步长的关系。 （提示：判断“平衡”的标准是） 。选多大较合理? 22 x（16） 用 Metropolis 随机游走的方法计算积分。40,402xdxexx（17） Laplace 方程及其边界条件为， 1, 1, 0, 01 ,0 ,0,2yyxxyx用随机游走的蒙特卡洛方法数值求解正方形场域的势函数。11, 10yx第三章第三章 习题习题（1）利用蒙特卡洛方法计算三维、四维
5、、五维和六维空间的单位半径球的体积。（2）利用分布密度函数做重要抽样来求积分，并分析误差与投点数的关系。xAexf)( 02/5dxexIx（3）用事例证明蒙特卡洛求积分的标准误差为，NAA1222其中 A 为物理观测量，N 为蒙特卡洛投点个数。 （4）采用 Metropolis 方法产生一维分子速度分布密度函数为22)(veCvvf的游走样本点，并将其分布和上述分布函数曲线进行比较（上式中为常数） 。,C（5）写出采用 Metropolis 方法对高斯分布的抽样框图和程序（)2/exp()(2xAxf和为常数） 。A （6）编写采用 Metropolis 算法计算一维积分dxxx222/ex
6、p的程序。用该程序计算三维积分.dzrzdyydxx2/exp2222然后将得到的结果与解析计算得到的精确结果进行比较，并分析模拟游走点数与误差 的关系。 （7）对以下一维扩散方程， 1022), 1 (), 0()()0 ,(,UtUUtUxfxUtU xU可以通过在的矩形格点上的随机游走来求解（其中，分别为和 划分的格hhxt点长度） 。在方向向前和向后游走的几率为，而在 方x12 0301/2ht向向前和向后游走的几率为。试编程予以计算。12 0402/21h（8）编写程序，采用路径积分量子力学蒙特卡洛方法求液态的基态能量。4He（9）编写程序，采用变分量子力学蒙特卡洛方法求氢分子的基态
7、能量。其中两质子和两电 子应当按四体系统来处理。 （10）修改习题（8）的程序，采用格林函数量子蒙特卡洛方法求氢分子基态能量，并与习 题（8）的结果进行比较。第四章第四章 习题习题(1) 用有限差分法发展一个程序，数值求解正方形场域的拉普拉斯01,01xy方程，。 1, 1, 0, 01 ,0 ,0,2yyxxyx（2）用有限差分法发展一个程序，数值求解极坐标下的泊松方程。 212,4,VbVarr然后，选择，边界条件时，两个圆0,r1,0,2, 121VVba圈中间的势分布。 （3）第（2）题中若， 其它，0,/exp101 ,braaarr其余取值相同，数值求解两个圆圈中间的势分布。 （4
8、）在一个二维的反应堆中，中子的扩散方程为LL，0sinsin,2 2222 Ly Lxyxayx用的正方形网格离散化后，证明它的有限差分方程满足hhhyx。1 2 22 2sin8sinsin Lh haLhj Lhiij第五章第五章 习题习题（1）公式(5.2.20)是更一般的公式en ml jk inlknlkdxdyNNN!2!的特殊情况。试给出证明。（2）用有限元素法发展一个程序，数值求解正方形场域的拉普01,01xy拉斯方程， 1, 1, 0, 01 ,0 ,0,2yyxxyx （3）修改（2）题中的程序，仍然采用有限元素法数值求解三角形场域的拉普拉斯方程，xyxD1, 10 11,
9、 0, 00 ,0,2xyxyxyx第六章第六章 习题习题（1）写用 Verlet 速度算法求解三维分子运动方程的程序。（2）编写一个三维，元胞尺寸为的周期边界条件计算程序。3L（3）试做总能量固定的单原子系统的分子动力学模拟。元胞为，划10zyxLLL分为的正方形网格。元胞内原子数。原子质量。位势为10101064N1m Lenard-Jones 势，其中，边界条件为周期性边界条件，初始位置是随机1 分布在正则节点上，初始速度为按-1，1随机分布。分子动力学模拟步长取为，模拟 100-200 步后原子的速度分布和位置分布如何？02. 0t （4）试做二维单原子系统的分子动力学模拟。系统温度保
10、持固定,模拟参数及85. 0T 其他条件同上题。第七章第七章 习题习题（1）用 Mathematica 语言定义求解一元二次方程的函数，该函数还要02cbxax求能处理各种常数 a,b,c 的情况。 （2）用 Mathematica 语言定义一个能画出任意给定 n 值的正 n 边形的函数 （3）用 Mathematica 语言定义一个操作函数，它可以在给定二维矩形 x-y 平面区间，按给定步长划分矩形网格，并列出节点的坐标表。yxhh ,iiyx ,（4）接着上题，用 Mathematica 语言定义一个操作函数，对确定步数 n，生成一个在 x- y 平面格点上 n 步服从均匀分布的随机游走的
11、图。 （5）用 Mathematica 语言实现一个产生任意阶的勒让德多项式的 Mathematica 程序包。勒 让德多项式的递推公式为：， ， .1)(0xPxxP)(1 ) 1(12)(11nxnPxxPnxPnnn（6）用 Mathematica 语言实现一个产生任意阶的埃米尔特多项式的 Mathematica 程序包。 埃米尔特多项式的递推公式为： ， ， 1)(0xHxxH)(1 xnHxxHxHnnn1122)(（7）Mathematica 语言编写一个从某点出发求多元函数的局部极小或极大值的程序包。 （8）用 Mathematica 语言编写一个程序包，它能实现平面图形的（a）
12、平移， （b）旋转， (c)对 x 座标轴的反射。第八章第八章 习题习题（1）画出不同主量子数、轨道量子数下，氢原子径向部分波函数随的变化图形，并ln,r讨论原子序数变化的作用。Z（2）画出不同轨道量子数、磁量子数、下，氢原子部分波函数随和的变化ml,的三维图形。 （3）采用诺曼诺夫(Numerov)法（参见附录 C） ，编写一个程序求解电子的一维薛定格方程 的最低的两个能量本征值和波函数。该电子所在势阱的势函数为(我们选择原子单位a.u.，即)1cemeh。 其它,5050,10)(xxxV（4）编写求解不同势阱高度和宽度的一维势阱薛定格方程波函数和能量谱的程序包。 势阱的势函数为。 其它0)(0axVxV（5）试用 Schroedinger.m 程序计算时薛定格方程的基态，第一激发态的能量值，2)(rrV并与变分法求得的结果进行比较。
