清华大学2006数学分析真题参考答案 1．若数列{}n x 满足条件11221n n n n x x x x x x M ----+-++-≤g g g 则称{}n x 为有界变差数列，证：令10y =，11221n n n n n y x x x x x x ---=-+-++-g g g （n=2,3,….） 那么{}n y 单调递增，由条件知{}n y 有界， {}n y ∴收敛 ，从而0,0N ε?>?>，使当n m N >>时，有 n m y y ε-<，此即：11211n n n n m m x x x x x x ε---+--+-++-
7、设y0=28，按递推公式 y n=y n?1? 1 100 783，n=1,2,… 计算y100，若取≈27.982，试问计算y100将有多大误差？ 答：y100=y99?1 100783=y98?2 100 783=?=y0?100 100 783=28?783 若取783≈27.982，则y100≈28?27.982=0.018，只有2位有效数字，y100的最大误差位0.001 10、设f x=ln?(x? x2?1)，它等价于f x=?ln?(x+ x2?1)。分别计算f30，开方和对数取6位有效数字。试问哪一个公式计算结果可靠？为什么？ 答： x2?1≈29.9833 则对于f x=ln x?2?1，f30≈?4.09235 对于f x=?ln x+2?1，f30≈?4.09407 而f30= ln?(30?2?1) ，约为?4.09407，则f x=?ln?(x+ x2?1)计算结果更可靠。这是因为在公式f x=ln?(x? x2?1)中，存在两相近数相减（x? x2?1）的情况，导致算法数值不稳定。 11、求方程x2+62x+1=0的两个根，使它们具有四位有效数字。 答：x12=?62±622?4 2 =?31±312?1 则 x1=?31?312?1≈?31?30.98=?61.98 x2=?31+312?1= 1 31+312?1 ≈? 1 ≈?0.01613
12.（1）、计算101.1?101，要求具有4位有效数字 答：101.1?101= 101.1+101≈0.1 10.05+10.05 ≈0.004975 14、试导出计算积分I n=x n 4x+1dx 1 的一个递推公式，并讨论所得公式是否计算稳定。 答：I n=x n 4x+1dx 1 0= 1 4 4x+1x n?1?1 4 x n?1 4x+1 dx= 1 1 4 x n?1 1 dx?1 4 x n?1 4x+1 dx 1 = 1 4n ? 1 4 I n?1，n=1，2… I0= 1 dx= ln5 1 记εn为I n的误差，则由递推公式可得 εn=?1 εn?1=?=(? 1 )nε0 当n增大时，εn是减小的，故递推公式是计算稳定的。
第一章1-1： 已知：V=72cm3 m=129.1g m s =121.5g G s =2.70 则： 129.1121.5 6.3% 121.5 s s m m w m -- === 3 3 3 3 129.1 *1017.9/ 72 121.5 45 2.7 724527 1.0*27121.5 *1020.6/ 72 s s s V s sat w V s sat sat m g g KN m v m V cm V V V cm m V m g g g KN m V V γρ ρ ρ γρ ==== === =-=-= ++ ===== 3 3 20.61010.6/ 121.5 *1016.9/ 72 sat w s d sat d KN m m g KN m V γγγ γ γγγγ '=-=-= === ' >>> 则 1-2： 已知：G s =2.72 设V s =1cm3 则 3 3 3 3 2.72/ 2.72 2.72 *1016/ 1.7 2.720.7*1 *1020.1/ 1.7 20.11010.1/ 75% 1.0*0.7*75%0.525 0.525 19.3% 2.72 0.525 2.72 1. s s s d d s V w w r w w V r w s w s g cm m g m g g KN m V m V g g KN m V KN m m V S g m w m m m g g V ρ γρ ρ γρ γγγ ρ γρ = = ==== ++ ==== '=-=-= = === === ++ === 当S时， 3 *1019.1/ 7 KN m =
高等土力学 第一章土的物质构成及分类 1蒙脱石和伊利石晶胞结构相同，但蒙脱石具有较大的胀缩性，为什么？ 2用土的结构说明为什么软粘土具有较大流变特性，原生黄土具湿陷性？ 3试述非饱和土中水的迁移特征及控制迁移速率的主要因素？ 4非饱和土中水的运移规律与饱和土中水的渗透规律有什么不同？ 试述非饱和土和饱和土中孔隙水迁移规律的异同点？ 5X射线衍射法是怎样分析粘土矿物成份的？ 6粘土表面电荷来源有哪几方面？利用粘粒表面带电性解释吸着水（结合水）形成机理？ 7非饱和土中土水势以哪种为主？如何测定非饱和土的土水势大小？ 8非饱和土中的土水势主要由哪个几个部分组成？非饱和土中水的迁移速率主要与哪几种因素有关？ 9请用粘性土的结构解释粘性土具有可塑性而砂土没有可塑性的机理。 10试简明解说土水势的各分量？ 11土的结构有哪些基本类型？各有何特征？ 12分散土的主要特征是什么？为什么有些粘性土具有分散性？ 13粘性土主要有哪些性质，它们是如何影响土的力学性质的？ 14为什么粘土颗粒具有可塑性、凝聚性等性质，而砂土颗粒却没有这些性质？ 15非饱和粘性土和饱和的同种粘性土（初始孔隙比相同）在相同的法向应力作用下压缩，达到稳定的压缩量和需要的时间哪个大，哪个小，为什么？ 16粘土的典型结构有哪几种，它们与沉积环境有什么联系，工程性质方面各有何特点？
17粘性土的结构与砂土的结构有什么不同？ 18为什么粘性土在外力作用下具有较大流变特性？ 19粘土矿物颗粒形状为什么大都为片状或针状，试以蒙脱石的晶体结构为例解释之。 第二章土的本构关系及土工有限元分析 1中主应力对土体强度和变形有什么影响？分别在普通三轴仪上和平面应变仪上做 试验，保持σ3为常量，增加σ1－σ3所得应力应变关系曲线有何不同？所得强度指标是否相同？ 2屈服面和硬化规律有何关系？ 3弹塑性柔度矩阵[C]中的元素应有哪三点特征？ 4剑桥弹塑性模型应用了哪些假定？欲得到模型参数应做哪些试验？ 5广义的“硬化”概念是什么？什么叫硬化参数？ 6什么是流动规则？什么叫塑性势？流动规则有哪两种假定？ 7弹塑性模型中，为什么要假定某种型式的流动法则，它在确定塑性应变中有何作用？ 8根据相适应的流动规则，屈服面和塑性应变增量的方向有何特征？ 9试解释为什么球应力影响塑性剪应变？ 10什么叫土的变形“交叉效应”？“交叉效应”对土的刚度矩阵[D]或柔度矩阵[C]有何影响？ 11什么叫应力路径？什么叫应力历史？试结合图示说明它们对土的变形的影响？ 12什么叫土的“各向异性”？考虑“各向异性”对土的刚度矩阵[D]或柔度矩阵[C]有何影响？ 13哪些因素影响土的变形？或土体变形有哪些特征？ 14什么叫剪缩？什么叫剪胀？什么样的土表现为剪胀，怎样的土表现为剪缩？邓肯双曲线模型能否反映剪胀，剪缩？为什么？修正剑桥模型能否反映？
清华大学大学物理习题库：量子物理 一、选择题 1．4185：已知一单色光照射在钠表面上，测得光电子的最大动能是1.2 eV ，而钠的红限波长是5400 ?，那么入射光的波长是 (A) 5350 ? (B) 5000 ? (C) 4350 ? (D) 3550 ? ［ ］ 2．4244：在均匀磁场B 内放置一极薄的金属片，其红限波长为??。今用单色光照射，发现有电子放出，有些放出的电子(质量为m ，电荷的绝对值为e )在垂直于磁场的平面内作半径为R 的圆周运动，那末此照射光光子的能量是： (A) 0λhc (B) 0λhc m eRB 2)(2+ (C) 0λhc m eRB + (D) 0λhc eRB 2+ ［ ］ 3．4383：用频率为??的单色光照射某种金属时，逸出光电子的最大动能为E K ；若改用 频率为2??的单色光照射此种金属时，则逸出光电子的最大动能为： (A) 2 E K (B) 2h ??- E K (C) h ??- E K (D) h ??+ E K ［ ］ 4．4737： 在康普顿效应实验中，若散射光波长是入射光波长的1.2倍，则散射光光子能量?与反冲电子动能E K 之比??/ E K 为 (A) 2 (B) 3 (C) 4 (D) 5 ［ ］ 5．4190：要使处于基态的氢原子受激发后能发射赖曼系(由激发态跃迁到基态发射的各谱线组成的谱线系)的最长波长的谱线，至少应向基态氢原子提供的能量是 (A) 1.5 eV (B) 3.4 eV (C) 10.2 eV (D) 13.6 eV ［ ］ 6．4197：由氢原子理论知，当大量氢原子处于n =3的激发态时，原子跃迁将发出： (A) 一种波长的光 (B) 两种波长的光 (C) 三种波长的光 (D) 连续光谱 ［ ］ 7．4748：已知氢原子从基态激发到某一定态所需能量为10.19 eV ，当氢原子从能量为－0.85 eV 的状态跃迁到上述定态时，所发射的光子的能量为 (A) 2.56 eV (B) 3.41 eV (C) 4.25 eV (D) 9.95 eV ［ ］ 8．4750：在气体放电管中，用能量为12.1 eV 的电子去轰击处于基态的氢原子，此时氢原子所能发射的光子的能量只能是 (A) 12.1 eV (B) 10.2 eV (C) 12.1 eV ，10.2 eV 和 1.9 eV (D) 12.1 eV ，10.2 eV 和 3.4 eV ［ ］ 9．4241： 若?粒子(电荷为2e )在磁感应强度为B 均匀磁场中沿半径为R 的圆形轨道运动，则?粒子的德布罗意波长是 (A) )2/(eRB h (B) )/(eRB h (C) )2/(1eRBh (D) )/(1eRBh ［ ］ 10．4770：如果两种不同质量的粒子，其德布罗意波长相同，则这两种粒子的 (A) 动量相同 (B) 能量相同 (C) 速度相同 (D) 动能相同 ［ ］
清华大学《大学物理》习题库试题及答案--08-电学习 题答案 本页仅作为文档页封面，使用时可以删除 This document is for reference only-rar21year.March
一、选择题 1．1003：下列几个说法中哪一个是正确的？ (A) 电场中某点场强的方向，就是将点电荷放在该点所受电场力的方向 (B) 在以点电荷为中心的球面上，由该点电荷所产生的场强处处相同 (C) 场强可由定出，其中q 为试验电荷，q 可正、可负，为试验电荷所受的电场力 (D) 以上说法都不正确 ［ ］ 2．1405：设有一“无限大”均匀带正电荷的平面。取x 轴垂直带电平面， 坐标原点在带电平面上，则其周围空间各点的电场强度随距离平面的位置坐 标x 变化的关系曲线为(规定场强方向沿x 轴正向为正、反之为负)： ［ ］ 3．1551：关于电场强度定义式，下列说法中哪个是正确的？ (A) 场强的大小与试探电荷q 0的大小成反比 (B) 对场中某点，试探电荷受力与q 0的比值不因q 0而变 (C) 试探电荷受力的方向就是场强的方向 (D) 若场中某点不放试探电荷q 0，则＝0，从而＝0 ［ ］ 4．1558：下面列出的真空中静电场的场强公式，其中哪个是正确的？ ［ ］ q F E / =F E /q F E =E F F E F E ( x
(A)点电荷q 的电场：(r 为点电荷到场点的距离) (B)“无限长”均匀带电直线(电荷线密度)的电场：(为带电直线到场点的垂直于直线的矢量) (C)“无限大”均匀带电平面(电荷面密度)的电场： (D) 半径为R 的均匀带电球面(电荷面密度)外的电场：(为球心到场点的矢量) 5．1035：有一边长为 a 的正方形平面，在其中垂线上距中心O 点a /2处，有一电荷为q 的正点电荷，如图所示，则通过该平面的电场强度通量为 (A) (B) (C) (D) ［ ］ 6．1056：点电荷 Q 被曲面S 所包围，从无穷远处引入另一点电荷q 至曲面外一点，如图所示，则引入前后： (A) 曲面S 的电场强度通量不变，曲面上各点场强不变 (B) 曲面S 的电场强度通量变化，曲面上各点场强不变 (C) 曲面S 的电场强度通量变化，曲面上各点场强变化 (D) 曲面S 的电场强度通量不变，曲面上各点场强变化 ［ ］ 7．1255：图示为一具有球对称性分布的静电场的E ～r 关系曲线。请指出该静电场是由下列哪种带电体产生的 (A) 半径为R 的均匀带电球面 (B) 半径为R 的均匀带电球体 (C) 半径为R 的、电荷体密度为的非均匀带电球体 2 04r q E επ= λr r E 302ελπ= r σ02εσ = E σr r R E 3 02εσ=r 0 3εq 4επq 0 3επq 0 6εq Ar =ρ q 1035图 q
9.令*()(21),[0,1]n n T x T x x =-∈，试证*{()}n T x 是在[0,1] 上带权()x ρ=的正交多项式，并求****0123(),(),(),()T x T x T x T x . 证明： 1 1 * *0 1 1 * *011**0 ()()()(21)(21)211()()()()()2()()()()()()()()n m n m n m n m n m n n m n m x T x T x dx x T x dx t x x T x T x dx t T t dt t T t dt T x x T x T x dx t T t ρρρ---=--=-== = ???? ?令，则 由切比雪夫多项式1 01=02 m n dt m n m n ππ ≠??? =≠??==??? 所以*{()}n T x 是在[0,1] 上带权()x ρ= *00*11* 22 2 2*33233()(21)1()(21)21 ()(21)2(21)188()(21)4(21)3(21)3248181 T x T x T x T x x T x T x x x x T x T x x x x x x =-==-=-=-=--=-=-=---=-+- 14.已知实验数据如下： 用最小二乘法求形如2y a bx =+的经验公式，并求均方误差 解： 法方程为
22222(1,)(1,1)(1,)(,)(,1)(,)a y x b x y x x x ?????? =???? ?????? ?? 即 5 5327271.453277277699369321.5a b ??????=???????????? 解得 0.972579 0.050035a b =?? =? 拟合公式为20.9725790.050035y x =+ 均方误差 2 4 2 2 0[]0.015023i i i y a bx σ==--=∑ 21.给出()ln f x x =的函数表如下： 用拉格朗日插值求ln 0.54的近似值并估计误差（计算取1n =及2n =） 解：1n =时，取010.5,0.6x x == 由拉格朗日插值定理有 1 100.60.5 0.693147 0.510826 0.50.(60.60.51.82321)0 1.()6047()52 j j j x x x L x f x l x ==------=-=∑ 所以1ln0.54(0.54)0.620219L ≈=- 误差为ln 0.54(0.620219)= 0.004032ε=-- 2n =时，取0120.4,0.5,0.6x x x === 由拉格朗日插值定理有
2．定义映射22:B R R →，()B x y =，满足y Ax =，其中 0.80.40.10.4A ??=????，2,x y R ∈ 则对任意的2 ,u v R ∈ 1111119 ||()()||||||||()||||||||||||||10B u B v Au Av A u v A u v u v -=-=-≤-=- 故映射B 对一范数是压缩的 由范数定义 ||||1 ||||max |||| 1.2 x A Ax ∞∞∞===，知必然存在0 x ， 0||||1 x ∞= 使得0|||||||| 1.2 Ax A ∞∞== 设012(,)T x x x = 取 12(,0),(0,)T T u x v x ==-，则 u v x -=，有 00||()()||||||||()|||||||||| 1.21||||||||B u B v Au Av A u v Ax A x u v ∞∞∞∞∞∞∞ -=-=-===>==- 故有||()()||B u B v ∞->||||u v ∞ -，从而映射B 对无穷范数不是压缩的 4. 证明：对任意的,[,]x y a b ∈ 由拉格朗日中值定理，有 ()()'()()() 1e G x G y G x y x y e ξ ξξ-=-=-+ 其中0111b b e e e e ξξ<≤<++ 所以 |()()||()||| 11b b e e G x G y x y x y e e ξξ-=-≤-++ 故G 为[,]a b 上的压缩映射 而 ()ln(1)ln x x G x e e x =+>= 即()G x x =无根
热学部分 一、选择题 1．4251：一定量的理想气体贮于某一容器中，温度为T ，气体分子的质量为m 。根据理想气体的分子模型和统计假设，分子速度在x 方向的分量平方的平均值 (A) (B) (C) (D) ［ ］ 2．4252：一定量的理想气体贮于某一容器中，温度为T ，气体分子的质量为m 。根据理想气体分子模型和统计假设，分子速度在x 方向的分量的平均值 (A) (B) (C) (D) 0 ［ ］ 3．4014：温度、压强相同的氦气和氧气，它们分子的平均动能和平均平动动能 有如下关系：(A) 和都相等 (B) 相等，而不相等 (C) 相等，而不相等 (D) 和都不相等 ［ ］ 4．4022：在标准状态下，若氧气(视为刚性双原子分子的理想气体)和氦气的体积比V 1 / V 2=1 / 2 ，则其内能之比E 1 / E 2为： (A) 3 / 10 (B) 1 / 2 (C) 5 / 6 (D) 5 / 3 ［ ］ 5．4023：水蒸气分解成同温度的氢气和氧气，内能增加了百分之几(不计振动自由度和化学能)？ (A) 66.7％ (B) 50％ (C) 25％ (D) 0 ［ ］ 6．4058：两瓶不同种类的理想气体，它们的温度和压强都相同，但体积不同，则单位体积内的气体分子数n ，单位体积内的气体分子的总平动动能(EK /V )，单位体积内的气体质量，分别有如下关系：(A) n 不同，(EK /V )不同，不同 (B) n 不同，(EK /V )不同，相同 (C) n 相同，(EK /V )相同，不同 (D) n 相同，(EK /V )相同，相同 ［ ］ 7．4013：一瓶氦气和一瓶氮气密度相同，分子平均平动动能相同，而且它们都处于平衡状态，则它们 (A) 温度相同、压强相同 (B) 温度、压强都不相同 (C) 温度相同，但氦气的压强大于氮气的压强 (D) 温度相同，但氦气的压强小于氮气的压强 ［ ］ 8．4012：关于温度的意义，有下列几种说法：(1) 气体的温度是分子平均平动动能的量度；(2) 气体的温度是大量气体分子热运动的集体表现，具有统计意义；(3) 温度的高低反映物质内部分子运动剧烈程度的不同；(4) 从微观上看，气体的温度表示每个气体分子的冷热程度。这些说法中正确的是 (A) (1)、(2)、(4)；(B) (1)、(2)、(3)；(C) (2)、(3)、(4)；(D) (1)、(3) 、(4)； ［ ］ 9．4039：设声波通过理想气体的速率正比于气体分子的热运动平均速率，则声波通过具有相同 温度的氧气和氢气的速率之比为 (A) 1 (B) 1/2 (C) 1/3 (D) 1/4 ［ ］ 10．4041：设图示的两条曲线分别表示在相同温度下氧气和氢气分子的速率分布曲线；令 和分别表示氧气和氢气的最概然速率，则： (A) 图中ａ表示氧气分子的速率分布曲线； /=4 (B) 图中ａ表示氧气分子的速率分布曲线； /＝1/4 (C) 图中ｂ表示氧气分子的速率分布曲线； /＝1/4 (D) 图中ｂ表示氧气分子的速率分布曲线； /＝ 4 ［ ］ m kT x 32= v m kT x 3312 =v m kT x /32=v m kT x /2 =v m kT π8= x v m kT π831=x v m kT π38= x v =x v εw εw εw w εεw ρρρρρ2 2H O /v v ()2 O p v ()2 H p v ()2 O p v ()2 H p v ()2O p v ()2H p v ()2 O p v ()2 H p v ()2 O p v ()2 H p v
高等数值分析实验一 工物研13 成彬彬2004310559 一．用CG，Lanczos和MINRES方法求解大型稀疏对称正定矩阵Ax=b 作实验中，A是利用A= sprandsym(S,[],rc,3)随机生成的一个对称正定阵，S是1043阶的一个稀疏阵 A= sprandsym(S,[],0.01,3); 检验所生成的矩阵A的特征如下： rank(A-A')=0 %即A=A’，A是对称的； rank(A)=1043 %A满秩 cond(A)= 28.5908 %A是一个“好”阵 1．CG方法 利用CG方法解上面的线性方程组 [x,flag,relres,iter,resvec] = pcg(A,b,1e-6,1043); 结果如下： Iter=35，表示在35步时已经收敛到接近真实x relres= norm(b-A*x)/norm(b)= 5.8907e-007为最终相对残差 绘出A的特征值分布图和收敛曲线： S=svd(A); %绘制特征值分布 subplot(211) plot(S); title('Distribution of A''s singular values');; xlabel('n') ylabel('singular values') subplot(212); %绘制收敛曲线 semilogy(0:iter,resvec/norm(b),'-o'); title('Convergence curve'); xlabel('iteration number'); ylabel('relative residual'); 得到如下图象：
为了观察CG方法的收敛速度和A的特征值分布的关系，需要改变A的特征值： (1).研究A的最大最小特征值的变化对收敛速度的影响 在A的构造过程中，通过改变A= sprandsym(S,[],rc,3)中的参数rc(1/rc为A的条件数)，可以达到改变A的特征值分布的目的： 通过改变rc=0.1，0.0001得到如下两幅图 以上三种情况下，由收敛定理2.2.2计算得到的至多叠代次数分别为：48，14和486，由于上实验结果可以看出实际叠代次数都比上限值要小较多。 由以上三图比较可以看出，A的条件数越大，即A的最大最小特征值的差别越大，叠代所需要的步骤就越多，收敛越慢。 (2)研究A的中间特征值的分布对于收敛特性的影响： 为了研究A的中间特征值的分布对收敛速度的影响，进行了如下实验： 固定A的条件数，即给定A的最大最小特征值，改变中间特征值得分布，再来生成A，具体的实现方法是，先将原来的生成A进行特征值分解： [U,S]=svd(A);
清华大学版土力学课后 答案修订版 IBMT standardization office【IBMT5AB-IBMT08-IBMT2C-ZZT18】
第一章1-1： 已知：V=72cm3 m=129.1g m s =121.5g G s =2.70 则： 129.1121.5 6.3% 121.5 s s m m w m -- === 1-2： 已知：G s =2.72 设V s =1cm3 则 3 3 3 3 2.72/ 2.72 2.72 *1016/ 1.7 2.720.7*1 *1020.1/ 1.7 20.11010.1/ 75% 1.0*0.7*75%0.525 0.525 19.3% 2.72 0.525 2.72 1. s s s d d s V w w r w w V r w s w s g cm m g m g g KN m V m V g g KN m V KN m m V S g m w m m m g g V ρ γρ ρ γρ γγγ ρ γρ = = ==== ++ ==== '=-=-= = === === ++ === 当S时， 3 *1019.1/ 7 KN m = 1-3：1-4：甲：乙：
则（1）、（4）正确 1-5： 1s w d G e ρρ= + 则 所以该料场的土料不适合筑坝，建议翻晒，使其含水率降低。 1-6： 式中D r =0.7 3max 1.96/d g cm ρ= 3min 1.46/d g cm ρ= 则可得：31.78/d g cm ρ= 1-7： 设 S=1， 则s V Sh h == 则压缩后： 2.7s s s m V G h == 2.7*28%w s m m w h == 则 2.7*28%w w w m V h ρ= = 2.7*28% 1.95s w V V h h +=+= 则 1.11h cm = 1-8： 甲：4525 1.334025 p L L p w w I w w --= = =-- 流塑状态 乙：2025 0.334025 p L L p w w I w w --= = =--- 坚硬（半固态）
清华大学杨顶辉数值分析第6次作业
9.令*()(21),[0,1]n n T x T x x =-∈，试证*{()}n T x 是在[0,1]上带权 2 ()x x x ρ= -****0123(),(),(),()T x T x T x T x . 证明： 1 1 **2 1 1 * *20 12 2 1**20 ()()()(21)(21)211()()()()()211()22 ()()1()1()()()()()1n m n m n m n m n m n n m n m x T x T x dx x T x dx x x t x x T x T x dx t T t dt t t t T t dt t T x x x T x T x dx t T t t ρρρ---=---=-=++-= --= -???? ?令，则 由切比雪夫多项式1 01=02 m n dt m n m n ππ ≠??? =≠??==??? 所以*{()}n T x 是在[0,1]上带权2 ()x x x ρ= - *00*11* 2 2 2 2*33233()(21)1()(21)21 ()(21)2(21)188()(21)4(21)3(21)3248181 T x T x T x T x x T x T x x x x T x T x x x x x x =-==-=-=-=--=-=-=---=-+- 14.已知实验数据如下： i x 19 25 31 38 44 i y 19.0 32.3 49.0 73.3 97.8 用最小二乘法求形如2y a bx =+的经验公式，并求均方误差 解： 法方程为
一、选择题 1．0018：某质点作直线运动的运动学方程为x ＝3t -5t 3 + 6 (SI)，则该质点作 (A) 匀加速直线运动，加速度沿x 轴正方向 (B) 匀加速直线运动，加速度沿x 轴负方向 (C) 变加速直线运动，加速度沿x 轴正方向 (D) 变加速直线运动，加速度沿x 轴负方向 ［ ］ 2．5003：一质点在平面上运动，已知质点位置矢量的表示式为 （其中a 、b 为常量），则该质点作 (A) 匀速直线运动 (B) 变速直线运动 (C) 抛物线运动 (D)一般曲线运动 ［ ］ 3．0015：一运动质点在某瞬时位于矢径的端点处， 其速度大小为 (A) (B) (C) (D) 4．0508：质点沿半径为R 的圆周作匀速率运动，每T 秒转一圈。在2T 时间间隔中，其平均速度大小与平均速率大小分别为 (A) 2p R /T ， 2p R/T (B) 0 ， 2πR /T (C) 0 ， 0 (D) 2πR /T ， 0. ［ ］ 5．0518：以下五种运动形式中，保持不变的运动是 (A) 单摆的运动 (B) 匀速率圆周运动 (C) 行星的椭圆轨道运动 (D) 抛体运动 (E) 圆锥摆运动 ［ ］ 6．0519：对于沿曲线运动的物体，以下几种说法中哪一种是正确的： (A) 切向加速度必不为零 (B) 法向加速度必不为零（拐点处除外） (C) 由于速度沿切线方向，法向分速度必为零，因此法向加速度必为零 (D) 若物体作匀速率运动，其总加速度必为零 (E) 若物体的加速度为恒矢量，它一定作匀变速率运动 ［ ］ 7．0602：质点作曲线运动，表示位置矢量，表示速度，表示加速度，S 表示路 程，a 表示切向加速度，下列表达式中， (1) ， (2) ， (3) ， (4) (A) 只有(1)、(4)是对的 (B) 只有(2)、(4)是对的 (C) 只有(2)是对的 (D) 只有(3)是对的 ［ ］ 8．0604：某物体的运动规律为，式中的k 为大于零的常量。当时，初速为v 0，则速度与时间t 的函数关系是 (A) ， (B) ， (C) ， (D) ［ ］ 9．0014：在相对地面静止的坐标系内，A 、B 二船都以2 m/s 速率匀速行驶，A 船沿x 轴正向，B 船沿y 轴正向。今在A 船上设置与静止坐标系方向相同的坐标系(x 、y 方向单位 矢用、表示)，那么在A 船上的坐标系中，B 船的速度（以m/s 为单位）为 j bt i at r 2 2+=()y x r , t r d d t r d d t r d d 2 2d d d d ?? ? ??+??? ??t y t x a a r v a a t = d /d v v =t r d /d v =t S d /d t a t =d /d v t k t 2 d /d v v -=0=t v 0 2 2 1v v += kt 2 2 1v v +- =kt 02 12 1v kt v += 2 12 1v kt v + - =i j
线性代数方程组的数值解法 实验1.主元的选取与算法的稳定性 问题提出：Gauss 消去法是我们在线性代数中已经熟悉的。但由于计算机的数值运算是在一个有限的浮点数集合上进行的，如何才能确保Gauss 消去法作为数值算法的稳定性呢？Gauss 消去法从理论算法到数值算法，其关键是主元的选择。主元的选择从数学理论上看起来平凡，它却是数值分析中十分典型的问题。 实验内容：考虑线性方程组 n n n R b R A b Ax ∈∈=?,, 编制一个能自动选取主元，又能手动选取主元的求解线性方程组的Gauss 消去过程。 实验要求： （1）取矩阵?? ???? ? ?????????=???????????? ? ?? ?=141515 7,68 168 16816 b A ，则方程有解T x )1,,1,1(* =。取n=10 计算矩阵的条件数。让程序自动选取主元，结果如何？ （2）现选择程序中手动选取主元的功能。每步消去过程总选取按模最小或按模尽可能小的元素作为主元，观察并记录计算结果。若每步消去过程总选取按模最大的元素作为主元，结果又如何？分析实验的结果。 （3）取矩阵阶数n=20或者更大，重复上述实验过程，观察记录并分析不同的问题及消去过程中选择不同的主元时计算结果的差异，说明主元素的选取在消去过程中的作用。 （4）选取其他你感兴趣的问题或者随机生成矩阵，计算其条件数。重复上述实验，观察记录并分析实验结果。 1.1程序清单 n=input('矩阵A 的阶数:n='); A=6*diag(ones(1,n))+diag(ones(1,n-1),1)+8*diag(ones(1,n-1),-1); b=A*ones(n,1); p=input('计算条件数使用p-范数,p='); cond_A=cond(A,p) [m,n]=size(A); Ab=[A b]; r=input('选主元方式(0:自动;1:手动)，r=');
书上18页表1——5的公式的推导过程（对后面的章节的知识补充）： 1、土的密度：即表示单位土体的质量，单位：g/cm 3或kg/m 3 ρ= m V = m s +m w V a +V w +V s 2、土的容重：即表示单位土体的重量，单位：N/cm 3或KN/m 3 γ=ρg 3、土的比重：土粒质量与同体积纯水在4摄氏度下的质量之比（无量纲） G s = m s V s .ρw = ρs ρw 4、土壤含水量：土中水的质量与土的质量的比值 ω % = m w m s = m ?m s m s 5、干密度：单位体积土的质量，单位：g/cm 3或kg/m 3 基本公式：ρd =m s V (1) 又∵G s =m s V s .ρw →m s =V s .ρw .G s ；∴ρd =m s V = V s .ρw .G s V 又∵ V s V = V s V s +V v = 11+e (或者1+e=1+ V v V s = V V s , ∴ V s V = 1 1+e ) ∴ρd =ρw ×G s (1+e) (2) ρd = m s V = m s V .m m =m V .m s m = ρ1+ω (1+ω=1+ m w m s = m m s ,∴ m s m = 11+ω ) (3) 6、孔隙比：孔隙体积与固体颗粒实体体积之比 基本公式：e=V v V s (1) 又∵e= V v V s = V ?V s V s =V V s ?1， 又∵V =m s ρd ，V s = m s ρs ∴e= V V s ?1= m s ρd / m s ρs ?1=ρs ρd ?1 (2) ∵ρd = ρ1+ω ∴e= ρs ρd ?1= ρs 1+ω ρ?1 (3) 又因为n= V v V , e= V v V s =V v V ?V v , 1e = V ?V v V v = V V v ?1=1 n ?1 ∴e= n 1?n (4) 7、孔隙率：孔隙体积与土样总体积之比 基本公式：n= V v V (1)
第一章 绪论 1．设0x >,x 的相对误差为δ，求ln x 的误差。 解：近似值*x 的相对误差为* **** r e x x e x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈ 进而有(ln *)x εδ≈ 2．设x 的相对误差为2%，求n x 的相对误差。 解：设()n f x x =，则函数的条件数为'()||() p xf x C f x = 又1'()n f x nx -=, 1 ||n p x nx C n n -?∴== 又((*))(*)r p r x n C x εε≈? 且(*)r e x 为2 3．下列各数都是经过四舍五入得到的近似数，即误差限不超过最后一位的半个单位，试指 出它们是几位有效数字：*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =,*57 1.0.x =? 解：*1 1.1021x =是五位有效数字； *20.031x =是二位有效数字； *3385.6x =是四位有效数字； *456.430x =是五位有效数字； *57 1.0.x =?是二位有效数字。 4．利用公式(2.3)求下列各近似值的误差限：(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x . 其中****1234,,,x x x x 均为第3题所给的数。 解： 5计算球体积要使相对误差限为1，问度量半径R 时允许的相对误差限是多少？ 解：球体体积为343 V R π= 则何种函数的条件数为 又(*)1r V ε=
故度量半径R 时允许的相对误差限为1(*)10.333r R ε= ?≈ 6．设028Y =，按递推公式1n n Y Y -= （n=1,2,…） 计算到100Y 27.982≈（5位有效数字），试问计算100Y 将有多大误差？ 解：1n n Y Y -= …… 依次代入后，有1000100Y Y =- 即1000Y Y =， 27.982≈, 100027.982Y Y ∴=- 100Y ∴的误差限为31102 -?。 7．求方程25610x x -+=的两个根，使它至少具有427.982=）。 解：2 5610x x -+=， 故方程的根应为1,228x = 故 1282827.98255.982x =≈+= 1x ∴具有5位有效数字 2x 具有5位有效数字 8．当N 充分大时，怎样求 1211N N dx x ++?？ 解 1 21arctan(1)arctan 1N N dx N N x +=+-+? 设arctan(1),arctan N N αβ=+=。 则tan 1,tan .N N αβ=+= 9．正方形的边长大约为了100cm ，应怎样测量才能使其面积误差不超过2 1cm ？ 解：正方形的面积函数为2()A x x = (*)2*(*)A A x εε∴=. 当*100x =时，若(*)1A ε≤,
11. 解：计算中保留5位有效数字，第一步，选取作为主元，则 消去，得 第二步，选择 作为主元，则 消去，得 回代计算得到方程的解为 12. （1）证明： 先证明 的对称性，易得 再证的正定性，只要证明的顺序主子式 21311 0.250253.9960.0020.0005005 3.996l l = =-==-(2)(2) 3.9960 5.562547.4178(|)00.61077 1.00100.474700 2.0028 2.00200.40371A b ?? ? =--- ? ??? 320.61077 0.30496 2.0028 l -= =-(3)(3) 3.9960 5.562547.4178(|)0 2.0028 2.00200.40371000.390470.35158A b ?? ? = ? ?--?? 1231.92729,0.69847,0.90040 x x x ==-=(2)(2)11111111 11111111(2)(2)22 ,,2,3,...,,2,3,...,A ===A ij ij i j ji ji j i i j i j ij ji i j j i T ij ji a a l a a a l a i n j n a a a a a a l a l a a a a a A =-=-=====由于对称正定，则 ，则 ，即(2) 22()2()0i ii n nn a A a a ?? ? ? ?→ ? ? ?? ? O O
易得将作Gauss 消去，最终得到 由于这种变换不改变矩阵的行列式，则 由于A 对称正定，则，因此，即的顺序主子式大于零 综上， 对称正定。 （2）证明： 只需证明 由于 则 由于A 严格对角占优，则 则严格对角占优。 13. 解：显然A 对称，，则A 为对称正定矩阵，用平方 根法求得下三角矩阵L 为 由得 ，再由 得 (2)(3)() 22313...,2,3,...,i i ii a a a i n -=?=(2)(2) 2 ||||0,2,3,...,n ii ij j j i a a i n =≠->=∑(2)11ij ij i j a a l a =-(2)(2)11112211112222 11|||||||||||||||||||||| n n ii ij ii i i ij i j j j j i j i n n n n i ii ij i j ii ij j j j j j j i j i a a a l a a l a a a a l a a a a a ==≠≠====≠≠-=---≥ --=--∑∑∑∑∑∑12 1111 2 11(2)(2)21||||||,||||,1 || ||||||||0 n j n n j ii ij j j j j i n n ii ij ii ij j j j i j i a a a a a a a a a a ===≠==≠≠>><->->∑∑∑∑∑则 4L=12233?? ? ? ?-??