土木工程数值分析及应用案例2015.11
《土木工程数值分析及应用》案例讨论
一、验证圣维南原理
用数值方法验证圣维南原理
二、简支梁3种类型单元分析及单元性能比较
梁单元、3节点、4节点平面单元。
三、悬臂柱具有n个质量，比较基底剪力法和振型分解反应谱法的结果
编程求悬臂柱的周期和振型，用抗震规范的方法计算地震作用。
四、大型稀疏矩阵的求解算法比较
稀疏矩阵的求解算法与常规的LU分解法比较，分别从不同的稀疏率、矩阵大小等方面比较
五、框架的位移分析
编程计算小型框架的位移
六、基于Mindlin公式的群桩基础柔度矩阵及刚度矩阵
先用Mindlin公式计算各桩沉降，得到柔度矩阵，然后用数值方法将柔度矩阵转换成刚度矩阵。
七、带柱帽无梁楼盖近似计算方法的精确度分析
近似计算方法：等代框架法、经验系数法。将近似方法与有限元分析结果进行对比。
八、用高斯数值积分方法求柱下独立基础的沉降
规范采用分层总和法计算沉降，其中平均附加应力系数查表得到。根据e-p 曲线沿土层厚度方向积分计算沉降，并与规范方法进行对比。
九、消防车等效均布荷载的函数曲线拟合
等效均布荷载与柱网大小有关，根据《荷载规范》的方法分别计算不同柱网尺寸下的等效均布荷载，以柱网尺寸为变量，拟合出等效均布荷载的函数表达式。
十、边坡稳定性分析的数值模拟
由学生提出
十一、Gauss和LU分解法
由学生提出
十二、多质点隔震结构分析
由学生提出
十三、基坑围护结构的数值模拟由学生提出
十四、水流对桥墩影响的数值模拟由学生提出
数值分析试题及答案汇总
数值分析试题 一、 填空题（2 0×2′） 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值，则x 有 2 位 有效数字。 2. 若f (x )=x 7－x 3＋1，则f [20,21,22,23,24,25,26,27]= 1 ， f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设，‖A ‖∞＝___5 ____，‖X ‖∞＝__ 3_____， ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =?(x )在有解区间满足 |?’(x )| <1 ，则使用该迭代 函数的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差商 公式的 前插公式 ，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的 后插公式 ；如果要估计结果的舍入误差，应该选用插值公式中的 拉格朗日插值公式 。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是：=∑=n i i x a 0)( 1 ；所以当 系数a i (x )满足 a i (x )>1 ，计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%，至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ，线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收 敛于方程组的精确解x *的充分必要条件是 ?(B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。 11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。 12. 线性方程组的松弛迭代法是通过逐渐减少残差r i (i =0,1,…,n )来实现的，其中的残差 r i ＝ (b i -a i1x 1-a i2x 2-…-a in x n )/a ii ，(i =0,1,…,n )。 13. 在非线性方程f (x )=0使用各种切线法迭代求解时，若在迭代区间存在唯一解，且f (x )
数值分析学期期末考试试题与答案(A)
期末考试试卷（A 卷） 2007学年第二学期 考试科目： 数值分析 考试时间：120 分钟 学号 姓名 年级专业 一、判断题（每小题2分，共10分） 1. 用计算机求 1000 1000 1 1 n n =∑时，应按照n 从小到大的顺序相加。 （ ） 2. 为了减少误差,进行计算。 （ ） 3. 用数值微分公式中求导数值时，步长越小计算就越精确。 （ ） 4. 采用龙格－库塔法求解常微分方程的初值问题时，公式阶数越高，数值解越精确。（ ） 5. 用迭代法解线性方程组时，迭代能否收敛与初始向量的选择、系数矩阵及其演变方式有 关，与常数项无关。 （ ） 二、填空题（每空2分，共36分） 1. 已知数a 的有效数为0.01，则它的绝对误差限为________，相对误差限为_________. 2. 设1010021,5,1301A x -????????=-=-????????-???? 则1A =_____，2x =______，Ax ∞ =_____. 3. 已知5 3 ()245,f x x x x =+-则[1,1,0]f -= ,[3,2,1,1,2,3]f ---= . 4. 为使求积公式 1 1231 ()()(0)33 f x dx A f A f A f -≈- ++? 的代数精度尽量高，应使1A = ，2A = ，3A = ，此时公式具有 次的代数精度。 5. n 阶方阵A 的谱半径()A ρ与它的任意一种范数A 的关系是 . 6. 用迭代法解线性方程组AX B =时，使迭代公式(1) ()(0,1,2,)k k X MX N k +=+=产 生的向量序列{ }() k X 收敛的充分必要条件是 . 7. 使用消元法解线性方程组AX B =时，系数矩阵A 可以分解为下三角矩阵L 和上三角矩
数值分析最佳习题(含答案)
第一章 绪论 姓名 学号 班级 习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。 1 若误差限为5 105.0-?,那么近似数有几位有效数字？（有效数字的计算） 解：2*103400.0-?=x ，325* 102 1 1021---?=?≤-x x 故具有3位有效数字。 2 14159.3=π具有4位有效数字的近似值是多少？（有效数字的计算） 解：10314159.0?= π，欲使其近似值* π具有4位有效数字，必需 41*1021 -?≤-ππ，3*3102 11021--?+≤≤?-πππ，即14209.314109.3*≤≤π 3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +，b a ?有几位有效数字？（有效数字的计算） 解：3* 1021-?≤ -a a ，2*102 1 -?≤-b b ，而1811.2=+b a ，1766.1=?b a 2123****102 1 10211021)()(---?≤?+?≤-+-≤+-+b b a a b a b a 故b a +至少具有2位有效数字。 2 123*****102 1 0065.01022031.1102978.0)()(---?≤=?+?≤-+-≤-b b a a a b b a ab 故b a ?至少具有2位有效数字。 4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差？（误差的计算） 解：已知 δ=-* *x x x ，则误差为 δ=-= -* **ln ln x x x x x 则相对误差为 * * ** * * ln ln 1ln ln ln x x x x x x x x δ = -= - 5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5* =，已知 cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v 2π=的绝对误差限与相对误差 限。（误差限的计算） 解： * 2******2),(),(h h r r r h r r h v r h v -+-≤-ππ 绝对误差限为 π ππ252.051.02052)5,20(),(2=??+????≤-v r h v
数值分析试卷及其答案
1、（本题5分）试确定7 22 作为π的近似值具有几位有效数字，并确定其相对误差限。 解 因为 7 22 =3.142857…=1103142857 .0-? π=3.141592… 所以 312102 11021005.0001264.0722--?=?=<=- π （2分） 这里，3,21,0=-=+-=n n m m 由有效数字的定义可知7 22 作为π的近似值具有3位有效数字。 （1分） 而相对误差限 3102 1 0005.00004138.0001264.07 22-?= <≈= -= π π πε r （2分） 2、（本题6分）用改进平方根法解方程组：??? ?? ??=????? ??????? ??--654131*********x x x ； 解 设???? ? ??????? ? ?????? ??===????? ??--11111 1 131321112323121 32 132 31 21 l l l d d d l l l LDL A T 由矩阵乘法得： 5 7,21,215 27 ,25,2323121321- ==-== -==l l l d d d （3分） 由y D x L b Ly T 1 ,-==解得 T T x y )9 23 ,97,910(,)563, 7,4(== （3分） 3、（本题6分）给定线性方程组???????=++-=+-+=-+-=-+17 7222382311387 510432143213 21431x x x x x x x x x x x x x x 1）写出Jacoib 迭代格式和Gauss-Seidel 迭代格式； 2）考查Jacoib 迭代格式和Gauss-Seidel 迭代格式的敛散性； 解 1）Jacoib 迭代格式为
数值分析期末考试复习题及其答案.doc
数值分析期末考试复习题及其答案 1. 已知325413.0,325413* 2* 1==X X 都有6位有效数字，求绝对误差限。（4分） 解： 由已知可知,n=6 5.01021 ,0,6,10325413.0016*1=?= =-=?=ε绝对误差限n k k X 2分 620* 21021,6,0,10325413.0-?=-=-=?=ε绝对误差限n k k X 2分 2. 已知?????=001A 220 - ???? ?440求21,,A A A ∞ （6分） 解： {},88,4,1max 1==A 1分 {},66,6,1max ==∞A 1分 () A A A T max 2λ= 1分 ?????=001A A T 420 ?? ?? ? -420?????001 220 - ?????440=?????001 080 ???? ?3200 2分 {}3232,8,1max )(max ==A A T λ 1分 24322==A 3. 设3 2 )()(a x x f -= （6分） ① 写出f(x)=0解的Newton 迭代格式 ② 当a 为何值时，)(1k k x x ?=+ （k=0,1……）产生的序列{}k x 收敛于2 解： ①Newton 迭代格式为： x a x x x a x a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(2 2 32 1 += +=---=-=+? 3分
②时迭代收敛即当222,112 10)2(',665)('2<<-<-=-=a a x a x ?? 3分 4. 给定线性方程组Ax=b ，其中：? ??=1 3A ??? 22，??????-=13b 用迭代公式)()()()1(k k k Ax b x x -+=+α（k=0,1……）求解Ax=b ，问取什么实数α，可使迭代收 敛 （8分） 解： 所给迭代公式的迭代矩阵为?? ? --? ??--=-=ααααα21231A I B 2分 其特征方程为 0) 21(2)31(=----= -αλα ααλλB I 2分 即，解得αλαλ41,121-=-= 2分 要使其满足题意，须使1)(岩土工程计算原理和方法
岩土工程数值计算原理与方法 随着计算机的计算速度和存储能力的飞速发展以及计算方法的日益完善，数值模拟方法已经成为研究未知领域的强有力的工具。在岩土工程计算与分析中数值计算原理与方法也发展很快。特别是有限元的发展，促进了岩土工程研究、工程预测、优化设计和计算机辅助设计等的发展。但在工程实际中使用数值计算原理与方法却存在一些问题：例如有些人因缺乏对有限元和工程性质的深入了解，而有限元的迅速发展给他们造成一种假象，认为它是万能的，可以处理几乎所有的岩土工程问题；同时他们又被有限元计算结果的精度所迷惑，不了解这些精确结果后面所隐藏的不确定性，也不了解这些数值方法所采用本构模型的局限性以及相应参数的不确定性；因这些不确定性导致数值计算原理与方法的预测结果与实际情况和实际经验相差很大，又由于部分人计算偏于保守，使得岩土工程师难以接受现代数值计算原理与方法。 1. 岩土工程数值计算原理与方法也具有两面性。 有些人偏向于用其进行岩土工程的分析计算的原因在于： (1)数值计算原理与方法能够做任何传统的分析方法所能做到的分析与计算，而且做得更多、更好。 (2)数值计算原理与方法能够给出复杂数学模型的解。因而能够从机理上预测工程性质，而不是统计和经验性的描述，这是一大优点；而简化或经验分析方法有时只能描述其表面或形式上(统计)的关系，缺乏物理机制的描述和探讨。 (3)该方法既能处理简单问题，也能处理复杂问题。 数值计算原理与方法难以被其他人接受的原因在于： (1)使用复杂，难以被很好的掌握。 (2)数值计算原理与方法本身的不确定性(指与精确的解析方法相比所产生的不确定性，特别是在岩土动力非线性问题中这种不确定性会很大)导致预测结果与工程实际不符。 (3)数值计算原理与方法所使用的物理模型或本构模型有局限性，难以反映实际情况，导致预测结果与工程实际不符。 (4)采用复杂模型要求较多的参数，而这些参数难以用简单试验获得。 (5)既然数值计算原理与方法和传统的分析方法都具有很大的不确定性，还不如采用传统的分析方法，因为传统的方法简单、实用。 (6)精确的数值分析结果会误导使用者迷信这些结果的精确性，而没有认识到其后面隐
数值分析整理版试题及答案
数值分析整理版试题及答案
例1、 已知函数表 x -1 1 2 ()f x -3 0 4 求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解： （1）k x -1 1 2 k y -3 0 4 插值基函数分别为 ()()()()()()()()()() 1200102121()1211126 x x x x x x l x x x x x x x ----= ==-------- ()()()()()()()() ()()021******* ()1211122x x x x x x l x x x x x x x --+-= ==-+---+- ()()()()()()()()()()0122021111 ()1121213 x x x x x x l x x x x x x x --+-= ==-+--+- 故所求二次拉格朗日插值多项式为 () ()()()()()()()()()()2 20 2()11131201241162314 121123537623k k k L x y l x x x x x x x x x x x x x ==?? =-? --+?-+-+?+-????=---++-=+-∑ （2）一阶均差、二阶均差分别为
[]()()[]()()[][][]010********* 011201202303 ,11204 ,412 3 4,,5 2,,126 f x f x f x x x x f x f x f x x x x f x x f x x f x x x x x ---===-----= = =----=== --- k x ()k f x 一阶 二阶 -1 -3 1 0 3/ 2 2 4 4 5/6 故所求Newton 二次插值多项式为 ()()[]()[]()() ()()()20010012012,,,35 311126537623P x f x f x x x x f x x x x x x x x x x x x =+-+--=-+ +++-=+- 例2、 设2 ()32f x x x =++，[0,1]x ∈，试求()f x 在[0, 1]上关于()1x ρ=，{} span 1,x Φ=的最佳平方逼近多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，且()1x ρ=，这样，有
数值分析判断题及答案
} 判断题及答案 (认为正确的在题后的括号中打√，认为错误的打×) 1. 按四舍五入的原则，的具有5位有效数字的近似数是. （对） 2. 在做数值计算时，为减少误差，应该尽可能的避免大数做分母。 （错） 3. 计算机上将1000个数量级不同的数相加，不管次序如何结果都是一样的。 （错） 4. 高精度的运算可以改善问题的病态性. （错） 5. 在插值条件相同的情况下，使用Lagrange 插值法和Newton 插值法，所得到 的插值多项式相同。 （对） ！ 6. 假设()(0,,)i l x i n =是Lagrange 插值基函数，则11()0,()1n n l x l x ==。 （对） 7. 高次插值多项式不能令人满意的主要原因是不会出现龙格现象。 （错） 8. Newton 插值方法的一个优点是在增加新的插值节点后，原来计算结果还可以 使用。 （对） 9.曲线拟合和插值是一回事。 （错） 10.二次拟合曲线过给定的所有数据点。 （错） 11.矛盾方程组的法方程组的解就是该矛盾方程组的精确解。 （错） 多项式()n P x 当n 是偶数时是偶函数，当n 是奇数时是奇函数。（对） 13.切比晓夫多项式所满足的递推关系是11()2()(),(1,2,)n n n T x xT x T x n +-=-=。 （对） 《 14.假设()n T x 是[-1,1]上首项系数为1的切比晓夫多项式，()n n Q x H ∈是任一个 首项系数为1的多项式，则1111 max |()|max |()|n n x x T x Q x -≤≤-≤≤≤。 （对） 15.梯形公式和两点高斯公式的代数精度是一样的。 （错） 16.假设n x R ∈，则1||||||||x x ∞<。 （错） 17. 假设n n x R ?∈，则1||||||||x x ∞<。 （错） 18. 假设n x R ∈，则1||||||||x n x ∞<。 （错） 19.只要矩阵n n A R ?∈非奇异，则求解线性方程组Ax b =的直接顺序消去法或直接 LU 分解法可以得到方程组的解。 （错） 20.对称正定的方程组总是良态的。 （错） 21.奇异矩阵的范数一定是零。 （错） &
数值分析试卷及其答案1
1. 已知325413.0,325413*2*1==X X 都有6位有效数字，求绝对误差限。（4分） 解： 由已知可知6 5.0102 1 ,0,6,10325413.0016*1=?= =-=?=ε绝对误差限n k k X 2分 620*2102 1 ,6,0,10325413.0-?= -=-=?=ε绝对误差限n k k X 2分 2. 已知?? ???=0 01 A 220- ?????440求21,,A A A ∞ （6分） 解： {}, 88,4,1max 1==A 1分 {}, 66,6,1max ==∞A 1分 () A A A T max 2λ= 1分 ?????=0 1 A A T 4 2 ???? ? -420?????0 01 2 20 - ???? ?440= ?????0 01 80 ???? ?3200 2分 {}32 32,8,1max )(max ==A A T λ
1分 24322==A 3. 设32)()(a x x f -= （6分） ① 写出f(x)=0解的迭代格式 ② 当a 为何值时，)(1k k x x ?=+ （0,1……）产生的序列{}k x 收敛于 2 解： ①迭代格式为： x a x x x a x a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(2 2 32 1 += +=---=-=+? 3 分 ②时迭代收敛即当222,112 10)2(',665)('2<<-<-=-= a a x a x ?? 3分 4. 给定线性方程组，其中：?? ?=13A ?? ?2 2，?? ? ???-=13b 用迭代公式 )()()()1(k k k Ax b x x -+=+α（0,1……）求解，问取什么实数α ，可使 迭代收敛 （8分） 解： 所给迭代公式的迭代矩阵为?? ? --???--=-=ααααα21231A I B 2分
数值分析整理版试题及答案
例1、 已知函数表 求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解： （1） 故所求二次拉格朗日插值多项式为 （2）一阶均差、二阶均差分别为 例2、 设2 ()32f x x x =++，[0,1]x ∈，试求()f x 在[0,1]上关于()1x ρ=，{}span 1,x Φ=的最佳平 方逼近多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，且()1x ρ=，这样，有 所以，法方程为
011231261192 34a a ??????????=?????????? ?????????? ，经过消元得012311 62110123a a ??? ???????=???????????????????? 再回代解该方程，得到14a =，011 6 a = 故，所求最佳平方逼近多项式为* 111()46S x x =+ 例3、 设()x f x e =，[0,1]x ∈，试求()f x 在[0,1]上关于()1x ρ=，{}span 1,x Φ=的最佳平方逼近 多项式。 解： 若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，这样，有 所以，法方程为 解法方程，得到00.8732a =，1 1.6902a =， 故，所求最佳平方逼近多项式为 例4、 用4n = 的复合梯形和复合辛普森公式计算积分1 ? 。 解： （1）用4n =的复合梯形公式 由于2h =，( )f x =()121,2,3k x k k =+=，所以，有 （2）用4n =的复合辛普森公式 由于2h =，( )f x =()121,2,3k x k k =+=，()12 220,1,2,3k x k k + =+=，所以，有 例5、 用列主元消去法求解下列线性方程组的解。 解：先消元 再回代，得到33x =，22x =，11x = 所以，线性方程组的解为11x =，22x =，33x = 例6、 用直接三角分解法求下列线性方程组的解。 解： 设 则由A LU =的对应元素相等，有 1114u = ，1215u =，1316u =， 2111211433l u l =?=，3111311 22 l u l =?=， 2112222211460l u u u +=?=-，2113232311 545l u u u +=?=-，
北航2010-2011年研究生数值分析期末模拟试卷1-3
数值分析模拟试卷1 一、填空（共30分，每空3分） 1 设??? ? ??-=1511A ，则A 的谱半径=)(a ρ______，A 的条件数=________. 2 设 ,2,1,0,,53)(2==+=k kh x x x f k ，则],,[21++n n n x x x f ＝________, ],,[321+++n n n n x x x x f ，＝________. 3 设?????≤≤-++≤≤+=2 1,121 0,)(2 323x cx bx x x x x x S ，是以0，1，2为节点的三次样条函数，则b=________,c=________. 4 设∞=0)]([k k x q 是区间［0，1］上权函数为x x =)(ρ的最高项系数为1的正交多项式族，其中1)(0=x q ，则 ?=1 )(dx x xq k ________，=)(2 x q ________. 5 设???? ??????=11001a a a a A ，当∈a ________时，必有分解式，其中L 为下三角阵，当 其对角线元素)3,2,1(=i L ii 满足条件________时，这种分解是唯一的. 二、（14分）设4 9,1,41,)(2102 3 === =x x x x x f , （1）试求)(x f 在]4 9,41[上的三次Hermite 插值多项式)(x H 使满足 2,1,0),()(==i x f x H i i ，)()(11x f x H '='. （2）写出余项)()()(x H x f x R -=的表达式. 三、（14分）设有解方程0cos 2312=+-x x 的迭代公式为n n x x cos 3 2 41+ =+， （1） 证明R x ∈?0均有? ∞ →=x x n x lim （? x 为方程的根）； （2） 取40=x ，用此迭代法求方程根的近似值，误差不超过，列出各次迭代值； （3）此迭代的收敛阶是多少？证明你的结论. 四、(16分) 试确定常数A ，B ，C 和，使得数值积分公式 有尽可能高的代数精度. 试问所得的数值积分公式代数精度是多少？它是否为Gauss 型的？
岩土工程数值分析学习笔记(DOC)
岩土工程数值分析读书笔记 摘要：阅读笔记分为两部分：理论学习和plaxis模拟相关问题。 理论部分 0岩土工程数值分析简介 岩土工程问题解析分析是以弹塑性力学理论和结构力学作为理论依据，适用于解决连续介质、各向同性材料、未知量少、边界条件简单的工程问题，存在很大的局限性。 岩土工程问题数值分析是借助于计算机的计算能力，适用于解决材料复杂、边界条件复杂、任意荷载、任意几何形状，适用范围广。 岩土工程数值分析发展过程： 20世纪40年代，使用差分法解决了土工中的渗流及固结问题，如土坝渗流及浸润线的求法、土坝及地基的固结等。 20世纪60年代，使用有限元法成解决了土石坝的静力问题的求解。 20世纪70年代，使用有限元法解决了土石坝及高楼(包括地基)的抗震分析。 20世纪80年代，边界元法异军突起，解决了半无限域的边界问题；地基的静力及动力问题都使用边界元法得到了有效地解决。 岩土工程数值分析的方法有两类，一类方法是将土视为连续介质，随后又将其离散化，如有限单元法、有限差分法、边界单元法、有限元线法、无单元法以及各种方法的耦合。另一类计算方法是考虑岩土材料本身的不连续性，如裂缝及不同材料间界面的界面模型和界面单元的使用，离散元法，不连续变形分析，流形元法，颗粒流等数值计算方法。 1数值分析过程中存在的问题及解决措施 问题：(1)对岩土工程数值分析方法缺乏系统的知识和深入的理解，出现问题时不知道在什么情况下属于理论问题或数学模型问题；在什么情况下是属于计算方法问题或本构模型问题；在什么情况下是参数的确定问题或计算本身的问题等。 (2)各种本构模型固有的局限性。具有多相性土的物理力学性质太复杂，难以准确地用数学模型和本构模型描述。例如邓肯一张模型不能反映剪胀性，不能反映压缩与剪切的交叉影响； (3)现有的试验手段和设备不能提供适当、合理和精确的参数。靠少数样本点所获得的参数难以准确地描述整个空间场地的物理力学性能；土的参数因土样扰动难以高质量的获取，其精度很差。 (4)数学模型还会给人造成一种错觉，让人觉得其计算结果也一定会更好、更可靠。这样可能使人们忽略了精确的数学公式也照样会有出错的可能性。只有当输入参数的质量和精度很高，并能与数学模型的精度相匹配时，才有可能得到较为准确的计算结果。 措施：(1)加强对土的本构模型的教学与培训，了解和掌握各种土的本构模型的优点和局限性以及模型参数的离散性。 (2)在使用数值分析方法的同时，不断地积累使用经验，包括他人的经验。
数值分析期末试卷
数值分析2006 — 2007学年第学期考试 课程名称：计算方法 A 卷 考试方式：开卷[] 闭卷[V ] 半开卷[] IV 类 充要条件是a 满足 二、（18分）已知函数表如下 1?设 f(0) = 0, f (1) =16 , f( 2) =46，则 f [0,1]= ,f[0,1,2]二 2 ?设 AJ <2 -3 -1 ，则X ，A := A 1 1 j — 3 ?计算积分 xdx ，取4位有效数字。用梯形公式求得的近似值为 "0.5 （辛普森）公式求得的近似值为 ，用 Spsn 4?设f （x ）二xe x -3，求方程f （x ） =0近似根的牛顿迭代公式是 ，它的收 敛阶是 5 ?要使求积公式 1 1 [f (x)dx 拓一(0) + A , f (x 1)具有2次代数精度，则 捲= _________________ , 0 4 6 ?求解线性方程组 x 1 ax 2 = 4 , 12_3 （其中a 为实数）的高斯一赛德尔迭代格式收敛的 10 11 12 13 In x 2.3026 2.3979 2.4849 2.5649
三、(20分)构造如下插值型求积公式，确定其中的待定系数，使其代数精度尽可能高, 并指出所得公式的代数精度。 2 f (x)dx ： A o f (0) A f (1) A2f(2) o
X 2 4 6 8 y 2 11 28 40 五、（14分）为求方程X ’ -X 2 -1 =0在X o =1.5附近的一个根，将方程改写为下列等价 形式，并建立相应的迭代公式: 试问上述两种迭代公式在 x 0 =1.5附近都收敛吗？为什么？说明理由。 （1）X =1 ?丄，迭代公式 X 1 X k 1 = 1 - X k （2） X 2二1 ，迭代公式 X —1 2 （X k ）； X k 1
《数值计算方法》试题集及答案
《计算方法》期中复习试题 一、填空题： 1、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案：2.367，0.25 2、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2 x 的系数为 ，拉 格朗日插值多项式为 。 答案：-1， )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 3、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字； 4、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( )； 答案 )(1)(1n n n n n x f x f x x x '--- =+ 5、对1)(3 ++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 )； 6、计算方法主要研究( 截断 )误差和( 舍入 )误差； 7、用二分法求非线性方程f (x )=0在区间(a ,b )内的根时，二分n 次后的误差限为 ( 1 2+-n a b )； 8、已知f (1)＝2，f (2)＝3，f (4)＝5.9，则二次Newton 插值多项式中x 2系数为( 0.15 )； 11、 两点式高斯型求积公式?1 d )(x x f ≈( ?++-≈1 )] 321 3()3213([21d )(f f x x f )，代数精度 为( 5 )；
12、 为了使计算 32)1(6)1(41310-- -+-+ =x x x y 的乘除法次数尽量地少，应将该表达 式改写为 11 ,))64(3(10-= -++=x t t t t y ，为了减少舍入误差，应将表达式1999 2001-改写为 199920012 + 。 13、 用二分法求方程01)(3 =-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间为 0.5，1 ,进行两步后根的所在区间为 0.5，0.75 。 14、 计算积分?1 5 .0d x x ,取4位有效数字。用梯形公式计算求得的近似值为 0.4268 ，用 辛卜生公式计算求得的近似值为 0.4309 ，梯形公式的代数精度为 1 ，辛卜生公式的代数精度为 3 。 15、 设46)2(,16)1(,0)0(===f f f ,则=)(1x l )2()(1--=x x x l ，)(x f 的二次牛顿插值 多项式为 )1(716)(2-+=x x x x N 。 16、 求积公式?∑=≈b a k n k k x f A x x f )(d )(0的代数精度以( 高斯型 )求积公式为最高，具有 ( 12+n )次代数精度。 17、 已知f (1)=1,f (3)=5,f (5)=-3,用辛普生求积公式求?5 1 d )(x x f ≈( 12 )。 18、 设f (1)=1， f (2)=2，f (3)=0，用三点式求≈')1(f ( 2.5 )。 19、如果用二分法求方程043 =-+x x 在区间]2,1[内的根精确到三位小数，需对分（ 10 ）次。 20、已知?????≤≤+-+-+-≤≤=31)1()1()1(2110)(2 33x c x b x a x x x x S 是三次样条函数，则 a =( 3 )，b =（ 3 ），c =（ 1 ）。 21、)(,),(),(10x l x l x l n 是以整数点n x x x ,,,10 为节点的Lagrange 插值基函数，则 ∑== n k k x l 0)(( 1 )，∑== n k k j k x l x 0 )(( j x )，当2≥n 时 = ++∑=)()3(20 4x l x x k k n k k ( 32 4++x x )。 22、区间[]b a ,上的三次样条插值函数)(x S 在[]b a ,上具有直到_____2_____阶的连续导数。
数值分析试题及答案
数值分析试题及答案 一、单项选择题（每小题3分，共15分） 1. 3.142和3.141分别作为的近似数具有（）和（）位有效数字. A．4和3 B．3和2 C．3和4 D．4和4 2. 已知求积公式，则＝（） A． B．C．D． 3. 通过点的拉格朗日插值基函数满足（） A．＝0，B．＝0， C．＝1，D．＝1， 4. 设求方程的根的牛顿法收敛，则它具有（）敛速。 A．超线性B．平方C．线性D．三次 5. 用列主元消元法解线性方程组作第一次消元后得到的第3个方程（）. A．B． C．D． 单项选择题答案 1.A 2.D 3.D 4.C 5.B 得分评卷 人 二、填空题（每小题3分，共15分） 1. 设, 则， . 2. 一阶均差 3. 已知时，科茨系数，那么 4. 因为方程在区间上满足，所以在区间内有根。 5. 取步长，用欧拉法解初值问题的计算公式.填空题答案
1. 9和 2. 3. 4. 5. 得分评卷 人 三、计算题（每题15分，共60分） 1. 已知函数的一组数据：求分段线性插值函数，并计算的近似值. 计算题1.答案 1. 解， ， 所以分段线性插值函数为 2. 已知线性方程组 （1）写出雅可比迭代公式、高斯－塞德尔迭代公式； （2）对于初始值，应用雅可比迭代公式、高斯－塞德尔迭代公式分别计算（保留小数点后五位数字）. 计算题2.答案 1.解原方程组同解变形为 雅可比迭代公式为 高斯－塞德尔迭代法公式 用雅可比迭代公式得 用高斯－塞德尔迭代公式得 3. 用牛顿法求方程在之间的近似根 （1）请指出为什么初值应取2？ （2）请用牛顿法求出近似根，精确到0.0001. 计算题3.答案
数值分析期末试题
数值分析期末试题 一、填空题（20102=?分） （1)设??? ? ? ??? ??---=28 3 012 251A ，则=∞ A ______13_______。 （2)对于方程组?? ?=-=-3 4101522121x x x x ，Jacobi 迭代法的迭代矩阵是=J B ?? ? ? ??05.25.20。 （3)3*x 的相对误差约是*x 的相对误差的 3 1倍。 （4）求方程)(x f x =根的牛顿迭代公式是) ('1)(1n n n n n x f x f x x x +-- =+。 （5）设1)(3 -+=x x x f ，则差商=]3,2,1,0[f 1 。 （6）设n n ?矩阵G 的特征值是n λλλ,,,21 ，则矩阵G 的谱半径=)(G ρi n i λ≤≤1max 。 （7）已知?? ? ? ??=1021 A ，则条件数=∞ )(A Cond 9 （8）为了提高数值计算精度，当正数x 充分大时,应将)1ln(2 -- x x 改写为 )1ln(2 ++ -x x 。 （9）n 个求积节点的插值型求积公式的代数精确度至少为1-n 次。 （10）拟合三点))(,(11x f x ，))(,(22x f x ，))(,(33x f x 的水平直线是)(3 1 3 1 ∑== i i x f y 。 二、（10分）证明：方程组? ?? ??=-+=++=+-1 211 2321321321x x x x x x x x x 使用Jacobi 迭代法求解不收敛性。 证明：Jacobi 迭代法的迭代矩阵为 ???? ? ?????---=05 .05 .01015.05.00J B J B 的特征多项式为
2012研究生数值分析课期末考试复习题及答案
一、填空 1. 设 2.3149541...x * =，取5位有效数字，则所得的近似值x= 2.3150 . 2.设一阶差商 ()()()21122114 ,321f x f x f x x x x --= = =---， ()()()322332 615 ,422f x f x f x x x x --= = =-- 则二阶差商 ()123,,______ f x x x =11/6 3. 设(2,3,1)T X =--, 则2||||X = 14 ，=∞||||X 3 。p49 4. 4．求方程 2 1.250x x --= 的近似根，用迭代公式 1.25x x =+，取初始值 01 x =， 那么 1______x =。 1.5 5．解初始值问题 00 '(,)()y f x y y x y =?? =?近似解的梯形公式是 1______k y +≈。 ()()[]11,,2 ++++k k k k k y x f y x f h y 6、 1151A ??= ? -??，则A 的谱半径 ＝ 6 。 7、设 2()35, , 0,1,2,... , k f x x x kh k =+== ，则 []12,,n n n f x x x ++= —————— ————3 和 []123,,,n n n n f x x x x +++= _______________0_____ 。 8、 若线性代数方程组AX=b 的系数矩阵A 为严格对角占优阵，则雅可比迭代和高斯-塞德尔迭代都 收敛 。 9、解常微分方程初值问题的欧拉（Euler ）方法的局部截断误差为_______O(h ） ___。
关于岩土工程的数值计算方法的综述
题目：关于岩土工程的数值计算方法的综述学院：资源与土木工程学院 专业：岩土工程 学号： 姓名：
关于岩土工程的数值计算方法的综述 我通过学习和查阅相关资料文献了解到，近年来，数值计算模拟分析在岩土工程中越来越受欢迎，随着城市的建设，地下工程所处的环境越来越复杂，影响的因素也是越来越多，所以依靠传统的解析计算难以实现，计算机的数值模拟恰恰解决的了岩土的计算的问题，它可以模拟各种复杂情况下岩土问题。就岩土工程而言，由于岩土介质涉及本构关系、力学参数、自身构造以及边界条件等的复杂多变性，在未采用计算机数值方法以前，对于复杂、重要的岩土工程，如果用传统的弹性力学或弹塑性力学的解析法难以求解时，只好采用物理模拟或其他方法从宏观上把握工程的受力和变形特征。随着计算机数值分析方法的出现和发展，情况发生了巨大的变化。计算机数值方法已经能够较好的模拟非均匀质体、各向异性介质面临的复杂边界条件问题，也能处理岩土工程中不连续性界面、渗流问题、岩土损伤断裂问题以及复杂的岩土工程结构分析问题，对于涉及时间因素的动力问题、蠕变问题，特别是耦合问题，数值模拟计算方法极大的加强了解决岩土工程的能力。 数值计算方法其主要有有限单元法、有限差分法、边界元法、离散元法和流形元法等。 有限单元法：有限单元法发展非常迅速，至今已经成为求解复杂工程问题的有力工具，并在岩土工程领域广泛的采用，主要的分析软件ANSYS。 有限单元法的最基本的元素是单元和节点，基本计算步骤的第一步为离散化，问题域的连续体被离散为单元与节点的组合，连续体内部分的应力及位移通过节点传递，每个单元可以具有不同的物理特征，这样，便可以得到在物理意义上与原来的连续体相近似的模型。第二步为单元分析，一般以位移法为基本方法，建立单元的刚度矩阵。第三步由单元的刚度矩阵集合成总体刚度矩阵，并由此建立系统的整体方程组。第四步进入计算模型的边界条件，求解方程组，求得节点位移。第五步求出各单元的应变、应力及主应力。 有限差分法：有限差分法在岩土工程中是应用非常广泛的方法，在数值计算模拟上有很大的贡献，主要的应用软件为FLAC3D。 基本思想是把连续的定解区域用有限个离散点构成的网格来代替，这些离散点称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似；把原方程和定解条件中的微商用差商来近似，积分用积分和来近似，于是原微分方
数值分析期末试题
一、（8分）用列主元素消去法解下列方程组： ??? ??=++-=+--=+-11 2123454 321321321x x x x x x x x x 二、（10分）依据下列数据构造插值多项式：y(0)=1,y(1)= —2,y '(0)=1, y '(1)=—4 三、（12分）分别用梯形公式和辛普生公式构造 复化的梯形公式、复化的辛普生公式并利用复化的梯形公式、复化的辛普生公式计算下列积分： ? 9 1dx x n=4 四、（10分）证明对任意参数t ，下列龙格－库塔方法是二阶的。 五、（14分）用牛顿法构造求c 公式，并利用牛顿法求115。保留有效数字五位。 六、（10分）方程组AX=B 其中A=????????? ?10101a a a a 试就AX=B 建立雅可比迭代法和高斯-赛德尔迭代法，并讨论a 取何值时 迭代收斂。 七、（10分）试确定常数A,B,C,a,使得数值积分公式?-++-≈2 2 ) (}0{)()(a Cf Bf a Af dx x f 有尽可能多的 代数精确度。并求该公式的代数精确度。 八、{6分} 证明： A ≤ 其中A 为矩阵，V 为向量. 第二套 一、（8分）用列主元素消去法解下列方程组： ??? ??=++=+-=+3 2221 43321 32132x x x x x x x x 二、（12分）依据下列数据构造插值多项式：y(0)=y '(0)=0, y(1)=y '(1)= 1,y(2)=1 三、（14分）分别用梯形公式和辛普生公式构造 复化的梯形公式、复化的辛普生公式,并利用复化的梯形公式、 复化的辛普生公式及其下表计算下列积分： ?2 /0 sin πxdx ????? ? ? -+-+=++==++=+1 3121231)1(,)1(() ,(),()(2 hk t y h t x f k thk y th x f k y x f k k k h y y n n n n n n n n