材料物理习题集
第一章固体中电子能量结构和状态（量子力学基础）
1.一电子通过5400V电位差的电场，（1）计算它的德布罗意波长；（2）计算它的波数；
（3）计算它对Ni晶体（111）面（面间距d=2.04×10-10m）的布拉格衍射角。（P5）
1
2
34
1
31192
11
11
o'
(2)
6.610
=
(29.1105400 1.610)
=1.6710
2
K 3.7610
sin
sin218
2
h h
p
mE
m
d
d
λ
π
λ
θλ
λ
θθ
-
--
-
=
?
?????
?
=?
=
=?=
解：（1）=
（2）波数=
（3）2
2.有两种原子，基态电子壳层是这样填充的
;
;
s s s
s s s s
22623
22626102610
（1）1、22p、33p
（2）1、22p、33p3d、44p4d
，请分别写出n=3的所有电子的四个量子数的可能组态。（非书上内容）
3. 如电子占据某一能级的几率是1/4，另一能级被占据的几率为3/4，分别计算两个能级
的能量比费米能级高出多少k T ？（P15）
1
()exp[]11ln[1]()()1/4ln 3()3/4ln 3F F F F f E E E kT
E E kT f E f E E E kT
f E E E kT =
-+?-=-=-=?=-=-?解：由将代入得将代入得
4. 已知Cu 的密度为8.5×103kg/m 3，计算其E 0
F 。（P16） 22
03
23426233
3118(3/8)2(6.6310)8.510 =(3 6.0210/8)291063.5
=1.0910 6.83F
h E n m J eV ππ---=????????=解：
由
5. 计算Na 在0K 时自由电子的平均动能。（Na 的摩尔质量M=22.99，
.0ρ?33=11310kg/m ）（P16）
22
03
23426
233
311900(3/8)2(6.6310) 1.01310 =(3 6.0210/8)291022.99
=5.2110 3.253 1.085F
F h E n m
J eV
E E eV ππ---=????????===解：由由 6. 若自由电子矢量K 满足以为晶格周期性边界条件x x L ψψ+()=()和定态薛定谔方程。
试证明下式成立：e iKL =1
()()()()1
iKx
iK x L iKx iKL iKx
iKL x Ae x L Ae Ae e x Ae e ???+∴=∴+====?=Q g 解：由于满足薛定谔定态方程
又满足周期性边界条件
7. d h r K K cos r /2θ?=*hkl *hkl 已知晶面间距为，晶面指数为（ k l ）的平行晶面
的倒易矢量为，一电子波与该晶面系成角入射，试证明产生布拉格反射的临界波矢量的轨迹满足方程。
8. 试用布拉格反射定律说明晶体电子能谱中禁带产生的原因。
（P20）
9. 试用晶体能带理论说明元素的导体、半导体、绝缘体的导电性质。
答： （画出典型的能带结构图，然后分别说明）
10. 过渡族金属物理性质的特殊性与电子能带结构有何联系？（P28）
答：过渡族金属的d 带不满，且能级低而密，可容纳较多的电子，夺取较高的s 带中的电子，降低费米能级。
补充习题
1. 为什么镜子颠倒了左右而没有颠倒上下？
2. 只考虑牛顿力学，试计算在不损害人体安全的情况下，加速到光速需要多少时间？
3. 已知下列条件，试计算空间两个电子的电斥力和万有引力的比值
1122
3119922
G 6.6710 9.1110 1.6010 8.9910e e N m kg m kg
q C N m C ε------=?=?=?=?g g g g 万有引力常数电子质量电子电量介电常数
12212271
1228
1243
5.510/ 2.3102.4110m m F G
r q q F k r Gm m F F kq q ---=??==?=?引斥引斥解：=
4. 画出原子间引力、斥力、能量随原子间距变化的关系图。
5. 面心立方晶体，晶格常数a=0.5nm ，求其原子体密度。
223-7344 3.210/0.510cm cm =??解：由于每个面心立方晶胞含个原子，所以原子体密度为：
原子原子（）
6. 简单立方的原子体密度是223310cm -?。假定原子是钢球并与最近的相邻原子相
切。确定晶格常数和原子半径。
22-331310cm a 0.3221r 0.1612
nm a
a nm =??===解：每个简单立方晶胞含有一个原子：
第二章 材料的电性能
1. 铂线300K 时电阻率为1×10-7Ω·m ，假设铂线成分为理想纯。试求1000K 时的电阻率。
（P38）
T 0772*******(1)
1+T 1+T 5110 2.27101+1+ 2.2
T m T T ρραρααρρραα--=+=?==??=?Ωg g 解：
2. 镍铬丝电阻率（300K ）为1×10-6Ω·m ，加热到4000K 时电阻率增加5%，假定在此温
度区间内马西森定则成立。试计算由于晶格缺陷和杂质引起的电阻率。（P38）
3. 为什么金属的电阻温度系数为正的？ （P37-38）
答：当电子波通过一个理想晶体点阵时（0K ），它将不受散射；只有在晶体点阵完整性遭到破坏的地方，电子波才受到散射（不相干散射），这就是金属产生电阻的根本原因，因此随着温度升高，电阻增大，所以金属的电阻温度系数为正。
4. 试说明接触电阻产生的原因和减小这个电阻的措施。（P86）
接触电阻产生的原因有两个：一是因为接触面不平，真正接触面比看到的要小，电流通过小的截面必然产生电阻，称为会聚电阻。二是无论金属表面怎样干净，总是有异物形成的膜，可能是周围气体、水分的吸附层。因此，一般情况下，接触金属时首先接触到的是异物薄膜，这种由于膜的存在而引起的电阻称为过渡电阻。
5. 镍铬薄膜电阻沉积在玻璃基片上其形状为矩形1mm ×5mm ，镍铬薄膜电阻率为1×
10-6Ω·m ，两电极间的电阻为1K Ω，计算表面电阻和估计膜厚。
工程硕士《应用概率统计》复习题 考试要求：开一页；题目类型：简答题和大题；考试时间：100分钟。 1. 已知 0.5,)( 0.4,)( 0.3,)(===B A P B P A P 求)(B A P ?。 解：因为 0.7,0.3-1)(-1(A)===A P P 又因为, ,-- A B A B A A B A AB ?== 所以 0.2,0.5-7.0)( -(A))(A ===B A P P B P 故 0.9.0.2-0.40.7P(AB)-P(B)(A))(A =+=+=?P B P 2.设随机变量)1(,9 5 )1(),,4(~),,2(~≥=≥Y P X P p b Y p b X 求并且。 解： . 8165 31-1-10)(Y -11)(Y ),3 1,4(~,31,94-1-1-10)(X -1)1(,9 5)1(),,2(~422 ====≥=====≥=≥）（故从而解得）所以（） （而且P P b Y p p p P X P X P p b X 3．随机变量X 与Y 相互独立，下表中给出了X 与Y 的联合分布的部分数值，请将表中其
4．设随机变量Y 服从参数2 1=λ的指数分布，求关于x 的方程0322 =-++Y Yx x 没有实根的概率。 解：因为当时没有实根时，即0128Y -Y 03)-4(2Y -Y 2 2 <+<=?，故所求的概率为}6Y P{20}128Y -P{Y 2 <<=<+，而Y 的概率密度 ?? ???≤>=0,00 ,21f(y)21-y y e y ，从而36221 -621-1dy 21f(y)dy 6}Y {2e e e P y ===<
数值分析模拟试卷1 一、填空（共30分，每空3分） 1 设??? ? ??-=1511A ，则A 的谱半径=)(a ρ______，A 的条件数)(1A cond =________. 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 型的？
材料物理性能习题与解答
目录 1 材料的力学性能 (2) 2 材料的热学性能 (12) 3 材料的光学性能 (17) 4 材料的电导性能 (20) 5 材料的磁学性能 (29) 6 材料的功能转换性能 (37)
1材料的力学性能 1-1一圆杆的直径为2.5 mm、长度为25cm并受到4500N的轴向拉力，若直径拉细至 2.4mm，且拉伸变形后圆杆的体积不变,求在此拉力下的真应力、真应变、名义应力和名义应变，并比较讨论这些计算结果。 解：根据题意可得下表 由计算结果可知：真应力大于名义应力，真应变小于名义应变。 1-2一试样长40cm,宽10cm,厚1cm，受到应力为1000N拉力，其氏模量为3.5×109 N/m2，能伸长多少厘米? 解： 拉伸前后圆杆相关参数表 ) ( 0114 .0 10 5.3 10 10 1 40 1000 9 4 0cm E A l F l E l l= ? ? ? ? ? = ? ? = ? = ? = ? - σ ε 10 909 .4 0? 0851 .0 1 = - = ? = A A l l ε 名义应变
1-3一材料在室温时的氏模量为3.5×108 N/m 2,泊松比为0.35，计算其剪切模量和体积模量。 解：根据 可知： 1-4试证明应力-应变曲线下的面积正比于拉伸试样所做的功。 证： 1-5一陶瓷含体积百分比为95%的Al 2O 3 (E = 380 GPa)和5%的玻璃相(E = 84 GPa)，试计算其上限和下限弹性模量。若该陶瓷含有5 %的气孔，再估算其上限和下限弹性模量。 解：令E 1=380GPa,E 2=84GPa,V 1=0.95,V 2=0.05。则有 当该陶瓷含有5%的气孔时，将P=0.05代入经验计算公式E=E 0(1-1.9P+0.9P 2)可得，其上、下限弹性模量分别变为331.3 GPa 和293.1 GPa 。 1-6试分别画出应力松弛和应变蠕变与时间的关系示意图，并算出t = 0,t = ∞ 和t = τ时的纵坐标表达式。 解：Maxwell 模型可以较好地模拟应力松弛过程： Voigt 模型可以较好地模拟应变蠕变过程： )21(3)1(2μμ-=+=B G E ) (130)(103.1)35.01(2105.3)1(288MPa Pa E G ≈?=+?=+=μ剪切模量) (390)(109.3) 7.01(3105.3)21(388 MPa Pa E B ≈?=-?=-=μ体积模量. ,.,1 1 2 1 212 12 1 2 1 21 S W VS d V ld A Fdl W W S W V Fdl V l dl A F d S l l l l l l ∝====∝= ===???? ? ?亦即做功或者： 亦即面积εεεεεεεσεσεσ)(2.36505.08495.03802211GPa V E V E E H =?+?=+=上限弹性模量) (1.323)84 05.038095.0()(1 12211GPa E V E V E L =+=+=--下限弹性模量). 1()()(0)0() 1)(()1()(10 //0 ----= = ∞=-∞=-=e e e E t t t στεσεεεσεττ；；则有：其蠕变曲线方程为：. /)0()(;0)();0()0((0)e (t)-t/e στσσσσσστ==∞==则有：：其应力松弛曲线方程为
习 题 一 解 答 １. 设Ａ、Ｂ、Ｃ表示三个随机事件，试将下列事件用Ａ、Ｂ、Ｃ及其运算符号表示出来： (1) Ａ发生，Ｂ、Ｃ不发生； (2) Ａ、Ｂ不都发生，Ｃ发生； (3) Ａ、Ｂ中至少有一个事件发生，但Ｃ不发生； (4) 三个事件中至少有两个事件发生； (5) 三个事件中最多有两个事件发生； (6) 三个事件中只有一个事件发生． 解：（1）C B A (2)C AB (3)()C B A ? (4)BC A C AB ABC ?? (5)ABC (6)C B A C B A C B A ?? ――――――――――――――――――――――――――――――――――――――― ２. 袋中有15只白球 5 只黑球，从中有放回地抽取四次，每次一只．设Ａi 表示“第i 次取到白球”（i ＝1，2，3，4 ），Ｂ表示“至少有 3 次取到白球”． 试用文字叙述下列事件： (1) 41 ==i i A A ， (2) A ,(3) B ， (4) 32A A ． 解：（1）至少有一次取得白球 （2）没有一次取得白球 （3）最多有2次取得白球 （4）第2次和第3次至少有一次取得白球 ――――――――――――――――――――――――――――――――――――――― ３. 设Ａ、Ｂ为随机事件，说明以下式子中Ａ、Ｂ之间的关系． (1) Ａ Ｂ＝Ａ （2）ＡＢ＝Ａ 解：（1）A B ? (2)A B ? ――――――――――――――――――――――――――――――――――――――― ４. 设Ａ表示粮食产量不超过500公斤，Ｂ表示产量为200-400公斤 ，Ｃ表示产量低于300公斤，Ｄ表示产量为250-500公斤，用区间表示下列事 件： (1) AB ， (2) BC ，(3) C B ，(4)C D B )( ，(5)C B A ． 解：(1)[]450,200; (2)[]300,200 (3)[]450,0 (4)[]300,200 (5)[]200,0 ――――――――――――――――――――――――――――――――――――――― ５. 在图书馆中任选一本书，设事件Ａ表示“数学书”，Ｂ表示“中文版”， Ｃ表示“ 1970 年后出版”．问： (1) ＡＢＣ表示什么事件？ (2) 在什么条件下，有ＡＢＣ＝Ａ成立？ (3) C ?Ｂ表示什么意思？ (4) 如果A ＝Ｂ，说明什么问题？ 解：（1）选了一本1970年或以前出版的中文版数学书 （2）图书馆的数学书都是1970年后出版的中文书 （3）表示1970年或以前出版的书都是中文版的书 （4）说明所有的非数学书都是中文版的，而且所有的中文版的书都不是数学书 ――――――――――――――――――――――――――――――――――――――― ６. 互斥事件与对立事件有什么区别？试比较下列事件间的关系． (1) X ＜ 20 与X ≥ 20 ； (2) X ＞ 20与X ＜ 18 ；
1-1一圆杆的直径为2.5 mm 、长度为25cm 并受到4500N 的轴向拉力，若直径拉细至2.4mm ，且拉伸变形后圆杆的体积不变,求在此拉力下的真应力、真应变、名义应力和名义应变，并比较讨论这些计算结果。 解： 由计算结果可知：真应力大于名义应力，真应变小于名义应变。 1-5一陶瓷含体积百分比为95%的Al 2O 3 (E = 380 GPa)和5%的玻璃相(E = 84 GPa)，试计算其上限和下限弹性模量。若该陶瓷含有5 %的气孔，再估算其上限和下限弹性模量。 解：令E 1=380GPa,E 2=84GPa,V 1=0.95,V 2=0.05。则有 当该陶瓷含有5%的气孔时，将P=0.05代入经验计算公式E=E 0(1-1.9P+0.9P 2)可得，其上、下限弹性模量分别变为331.3 GPa 和293.1 GPa 。 1-6试分别画出应力松弛和应变蠕变与时间的关系示意图，并算出t = 0,t = ∞ 和t = τ时的纵坐标表达式。 解：Maxwell 模型可以较好地模拟应力松弛过程： V oigt 模型可以较好地模拟应变蠕变过程： ) (2.36505.08495.03802211GPa V E V E E H =?+?=+=上限弹性模量 ) (1.323)84 05.038095.0()(112211GPa E V E V E L =+=+=--下限弹性模量 ). 1()()(0)0() 1)(()1()(1 //0 ----= = ∞=-∞=-=e E E e e E t t t στεσεεεσετ τ ；；则有：其蠕变曲线方程为：. /)0()(;0)();0()0((0)e (t)-t/e στσσσσσστ ==∞==则有：：其应力松弛曲线方程为1.0 1.0 0816.04.25 .2ln ln ln 2 2 001====A A l l T ε真应变)(91710 909.44500 60MPa A F =?==-σ名义应力0851 .0100 =-=?=A A l l ε名义应变)(99510 524.445006MPa A F T =?==-σ真应力
《材料物理性能》 第一章材料的力学性能 1-1一圆杆的直径为2.5 mm 、长度为25cm 并受到4500N 的轴向拉力，若直径拉细至 2.4mm ，且拉伸变形后圆杆的体积不变,求在此拉力下的真应力、真应变、名义应力和名义应变，并比较讨论这些计算结果。 解： 由计算结果可知：真应力大于名义应力，真应变小于名义应变。 1-5一陶瓷含体积百分比为95%的Al 2O 3 (E = 380 GPa)和5%的玻璃相(E = 84 GPa)，试计算其上限和下限弹性模量。若该陶瓷含有5 %的气孔，再估算其上限和下限弹性模量。 解：令E 1=380GPa,E 2=84GPa,V 1=,V 2=。则有 当该陶瓷含有5%的气孔时，将P=代入经验计算公式E=E 0+可得，其上、下限弹性模量分别变为 GPa 和 GPa 。 1-11一圆柱形Al 2O 3晶体受轴向拉力F ，若其临界抗剪强度 τf 为135 MPa,求沿图中所示之方向的滑移系统产生滑移时需要的最小拉力值，并求滑移面的法向应力。 0816 .04.25.2ln ln ln 22 001====A A l l T ε真应变) (91710909.44500 60MPa A F =?==-σ名义应力0851 .010 0=-=?=A A l l ε名义应变) (99510524.445006MPa A F T =?== -σ真应力)(2.36505.08495.03802211GPa V E V E E H =?+?=+=上限弹性模量) (1.323)84 05.038095.0()(1 12211GPa E V E V E L =+=+=--下限弹性模量
1-1一圆杆的直径为2.5 mm 、长度为25cm 并受到4500N 的轴向拉力，若直径拉细至2.4mm ，且拉伸变形后圆杆的体积不变,求在此拉力下的真应力、真应变、名义应力和名义应变，并比较讨论这些计算结果。 解： 由计算结果可知：真应力大于名义应力，真应变小于名义应变。 1-5一瓷含体积百分比为95%的Al 2O 3 (E = 380 GPa)和5%的玻璃相(E = 84 GPa)，试计算其上限和下限弹性模量。若该瓷含有5 %的气孔，再估算其上限和下限弹性模量。 解：令E 1=380GPa,E 2=84GPa,V 1=0.95,V 2=0.05。则有 当该瓷含有5%的气孔时，将P=0.05代入经验计算公式E=E 0(1-1.9P+0.9P 2)可得，其上、下限弹性模量分别变为331.3 GPa 和293.1 GPa 。 1-6试分别画出应力松弛和应变蠕变与时间的关系示意图，并算出t = 0,t = ∞ 和t = τ时的纵坐标表达式。 解：Maxwell 模型可以较好地模拟应力松弛过程： ) (2.36505.08495.03802211GPa V E V E E H =?+?=+=上限弹性模量) (1.323)84 05.038095.0()(1 12211GPa E V E V E L =+=+=--下限弹性模量. /)0()(;0)();0()0((0)e (t)-t/e στσσσσσστ==∞==则有：：其应力松弛曲线方程为0816.04.25.2ln ln ln 2 2 001====A A l l T ε真应变) (91710 909.44500 60MPa A F =?==-σ名义应力0851 .010 0=-=?=A A l l ε名义应变) (99510524.44500 6 MPa A F T =?= =-σ真应力
数值分析第三次大作业 一、算法的设计方案： （一）、总体方案设计： x y当作已知量代入题目给定的非线性方程组，求（1）解非线性方程组。将给定的(,) i i
得与(,)i i x y 相对应的数组t[i][j],u[i][j]。 （2）分片二次代数插值。通过分片二次代数插值运算，得到与数组t[11][21],u[11][21]]对应的数组z[11][21]，得到二元函数z=(,)i i f x y 。 （3）曲面拟合。利用x[i],y[j],z[11][21]建立二维函数表，再根据精度的要求选择适当k 值，并得到曲面拟合的系数矩阵C[r][s]。 （4）观察和(,)i i p x y 的逼近效果。观察逼近效果只需要重复上面（1）和（2）的过程，得到与新的插值节点(,)i i x y 对应的(,)i i f x y ，再与对应的(,)i i p x y 比较即可，这里求解 (,)i i p x y 可以直接使用（3）中的C[r][s]和k 。 （二）具体算法设计： （1）解非线性方程组 牛顿法解方程组()0F x =的解* x ，可采用如下算法： 1）在* x 附近选取(0) x D ∈，给定精度水平0ε>和最大迭代次数M 。 2）对于0,1, k M =执行 ① 计算() ()k F x 和()()k F x '。 ② 求解关于() k x ?的线性方程组 () ()()()()k k k F x x F x '?=- ③ 若() () k k x x ε∞∞ ?≤，则取*()k x x ≈，并停止计算；否则转④。 ④ 计算(1) ()()k k k x x x +=+?。 ⑤ 若k M <，则继续，否则，输出M 次迭代不成功的信息，并停止计算。 （2）分片双二次插值 给定已知数表以及需要插值的节点，进行分片二次插值的算法： 设已知数表中的点为： 00(0,1,,) (0,1,,)i j x x ih i n y y j j m τ=+=???=+=?? ，需要插值的节点为(,)x y 。 1) 根据(,)x y 选择插值节点(,)i j x y ： 若12h x x ≤+ 或12 n h x x ->-，插值节点对应取1i =或1i n =-，
第七章课后习题答案 7.2 设总体X ~ N(12,4), X^XzJII’X n 为简单随机样本,求样本均值与总体均值之 差的绝对 值大于1的概率. X 解：由于 X ~ N(12,4),故 X 一 ~ N(0,1) /V n 1 ( 2 0.8686 1) 0.2628 10 7.3 设总体X ?N(0,0.09),从中抽取n 10的简单随机样本，求P X ： 1.44 i 1 X i 0 X i 0 X i ~N(0,°.09),故亠-X0r~N(0,1) X 所以 ~ N(0,1)，故U n P{ X 1} 1 P{ X 1} 解: 由于X ~ N (0,0.09),所以 10 所以 X i 2 2 是)?(10) 所以 10 10 X ： 1.44 P i 1 i 1 X i 2 (倉 1.44 P 0.09 2 16 0.1 7.4 设总体 X ~ N( , 2), X 1,X 2,|||,X n 为简单随机样本 2 ,X 为样本均值，S 为样 本方差，问U n X 2 服从什么分布? 解: (X_)2 2 ( n )2 X __ /V n ,由于 X ~ N( , 2)， 2 ~ 2(1)。 1 —n
7.6 设总体X ~ N( , 2), Y?N( , 2)且相互独立，从X,Y中分别抽取 m 10, n215的简单随机样本，它们的样本方差分别为S2,M，求P(S2 4S ； 0)。 解： S2 P(S24S2 0) P(S24S；) P 12 4 由于X ~ N( , 2), Y~ N( , 2)且相互独立S2 所以S12~ F(10 1,15 1)，又由于F°oi(9,14) 4.03 S2 即P F 4 0.01
数值分析第二次大作业 史立峰 SY1505327
一、 方案 （1）利用循环结构将sin(0.50.2)() 1.5cos( 1.2)() {i j i j ij i j i j a +≠+==(i,j=1,2,……,10)进行赋值，得到需要变换的 矩阵A ； （2）然后，对矩阵A 利用Householder 矩阵进行相似变换，把A 化为上三角矩阵A (n-1)。 对A 拟上三角化，得到拟上三角矩阵A (n-1)，具体算法如下： 记A(1)=A ，并记A(r)的第r 列至第n 列的元素为()n r r j n i a r ij ,,1,;,,2,1) ( +==。 对于2,,2,1-=n r 执行 1. 若 ()n r r i a r ir ,,3,2) ( ++=全为零，则令A(r+1) =A(r),转5；否则转2。 2. 计算 () ∑+== n r i r ir r a d 1 2 )( ()( )r r r r r r r r r r d c a d a c ==-=++则取,0sgn ) (,1)(,1若 )(,12r r r r r r a c c h +-= 3. 令 () n T r nr r r r r r r r r R a a c a u ∈-=++) ()(,2)(,1,,,,0,,0 。 4. 计算 r r T r r h u A p /)(= r r r r h u A q /)(= r r T r r h u p t /= r r r r u t q -=ω T r r T r r r r p u u A A --=+ω)()1( 5. 继续。 （3）使用带双步位移的QR 方法计算矩阵A (n-1)的全部特征值，也是A 的全部特征值，具体算法如下： 1. 给定精度水平0>ε和迭代最大次数L 。 2. 记n n ij n a A A ?-==][) 1()1()1(，令n m k ==,1。
1、试说明下列磁学参量的定义和概念：磁化强度、矫顽力、饱和磁化强度、磁导率、磁化率、剩余磁感应强度、磁各向异性常数、饱和磁致伸缩系数。 a、磁化强度：一个物体在外磁场中被磁化的程度，用单位体积内磁矩的多少来衡量，成为磁化强度M b、矫顽力Hc：一个试样磁化至饱和，如果要μ=0或B=0，则必须加上一个反向磁场Hc，成为矫顽力。 c、饱和磁化强度：磁化曲线中随着磁化场的增加，磁化强度M或磁感强度B开始增加较缓慢，然后迅速增加，再转而缓慢地增加，最后磁化至饱和。Ms成为饱和磁化强度，Bs成为饱和磁感应强度。 d、磁导率：μ=B/H，表征磁性介质的物理量，μ称为磁导率。 e、磁化率：从宏观上来看，物体在磁场中被磁化的程度与磁化场的磁场强度有关。 M=χ·H，χ称为单位体积磁化率。 f、剩余磁感应强度：将一个试样磁化至饱和，然后慢慢地减少H，则M也将减少，但M并不按照磁化曲线反方向进行，而是按另一条曲线改变，当H减少到零时，M=Mr或Br=4πMr。（Mr、Br分别为剩余磁化强度和剩余磁感应强度） g、磁滞消耗：磁滞回线所包围的面积表征磁化一周时所消耗的功，称为磁滞损耗Q（ J/m3） h、磁晶各向异性常数：磁化强度矢量沿不同晶轴方向的能量差代表磁晶各向异性能，用Ek表示。磁晶各向异性能是磁化矢量方向的函数。 i、饱和磁致伸缩系数：随着外磁场的增强，致磁体的磁化强度增强，这时|λ|也随之增大。当H=Hs时，磁化强度M达到饱和值，此时λ=λs，称为饱和磁致伸缩所致。 2、计算Gd3+和Cr3+的自由离子磁矩Gd3+的离子磁矩比Cr3+离子磁矩高的原因是什么 Gd3+有7个未成对电子,Cr3+ 3个未成对电子. 所以, Gd3+的离子磁矩为7μB, Cr3+的离子磁矩为3μB. 3、过渡族金属晶体中的原子（或离子）磁矩比它们各自的自由离子 磁矩低的原因是什么 4、试绘图说明抗磁性、顺磁性、铁磁性物质在外场B=0的磁行为。
第七章课后习题答案 7.2 设总体12~(12,4),,,,n X N X X X L 为简单随机样本,求样本均值与总体均值之 差的绝对值大于1的概率. 解:由于~(12,4)X N , ~(0,1)X N {1}1{1}1P X P X P μμ?->=--≤=-≤ 112(11(20.86861)0.262822P ??=-≤=-Φ-=-?-=?????? 7.3 设总体~(0,0.09),X N 从中抽取10n =的简单随机样本，求1021 1.44i i P X =?? >???? ∑. 解：由于~(0,0.09),X N 所以~(0,0.09),i X N 故 ~(0,1)0.3 i i X X N σ --= 所以 10 2 21 ( )~(10)0.3 i i X χ=∑ 所以{}1010222 11 1.441.44()160.10.3 0.09i i i i X P X P P χ==????>=>=>=????????∑∑ 7.4 设总体2 ~(,),X N μσ12,,,n X X X L 为简单随机样本, X 为样本均值，2 S 为样 本方差，问2 X U n μσ?? -= ??? 服从什么分布？ 解： 2 2 2 X X X U n μσ????-=== ???,由于2 ~(,)X N μσ， ~(0,1)N ，故2 2 ~(1)X U χ??=。
7.6 设总体2 ~(,),X N μσ2 ~(,)Y N μσ且相互独立，从,X Y 中分别抽取1210,15n n ==的简单随机样本，它们的样本方差分别为22 12,S S ，求2212(40)P S S ->。 解： 22 22211 2 1 2 22(40)(4)4S P S S P S S P S ?? ->=>=> ??? 由于2 ~(,),X N μσ2 ~(,)Y N μσ且相互独立 所以2 122 ~(101,151)S F S --，又由于0.01(9,14) 4.03F = 即()40.01P F >=
北京航空航天大学 数值分析大作业八 学院名称自动化 专业方向控制工程 学号 学生姓名许阳 教师孙玉泉 日期2014 年11月26 日
一．题目 关于x , y , t , u , v , w 的方程组(A.3) ???? ?? ?=-+++=-+++=-+++=-+++79 .0sin 5.074.3cos 5.007.1cos sin 5.067.2cos 5.0y w v u t x w v u t y w v u t x w v u t (A.3) 以及关于z , t , u 的二维数表（见表A-1）确定了一个二元函数z =f (x , y )。 表A-1 二维数表 t z u 0 0.4 0.8 1.2 1.6 2 0 -0.5 -0.34 0.14 0.94 2.06 3.5 0.2 -0.42 -0.5 -0.26 0.3 1.18 2.38 0.4 -0.18 -0.5 -0.5 -0.18 0.46 1.42 0.6 0.22 -0.34 -0.58 -0.5 -0.1 0.62 0.8 0.78 -0.02 -0.5 -0.66 -0.5 -0.02 1.0 1.5 0.46 -0.26 -0.66 -0.74 -0.5 1. 试用数值方法求出f (x , y ) 在区域}5.15.0,8.00|), {≤≤≤≤=y x y x D （上的近似表达式 ∑∑===k i k j s r rs y x c y x p 00 ),( 要求p (x , y )以最小的k 值达到以下的精度 ∑∑==-≤-=10020 7210)],(),([i j i i i i y x p y x f σ 其中j y i x i i 05.05.0,08.0+==。 2. 计算),(),,(* ***j i j i y x p y x f (i =1,2,…,8 ; j =1,2,…,5) 的值，以观察p (x , y ) 逼 近f (x , y )的效果，其中j y i x j i 2.05.0,1.0**+==。
有效电子数：不是所有的自由电子都能参与导电，在外电场的作用下，只有能量接近费密能的少部分电子，方有可能被激发到空能级上去而参与导电。这种真正参加导电的自由电子数被称为有效电子数。 K状态：一般与纯金属一样，冷加工使固溶体电阻升高，退火则降低。但对某些成分中含有过渡族金属的合金，尽管金相分析和Ｘ射线分析的结果认为其组织仍是单相的，但在回火中发现合金电阻有反常升高，而在冷加工时发现合金的电阻明显降低，这种合金组织出现的反常状态称为K状态。X射线分析发现，组元原子在晶体中不均匀分布，使原子间距的大小显著波动，所以也把K状态称为“不均匀固溶体”。 能带：晶体中大量的原子集合在一起，而且原子之间距离很近，致使离原子核较远的壳层发生交叠，壳层交叠使电子不再局限于某个原子上，有可能转移到相邻原子的相似壳层上去，也可能从相邻原子运动到更远的原子壳层上去，从而使本来处于同一能量状态的电子产生微小的能量差异，与此相对应的能级扩展为能带。 禁带：允许被电子占据的能带称为允许带，允许带之间的范围是不允许电子占据的，此范围称为禁带。 价带：原子中最外层的电子称为价电子，与价电子能级相对应的能带称为价带。 导带：价带以上能量最低的允许带称为导带。 金属材料的基本电阻：理想金属的电阻只与电子散射和声子散射两种机制有关，可以看成为基本电阻，基本电阻在绝对零度时为零。 残余电阻(剩余电阻)：电子在杂质和缺陷上的散射发生在有缺陷的晶体中，绝对零度下金属呈现剩余电阻。这个电阻反映了金属纯度和不完整性。 相对电阻率：ρ (300K)／ρ (4.2K)是衡量金属纯度的重要指标。 剩余电阻率ρ’：金属在绝对零度时的电阻率。实用中常把液氦温度(4.2K)下的电阻率视为剩余电阻率。 相对电导率：工程中用相对电导率( IACS%) 表征导体材料的导电性能。把国际标准软纯铜(在室温20 ℃下电阻率ρ= 0 .017 24Ω·mm2/ m)的电导率作为100% , 其他导体材料的电导率与之相比的百分数即为该导体材料的相对电导率。 马基申定则(马西森定则)：ρ＝ρ’＋ρ(T)在一级近似下，不同散射机制对电阻率的贡献可以加法求和。ρ’：决定于化学缺陷和物理缺陷而与温度无关的剩余电阻率。ρ(T)：取决于晶格热振动的电阻率(声子电阻率)，反映了电子对热振动原子的碰撞。 晶格热振动：点阵中的质点(原子、离子)围绕其平衡位置附近的微小振动。 格波：晶格振动以弹性波的形式在晶格中传播，这种波称为格波，它是多频率振动的组合波。 热容：物体温度升高1K时所需要的热量(J/K)表征物体在变温过程中与外界热量交换特性的物理量，直接与物质内部原子和电子无规则热运动相联系。 比定压热容：压力不变时求出的比热容。 比定容热容：体积不变时求出的比热容。 热导率：表征物质热传导能力的物理量为热导率。 热阻率：定义热导率的倒数为热阻率ω，它可以分解为两部分，晶格热振动形成的热阻(ωp)和杂质缺陷形成的热阻(ω0)。导温系数或热扩散率：它表示在单位温度梯度下、单位时间内通过单位横截面积的热量。热导率的单位：W/(m·K) 热分析：通过热效应来研究物质内部物理和化学过程的实验技术。原理是金属材料发生相变时，伴随热函的突变。 反常膨胀：对于铁磁性金属和合金如铁、钴、镍及其某些合金，在正常的膨胀曲线上出现附加的膨胀峰，这些变化称为反常膨胀。其中镍和钴的热膨胀峰向上为正，称为正反常；而铁和铁镍合金具有负反常的膨胀特性。 交换能：交换能E ex=－2Aσ1σ2cosφA—交换积分常数。当A＞0，φ＝0时，E ex最小，自旋磁矩自发排列同一方向，即产生自发磁化。当A＜0，φ＝180°时，E ex也最小，自旋磁矩呈反向平行排列，即产生反铁磁性。交换能是近邻原子间静电相互作用能，各向同性，比其它各项磁自由能大102～104数量级。它使强磁性物质相邻原子磁矩有序排列，即自发磁化。 磁滞损耗：铁磁体在交变磁场作用下，磁场交变一周，B-H曲线所描绘的曲线称磁滞回线。磁滞回线所围成的面积为铁 =? 磁体所消耗的能量，称为磁滞损耗，通常以热的形式而释放。磁滞损耗Q HdB 技术磁化：技术磁化的本质是外加磁场对磁畴的作用过程即外加磁场把各个磁畴的磁矩方向转到外磁场方向(和)或近似外磁场方向的过程。技术磁化的两种实现方式是的磁畴壁迁移和磁矩的转动。 请画出纯金属无相变时电阻率—温度关系曲线，它们分为几个阶段，各阶段电阻产生的机制是什么？为什么高温下电阻率与温度成正比？ 1—ρ电-声∝T( T > 2/ 3ΘD ) ； 2—ρ电-声∝T5 ( T< <ΘD )；
课后习题 第一章 1.德拜热容的成功之处是什么？ 答：德拜热容的成功之处是在低温下，德拜热容理论很好的描述了晶体热容，CV.M∝T的三次方 2.何为德拜温度？有什么物理意义？ 答：HD=hνMAX/k 德拜温度是反映晶体点阵内原子间结合力的一个物理量 德拜温度反映了原子间结合力，德拜温度越高，原子间结合力越强 3.试用双原子模型说明固体热膨胀的物理本质 答：如图，U1（T1）、U2（T2）、U3(T3）为不同温度时的能量，当原子热振动通过平衡位置r0时，全部能量转化为动能，偏离平衡位置时，动能又逐渐转化为势能；到达振幅最大值时动能降为零，势能打到最大。由势能曲线的不对称可以看到，随温度升高，势能由U1(T1）、U2（T2）向U3（T3）变化，振幅增加，振动中心就由r0'，r0''向r0'''右移，导致双原子间距增大，产生热膨胀 第二章 1.镍铬丝电阻率300K为1×10-6Ω·m加热到4000K时电阻率增加5%假定在此温度区间内马西森定则成立。试计算由于晶格缺陷和杂质引起的电阻率。 解：按题意:p(300k) = 10∧-6 则: p(400k) = (10∧-6)* (1+0.05) ----(1) 在400K温度下马西森法则成立，则: p(400k) = p(镍400k) + p(杂400k)
----(2) 又: p(镍400k) = p(镍300k) * [1+ α * 100] ----(3) 其中参数: α为镍的温度系数约 = 0.007 ; p(镍300k)(室温) = 7*10∧-6 Ω.cm) 将(1)和(3)代入(2)可算出杂质引起的电阻率 p(杂400k)。 2.为什么金属的电阻因温度升高而增大，而半导体的电阻却因温度的升高而减小？ 对金属材料，尽管温度对有效电子数和电子平均速率几乎没有影响，然而温度升高会使离子振动加剧，热振动振幅加大，原子的无序度增加，周期势场的涨落也加大。这些因素都使电子运动的自由称减小，散射几率增加而导致电阻率增大 而对半导体当温度升高时，满带中有少量电子有可能被激发到上面的空带中去，在外电场作用下，这些电子将参与导电。同时，满带中由于少了一些电子，在满带顶部附近出现了一些空的量子状态，满带变成了部分占满的能带，在外电场作用下，仍留在满带中的电子也能够起导电作用。 3.表征超导体性能的3个主要指标是什么?(P80) （表征超导体的两个基本特性完全的导电性和完全的抗磁性） 1)，临界转变温度TC，即成为超导态的最高温度 2)。临界磁场HC，即能破坏超导态的最小磁场，HC的大小与超导材料的性质有关 3)，临界电流密度JC，即材料保持超导状态的最大输入电流 第三章 1.什么是自发磁化？（P142) 在铁磁质内部存在着很强的“分子场”,在这种“分子场”的作用下，原
习 题 1. 指出有效数49×102，0.0490，490.00的绝对误差限、相对误差限和有效数字位数． 2. 将 3.142作为π的近似值，它有几位有效数字，相对误差限和绝对误差限各为多少？ 3. 要使101的近似值x * 的相对误差限不超过4102 1?×，问查开方表时x * 需要保留几位有效数字？ 4. 已知近似数x * 有两位有效数字，试估计其相对误差限． 5. 设x * 为x 的近似数, 证明n x * 的相对误差大约为x * 相对误差的n 1倍． 6. 某矩形的长和宽大约为100cm 和50cm, 应该选用最小刻度为多少cm 的测量工具, 才能保证计算出的面积误差（绝对值）不超过0.15cm 2． 7. 已知三角形面积c ab S sin 2 1=，测量a , b , c 时产生的相对误差为)(*a e r ，)(*b e r ，)(*c e r ，其中2 ,0*π<
课后习题 《材料物理性能》 第一章材料的力学性能 1-1一圆杆的直径为2.5 mm 、长度为25cm 并受到4500N 的轴向拉力，若直径拉细至2.4mm ，且拉伸变形后圆杆的体积不变,求在此拉力下的真应力、真应变、名义应力和名义应变，并比较讨论这些计算结果。 解： 由计算结果可知：真应力大于名义应力，真应变小于名义应变。 1-5一陶瓷含体积百分比为95%的Al 2O 3 (E = 380 GPa)和5%的玻璃相(E = 84 GPa)，试计算其上限和下限弹性模量。若该陶瓷含有5 %的气孔，再估算其上限和下限弹性模量。 解：令E 1=380GPa,E 2=84GPa,V 1=0.95,V 2=0.05。则有 当该陶瓷含有5%的气孔时，将P=0.05代入经验计算公式E=E 0(1-1.9P+0.9P 2) 0816.04.25.2ln ln ln 2 2 001====A A l l T ε真应变) (91710 909.44500 60MPa A F =?==-σ名义应力0851 .010 0=-=?=A A l l ε名义应变)(99510524.445006 MPa A F T =?== -σ真应力)(2.36505.08495.03802211GPa V E V E E H =?+?=+=上限弹性模量) (1.323)84 05.038095.0()(1 12211GPa E V E V E L =+=+=--下限弹性模量
可得，其上、下限弹性模量分别变为331.3 GPa 和293.1 GPa 。 1-11一圆柱形Al 2O 3晶体受轴向拉力F ，若其临界抗剪强度τf 为135 MPa,求沿图中所示之方向的滑移系统产生滑移时需要的最小拉力值，并求滑移面的法向应力。 解： 1-6试分别画出应力松弛和应变蠕变与时间的关系示意图，并算出t = 0,t = ∞ 和t = τ时的纵坐标表达式。 解：Maxwell 模型可以较好地模拟应力松弛过程： Voigt 模型可以较好地模拟应变蠕变过程： ). 1()()(0)0() 1)(()1()(10 //0 ----= = ∞=-∞=-=e E E e e E t t t στεσεεεσεττ；；则有：其蠕变曲线方程为：. /)0()(;0)();0()0((0)e (t)-t/e στσσσσσστ==∞==则有：：其应力松弛曲线方程为 ) (112)(1012.160cos /0015.060cos 1017.3) (1017.360cos 53cos 0015.060cos 0015.053cos 8 2332min 2MPa Pa N F F f =?=? ???=?=? ???=?? ?? = πσπ τπ τ：此拉力下的法向应力为为：系统的剪切强度可表示由题意得图示方向滑移
材料物理习题集 第一章 固体中电子能量结构和状态（量子力学基础） 1. 一电子通过5400V 电位差的电场，（1）计算它的德布罗意波长；（2）计算它的波数；（3） 计算它对Ni 晶体（111）面（面间距d =2.04×10-10m ）的布拉格衍射角。（P5） 12 34 131 192 1111 o ' (2) 6.610 = (29.110 5400 1.610 ) =1.67102K 3.7610sin sin 2182h h p mE m d d λπ λ θλ λ θθ----=???????=?==?=解：（1）= （2）波数= （3）2 2. 有两种原子，基态电子壳层是这样填充的 ; ; s s s s s s s 226232 2 6 2 6 10 2 6 10 （1）1、22p 、33p （2）1、22p 、33p 3d 、44p 4d ，请分别写出n=3的所有电子的四个量 子数的可能组态。（非书上内容）
3. 如电子占据某一能级的几率是1/4，另一能级被占据的几率为3/4，分别计算两个能级 的能量比费米能级高出多少k T ？（P15） 1()exp[]1 1 ln[1] ()()1/4ln 3()3/4ln 3F F F F f E E E kT E E kT f E f E E E kT f E E E kT = -+?-=-=-=?=-=-?解：由将代入得将代入得 4. 已知Cu 的密度为8.5×103kg/m 3，计算其E 0F 。 （P16） 2 2 03 23426 23 3 31 18(3/8)2(6.6310)8.510 =(3 6.0210/8)291063.5 =1.0910 6.83F h E n m J eV ππ---=????????=解： 由 5. 计算Na 在0K 时自由电子的平均动能。（Na 的摩尔质量M=22.99， .0ρ?33 =11310kg/m ）（P16）