(A)1. 3.142和3.141分别作为π的近似数具有（）和（）为有效数字(有效数字) A. 4和3 B. 3和2 C. 3和4 D. 4和4 解..14159.3==*πx ,1103142.0?=a 时,1=m ,3102 1...00041.0)(-*?≤ =-=a x a E m-n= -3,所以n=4,即有4位有效数字。当1103141.0?=a 时,1=m , 2102 1005.0...00059.0)(-*?=≤=-=a x a E ，m-n= -2,所以n=3,即有3位有效数字。 (A)2. 为了减少误差，在计算表达式19992001-时，应该改为 199920012+计算，是属于（）来避免误差。（避免误差危害原则） A.避免两相近数相减； B.化简步骤，减少运算次数； C.避免绝对值很小的数做除数； D.防止大数吃小数 解：由于2001和1999相近，两数相减会使误差大，因此化加法为减法，用的方法是避免误差危害原则。 (B)3.下列算式中哪一个没有违背避免误差危害原则（避免误差危害原则） A.计算123460.60.612345++- B.计算 25612520000450?- C.计算10.99994- D.计算11x x +- 解:A 会有大数吃掉小数的情况C 中两个相近的数相减，D 中两个相近的数相减也会增大误差 (D)4.若误差限为5105.0-?,那么近似数0.003400有()位有效数字。(有效数字) A. 5 B. 4 C. 7 D. 3 解：51021)(-?= a E 即m-n= -5,2103400.0-?=a ，m= -2，所以n=3，即有3位有效数字 (A)5.设*x 的近似数为40.32710a =?，如果a 具有3位有效数字，则a 的相对误差限为 ()（有效数字与相对误差的关系） A ． 35103-g B. 33105-g C. 53105-g D. 5103 g -2 解：因为40.32710a =?所以31=a ,因为a 有3位有效数字，所以n=3,由相对误差和有效数字的关系可得a 的相对误差限为 31103510.5--?== n r a δ
大作业 三 1. 给定初值 0x 及容许误差 ，编制牛顿法解方程f (x )=0的通用 程序. 解：Matlab 程序如下： 函数m 文件：fu.m function Fu=fu(x) Fu=x^3/3-x; end 函数m 文件：dfu.m function Fu=dfu(x) Fu=x^2-1; end 用Newton 法求根的通用程序Newton.m clear; x0=input('请输入初值x0:'); ep=input('请输入容许误差:'); flag=1; while flag==1 x1=x0-fu(x0)/dfu(x0); if abs(x1-x0) while flag1==1 && m<=10^3 x1=x0-fu(x0)/dfu(x0); if abs(x1-x0) 中考物理力学综合计算题 1.如图24所示,质量为60kg 的工人在水平地面上,用滑轮组把货物运到高处。第一次运送货 物时,货物质量为130kg,工人用力F 1匀速拉绳,地面对工人的支持力为N 1,滑轮组的机械效率为η1;第二次运送货物时,货物质量为90 kg,工人用力F 2匀速拉绳的功率为P 2,货箱以0.1m/s 的速度匀速上升,地面对人的支持力为N 2, N 1与 N 2之比为2:3。(不计绳重及滑轮摩擦, g 取10N/kg) 求:(1)动滑轮重和力F 1的大小; (2)机械效率η1; (3) 功率P 2。 2.火车与公路交叉处设置人工控制的栏杆，图22是栏杆的示意图。栏杆全长AB =6m ，在栏杆的左端安装配重，使栏杆和配重总体的重心位于O 点。栏杆的P 点安装转轴，转轴与支架C 连结，使栏杆能绕P 在竖直平面无摩擦转动，支架C 用两块木板做成，中间空隙可以容纳栏杆。栏杆的B 端搁置在支架D 上，当支架D 上受到压力为F D 时，栏杆恰好在水平位置平衡。当体重为G 人的管理人员双脚站在水平地面时，他对地面的压强是p 1；当他用力F 1竖直向下压A 端，使栏杆的B 端刚好离开支架，此时人双脚对地面的压强是p 2。管理人员继续用力可使栏杆逆时针转动至竖直位置，并靠在支架C 上。火车要通过时，他要在A 端用力F 2使栏杆由竖直位置开始离开支架C ，使栏杆能顺时针转动直至栏杆B 端又搁置在支架D 上。已知AP =OP =1m ，PE = 2 3 m ，O 点到栏杆下边缘的距离OE =0.5m ，p 1∶p 2=2∶1，栏杆与配重的总重G 杆=2403N 。 求：（1）F D （2）G 人 （3）F 2的最小值，此时F 2的方向。（计算和结果可带根号）（6分） 图24 《计算方法》习题答案 第一章 数值计算中的误差 1．什么是计算方法？（狭义解释） 答：计算方法就是将所求的的数学问题简化为一系列的算术运算和逻辑运算，以便在计算机上编程上机，求出问题的数值解，并对算法的收敛性、稳定性和误差进行分析、计算。 2．一个实际问题利用计算机解决所采取的五个步骤是什么？ 答：一个实际问题当利用计算机来解决时，应采取以下五个步骤： 实际问题→建立数学模型→构造数值算法→编程上机→获得近似结果 4．利用秦九韶算法计算多项式4)(5 3 -+-=x x x x P 在3-=x 处的值，并编程获得解。 解：400)(2 3 4 5 -+?+-?+=x x x x x x P ，从而 所以，多项式4)(5 3 -+-=x x x x P 在3-=x 处的值223)3(-=-P 。 5．叙述误差的种类及来源。 答：误差的种类及来源有如下四个方面： （1）模型误差：数学模型是对实际问题进行抽象，忽略一些次要因素简化得到的，它是原始问题的近似，即使数学模型能求出准确解，也与实际问题的真解不同，我们把数学模型与实际问题之间存在的误差称为模型误差。 （2）观测误差：在建模和具体运算过程中所用的一些原始数据往往都是通过观测、实验得来的，由于仪器的精密性，实验手段的局限性，周围环境的变化以及人们的工作态度和能力等因素，而使数据必然带有误差，这种误差称为观测误差。 （3）截断误差：理论上的精确值往往要求用无限次的运算才能得到，而实际运算时只能用有限次运算的结果来近似，这样引起的误差称为截断误差（或方法误差）。 （4）舍入误差：在数值计算过程中还会用到一些无穷小数，而计算机受机器字长的限制，它所能表示的数据只能是一定的有限数位，需要把数据按四舍五入成一定位数的近似的有理数来代替。这样引起的误差称为舍入误差。 6．掌握绝对误差（限）和相对误差（限）的定义公式。 答：设* x 是某个量的精确值，x 是其近似值，则称差x x e -=* 为近似值x 的绝对误差（简称误差）。若存在一个正数ε使ε≤-=x x e * ，称这个数ε为近似值x 的绝对误差限（简称误差限或精度）。 把绝对误差e 与精确值* x 之比* **x x x x e e r -==称为近似值x 的相对误差，称 第一章 1、ln2＝0.69314718…，精确到 10－3 的近似值是多少？ 解 精确到 10－3＝0.001，即绝对误差限是 e ＝0.05％，故至少要保留小数点后三位才可以。 ln2≈0.693。 2、设115.80,1025.621≈≈x x 均具有5位有效数字，试估计由这些数据计算21x x ， 21x x +的绝对误差限 解：记126.1025, 80.115x x == 则有11232411 10, | 102|||2 x x x x --≤?-≤?- 所以 121212121212211122||||||||||||x x x x x x x x x x x x x x x x x x -=-+-+≤-- 3411 80.11610 6.10102522 0.007057-==??+≤?? 1212112243|()|||11 |10100.0005522 |x x x x x x x x --≤≤?+?=+-+-+- 3、一个园柱体的工件,直径d 为10.250.25mm,高h 为40.00 1.00mm,则它的体 积V 的近似值、误差和相对误差为多少。 解： ()() 22222222 4 314210254000000330064 221025400002510251002436444 3300624362436 0073873833006 , .....; ()()()......, ..().()..% .r d h V d h V mm d h V dh d d h V mm V V V πππππεεεεε= ≈=??===+=???+?==±====第二章： 1、分别利用下面四个点的Lagrange 插值多项式和Newton 插值多项式N 3(x )， 计算L 3(0.5)及N 3(-0.5) x －2 －1 0 1 f (x ) －1 1 2 大作业 三 1. 给定初值 0x 及容许误差 ，编制牛顿法解方程f (x )=0的通用程序. 解：Matlab 程序如下： 函数m 文件：fu.m function Fu=fu(x) Fu=x^3/3-x; end 函数m 文件：dfu.m function Fu=dfu(x) Fu=x^2-1; end 用Newton 法求根的通用程序Newton.m clear; x0=input('请输入初值x0:'); ep=input('请输入容许误差:'); flag=1; while flag==1 x1=x0-fu(x0)/dfu(x0); if abs(x1-x0) while flag==1 sigma=k*eps; x0=sigma; k=k+1; m=0; flag1=1; while flag1==1 && m<=10^3 x1=x0-fu(x0)/dfu(x0); if abs(x1-x0) end end fprintf('最大的sigma 值为：%f\n',sigma); 2．求下列方程的非零根 5130.6651()ln 05130.665114000.0918 x x f x x +?? =-= ?-???解： Matlab 程序为： （1）主程序 clear clc format long x0=765; N=100; errorlim=10^(-5); x=x0-f(x0)/subs(df(),x0); n=1; 24．（20分） 雨滴落到地面的速度通常仅为几米每秒，这与雨滴下落过程中受到空气阻力有关。雨滴间无相互作用且雨滴质量不变，重力加速度为g 。 （1）质量为m 的雨滴由静止开始，下落高度h 时速度为u ，求这一过程中克服空气阻力 所做的功W 。 （2）将雨滴看作半径为r 的球体，设其竖直落向地面的过程中所受空气阻力 f = kr 2v 2， 其中v 是雨滴的速度，k 是比例系数。 a ．设雨滴的密度为ρ，推导雨滴下落趋近的最大速度v m 与半径r 的关系式； b ．示意图中画出了半径为r 1、r 2（r 1> r 2）的雨滴在空气中无初速下落的v -t 图线，其中______对应半径为r 1的雨滴（选填①、②）；若不计空气阻力，请在图中画出雨滴无初速下落的v -t 图线。 （3）由于大量气体分子在各方向运动的几率相等，其对静止雨滴的作用力为零。将雨滴 简化为垂直于运动方向面积为S 的圆盘，证明：圆盘以速度v 下落时受到的空气阻力f ∝2v （提示：设单位体积内空气分子数为n ，空气分子质量为m 0）。 22．（16分） 2022年将在我国举办第二十四届冬奥会，跳台滑雪是其中最具观赏性的项目之一。某滑道示意图如下，长直助滑道AB 与弯曲滑道BC 平滑衔接，滑道BC 高h = 10 m ，C 是半径R = 20 m 圆弧的最低点。质量m = 60 kg 的运动员从A 处由静止开始匀加速下滑，加速度a = 4.5 m/s 2，到达B 点时速度v B = 30 m/s 。取重力加速度210m/s g =。 （1）求长直助滑道AB 的长度L ； （2）求运动员在AB 段所受合外力的冲量I 的大小； （3）若不计BC 段的阻力，画出运动员经过C 点时的受力 图，并求其所受支持力N F 的大小。 B h C A 第一章误差 1问3.142, 3.141, 22分别作为n 的近似值各具有几位有效数字? 7 分析 利用有效数字的概念可直接得出。 解 n =3.141 592 65 … =3.141 59 …-3.142 85…=-0.001 26 …知 2 10 3 | 22| 1 10 2 因而X 3具有3位有效数字。 2 已知近似数X*有两位有效数字，试求其相对误差限。 分析本题显然应利用有效数字与相对误差的关系。 解 利用有效数字与相对误差的关系。这里 n=2,a 1是1到9之间的数字。 分析本题利用有效数字与相对误差的关系。 解 a 1是1到9间的数字。 * (X)0.3 % 1000 2 102 2 (9 1) ? 2(a? 1) 10' 设x*具有n 位有效数字，令-n+仁-1，则n=2,从而x*至少具有2位有效数字。 4计算sin 1.2，问要取几位有效数字才能保证相对误差限不大于 0.01%。 分析本题应利用有效数字与相对误差的关系。 解 设取n 位有效数字，由sin 1.2=0.93…，故a 1=9。 *(x) | 悩盍 10n1 o.。1% 104 解不等式丄 10 n 1 10 4知取n=4即可满足要求。 2ai 5 计算盂盘，视已知数为精确值，用4位浮点数计算。 因而 因而 记 X 1=3.142, X 2=3.141 , X 3= ^2 . 由 n - X 1=3.141 59 …-3.142=-0.000 40 …知 1 103 | 2 X 1具有4位有效数字。 由 n - X 2=3.141 59…-3.141=-0.000 59 …知 1 103 I X 2| X 2具有3位有效数字。 104 10 2 22 7 I *(x)| | X X* I |X*| 1 1 10n1 1021 5 % 已知近似数的相对误差限为 0.3%，问X*至少有几位有效数字? （适合课程《数值方法A 》和《数值方法B 》） 第一章 绪 论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2％,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位 有效数字: ***** 123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====? 4. 利用公式(3.3)求下列各近似值的误差限: * * * * * * * * 12412324(),(),()/,i x x x ii x x x iii x x ++其中* * * * 1234,,,x x x x 均为第3题所给的数. 5. 计算球体积要使相对误差限为1％,问度量半径R 时允许的相对误差限是多少? 6. 设028,Y =按递推公式 11783 100 n n Y Y -=- ( n=1,2,…) 计算到100Y .若取783≈27.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字(783≈27.982). 8. 当N 充分大时,怎样求 2 11N dx x +∞+?? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设2 12S gt = 假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对误差增加, 而相对误差却减小. 11. 序列{}n y 满足递推关系1101 n n y y -=-(n=1,2,…),若02 1.41y =≈(三位有效数字),计算到10 y 时误差有多大?这个计算过程稳定吗? 12. 计算6 (21)f =-,取 2 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 6 3 11,(322), ,9970 2. (21) (322) --++ 13. 2 ()ln(1)f x x x =- -,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等 价公式 2 2 ln(1)ln(1)x x x x - -=-+ + 计算,求对数时误差有多大? 14. 试用消元法解方程组{ 10 10 12121010; 2. x x x x +=+=假定只用三位数计算,问结果是否可靠? 15. 已知三角形面积 1sin , 2 s ab c = 其中c 为弧度, 02c π << ,且测量a ,b ,c 的误差分别为,,.a b c ???证 明面积的误差s ?满足 . s a b c s a b c ????≤ ++ 第二章 插值法 1. 根据( 2.2)定义的范德蒙行列式,令 目录 第一章非线性方程求根 (3) 1.1迭代法 (3) 1.2牛顿法 (4) 1.3弦截法 (5) 1.4二分法 (6) 第二章插值 (7) 2.1线性插值 (7) 2.2二次插值 (8) 2.3拉格朗日插值 (9) 2.4分段线性插值 (10) 2.5分段二次插值 (11) 第三章数值积分 (13) 3.1复化矩形积分法 (13) 3.2复化梯形积分法 (14) 3.3辛普森积分法 (15) 3.4变步长梯形积分法 (16) 第四章线性方程组数值法 (17) 4.1约当消去法 (17) 4.2高斯消去法 (18) 4.3三角分解法 (20) 4.4雅可比迭代法 (21) 4.5高斯—赛德尔迭代法 (23) 第五章常积分方程数值法 (25) 5.1显示欧拉公式法 (25) 5.2欧拉公式预测校正法 (26) 5.3改进欧拉公式法 (27) 5.4四阶龙格—库塔法 (28) 数值计算方法 第一章非线性方程求根 1.1迭代法 程序代码： Private Sub Command1_Click() x0 = Val(InputBox("请输入初始值x0")) ep = Val(InputBox(请输入误差限ep)) f = 0 While f = 0 X1 = (Exp(2 * x0) - x0) / 5 If Abs(X1 - x0) < ep Then Print X1 f = 1 Else x0 = X1 End If Wend End Sub 例：求f(x)=e2x-6x=0在x=0.5附近的根（ep=10-10） 1.2牛顿法 程序代码： Private Sub Command1_Click() b = Val(InputBox("请输入被开方数x0")) ep = Val(InputBox(请输入误差限ep)) f = 0 While f = 0 X1 = x0 - (x0 ^ 2 - b) / (2 * b) If Abs(X1 - x0) < ep Then Print X1 f = 1 Else x0 = X1 End If Wend End Sub 例：求56的值。（ep=10-10） 1、 0.1％,要取几位有效数字？ （ c ） (a) 2 (b) 3 (c) 4 (d) 5 2、若* 12.30x =是经过四舍五入得到的近似数，则它有几位有效数字？ （ c ） (a) 2 (b) 3 (c) 4 (d) 5 3、已知n ＋1个互异节点(x 0,y 0), (x 1,y 1),…, (x n ,y n )和过这些点的拉格朗日插值基函数l k (x )(k =0,1,2,…,n )，且ω(x )=(x －x 0) (x －x 1)… (x －x n )．则n 阶差商f (x 0,x 1,…, x n )= ( ) (a) ∑=n k k k y x l 0 )( (b) ∑='n k k k k x l y 0)( (c) ∑=n k k k x y 0)(ω (d) ∑='n k k k x y 0)(ω 4、已知由数据（0，0），（0.5，y ），（1，3），（2，2）构造出的三次插值多项式 33()6 P x x y 的 的系数是，则 等于 （ ） (a) -1.5 (b) 1 (c) 5.5 (d) 4.25 5、设(0,1,2,3,4)i x i =为互异结点，()i l x 为拉格朗日插值基函数,则 4 2 () ()i i i x x l x =-∑等于 （ a ） (a) 0 (b) 1 (c) 2 (d) 4 4()[,],()()(),()(),( )(), ' () ' (),22 ()()_________________________f x C a b H x a b a b H a f a H b f b H f H a f a f x H x ∈++====-=设是满足下列插值条件的三次多项式：则插值余项 1、 是以0，1，2为节点的三次样条函数，则b=-2,c=3 2、 已知(1)0,(1)3,(2)4,f f f =-=-=写出()f x 的牛顿插值多项式 2()P x =___2537 623x x +-__,其余项表达式 R(x)=__() (1)(1)(4) [1,4]6 f x x x ξξ'''-+-∈-_______________________ 3、 确定求积公式1 0121 ()(1)(0)'(1)f x dx A f A f A f -≈-++? 中的待定参数，使其代数精度 尽量高,则A 0=_ 29__________, A 1=__169________, A 2=_29 _______,代数精度=__2_________。 2020年中考物理力学综合计算题 1.如图24所示,质量为60kg 的工人在水平地面上,用滑轮组把货物运到高处。第一次运送货 物时,货物质量为130kg,工人用力F 1匀速拉绳,地面对工人的支持力为N 1,滑轮组的机械效率为η1;第二次运送货物时,货物质量为90 kg,工人用力F 2匀速拉绳的功率为P 2,货箱以0.1m/s 的速度匀速上升,地面对人的支持力为N 2, N 1与 N 2之比为2:3。(不计绳重及滑轮摩擦, g 取10N/kg) 求:(1)动滑轮重和力F 1的大小; (2)机械效率η1; (3) 功率P 2。 2.火车与公路交叉处设置人工控制的栏杆，图22是栏杆的示意图。栏杆全长AB =6m ，在栏杆的左端安装配重，使栏杆和配重总体的重心位于O 点。栏杆的P 点安装转轴，转轴与支架C 连结，使栏杆能绕P 在竖直平面无摩擦转动，支架C 用两块木板做成，中间空隙可以容纳栏杆。栏杆的B 端搁置在支架D 上，当支架D 上受到压力为F D 时，栏杆恰好在水平位置平衡。当体重为G 人的管理人员双脚站在水平地面时，他对地面的压强是p 1；当他用力F 1竖直向下压A 端，使栏杆的B 端刚好离开支架，此时人双脚对地面的压强是p 2。管理人员继续用力可使栏杆逆时针转动至竖直位置，并靠在支架C 上。火车要通过时，他要在A 端用力F 2使栏杆由竖直位置开始离开支架C ，使栏杆能顺时针转动直至栏杆B 端又搁置在支架D 上。已知AP =OP =1m ，PE = 2 3 m ，O 点到栏杆下边缘的距离OE =0.5m ，p 1∶p 2=2∶1，栏杆与配重的总重G 杆=2403N 。 求：（1）F D （2）G 人 （3）F 2的最小值，此时F 2的方向。（计算和结果可带根号）（6分） 图24 第一章 误差 1 问3.142，3.141，7 22分别作为π的近似值各具有几位有效数字？ 分析 利用有效数字的概念可直接得出。 解 π=3.141 592 65… 记x 1=3.142，x 2=3.141，x 3=7 22. 由π- x 1=3.141 59…-3.142=-0.000 40…知 3411110||1022 x π--?<-≤? 因而x 1具有4位有效数字。 由π- x 2=3.141 59…-3.141=-0.000 59…知 223102 1||1021--?≤- 习题1 1． 填空题 (1) 为便于算法在计算机上实现，必须将一个数学问题分解为 的 运算； (2) 在数值计算中为避免损失有效数字，尽量避免两个 数作减法运算；为避免 误差的扩大，也尽量避免分母的绝对值 分子的绝对值； (3) 误差有四大来源,数值分析主要处理其中的 和 ； (4) 有效数字越多,相对误差越 ； 2. 用例1.4的算法计算10,迭代3次,计算结果保留4位有效数字. 3. 推导开平方运算的误差限公式,并说明什么情况下结果误差不大于自变量误差. 4. 以下各数都是对准确值进行四舍五入得到的近似数,指出它们的有效数位、误差限和相对误差限. 95123450304051104000003346087510., ., , ., .x x x x x -==?===? 5. 证明1.2.3之定理1.1. 6. 若钢珠的的直径d 的相对误差为1.0%,则它的体积V 的相对误差将为多少。（假定钢珠为标准的球形） 7. 若跑道长的测量有0.1%的误差,对400m 成绩为60s 的运动员的成绩将会带来多大的误差和相对误差. 8. 为使20的近似数相对误差小于0.05%,试问该保留几位有效数字. 9. 一个园柱体的工件,直径d 为10.25±0.25mm,高h 为40.00±1.00mm,则它的体积V 的近似值、误差和相对误差为多少． 10 证明对一元函数运算有 r r xf x f x k x k f x εε'≈= () (())(),() 其中 并求出157f x x x ==()tan ,.时的k 值，从而说明f x x =()tan 在2 x π ≈时是病态问题． 11. 定义多元函数运算 1 1 1,,(),n n i i i i i i S c x c x εε====≤∑∑其中 求出S ε()的表达式，并说明i c 全为正数时，计算是稳定的，i c 有正有负时，误差难以控制． 12. 下列各式应如何改进,使计算更准确： 数值分析报大作业 班级：铁道2班 专业：道路与铁道工程 姓名：蔡敦锦 学号：13011260 一、序言 该数值分析大作业是通过C语言程序编程在Microsoft Visual C++ 6.0编程软件上运行实现的。本来是打算用Matlab软间来计算非线性方程的根的。学习Matlab也差不多有一个多月了，感觉自己编程做题应该没什么问题了；但是当自己真心的去编程、运行时才发现有很多错误，花了一天时间修改、调试程序都没能得到自己满意的结果。所以，我选择了自己比较熟悉的C程序语言来编程解决非线性的求值问题，由于本作业是为了比较几种方法求值问题的收敛速度和精度的差异，选择了一个相对常见的非线性函数来反映其差异，程序运行所得结果我个人比较满意。编写C语言，感觉比较上手，程序出现问题也能比较熟练的解决。最终就决定上交一份C程序语言编程的求值程序了！ 二、选题 本作业的目的是为了加深对非线性方程求根方法的二分法、简单迭代法、、牛顿迭代法弦截法等的构造过程的理解；能将各种方法的算法描述正确并且能够改编为程序并在计算机上实现程序的正确合理的运行，能得到自己满意的结果，并且能调试修改程序中可能出现的问题和程序功能的增减修改。本次程序是为了比较各种方法在求解同一非线性方程根时，在收敛情况上的差异。 为了达到上面的条件我选择自己比较熟悉的语言—C语言来编程，所选题目为计算方程f(x)=x3-2x-5=0在区间[2,3]内其最后两近似值的差的绝对值小于等于5 ?的根的几种方法的比较。 110- 本文将二分法、牛顿法、简单迭代法、弦截法及加速收敛法这五种方法在同一个程序中以函数调用的方式来实现，比较简洁明了，所得结果能很好的比较，便于分析；发现问题和得出结论。 第一章 误差 1 问，，7 22分别作为π的近似值各具有几位有效数字？ 分析 利用有效数字的概念可直接得出。 解 π= 592 65… 记x 1=，x 2=，x 3=7 22. 由π- x 1= 59…= 40…知 34111 10||1022 x π--?<-≤? 因而x 1具有4位有效数字。 由π- x 2= 59…= 59…知 223102 1||1021--?≤- 第一章 绪论（12） 1、设0>x ，x 的相对误差为δ，求x ln 的误差。 [解]设0*>x 为x 的近似值，则有相对误差为δε=)(*x r ，绝对误差为**)(x x δε=，从而x ln 的误差为δδεε=='=* ****1)()(ln )(ln x x x x x ， 相对误差为* * ** ln ln ) (ln )(ln x x x x r δ εε= = 。 2、设x 的相对误差为2%，求n x 的相对误差。 [解]设*x 为x 的近似值，则有相对误差为%2)(*=x r ε，绝对误差为**%2)(x x =ε，从而n x 的误差为n n x x n x n x x n x x x ** 1 *** %2%2) ()()()(ln * ?=='=-=εε， 相对误差为%2) () (ln )(ln *** n x x x n r == εε。 3、下列各数都是经过四舍五入得到的近似数，即误差不超过最后一位的半个单位，试指出它们是几位有效数字： 1021.1*1=x ，031.0*2=x ，6.385*3=x ，430.56*4=x ，0.17*5 ?=x 。 [解]1021.1*1 =x 有5位有效数字；0031.0* 2=x 有2位有效数字；6.385*3=x 有4位有效数字；430.56* 4 =x 有5位有效数字；0.17*5?=x 有2位有效数字。 4、利用公式（3.3）求下列各近似值的误差限，其中* 4*3*2*1,,,x x x x 均为第3题所给 的数。 （1）* 4*2*1x x x ++； [解]3 334* 4*2*11** *4*2*1*1005.1102 1 10211021)()()()()(----=?=?+?+?=++=? ??? ????=++∑x x x x x f x x x e n k k k εεεε； （2）* 3*2 *1x x x ； 初中物理力学综合计算 题(有答案) -CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN 1、图31是某建筑工地利用滑轮组和卷扬机提起重物的示意图。当以速度v 1匀速提起质量为m 1的建筑材料时，滑轮组的机械效率为η1，卷扬机拉力的功率为P 1；当以速度v 2匀速提起质量为m 2的建筑材料时，滑轮组的机械效率为η2，卷扬机拉力的功率为P 2。若η2-η1=5%，P 1: P 2=2:3，m 1=90kg ，动滑轮受到的重力G 动=100N 。滑轮与轴的摩擦、细绳受到的重力忽略不计，g=10N/kg 。求： （1）提起质量为m 1的建筑材料时卷扬机对绳的拉力F 1； （2）两次工作过程中，建筑材料上升的速度v 1与v 2之比。 2．在生产玻璃过程中，常用位于天车上的卷扬机（其内部有电动机提供动力）通过滑轮组和真空吸盘提升玻璃，如图22甲所示。当卷扬机通过滑轮组提升质量为60kg 的玻璃并使玻璃以速度v 1匀速上升时，卷扬机对滑轮组绳端的拉力为F 1，天车对卷扬机的支持力为N 1，拉力为F 1的功率为P ，滑轮组的机 械效率为η；当卷扬机通过滑轮组提升质量为80kg 的玻璃并使玻璃以速度v 2匀 速上升时，卷扬机对滑轮组绳端的拉力为 F 2，天车对卷扬机的支持力为N 2。已 知拉力F 1所做功随时间变化的图像如图22乙所示，卷扬机的质量为120 kg ，图31 卷扬机 240 W/J 480 720 960 0 2.0 卷扬机 A B C 玻璃 吸盘 天车 滑轮A 、B 的质量均为4kg ，3v 1=5v 2，η=75%，吸盘和绳的质量及滑轮与轴的摩擦均可忽略不计，g 取10N/kg 。求： （1）P 的大小； （2）v 2的大小； （3）N 1与N 2的比值。 解：（1）由题中W －t 图像解得P =2s 960J = t W =480W …………………（2分） （2）根据η = %75W 480N 600111=?== = v P gv m P P W W 有总 有………………………（1分） 解得：v 1 =0.6m/s 已知：3v 1=5v 2 解得：v 2=0.36m/s ……………………………………………………………（1分） （3）设动滑轮C 所受重力为G 0，卷扬机提升60kg 玻璃时，滑轮组的机械效率为η=75% 所以有 η =%753B 011=++= g m G g m g m W W 总 有，代入数据解得G 0=80 N ………（1 分） 第一次提升60kg 玻璃的过程中，玻璃、动滑轮C 受力分析如答图5（1）所示，动滑轮B 受力分析如答图5（2）所示，卷扬机受力分析如答图5（3）所示。 第一章 题12 给定节点01x =-，11x =，23x =，34x =，试分别对下列函数导出拉格朗日插值余项： （1） （1） 3 ()432f x x x =-+ （2） （2） 4 3 ()2f x x x =- 解 （1）(4) ()0f x =， 由拉格朗日插值余项得(4)0123() ()()()()()()0 4!f f x p x x x x x x x x x ξ-=----=； （2）(4) ()4!f x = 由拉格朗日插值余项得 01234! ()()()()()() 4! f x p x x x x x x x x x -= ----(1)(1)(3)(4)x x x x =+---. 题15 证明：对于()f x 以0x ，1x 为节点的一次插值多项式()p x ，插值误差 012 10()()()max () 8x x x x x f x p x f x ≤≤-''-≤. 证 由拉格朗日插值余项得 01() ()()()()2!f f x p x x x x x ξ''-= --，其中01x x ξ≤≤， 01 0101max ()()()()()()()() 2!2!x x x f x f f x p x x x x x x x x x ξ≤≤''''-=--≤-- 01210()max () 8x x x x x f x ≤≤-''≤. 题22 采用下列方法构造满足条件(0)(0)0p p '==，(1)(1)1p p '==的插值多项式 ()p x ： （1） （1） 用待定系数法； （2） （2） 利用承袭性，先考察插值条件(0)(0)0p p '==，(1)1p =的插值多项式 ()p x . 解 （1）有四个插值条件，故设230123()p x a a x a x a x =+++，2 123()23p x a a x a x '=++， 代入得方程组001231123010231 a a a a a a a a a =? ?+++=?? =? ?++=? 解之，得01230 021 a a a a =??=?? =??=-? 精品学习资料--------------极力推荐中考物理力学综合计算题含答案
计算方法的课后答案
数值分析复习题要答案
数值分析大作业三 四 五 六 七
力学综合计算题
(完整版)计算方法习题第一、二章答案
数值分析习题集及答案
数值计算方法大作业
数值分析1-4习题及答案
2020年中考物理力学综合计算题含答案
计算方法-习题第一、二章答案
数值分析习题
数值分析大作业
计算方法习题第一、二章答案
数值分析课程第五版课后习题答案(李庆扬等)
初中物理力学综合计算题(有答案)
数值分析课后习题答案
初中物理力学综合计算题-(有答案)
1、图 31 是某建筑工地利用滑轮组和卷扬机提起重物的示意图。当
以速度 v1 匀速提起质量为 m1 的建筑材料时，滑轮组的机械效率为 η1，
卷扬机拉力的功率为 P1；当以速度 v2 匀速提起质量为 m2 的建筑材
料时，滑轮组的机械效率为 η2，卷扬机拉力的功率为 P2。若 η2-η1=5%，
P1: P2=2:3，m1=90kg，动滑轮受到的重力 G 动=100N。滑轮与轴的摩
擦、细绳受到的重力忽略不计，g=10N/kg。求：
（1）提起质量为 m1 的建筑材料时卷扬机对绳的拉力 F1；
（2）两次工作过程中，建筑材料上升的速度 v1 与 v2 之比。
卷扬机
建筑材料
图 31
2．在生产玻璃过程中，常用位于天车上的卷扬机（其内部有电动机提供动力）通过滑
轮组和真空吸盘提升玻璃，如图 22 甲所示。当卷扬机通过滑轮组提升质量为 60kg 的玻璃并
使玻璃以速度 v1 匀速上升时，卷扬机对滑轮组绳端的拉力为 F1，天车对卷扬机的支持力为
N1，拉力为 F1 的功率为 P，滑轮组的机械效率为 η；当卷扬机通过滑轮组提升质量为 80kg
的玻璃并使玻璃以速度 v2 匀速上升时，卷扬机对滑轮组绳端的拉力为 F2，天车对卷扬机的
支持力为 N2。已知拉力 F1 所做功随时间变化的图像如图 22 乙所示，卷扬机的质量为 120 kg，
滑轮 A、B 的质量均为 4kg，3v1=5v2，η=75%，吸盘和绳的质量及滑轮与轴的摩擦均可忽略
不计，g 取 10N/kg。求： （1）P 的大小； （2）v2 的大小；
卷扬机 天车
W/J 960
（3）N1 与 N2 的比值。
A
720
解：（1）由题中 W－t 图像解得B P= W 960J =480W ………480…………（2 分） t 2s
240
（2）根据 η= W有 W总
P有 P
m1gCv1 吸盘P
600N v1 75% ………………………（1 分）
480W玻璃
0
0.5 1.0 1.5 2.0
t/s
解得：v1 =0.6m/s
甲 图 22
乙
已知：3v1=5v2
解得：v2=0.36m/s……………………………………………………………（1 分） （3）设动滑轮 C 所受重力为 G0，卷扬机提升 60kg 玻璃时，滑轮组的机械效率为 η=75%
1/5
相关主题
文本预览