西北工业大学计算方法第二周作业答案
- 格式:pdf
- 大小:254.34 KB
- 文档页数:6
计算方法课后习题答案计算方法课后习题答案计算方法是一门重要的学科,它为我们提供了解决数学问题的方法和工具。
在学习这门课程时,我们经常会遇到一些习题,这些习题旨在帮助我们巩固所学的知识并提高我们的计算能力。
然而,习题的解答并非总是容易的,有时候我们可能会遇到困难。
因此,我将在本文中为大家提供一些计算方法课后习题的答案,希望能够帮助大家更好地理解和掌握这门学科。
1. 线性方程组的解法线性方程组是计算方法中的一个重要概念。
解决线性方程组的方法有很多种,其中最常用的方法是高斯消元法。
这种方法通过行变换将线性方程组转化为简化的行阶梯形矩阵,从而求得方程组的解。
下面是一个例子:2x + 3y = 84x - 5y = -7通过高斯消元法,我们可以得到方程组的解为x = 1,y = 2。
2. 数值积分的计算数值积分是计算方法中的另一个重要概念。
它可以用来计算曲线下的面积或者求解定积分。
常用的数值积分方法有梯形法则、辛普森法则等。
下面是一个例子:计算定积分∫[0, 1] (x^2 + 2x)dx。
通过梯形法则,我们可以得到定积分的近似值为1.5。
3. 插值和拟合插值和拟合是计算方法中的重要概念,它们可以用来估计未知数据点的值。
插值是通过已知数据点之间的连线或曲线来估计未知点的值,而拟合是通过已知数据点的函数来估计未知点的值。
下面是一个例子:已知数据点 (1, 3), (2, 5), (3, 8),通过插值和拟合方法来估计点 (4, ?) 的值。
通过线性插值,我们可以得到点 (4, 11) 的值。
通过多项式拟合,我们可以得到点 (4, 10.5) 的值。
4. 数值微分的计算数值微分是计算方法中的另一个重要概念,它可以用来估计函数在某一点的导数值。
常用的数值微分方法有前向差分法、后向差分法和中心差分法。
下面是一个例子:计算函数 f(x) = x^2 在点 x = 2 处的导数值。
通过中心差分法,我们可以得到导数的近似值为 4。
参考答案 第一章1 *1x =1.7; *2x =1.73; *3x =1.732 。
2.3. (1) ≤++)(*3*2*1x x x e r 0.00050; (注意:应该用相对误差的定义去求) (2) ≤)(*3*2*1x x x e r 0.50517; (3) ≤)/(*4*2x x e r 0.50002。
4.设6有n 位有效数字,由6≈2.4494……,知6的第一位有效数字1a =2。
令3)1()1(1*1021102211021)(-----⨯≤⨯⨯=⨯=n n r a x ε 可求得满足上述不等式的最小正整数n =4,即至少取四位有效数字,故满足精度要求可取6≈2.449。
5. 答:(1)*x (0>x )的相对误差约是*x 的相对误差的1/2倍;(2)nx )(* 的相对误差约是*x 的相对误差的n 倍。
6. 根据********************sin 21)(cos 21sin 21)(sin 21sin 21)(sin 21)(c b a c e c b a c b a b e c a c b a a e c b S e r ++≤=******)()()(tgcc e b b e a a e ++ 注意当20*π<<c 时,0**>>c tgc ,即1*1*)()(--<c tgc 。
则有)()()()(****c e b e a e S e r r r r ++<7.设20=y ,41.1*0=y ,δ=⨯≤--2*001021y y 由 δ1*001*111010--≤-=-y y y y ,δ2*111*221010--≤-=-y y y yMδ10*991*10101010--≤-=-y y y y即当0y 有初始误差δ时,10y 的绝对误差的绝对值将减小1010-倍。
而11010<<-δ,故计算过程稳定。
机械原理课后习题答案第二章 机构的结构分析题2-11 图a 所示为一简易冲床的初拟设计方案。
设计者的思路是:动力由齿轮1输入,使轴A 连续回转;而固装在轴A 上的凸轮2与杠杆3组成的凸轮机构使冲头4上下运动,以达到冲压的目的。
试绘出其机构运动简图(各尺寸由图上量取),分析是否能实现设计意图,并提出修改方案。
解:1)取比例尺,绘制机构运动简图。
(图2-11a)2)要分析是否能实现设计意图,首先要计算机构的自由度。
尽管此机构有4个活动件,但齿轮1和凸轮2是固装在轴A 上,只能作为一个活动件,故 3=n 3=l p 1=h p01423323=-⨯-⨯=--=h l p p n F原动件数不等于自由度数,此简易冲床不能运动,即不能实现设计意图。
分析:因构件3、4与机架5和运动副B 、C 、D 组成不能运动的刚性桁架。
故需增加构件的自由度。
3)提出修改方案:可以在机构的适当位置增加一个活动构件和一个低副,或用一个高副来代替一个低副。
(1) 在构件3、4之间加一连杆及一个转动副(图2-11b)。
(2) 在构件3、4之间加一滑块及一个移动副(图2-11c)。
(3) 在构件3、4之间加一滚子(局部自由度)及一个平面高副(图2-11d)。
11(c)题2-11(d)54364(a)5325215436426(b)321讨论:增加机构自由度的方法一般是在适当位置上添加一个构件(相当于增加3个自由度)和1个低副(相当于引入2个约束),如图2-1(b )(c )所示,这样就相当于给机构增加了一个自由度。
用一个高副代替一个低副也可以增加机构自由度,如图2-1(d )所示。
题2-12 图a 所示为一小型压力机。
图上,齿轮1与偏心轮1’为同一构件,绕固定轴心O 连续转动。
在齿轮5上开有凸轮轮凹槽,摆杆4上的滚子6嵌在凹槽中,从而使摆杆4绕C 轴上下摆动。
同时,又通过偏心轮1’、连杆2、滑杆3使C 轴上下移动。
最后通过在摆杆4的叉槽中的滑块7和铰链G 使冲头8实现冲压运动。
《计算方法》习题答案第一章 数值计算中的误差1.什么是计算方法?(狭义解释)答:计算方法就是将所求的的数学问题简化为一系列的算术运算和逻辑运算,以便在计算机上编程上机,求出问题的数值解,并对算法的收敛性、稳定性和误差进行分析、计算。
2.一个实际问题利用计算机解决所采取的五个步骤是什么?答:一个实际问题当利用计算机来解决时,应采取以下五个步骤: 实际问题→建立数学模型→构造数值算法→编程上机→获得近似结果 4.利用秦九韶算法计算多项式4)(53-+-=x x x x P 在3-=x 处的值,并编程获得解。
解:400)(2345-+⋅+-⋅+=x x x x x x P ,从而所以,多项式4)(53-+-=x x x x P 在3-=x 处的值223)3(-=-P 。
5.叙述误差的种类及来源。
答:误差的种类及来源有如下四个方面:(1)模型误差:数学模型是对实际问题进行抽象,忽略一些次要因素简化得到的,它是原始问题的近似,即使数学模型能求出准确解,也与实际问题的真解不同,我们把数学模型与实际问题之间存在的误差称为模型误差。
(2)观测误差:在建模和具体运算过程中所用的一些原始数据往往都是通过观测、实验得来的,由于仪器的精密性,实验手段的局限性,周围环境的变化以及人们的工作态度和能力等因素,而使数据必然带有误差,这种误差称为观测误差。
(3)截断误差:理论上的精确值往往要求用无限次的运算才能得到,而实际运算时只能用有限次运算的结果来近似,这样引起的误差称为截断误差(或方法误差)。
(4)舍入误差:在数值计算过程中还会用到一些无穷小数,而计算机受机器字长的限制,它所能表示的数据只能是一定的有限数位,需要把数据按四舍五入成一定位数的近似的有理数来代替。
这样引起的误差称为舍入误差。
6.掌握绝对误差(限)和相对误差(限)的定义公式。
答:设*x 是某个量的精确值,x 是其近似值,则称差x x e -=*为近似值x 的绝对误差(简称误差)。
计算方法习题答案在数学和工程领域,计算方法是指解决数学问题的一系列算法和程序。
以下是一些常见的计算方法习题及其答案。
习题1:求解线性方程组考虑线性方程组:\[ \begin{align*}3x + 2y &= 7, \\4x - y &= 5.\end{align*} \]答案:使用高斯消元法,我们首先将第二个方程乘以2,然后从第一个方程中减去得到:\[ \begin{align*}3x + 2y &= 7, \\0x + 9y &= 17.\end{align*} \]解得 \( y = \frac{17}{9} \)。
将 \( y \) 的值代入第一个方程,解得 \( x = 1 \)。
因此,解为 \( x = 1, y = \frac{17}{9} \)。
习题2:数值积分给定函数 \( f(x) = x^2 \),求在区间 [0, 1] 上的积分。
答案:使用梯形法则进行数值积分,取两个子区间:\[ \int_{0}^{1} x^2 dx \approx \frac{1}{2} \left( f(0) + f(1) \right) = \frac{1}{2} \left( 0 + 1 \right) = 0.5. \]习题3:求解常微分方程的初值问题考虑初值问题:\[ y' = 3x^2 - 2y, \quad y(0) = 1. \]答案:使用欧拉方法,取步长 \( h = 0.1 \),计算 \( y \) 的值:\[ y_{n+1} = y_n + h f(x_n, y_n). \]从 \( y_0 = 1 \) 开始,计算得到:\[ y_1 = 1 + 0.1(0 - 2) = 1.2, \]\[ y_2 = 1.2 + 0.1(0.01 - 2.4) = 1.4, \]以此类推,可以得到 \( y \) 在区间 [0, 1] 上的近似值。
习题4:数值解非线性方程给定方程 \( f(x) = x^3 - x - 1 = 0 \),求根。
2002-2003第一学期一.计算及推导(5*8)1.已知* 3.141,x x π==,试确定*x 近似x 的有效数字位数。
2.有效数***1233.105,0.001,0.100x x x =-==,试确定***123x x x ++的相对误差限。
3.已知3()0.50.12f x x x =++,试计算差商[]0,1,2,3f 4.给出拟合三点(0,1),(1,0)A B ==和(1,1)C =的直线方程。
5.推导中矩形求积公式''31()()()()()224b aa b f x dx b a f f b a η+=-+-⎰ 6.试证明插值型求积公式()()nbi i ai f x dx A f x =≈∑⎰的代数精确度至少是n 次。
7.已知非线性方程()x f x =在区间[],a b内有一实根,试写出该实根的牛顿迭代公式。
8.用三角分解法求解线性方程组123121022331302x x x ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦要用二次插值多项式计算(0.63891)f 的近似值,试选择合适的插值节点进行计算,并说明所选用节点依据。
(保留5位有效数字)(12分) 三. 已知方程ln 0x x +=在(0,1)内有一实根α(1)给出求该实根的一个迭代公式,试之对任意的初始近似0(0,1)x ∈迭代法都收敛,并证明其收敛性。
(2)00.5x =试用构造的迭代公式计算α的近似值n x ,要求3110n n x x ---≤。
四. 设有方程组112233131232a x b a x b a x b ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦当参数a 满足什么条件时,雅可比方法对任意的初始向量都收敛。
写出与雅可比方法对应的高斯赛德尔迭代公式。
(12分) 五.用欧拉预估校正法求解初值问题 '2 (00.2)(0)1x y y x y y ⎧=-≤≤⎪⎨⎪=⎩ 取h=0.1,小数点后保留5位。