矩阵分析及其应用
3.1矩阵序列
定义3.1设矩阵序列｛应)｝,其中A?)=(#))￡Cms,当k—oo, 佝时，称矩阵序列｛A00｝收敛，并称矩阵A=(佝)为矩阵序列｛A00｝的极限,或称｛A00｝收敛于A,记为lim A a)= A或A,k)-> A ks
不收敛的矩阵序列称为发散的。
由定义，矩阵序列A(k)发散的充要条件为存在ij使
得数列站发散。
类似地，我们可以定义矩阵收敛的Cauchy定义
定义31矩阵序列｛A00｝收敛的充要条件为
对任给￡>0存在N(E),当k,l> N(E)时有
IIA(k)-A(/)ll < ￡
其中11.11为任意的广义矩阵范数。
例 1 A(n)
e~n
sin(-)
n y，sin(R) k=l K 7
如果直接按定义我们因为求不出A㈤的极限从
而很难应用定义3.1证明收敛。
相反，由于t^< t^<
v 1/m
从而只要/充分大，则当m, n > /时就有
n
z sin(A)
这样A")收
定理3.1 A(k)->A的充要条件为
HA'10-AII T O
证明：利用广义矩阵范数的等价性定理，仅对co范数可以证明。
即ci IIA(k) -AIL < IIA(k) -All< c2 IIA(k) -AIL
性质 1.设A(k,—> A mxn, B,k,—> B mxn>则
a- A(k)+P ? B(k) -> a- A+P B, V a,PeC
性质2.设A(k)—> A mxn, B,k)—> B nx/,则
A(k)由如一A B
证明：由于矩阵范数地等价性，我们E以只讨论相容的
矩阵范数。
IIA(k).B(k)-A-BII < II A(k) -B(k) -A-B(k)ll+IIAB(k)- A-BII
注意IIB(k)||_||BII,则结论可得。 特别地有 性质2，. A(k U A的充要条件为 A(k) x—Ax,对任意x成立 或者y H A fk) x-> yH Ax,对任意x,y成立. (在无穷维空间中称为弱收敛，但在有限维空间中 和一般收敛性定义是等价的) 对于Hermite(对称)矩阵我们有如下的定理: 设A?), k=l,2,?..,和A都为Hermite矩阵，那么 A(k?A的充要条件为 x”A时X—>x”Ax,对任意x成立 推论：设A如，k=l,2,...,为半正定的Hermite矩阵，且单调减少，即状和4J")为半正定Hermite矩阵，那么4的有极限. 性质3设泌幻和A都为可逆矩阵，且成则 (4伏 证明：因为Af(A如)所以存在K,当必K时有 III-AT?(A(*))II V]/2 我们有(A u))-,= A%( I- A-1- (A(k)? (A(k)r l 从而ll(A(k))-,ll 当k>K时，有 ll(A(k))_,ll 即ll(A(k))-1ll<2-IIA_,ll 因为A—、(A00)% A—】(A 从而II A-1- (A(k))~,ll (当k>K 时) (当krec 时) T O 由定理3.1有(A W A-I 定义3.2矩阵序列｛A00｝称为有界的，如果存在常数 M>0,使得对一切k都有 Il 定理：有界的矩阵序列(A 定义3.3设A为方阵，且当k->oo时有A k-?0,则称A为收敛矩阵。 定理3.2(迭代法基本定理)AkrO的充要条件为谱半径P(A)<1. 证明：必要性：设A*T O,证明p(A) 对4的任意特征值%和相应的特征向量x有 这样我们有4=4 从而 从而有IA.I k 这样有IRvl,由于尢为A的任意特征值， 所以p(A) 充分性。已知p(A) 取￡=(l-p(A))/2 >0,由定理2.10有，存在某种相容的 矩阵范数II.II M使得IIAII M< p(A)+ ￡<1 从而IIA k llM<(IIAII M)k<(P(A)+ e)k 所以当k—8有IIAK|I MT O,从而A5 定理3.3 AJ0的充分条件为存在矩阵范数II.II M使得 IIAII M vl 3.2矩阵级数 定义3.4设矩阵序列｛A%,其中A(k)=(t/^)eC nxn,由它们 形成的无穷和A%A⑴+...+A(k)+...称为矩阵级数, 记为￡"),即有k=0 00 ￡/')=A%A“)+...+A(k)+... N oo 定义3.5记S性，称其为矩阵级数Z A"')的部分和. Jt=() Jt=() 如果矩阵序列｛S(N)｝收敛，且有极限S,即有 S(N)->s 00 那么称矩阵级数闵收敛，且和为s,记为 *=0 S=N A"） k=0 不收敛的矩阵级数称为发散的。 显然^A（k） =S 是指￡#）二%?，V/J 虹0 虹0 即矩阵级数收敛是指它的每个分量所构成的数项级数收敛。 00 性质：矩阵级数闵收敛的充要条件为对任意向量X, k=0 00 向量级数Z A“少收敛。 虹0 定义3.6设矩阵级数￡人°）的每个分量。"所构成的数项k=0 00 S 级数绝对收敛，则称矩阵级数￡ A⑴绝对收敛。 k=Q k=0 关于绝对收敛，我们有如下的定理： 性质1.绝对收敛的￡A以）交换求和次序不改变其绝对 *=0 收敛性和极限值。 性质2.矩阵级数￡ 幻绝对收敛的充要条件为正项级数k=0 ￡11人⑴II收敛。 k=0 性质3.如果矩阵级数￡人伏）（绝对）收敛，那么 k=0 /t=（） 也是（绝对）收敛,且有 k=i）k=0 性质4.设C"11的两个矩阵级数 Si：A ⑴+A ⑵+..裁幻+... S2： B⑴+B⑵+...+B（k）+... 都绝对收敛，其和分别为A和B.则矩阵级数 S3： A⑴B⑴+ [A⑴B⑵+ A⑵B⑴]+... +[ A⑴ B很)+ A⑵ B(k-I) +.. .+A°° B⑴]+... 绝对收敛且和为AB. 证明：由于S「A⑴+A⑵+...+A°°+...绝对收敛的充要条件为正项级数IIA⑴II+IIA⑵ll+?..+IIA(k)||+...收敛且与排列无关。 我们证明的思路是证明正项级数： IIA ⑴ B ⑴ 11+ IIA ⑴ B ⑵+ A ⑵ B ⑴ II+... +IIA⑴ B(k)+ A⑵ B(k-1)+...+犬)B⑴II+... 收敛。引用魏氏定理，我们仅需验证下列正项级数： IIA⑴ll?IIB⑴ 11+ { IIA⑴ II.IIB⑵ 11+ IIA⑵ IIJIB⑴ 11}+... + {IIA⑴ IMIB(k)ll+ IIA⑵ 11?1 旧(k-1)ll+...+IIA(k) IIJIB⑴ 11}+... 收敛。这由题设 正项级数IIA⑴II+IIA⑵II+.. .+1山唧+... 和正项级数IIB⑴11+1旧⑵II+.. .+IIB唧+... 的收敛性可得。 定理3.4幕级数I+A+A2+...+A k+...ifc敛的充要条件为 A的谱半径p(A) 若有矩阵范数11.11使得IIAIIvl,则 ll(I-A)-1- (I+A+A2+.. .+A k)ll 证明：必要性.由于I+A+A2+...+AL...收敛，从而 S(k)= I+A+A2+...+A k收敛。记T(k)= I+A+A2+...+A k+1, Ak+l=T(k)_ S(k)收敛,且T、S(k)TO,这样我们有A k^0,从而p(A) 充分性：设p(A) I+A+A2+...+A k=(I-A)_1 -(I-A)-1 A k+1 因A J O,所以 I+A+A~+.. .+人杞+...—(I—A) I 又因为 (I-A)-1 - (I+A+A2+...+A k)= (I-A)-1 A k+1 从而 -(I+A+A2+...+A k)ll=ll (I-A)-1 A k+,ll 设B=(I-A)-1A k+1,Affi(I-A)B=A k+,即 B=AB+ A” 从而 IIBII< IIAII-IIBII+ IIA k+I ll< IIAII-IIBII+ IIAII k+1 因为矩阵范数ll.ll使得IIAIIvl,所以IIBII 定理3.6设帛级数f (z) = ￡c盘的收敛半径为尸，如果方阵A满足p(A)< r,则矩阵界级数 00 /（A） = Z Q A*是绝对收敛的；如果p（A）>尸，k=0 00 是发散的。 A=0 证明：利用绝对收敛的性质。 反之，设A的特征值人满足UI=/?（A）, x为入相应的 特征l'*J量￡以（A七）=￡以（#尤）=（￡ c k分）尤， k=0 A=0 A=0 由于p（A） >r，那么（￡c、/）x发散（注意x为非零仙景）k=0 〃s 从而Z Q（A%）发散，这样Z G**发散。 A=0 k=0 矩阵函数 定义：设一元函数Az）能展开为Z的蓦级数 3 /'(Z)= ?N、zlvr) k=0 其中r>0表示该帛级数的收敛半径。当〃阶矩阵A的 S 谱半径p(A) k=0 00 为f(A),即f(A)= . k=0 性质1(代入规则)：若/Xz)能展开为z的蒂级数，且f (z)=g(z), 对Izl < r成立，则当p(A)< r时，f(A)=g(A). 矩阵函数举例： sin⑵=々一//3!+广/5! -... 则sin(A)=Z-A3/3!+A5/5!-... cos⑵=1 一『/2!+z'/4! -... cos(A)= I-A2/2!+A4/4!-... e:-1 +Z+Z2/2 ! +//3!+... /=I+A+A2/2!+A3/3!+... sir?⑵ + cos2U)= 1 可得：sin2(A)+cos2(A)= I 性质2二元函数f(x,y)能展开为x,.y的蓦级数J\x,y)=g(x,y). 若AB=BA^\ f(A,B)=^A,B)(二元函数的代入规则). 矩阵函数值的求法 1.待定系数法 设〃阶矩阵4的特征多项式(p(2)=det(27M).如果首1 多项式w(A)=/+b|/l T+...+b m_n+bm 满足:⑴\|/(A)=0;(2)侦人)整除板人)(矩阵A的最小多项式与特征多项式均满足这些条件).那么,wQ)的零点都是A的特征值.记w(人)的互异零点为 入],人焉,相应的重数为「1,...小(1"1+「2+...+板=01),则伯 W(“Q)=O (/=o,l,..s?l;i=l,2,...,s) 这里,表示wQ)的I阶导数(下同). 00 设、人沪= W(z)g(z)+r(z).其中r⑵是次数低于刀的k=0 多项式，于是可由产(入i) = r(/)( 确定r(z). 利用f(A)= w(A)g(A)+r(A)=r(A). 因此我们的问题就是给定函数,/k),由约束条件 r (z ) ( Ai )=fXk) /=0,l,...,r 「l;i=l,2,...,s 确定r ⑵。 若知道函数Ax)在X 。的函数值和直到〃阶导数值， 则由Tayler 展开式可得多项式： <*) W)=AAo)tf\x o )(x-xo)+[/， \xo)/2!](x-x o )- +[仁3)(如/3 !](XTO )3+. . .+[严)(如加!](XF ))" 使得1/%0)寸％0),/=0,1,2,...,〃成立. 相应的若知道函数冷)在Xo ，XE..E 的函数值 则有相应的Newton 插值公式 N/U)=f(xo)±/>[XO ，X1](E ；O )±/'[XO ，X1，X2](X-XO )(E1)+... 甘1工0叫,.M 〃](X ?X ())(E]).. .(X*) 使得 成■)寸3)，i =0,1,2,...,〃 成立. 其中N/(x)的项为均差。 均差的定义和简单性质如下： 均差及性质 定义：,/[^o^k]=(A U ()))/(Xk-x 0) 为,f(X)关于点XoAk 的一阶均差； /Uo,X| 知=(/ [心 JkLf [Xo,X|])/(XrX]) 为二阶均差。一般地，称 亢心1,...阪]=(/*[心1,...由2 知mxoN ，?? N-1])/(XE1) 为.f 3)的上阶均差。 性质1.#阶均差可表为函数值的线性组合； 性质2.均差与节点的排列次序无关； 性质3(均差与导数)若八x)在[。四上存在〃阶导数， 且结点Wl,*2,...*"]，则存在&w[x°,x ，使得 』=")(&)/〃！ 成立, 其中 x°=min{xo,X|X2...x /f ), x Jmax{xo/il2,...x”}。 推论:若/3在[白力]上存在n 阶导 数，则对任何ce[a,b], f\c,c,... ,c]=f {n Xc)/n!成立。 证明：在性质3中令X (),X I X 2….X H 同时趋于c 可得结论. 由此可见Newton 插值公式可以看作Tayler 展开式的 推广，Tayler 展开式只不过是插值节点重合的情形， 相应的函数值变为导数值而已。 均差的计算: X o X] X2 *3 X4 ./Uo) . /(引 ./fe) 伯 ,/UoAil ./Ul 网] ,心3&] /Ix o ,xi,x 2l ./IX|,X 2,X 3] ./[X2,X 3,X 4] .心0乂,的/3] /[XU2,X3，X4] /1 心1,双,工3，罚 相应地， 若有X )=X2=X 3,则有 X() ./Uo ) ./U"] X\ ./Ul) ,r (xi ) ,gxui] 必1) 广⑴) f 〃W)/2! X4 A^4) .心两 yUi,x b x4J /[xoAbXi^iAd 其它的类似.因此我们可以使用Newton 插值公式 得到满足插值条件的多项式. 2 0 0 例1 设A= 1 I 2 ,求? 1 -1 3_ 解 (p(X)=det(XI-A)=a-2) (X-5) (X+1) z ./(z)=e z 1 2 e 2t 5 e 5t (/(5)-/(2))/(5-2)=(e 5t -e 2t )/3 -1 e-t (f(-l)顼5))/(-l-5)=("-eT)/6 (e 5t -2e 2t +e-l )/18 从而 /(z)=e z 1 =(p ⑵P(z)+r(z), 其中 r(z)= e 21 +(z-2) (e 51 -e 2 l )/3+(z-2)(z-5) (e 51 -2e 21 +e'l )/18 =e 2t (l +(z-2)/3-(z-2)(z-5)/9) 1 2 -1 1 0 1 0 0 1 0 0 -1 -1 -I + e 5l ((z-2)/3+(z-2)(z-5)/18) + e -t (z-2)(z-5)/18 XA)=e ' t =r(A)= e 2 l (I+(A-2-I)/3-(A-2-I)(A-5-I)/9) + e 5 t ((A-2-I)/3+(A-2-I)(A-5-I)/l 8) + e 」A -2-I)(A-5-I)/18 1 0 o - ■ 0 0 0 一 5/9 4/9 4/3 +任 2/9 -2/9 1/3 2/3 -2/3 16/9_ 1/6 -/6 1/9 0 1/9 1/9 一 1/3 1/6 1/6 -2/9 - 2 0 O - 例 3.5 设 A = 1 1 1 ,求 e Al 1 一 1 3 解(p(X)=det(XI-A)=(X-2)3 z ,/(z)=e" 2 e 2t 2 e 2t /v (2)= te 2t 2 e 2t .质2)=免" f f, (2)/21= t 2 e 2l /2 从而 /(z)=e" =(p(z)P(z)+r ⑵， 其中 r(z)= e 2' +t e 2 l (z-2)+t 2e 2 l (z-2)2/2 y (A)=e A,=r(A)= e 2z (I+r (A-2?I)+ ?(A-2-I)2/2 ) lo o _ = e 2r t l-t 1 t -t 1+， 求RA),其中f(z)= 解：由于 det(XI - A)=( A-1 尸(入 + l )=(p(A ) 12-11 1 -5 -3 0 0 10 0 0 10 0 ,A 2 = 0 0 10 0 0 10 0 -1 -1 -1 0 0 0 - 1 -7 -4 0 10 0 0 0 10 0 -1 -I -1 求 r ⑵使得f (沪P(z) z f (z) 1 例： A = 16 1() Z -z J/= 1 0 .厂(1)=6 1 0 f'(l)=6 广'(1)/2!=75 -10 (0-6)/(-l-l)=3 (3-75)/(-1-1)=36 因此 r(z)= 0+6(z — 1)+ 75(Z -1)2+36(Z -1)3 从而/(A)=r(A)= 6(A-1)+ 75(A 一 I )2+36(A -I)3 即 r(A)= —9 0 0 0 2. 数项级数求和法。 利用 A m=kol+k | A+l<2A~+...+k m _| A m ' 代入计算即可。 3. 对角阵法 若A 相似于对角阵，求出它的对角矩阵和相似矩阵 即 A=PDP -1,则 f(A)= Pf(D)P -1 4. Jordan 标准型法 设A 的若当标准形为J,则存在可逆矩阵P 使得 P A P= J= nii "4 ... PC /(J 沪 ￡《.J i = k=0 ，叫 fW … 1 ? f(4) (叫 -1)! 广一')(4) 即) /(A) 8 8 8 心 AA)= ￡4PJ*P-、P(Z%J*)pT = P k=0 k=0 7(J.) _ =p ?■.p_, /(J J 可见，求/(A)可转化为求A的若当标准形J及变换矩阵P的问题。 三、广义矩阵函数 利用前面的定理3.6及其推论的结果，对于可以展开成昴级数形式的解析函数.怎),它的矩阵函数的定义方式就是在它的幕级数形式中用矩阵A替代z,就得到了由矩阵羸级数形式的矩阵函数f(A)的定义。但是，对于任意给定的函数f(z)能够展开成收敛界级数形式的要求条件很强，一般不容易满足，例如/⑵=ln⑵就不满足。借助于前面讨论的矩阵函数若当标准形求法的公式(*)和(**),我们可以拓宽矩阵函数的定义如下. 定义3.8设Aer xn的Jordan标准形为J,即有可逆矩阵P,使得 1 . ?? ?? 儿1 2 如果函数,彻在％处具有直到四-1阶导数(，=1,2,…粉, AJ/)= /(A) fW -Jzn-xw. /(J.) M)= P .?. /(JJ_ (**) 那么，称为对应于木)的矩阵函数。 这样定义的矩阵函数在Rz)可展开成慕级数的形式和前面的定义一致，而当/(z)不能展开成幕级数的形式时仍然可以定义，因此确实拓广了矩阵函数的定义。 定义3.9设解析函数，定义矩阵函数/(A)为 AA)= Ja-/-A)-7aMz 密) 其中O为包含A的全部特征根的单连通区域，而8Q表示区域的边界。 4矩阵的微分和积分 定义3.9如果矩阵A(t)=(?y(O)W x W的每一个元素a^t) 是变量r的可微函数，则A⑺关于i的导数(微商)定义为d 4'(。= ( —),心" dt 性质I. — (A(t)+B(t))= — A(t)+ — B(t) dt dt dt 性质 2. — (A(t)B(t))= — A(t) B(t)+A(t)- — 5(t) dt dt dt r d d 性质 3. 一(a(t)A(t))= a'(t)A(t)+ a(t) — A(t) dt dt 性质4.如果A⑺和—A(r)可交换，々)和【无关的一元解析函数, dt 则有g/(A(t))=，(A(t)) 4 A(t) dt dt 特别注意若A(r)和—A⑺不可交换，则上式不一定成立。 dt 00 8 证明：将念)展开成昇级数形式/(z) = ￡cb ,得到/(A) = ^c,A" °k=Q k=0 根据假设A(f)和—AQ)可交换可得dt 心⑴)/出=kA k~\t) d(A(t))/dt 代入矩阵函数,AA(/))得￡c k kA k'{(。籍=d A(t)/dt, f'(A(t))成立。 k=\ 考虑到A。)和—A(t)uT交换的条件比较强，我们可以给出如下的性质。dt 性质5.如果A⑺和B(t河交换，为)是与I无关的一元解析函数， 则有 4 tr(/(A(t))B(t))=tr[B(tyXA(t))竺 A(t) ] + tr[/(A(t))竺 B(t)] dt dt dt oc oc 证明：将.怎)展开成幕级数形式f(z) = ￡cb ,得到/(A) = X C X X=0 Jl=0 根据假设AQ)和B⑺可交换可得 tr[d(A%)8(。)/次]=tr伙?A—(。d(A(t))/dt)? B(t)]+tr[A*(。d(B(g] mxn df 咨n 代入矩阵函数 8 ^(tr[AA(/))B(t)])/^= tr[(^c k kA k~l (f)籍)BQ)]也[心*。))d(B(t)/dt\] k=\ d d 顼BW(A(t))秋tEg(t))石B(t)] 成立。 定义3.10如果矩阵A。)的每个元素徇什)都是区间际川上的连续函数，贝U定义A(t)在gq］上的积分为分量积分构成的矩阵.即 f A(t) dt=( f' aij(t)dt) 性质：(积分算了仍为线性算了) f (A(t)+B(t))dt= f A(t)dt+ I P?A(t)?Qdt=P( A(t)dt) Q *0 *0 当aij(t)都在代展］上连续时，就称A⑴在R,q］上连续.且有 d f 一A(s)ds=A(t) dt土 当aij(t)都在［a,b］上连续时 A r(s)ds=A(b)-A(a) 其它微分概念 1函数对矩阵的导数(包括向量) 为了简单起见，在木节中我们记0为第，个分量:为1,而其余分量为0的单位列向量，其维数根据上下文确定，以后不一一说明。 定义:设X=(&j)mxn,mn元函数 /{X)=/(&11, Win, &mn) 定义RX)对矩阵X的导数为 8f ■ ■ ■ , 8f . . . 』I 咨〃 i=i j=i °Gij 特别注 意在上面的等式 中，们为m维列 向量而勺为〃维 列向量,在本节中 我们会经常使用 这种约定。 矩阵的元素和矩阵符号的关系 设A=(6/.. )eC mxn,那么 这个关系式实际上可以为我们把一个关于分量的等式转化为关于矩阵■向量的形式。最后一个式了可以由第二个式了结合vec(?)向量化算子的定义推得。 性质1.1若g(。项x), x=x(t)为向量，其中俭)为n元函数，则d^/dt=dx T/dt-d/7dx. 性质1.2若川)=f(X), X=X(t)为矩阵，其中.KX)为"叫元函数, 则dg/dt=tr(dX v/dt-df/dX)性质1.3：关于迹函数的导数 d(tr(X7A))/dX=A. d(tr(A T X))/dX=A. d(tr(AX))/dX=A T. 例/(X)=X T A X,则d/7dx=(A+A T)x, 特别地，若A为对称矩阵，g(x)=x7Ax/2, 则dg/dx=Ax. 证明：由于df/dxi=(dx/dxi)7 A x+x7 A (dx/dx,) =e： Ax+ x r Ae； =](A+A，)x “ df n 所以d/7dx = ￡旦y = ￡ 与(A + A' )x =(A+A1)x, j=l dXj /=! =IAI tr(A-LdA/dt) \)d(\A\)/dt=\A\tv[A~l 'dA/dt] .f f dlAI T 2).由于d(\A\)/dA=YY —— i=l j=l da 。 根据I)的结果我们知道 d(\A\)/daij= L4I tr[A -1- dA/da^ L4I ]=14lejA -1^ 因此 d(L4l)"A=Z￡l AkjA -l e. i=l j=\ = ￡￡l AI (AT )M 《=1 A I (A-y /=l j=l 从而 J(loglAI)/JA=(l/L4l)-t/(L4l)/JA =(A T )-! 例1)设非奇异矩阵A(f)的对，可导,则d(\A\)/dt =IAI tr(A - 2) 对于任给非奇异矩阵A, d(loglAI)/dA =(A r )-1 3) 对于任给 AcR'"", J(logl<5I+A 7'AI)/JA= 2A(3I+A 『A)- dAldt) " da 证明：1)由于"141/出=￡|(。］,。2，???,一 dt 〃 1(%,? =IAl￡ —— /=1 da i ?，-y ，…， 以 dt _______ \A\ 考虑线性方程组A ?y=dajdi,其中《?为矩阵4的第，列。 根据Gramer 法则可得 da X =1 (%a,…,一T 、?"〃)SAI at 另一方面，我们知道y=A^? dajdi,所以有 月=1 (丹,劣,...，华1/1 A I 1牛 at dt 因此有 d\A\/ 14|￡矿妃 /=! dcij dt K^Cp K^Cp K4^^4^ Cp Cp Cp K^Cp K^Cp Cp K^CpCJ) K4^C|) Cp K^C|) K^C|) K^Cp 或者det(A)=^6Z. A..,其中禹?为在矩阵A中；=i 去掉第i行第j列的代数余了式。排列缶得如下矩阵 Ai A2 …Aw 7 &] A22 (2) A = . . . ^ ? ? ?? ? ? ?? _A H A” 2 …A nn _ 则根据代数余了式的定义可得 A(A)T=L4I7 因此A=L4I(A T)_,. 注意禹显然和⑶无关，因此我们有d(HI)/dA=A =IAI (A7 )"1 3)由于d(\SI+A T M)/dA= i=l j=l dcijj 根据1)的结果我们知道 d(\Sl+\T\ l)/d佝=VI+A’AI tr^^I+A^A)-1-d(81+UA)/daij\ =ISI+A'AI trf(^l+A7A)-1- (eje^A+A‘W《)] =IJI+A r AI( trKJI+A^Ar1-^^ A]+ tr"I+A7A)-iAW《)]} =I^I+A7AI( A(3I+A7A)-W+《(SI+A'A^A W)]} =2-l^l+A7AI I A(^I+A7A)-1^ m . 因此d(匹I+A’A I)04=2?I5I+A7AI￡￡4A0?I + A，A)%/ ?q《 /=! j=l m 〃 =2?I5I+A7A I ￡Z e i e i A0 ? I + A r A)_, i=i j=i m n =2d^IM r AI(^e,^)A(^4 + A7A)_1^^J 1=1 j=l =2.l徂+A『AI ?A(5I+A7'A)T 从而d(logl母+A7AI)〃4=(1/ISI+A7'A l)-J(I^I+A r A 1)04 =2-\(Sl+\T^ df dX 前面讨论的函数对矩阵的微分定义没有考虑特殊矩阵的性质，所以有时在公式表达时针对某 些特殊矩阵并不能得到很简洁的表达性质。 举例来说，对于对称矩阵X 来说，我们考虑函数XX)=tr(A 7X),其中A 为任意〃阶矩阵, 如果仍然使用前面的关于函数对矩阵的微分定义我们有 91 Af ?牛；=A+A ，—D ,=1 j=\ gj 其中D 表不A 的对角部分构成的对角矩阵。这种表达形式看起来不是很简洁。 这是因为对于对称矩阵A 来说，我们有如下关系式 A = +EE +e j e i ) /=1 /=1 J=1+1 a u = e] Ae i =e T j Ae l J J J 由前面的讨论，在很多情形下，如果X 为,WI 对称矩阵, .f(X)为〃(〃-1)个变元的函数 .人X)=AS11，???，氟，色2...,如&〃) 那么定义 I 8f 2咨, . . . 这时对于对称矩阵的函数,/(X)=tr(A r X), 我们有#〃X=(A+A 『)/2, 特别地，当A 也为对称矩阵时有d (tr(AX))/t/X=A 0 矩阵函数对矩阵的导数 (主要考虑向量值函数对向量的导数) 为了进行下面的讨论我们引入A 与B 的直积定义. 定义5.9设A=(a 萨Cf,B=(b 瞧CT 。，则称 如下的分块矩阵 df 咨n 特别的, 当 F 为mx ］的向量和X 为HX I 的向量时，那么 其中 Of\ dF dX 7 对于前面的定义我们E 以进行形式化的表示为: A? B = a 7,B a \2B — a 2n B *5 (5.4.1) %B ?.. (Kronecker 积)。 为A 与B 的直积 性质：如果A 为nix 1的向量和B 为nxl 的向量时，那么 定义 3.13 设 X=（&j ）mxn 的 砰）明（京 mn 元函数 G 〃i,?… F(X)= /H rl frs 对矩阵x 的导数如下: dF dX dF dF /=! J=1 dX ?dF —19.? .9Z7Z, J — ， (2) 类似地，我们可以定义当然只是形式上的表示而已. e n e q e n e e e e e a e 性质i：$=i±i±%("WM) OA i=\ j=\ k=\ /=! Wij 证明：根据==￡￡W。就代入化的定义可得。 临.?。乌dx 实际上有所谓的转换定理 转换定理：设X=(勺)和Y=()为)分别是mxn和pxq的矩阵，A,B,C和D 分别为pxmjixq,pxR和mxq的矩阵，它们可以是X, Y的函数，则如下两条是等价的 dY T ^― = AEj.(m x 〃)B + CE；(m x n)D , 袅=A1 Ejj(p X q)B，+ DE；(p X q)C , i = 1,…,p, j = l,?..,q 其中Egi)表示第i行第j列元素为1,而其余元素为0的mxn矩阵. 证明：为了简单起见，我们记门为第i个分量为1,而其余分景为。的单位列向量，其维数根据上下文确定。 首先证明⑴T(2),实际上 .r SY €k€} =e[ Aepj Bq + e\ Ce^ De t =a ki b ji + c kj d u 由于 如V V代/? ax s=i r=1dx st 22(岫j+C’d sj)逐; =EE(e； AL* BW + e； De *CeMe； s=l r=l 〃i 〃m 〃 =V V* e e. A1 e.e1. B7 e e7 + VV e .e. De -e-Ce e1 ——S S I J t t —— $ $ J I f t 5=1 t=\S=1 t=\ =￡以(A，Wj B f疙以+芸(DeR C疙以 S=1 E=1 5=1 7=1 =A，Ejj(pxq)B，+ DE： (pxq)C 而证明⑵T(l)是一个类似的过程。 性质2： y=Wx,函数Ay)是|何量y的函数 贝U dy'7dx=W T, d顷y))/dW=d(/(y))/dy?x‘ 性质3：设Rx)是向量x的函数,而x又是11的函数,则d/7du=dx T/du -d/7dx 根据定义dx l/du=(dx'/dz/|, dx^dz^.., dx'/dzz；,)1 其中(hJ/d飨=(dri/d灿,dxjduj) 类似地，设.心)是的量X的函数，而X又是|何量U的 的向量值函数，U是V的向量值函数，贝IJ d/7dv=du7 /dv-dx1 /du d/7dx 证明:由题设危)是向量ii的函数,所以 d/7dv= du^dv-d/Zdu 而必du= dx */du d//dx 从而df/dv=du1 /dv-dx1 /du -d/7dx 成立。 性质4：设,心)是以矩阵4的元素为变量的nm元函数, 而4又是标量，的函数，则 d/7dr=tr(dA1/dt -df/dA ) 应用：考虑如下前馈神经网络模型模型的学习问题： u=Wx V=(p(w)(即V,=(p(w,-),其中，?(?)为一元可微函数) p=Uv g=e(p)(即0 as),其电巾(?)为—■元可微函数) o =Vq e=(l/2). Ily-olF 其中：x, u,v,p f q,o为向量,W,U,V为矩阵. 求de/dV, de/dU,de/dW. 解：”性质2： y=Wx,函数./(y)是向量y的函数 则dy T/dx=W T, d(Ay))/dW=d(/(y))/dy? 性质3： d/7dv=du'/dv-dx'/du -d//dx " de/do=o-y 中科院矩阵分析与应用大作业 实现LU分解 QR分解 Householder reduction、Givens reduction Matlab 代码： function [] =juzhendazuoye A=input('请输入一个矩阵A='); x=input('请输入序号 1 LU分解 2 Gram-Schmidt分解 3 Householder reduction 4 Givens reduction:' ); if(x==1) %%*************LU分解*****************%% disp('PA=LU') m=size(A,1); % m等于矩阵A的行数 n=size(A,2); % n等于矩阵A的列数 if(m==n) % 判断矩阵A是不是方阵 % 如果矩阵A不是方阵那么就输出“error” U=A; % 把矩阵A赋值给矩阵U L=zeros(n); % 先将L设为单位阵 P=eye(n); % 首先将交换矩阵P设为单位矩阵 for j=1:n-1 for i=j+1:n if (U(j,j)~=0) %判断主元元素是否不为0 L(i,j)=U(i,j)/U(j,j); U(i,:)=U(i,:)-U(j,:)*U(i,j)/U(j,j); % U(j,j)为主元元素 else a=j+1; % 令a等于j+1 while((U(a,j)==0)&&(a 1.简述模式的概念及其直观特性，模式识别的分类，有哪几种方法。（6’） 答（1）：什么是模式？广义地说，存在于时间和空间中可观察的物体，如果我们可以区别它们是否相同或是否相似，都可以称之为模式。 模式所指的不是事物本身，而是从事物获得的信息，因此，模式往往表现为具有时间和空间分布的信息。 模式的直观特性：可观察性；可区分性；相似性。 答（2）：模式识别的分类： 假说的两种获得方法（模式识别进行学习的两种方法）： ●监督学习、概念驱动或归纳假说； ●非监督学习、数据驱动或演绎假说。 模式分类的主要方法： ●数据聚类：用某种相似性度量的方法将原始数据组织成有意义的和有用的各种数据 集。是一种非监督学习的方法，解决方案是数据驱动的。 ●统计分类：基于概率统计模型得到各类别的特征向量的分布，以取得分类的方法。 特征向量分布的获得是基于一个类别已知的训练样本集。是一种监督分类的方法， 分类器是概念驱动的。 ●结构模式识别：该方法通过考虑识别对象的各部分之间的联系来达到识别分类的目 的。（句法模式识别） ●神经网络：由一系列互相联系的、相同的单元（神经元）组成。相互间的联系可以 在不同的神经元之间传递增强或抑制信号。增强或抑制是通过调整神经元相互间联 系的权重系数来（weight）实现。神经网络可以实现监督和非监督学习条件下的分 类。 2.什么是神经网络？有什么主要特点？选择神经网络模式应该考虑什么因素？ （8’） 答（1）：所谓人工神经网络就是基于模仿生物大脑的结构和功能而构成的一种信息处 理系统（计算机）。由于我们建立的信息处理系统实际上是模仿生理神经网络，因此称它为人工神经网络。这种网络依靠系统的复杂程度，通过调整内部大量节点之间相互连接的关系，从而达到处理信息的目的。 人工神经网络的两种操作过程：训练学习、正常操作（回忆操作）。 答（2）：人工神经网络的特点： ●固有的并行结构和并行处理； ●知识的分布存储； ●有较强的容错性； ●有一定的自适应性； 人工神经网络的局限性： ●人工神经网络不适于高精度的计算； ●人工神经网络不适于做类似顺序计数的工作； ●人工神经网络的学习和训练往往是一个艰难的过程； ●人工神经网络必须克服时间域顺序处理方面的困难； ●硬件限制； ●正确的训练数据的收集。 答（3）：选取人工神经网络模型，要基于应用的要求和人工神经网络模型的能力间的 匹配，主要考虑因素包括： 矩阵分析及其应用 3.1矩阵序列 定义3.1设矩阵序列｛A (k )｝,其中A(k)=( a (k ) ) C m n ,当k a j" a u 时，称矩阵序列｛A (k)｝收敛，并称矩阵 A=( a ij )为矩 阵序列｛A (k)｝的极限，或称｛A (k)｝收敛于A,记为 lim A (k) A 或 A (k) A k 不收敛的矩阵序列称为发散的。 由定义，矩阵序列 A (k )发散的充要条件为存在 j 使 得数列a (k)发散。 类似地，我们可以定义矩阵收敛的 Cauchy 定义 定义3.1'矩阵序列｛A (k)｝收敛的充要条件为 对任给>0存在N()，当k, l N()时有 ||A (k) A (l)|| < 其中||.|为任意的广义矩阵范数。 sin 』) n n sin(k) 如果直接按定义我们因为求不出 A (n)的极限从而 从而只要I 充分大，则当m, n > l 时就有 sin(k) k 2 这样A (l)收敛。 定理3.1 A (k) A 的充要条件为 ||A (k) A|| 0 证明：利用广义矩阵范数的等价性定理，仅对 范数可以证 明。 即c 1 IL A (k) A|| ||A (k) AII C 2 ||A (k) AII 性质 0 若 A (k) A ，则 ||A (k) II IIAII 成立。 性质 1. 设 A (k) A m n ， B (k) B m n , 则 A (k)+ B (k) A+ B ， ,C 性质 2. 设 A (k) A m n ， B (k ) B n l ,贝U A (k) B (k) A B 证明：由于矩阵范数地等价性，我们可以只讨论相容的 矩阵范数。 ||A (k )B (k) A B|| || A (k) B (k) A B (k)||+||AB (k) A B|| || A (k) A|| ||B (k)||+||A||||B (k) B|| 例 1 A (n) k m 1 k(k 1) 相反，由于 现场质量管理与突破性快速改善 课程目标： 建立品质稳定的现场及快速突破性改善的团队和机制，通过系统化改善大幅度降低运营成本。参加人员： 生产现场管理人员、质量经理/主管、项目/质量/工艺工程师等。 课程大纲： 第一天9：00~16：30 ●第1讲:标准化作业 ●现场质量管理的基本要义 ?问题提出与团队组建 自我介绍/分组 ?培训目标及要求 团队突破训练 ?制造业的微笑曲线与痛苦指数 ?现场控制的指标意义 ?中美日制造现场管理比较 ?现场管理四层修炼 观看优秀现场录像 ?卓越现场管理7要点 ?回答学员问题 提问及回答老师问题 ●标准作业的意义与基本要求 ?如何减少员工随意操作带来的浪费及品质风险 ?如何减少对员工经验的依赖及流动风险 ?如何科学设计作业程序以减少劳动强度并提升效率 ?标准化与持续改进 ●标准化作业 ?标准作业的SMART原则 ?标准作业三票 ?SIP的要点与要求 ●培训的标准化 ?作业培训标准化要求 第二天9：00~16：30 ●第3讲:制造过程监控与审核 ●制造过程监控 ?对KPC的预控制 ?X-R控制图 ?过程能力分析与改善 ?KCC的识别与确定 ?KCC监控与数据分析 提问及回答老师问题 ●制造过程审核 ?墨菲定律 ?分层审核的目的/策略及要点 ?过程审核方法及要点 ?分层审核与过程审核的综合运用 分组练习1 ●第4讲:不合格品控制(S.I.R) ●不合格品与可疑品 ?定义不合格品与可疑品的方法及 意义 ?不合格品与可疑品的发生与发现 ●不合格品隔离/评审与处臵 ?不合格品隔离方法 ?不合格品快速评审 ?不合格品处臵方式及跟踪 ●废品率降低与质量成本分析 ?废品率的统计/分类与分析 ?质量损失统计/分析方法及结果应 用 提问及回答老师问题 第三天9：00~16：30 ●第6讲:质量问题快速根除 ●线索生成 ?如何利用现场数据迅速发现问 题的时空规律并选用最佳解决工具 ?对于长期存在的老大难问题, 如何运用4步拆装法迅速判定问题是 由哪些零件引起的，还是装配过程中 的失误造成的 ?对于零件制造过程中的问题, 如何快速问题是由哪些过程参数造成 的,还是材料或加工方法造成的 ?对于铸造/热处理/焊接/喷涂/玻 璃制造等特殊过程，如何锁定重要过 程参数以及最佳作业条件(水平) 分组讨论2 ●线索确认 ?经过前轮筛选，无论存在多少 可疑因子，如何快速确定问题的真正 原因以及最佳解决方法 ?当问题尚处于开发(样品或试 生产)阶段，如何从根源上有效解决 ●效果验证 ?如何用最小样本及臵信度验证 问题的真因及改善效果 分组练习3 ●第7讲:现场质量管理经验库 ●经验库的建立 ?经验库是一切改善的起点和终 一、中国科学院数学与系统科学研究院简介 中国科学院数学与系统科学研究院由中科院数学研究所、应用数学研究所、系统科学研究所及计算数学与科学工程计算研究所四个研究所整合而成，此外还拥有科学与工程计算国家重点实验室、中科院管理决策与信息系统重点实验室、中科院系统控制重点实验室、中科院数学机械化重点实验室、华罗庚数学重点实验室、随机复杂结构与数据科学重点实验室，以及中科院晨兴数学中心和中科院预测科学研究中心等。2010年11月成立国家数学与交叉科学中心，旨在从国家层面搭建一个数学与其它学科交叉合作的高水平研究平台。数学与系统科学研究院拥有完整的学科布局，研究领域涵盖了数学与系统科学的主要研究方向。共有16个硕士点和13个博士点（二级学科），分布在经济学、数学、系统科学、统计学、计算机科学与技术、管理科学与工程六个一级学科中，可以在此范围内招收和培养硕士与博士研究生。在2006年全国学科评估中，我院数学学科的整体评估得分为本学科的最高分数。数学与系统科学研究院硕士招生类别为硕士研究生、硕博连读生和专业学位硕士研究生。2019年共计划招收122名。 二、中国科学院大学概率论与数理统计专业招生情况、考试科目 三、中国科学院大学概率论与数理统计专业分数线2018年硕士研究生招生复试分数线 2017年硕士研究生招生复试分数线 四、中国科学院大学概率论与数理统计专业考研参考书目 616数学分析 现行（公开发行）综合性大学（师范大学）数学系用数学分析教程。 801高等代数 [1] 北京大学编《高等代数》，高等教育出版社，1978年3月第1版，2003年7月第3版，2003年9月第2次印刷. [2] 复旦大学蒋尔雄等编《线性代数》，人民教育出版社，1988. [3] 张禾瑞，郝鈵新，《高等代数》，高等教育出版社， 1997. 五、中国科学院大学概率论与数理统计专业复试原则 在中国科学院数学与系统科学研究院招生工作小组领导下，按研究所成立招收硕士研究生复试小组，设组长1人、秘书1人。 复试总成绩按百分制计算，其中专业知识成绩占60％，英语听力及口语测试成绩占20％，综合素质成绩占20%。 在面试环节，每位考生有5分钟自述，考查内容主要包括专业知识、外语（口语）水平和综合素质等。 1、专业知识面试重点考查考生对专业基础知识掌握的深度和广度，对知识灵活运用的程度以及考生的实验技能和实际动手能力等，了解考生从事科研工作的潜力和创新能力。 2、外语面试主要考查考生的听、说能力及语言运用能力。 3、思想品德的面试包括考生的政治态度、思想品德、工作学习态度、团队合作精神、科研道德、遵纪守法以及心理素质等内容。 4、体检主要了解考生的身体健康状况，也包括体能、体质和心理素质等。 5、研究生部通过“政审表”向考生所在单位的人事、政工或考生管理部门了解考生的思想品德情况和现实表现。“政审表”将根据中国科学院大学时间部署与调档函一并寄发，需由考生本人档案所在单位的人事（政工）部门加盖公章，随档案一并寄回。政审合格方可寄发录取通知书。 六、中国科学院大学概率论与数理统计专业录取原则 复试小组对本学科参加复试的考生根据初试成绩和复试成绩的综合评定，得出拟录取考生名单，经数学与系统科学研究院招生工作领导小组审核通过。 最终录取成绩：将考生初试成绩和复试成绩按一定比例加权平均后，得出录取成绩。加权平均采用下列公式： 录取成绩＝(初试成绩÷5)×40%＋复试成绩×60%。复试成绩不合格者不予录取；政审不合格、体检不合格者不予录取。 拟录取名单确定后将在网站上公示10个工作日 七、中国科学院大学概率论与数理统计考研复习建议 1、零基础复习阶段（6月前) 矩阵分析及其应用 3.1 矩阵序列 定义3.1 设矩阵序列{A (k)},其中A (k)=() (k ij a )∈C m ?n ,当k →∞， )(k ij a →a ij 时，称矩阵序列{A (k)}收敛，并称矩阵A=(a ij )为矩 阵序列{A (k)}的极限，或称{A (k)}收敛于A, 记为 A A k k =∞ →)(lim 或 A (k)→ A 不收敛的矩阵序列称为发散的。 由定义，矩阵序列A (k) 发散的充要条件为存在ij 使 得数列) (k ij a 发散。 类似地，我们可以定义矩阵收敛的Cauchy 定义 定义3.1' 矩阵序列{A (k)}收敛的充要条件为 对任给ε>0 存在N(ε), 当 k , l ≥ N(ε) 时有 ||A (k)-A (l )|| < ε 其中||.||为任意的广义矩阵范数。 例1 ???? ? ? ??- =∑=-n k n n k k e n n 12) ()sin()1sin(11A 如果直接按定义我们因为求不出A (n )的极限从而 很难应用定义3.1证明收敛。 相反，由于∑∑∑+=+=+=-≤≤n m k n m k n m k k k k k k 112 1 2 ) 1(1 1 ) sin( < 1/m 从而只要l 充分大，则当m, n > l 时就有 ε≤∑ +=n m k k k 1 2 ) sin( 这样A (l ) 收敛。 定理3.1 A (k)→ A 的充要条件为 ||A (k) -A||→0 证明：利用广义矩阵范数的等价性定理，仅对∞范数可以证明。 即 c 1 ||A (k) -A||∞ ≤ ||A (k) -A||≤ c 2 ||A (k) -A||∞ 性质0 若A (k)→ A ， 则 ||A (k)|| → ||A|| 成立。 中科院矩阵分析与应用大作业 实现LU分解QR分解Householder reduction、Givens reduction Matlab 代码： function [] =j uzhendazuoye A=input ('请输入?个矩阵A='); 2 Gram-Schmidt 分解 3 Householder reduction 4 x=input (*请输入序号1 LU分解 Givens reduction: 1); if (>:==!) 壮mmm分解mm%% disp('PA=LU1) m=size(A,1); %nt等于矩阵A的行数 n=size(A,2); %n等于矩阵A的列数 if (m==n) % 刊斯NA是不足方阵 % 如果矩阵A不是方阵那么就输出"error" U=A; %把矩阵至賦值给矩阵u L=zeros(n); %先将L设为单位阵 P=eye(n); %首先将交换矩阵P设为单位矩阵 for j =1:n-1 for i=j +1:n if (U(j, j)-=0) %判断主元元素是否不为0 L(i z j)=U(i z j) /U(j z j)； U(i f :)=U(i, :)-U(j, j)/U(j z j); % U(j, j)为主元元素 else a=j+l；% 令 a 等于j + 1 while ( (U (a, j ) ==0) && (a 第4章 矩阵分解与表示 (I)高斯消去法 假设矩阵A 的顺序主子式i D ≠0 (i=1,…,n-1), 则我们可以进行以下的顺序消元过程 1．消元过程 n k k i b m b b n k k j i a m a a k k ik k i k i k kj ik k ij k ij ,,2,1,,,2,1,,) () () 1()()()1( ++=-=++=-=++ 等价于用初等矩阵T k k k e l I L -=分别 左乘)(k A 和)(k b ，即 )()1(k k k A L A =+ (1) 其中，T k n k k k k k m m m l ),,,,0,,0(,,2,1 ++=, n k i a a m k kk k ik ik ,,1,/)()( +== 我们称ik m 为消元因子，)(k kk a 为主元素； 消元过程的一个重要性质是：消元过程不改 变矩阵的顺序主子矩阵的行列式(顺序主子式)的值。 例 ???? ??????---=012131121A ,顺序主子式为，1，5，-10 ???? ??????--?????→?++250050121)1*(2)3(),1()2(,顺序主子式为，1，5，-10 ???? ??????--??→?-200050121)2()3(，顺序主子式为，1，5，-10 引理：约化的主元素)(i ii a ≠0的充要条件是 矩阵A 的顺序主子式i D ≠0 (i=1,…,k); 推论：若矩阵A 的顺序主子式i D ≠0 (i=1,…,k)，则 1)1(11D a =，k i D D a i i i ii ,,2,1,/1)( ==-； 由此有若A 对称正定或严格对角占优，而 它们的顺序主子矩阵也是对称正定或严格 对角占优，从而顺序主子式不为0，顺序高斯 消去过程可进行； 2.回代过程： ()() ()()()1/()/, 1,2,,1n n n n nn n k k k k k kj kk j k x b a x b a a k n n =+?=??=-???=--?∑ 二、经验类 [quote]1：考中科院科大完全攻略！ 普物类 力学科大出版社杨维宏很好的教材 电磁学高教社赵凯划经典教材（科大出版社的也不错） 热学高教社褚圣麟经典教材（科大出版社的也不错）已经出版了对照的习题解答 上述3门是普物a b的考试范围，弄清楚课后习题足够了！ 电动力学郭硕鸿高教社已经出版了对照的习题解答 理论力学高教社已经出版了对照的习题解答 光学赵凯华北大出版社 量子类 量子力学卷１曾谨言科学出版社，最好同时购买习题集的上下册非常好搞清楚就足够了！ 周世勋高教社《量子力学》入门型已经出版了对照的习题解答！ 考科大、中科院的用这些足够了。还有哪些？大家提出我补充。现在资料更新很快，很多抖出了专门的习题集建议大家看最新的，00年以前的老掉牙的东西没什么用处了。 引用 2、各位朋友大家好:也谈中国科大物理辅导班笔记，物理教材！ 我是科大研究生想告诉大家，不要太指望辅导班笔记。 看到不少人受到误导心痛不已，其实复习就是很简单的事情，很多教材的选择也就是基础常见的就足够了， 高教版的基本都是非常经典的还要习题集的选择电磁学力学等太多了，不过建议大家看一些比较新的资料。 老掉牙的就算了n年了，编这些书的老师估计早就退休了！ 下面几个常见问题: 中国科大物理辅导班笔记，物理教材！（我觉得这个帖子很好） 1 辅导班何时开办？ 每年的11月中旬，到12月20左右出来！ 1 考科大用什么教材？ 其实这个问题很简单了，当然最好是科大教材了，如果是科大习题集最好了，现在科大教材变化很快毫无疑问最好的教材就是最新的。多少年来变化很大的，但是科大教材不是好教材，力学其实复旦的比较好，科大yangweihong的觉得很一般，不过习题不错。电磁学毫无疑问是高教社的zhaokaihua的好啊，科大张玉民的也是很一般的教材。原子物理也是推荐高教社chushe nglin的很经典的教材。但是教材归教材，习题集最好还是选择科大这个道理很简单了 1 为什么考科大物理？ 2科大物理国内一流国际闻名科大全公费住宿免费补助待遇每月500以上设备先进值得你去努力 2 外校能否报名？ 不能，就是科大校内的学生也要凭借学生证，不是科大物理系的就很难接受。 3 辅导班笔记含金量多高？ 辅导班笔记其实就是串讲班不叫辅导班，所以就是科大物理各门的大复习。知识点几乎面面据到！ 4 市场上的辅导班笔记可信么？ 这个我觉得还是大家自己判断为好。你相信外部人有么？自己决定！ 5 给你辅导班笔记怎么判断真假？ 首先要考虑对方可能会有么？如果可能有，对比一下是不是往年的笔记可信度多大？科大官方部不提供这个咨询服务。 6 如果没有辅导班笔记怎么复习？ 扎扎实实的复习力学电磁学原子物理量子力学建议使用科大版本教材，道理很简单了。其实很多其他教材也不错，科大的很多教材很差不想想你想像的那么好！ 引用 3、关于中科大中科院量子力学和普通物理考研试题的若干说明 热烈欢迎2008年考中科大中科院的同学们！！！ 2008年考研的要提前准备才充分！！！ 中国科学院的一些招生单位（包括物理所和高能所在内），在06年研究生入学考试抛弃了以前科大的命题，改由中科院研究生院命题。实际上就是由以前中科大的老师出题，变为中科院的研究机构的那些导师出题。（据了解，一些导师接到出题任务都很烦，因为科研压力大啊，出题就让自己带的研究生随便在习题集上 矩阵分析及其应用 3.1矩阵序列 定义3.1设矩阵序列｛应)｝,其中A?)=(#))￡Cms,当k—oo, 佝时，称矩阵序列｛A00｝收敛，并称矩阵A=(佝)为矩阵序列｛A00｝的极限,或称｛A00｝收敛于A,记为lim A a)= A或A,k)-> A ks 不收敛的矩阵序列称为发散的。 由定义，矩阵序列A(k)发散的充要条件为存在ij使 得数列站发散。 类似地，我们可以定义矩阵收敛的Cauchy定义 定义31矩阵序列｛A00｝收敛的充要条件为 对任给￡>0存在N(E),当k,l> N(E)时有 IIA(k)-A(/)ll < ￡ 其中11.11为任意的广义矩阵范数。 例 1 A(n) e~n sin(-) n y，sin(R) k=l K 7 如果直接按定义我们因为求不出A㈤的极限从 而很难应用定义3.1证明收敛。 相反，由于t^< t^< v 1/m 从而只要/充分大，则当m, n > /时就有 n z sin(A) 这样A")收 定理3.1 A(k)->A的充要条件为 HA'10-AII T O 证明：利用广义矩阵范数的等价性定理，仅对co范数可以证明。 即ci IIA(k) -AIL < IIA(k) -All< c2 IIA(k) -AIL 性质 1.设A(k,—> A mxn, B,k,—> B mxn>则 a- A(k)+P ? B(k) -> a- A+P B, V a,PeC 性质2.设A(k)—> A mxn, B,k)—> B nx/,则 A(k)由如一A B 证明：由于矩阵范数地等价性，我们E以只讨论相容的 矩阵范数。 IIA(k).B(k)-A-BII < II A(k) -B(k) -A-B(k)ll+IIAB(k)- A-BII 中科院矩阵分析与应用大作业 按照上面建立函数文件，脚本文件。名字不要随便改，要改的话要同内部函数名一起修改。下面放代码，总共七段代码，第一个是主程序进行人机交互，其余为接口函数实现矩阵的各个分解功能。另外容我小小的收点下载券，毕竟花了大概8小时左右。 %%矩阵分解：LU分解；QR分解(schmidt)；Householder；Givens %%2016.11.22 clc,clear all; Matrix=input('请输入需要进行分解的矩阵：'); disp('请选择需要对矩阵进行的操作:'); disp('1.LU分解');disp('2.QR分解'); disp('3.Householder约简');disp('4.Givens约简'); n=input(''); clc;%清除命令窗口内容 Matrix%显示输入矩阵 %%分解操作 switch n case1 LUfactor(Matrix);%LU分解 case2 QR(Matrix);%QR分解 case3 houseHolder(Matrix);%householder约简case4 Givens(Matrix);%Givens约简 end function LUfactor(A) [m n]=size(A); %判断矩阵是否为方阵 if m~=n error('矩阵非方阵，请重新输入！'); end %判断矩阵是否奇异 if det(A)==0 error('矩阵奇异，请重新输入！'); end %判断矩阵顺序主子式是否全不为0 n=size(A,1); flagMat=zeros(n,1); for i=1:n if det(A(1:i,1:i))==0 flagMat(i)=1; end end %以顺序主子式是否含0来决定采用的LU分解方式if any(flagMat)==0 disp('顺序主子式均不为0，采用A=LU'); disp('***LU分解结果***') [L U]=LUFull(A) else disp('顺序主子式存在0，采用PA=LU'); disp('***LU分解结果***') [L U P]=LUPart(A) end 第五章 特征值的估计及对称矩阵的极性 本章主要讨论数值代数中的三个特殊理论， 即 特征值的估计 广义特征值问题 实对称矩阵（一般是Hermite 矩阵）特征值的 极小极大原理，其次也涉及到一些特征值 和奇异值的扰动问题,最后简要地介绍矩阵 直积的一些性质及其在线性矩阵方程求解 方面的应用。这几方面的内容，在矩阵的 理论研究与实际应用当中都有着相当重要 的作用。 5.1特征值的估计 一、特征值的界 首先给出直接估计矩阵特征值模的上界的 一些方法 定理5.1 设A=(a rs )∈R n×n ,令 M=||2 1 max ,1sr rs n s r a a -≤≤ λ若表示A 任一特征值，则λ的虚部Im(λ) 满足不等式 2 ) 1(|)Im(|-≤n n M λ |Im(λ)|≤||A -A T ||2 / 2 |Im(λ)|≤||A -A T ||1 ?/2. 证明：设x+i ?y 为对应于λ的A 的特征向量， 则 A(x+i ?y)=(α+β?i)(x+i ?y) 其中λ=α+β?i.显然x,y 为实向量，且x,y 为 线性无关的 向量。 经整理A(x,y)=(x,y)B, 其中B=??? ? ??-αββα 。 从而(x,y)T A(x,y)=(x,y)T (x,y)B 展开有 ???? ??Ay y Ax y Ay x Ax x T T T T =α????? ??y y y x y x x x T T T T + β???? ? ? ?--x y y y x x y x T T T T (求等式两边矩阵的对角元之和，可得 α(x T x +y T y )=x T Ax +y T Ay (1) 等式两边矩阵的左上角单元减去右下角单元 可得： β(x T x +y T y )=x T (A -A T )y 1). 记B=A -A T ,则 |x T By|≤||x||2 ?||B||2?||y||2 从而 |β|≤||x||2 ?||B||2?||y||2 /((||x ||2)2 +(||y ||2)2) 利用ab /(a 2+b 2)≤1/2 可得 |β|≤||B||2 /2. 2). 由于|x T By|≤||Bx||1 ?||y||∞≤||B||1?||x||1 ?||y||∞ 从而 |β|≤||B||1 ?||x||1 ?||y||∞ /((||x ||2)2 +(||y ||2)2) 易证明 ||x||1 ?||y||∞ /((||x ||2)2 +(||y ||2)2) /2. (显然，不妨假设(||x ||2)2 +(||y ||2)2=1, 设||y ||∞=t =cos(α), 则y 必为t ? e j 的形式（为什么？）， 从而极值转化为求解如下最大值问题： max ||x||1, 满足约束(||x ||2)2=1-t 2 这样有均值不等式||x||1 x ||2 = -t 2)1/2, 从而我们需要求解t (1-t 2)1/2的最大值，设t =cos(α) 可得t (1-t 2)1/2的最大值为1/2. 从而得证。) 因此 |β|≤||B||1 3). 由于b ii =0, i =1,2,…,n , b ij = -b ji , 因此 |x T By|2=| 1 1()n ij i j j i i j i b x y x y -=>??-∑∑|2 ≤(2M )2 2 1||n i j j i i j i x y x y =>??- ??? ∑∑ (利用(a 1+a 2+…+a n )2≤ n ((a 1)2+(a 2)2+…+(a n )2) ≤(2M )2 (n (n -1)/2) 21||n i j j i i j i x y x y =>??- ??? ∑∑ 第 2 章范数理论及其应用 2.1向量范数及I p范数 定义：如果V 是数域K 上的线性空间，且对于 V的任一向量x，对应一个实数值ixil,它满足以下三个条件： 1）非负性：||x|| 0,且||x||=0 x=0; 2）齐次性：iikxii=iki iixii，k K； 3）三角不等式：||x+y|| ||x||+||y||. 则称||x|为V上向量x的范数，简称为向量范数。 可以看出范数||||为将V映射为非负数的函数。注意：2）中|k|当K为实数时为绝对值， 当K 为复数域时为复数的模。 虽然向量范数是定义在一般的线性空间上的，但是由于前面的讨论，我们知道任何n 维线性空间在一个基下都代数同构于常用的n维复（或实）列向量空间, 因此下面我们仅仅讨论n 维复（或实）列向量空间就足够了下面讨论如下：1?设||||为线性空间V n的范数，任取它的一个 基X i,X2,…,X n,则对于任意向量X,它可以表示为 x= 1X1+ 2X2+ …+ n X n 其中，（1, 2,…,n）T为X的坐标。 由此定义C n（或R n）中的范数如下： || ||C = （） = || 1X1+ 2X2+ …+ n X n|| 则容易验证|| ||C确实为C n中的范数. 2?反之，若|| |C为C n中的范数，定义V n的范数如下：||X||= （X）=|| ||c 其中X= 1X1+ 2X2+ …+ n X n。 则容易验证（X）确实为V n的范数。 这个例子充分说明了一般线性空间的范数和n维 复（或实）列向量空间的范数之间的关系。这也是为我们只讨论n 维复（或实）列向量空间的范数的理由. 范数首先是一个函数，它将线性空间的任意向量映射为非负实数。 范数与函数 性质 1. 范数是凸函数， 即|| （1 ）X+ y|| （1 ）||X||+ ||y|| 其中0 第五章特征值的估计及对称矩阵的极性本章主要讨论数值代数中的三个特殊理论，即 特征值的估计 广义特征值问题 实对称矩阵(一般是Hermite矩阵)特征值的极小极大原理，其次也涉及到一些特征值和奇异值的扰动问题，最后简要地介绍矩阵直积的一些性质及其在线性矩阵方程求解方面的应用。这几方面的内容，在矩阵的理论研究与实际应用当中都有着相当重要的作用。 5.1特征值的估计 一、特征值的界 首先给出直接估计矩阵特征值模的上界的 一些方法 定理 5.1 设A=(a rs) R n X1,令 1 , , M= ma彷总a sr| 若表示A任一特征值，则的虚部Im() 满足不等式 |Im( )| M n(n21) |Im( )| ||A A T||2 / 2 |Im( )| ||A A T||1n /2. 证明：设x+i y为对应于的A的特征向量, 则A(x+i y)=( + i)(x+i y) 其中=+ i.显然x,y为实向量，且x,y为线性无关的向量。 经整理A(x,y)=(x,y)B, 其中B= 从而(x,y) T A(x,y)=(x,y) T(x,y)B 展开有 i 1 j i T T X y X X T T y y y X (求等式两边矩阵的对角元之和，可得 (x T x+y T y)=x T Ax+y T Ay (1) 等式两边矩阵的左上角单元减去右下角单元 可得： (x T x+y T y)=x T (A A T )y 1) . 记 B=A A T ,则 |x T By| ||x||2||B||2||y||2 从而 1 1 1凶|2 ||B||2||y||2 /((||x||2)2 +(||y|2)2) 利用 ab/(a 2+b 2) 1/2 可得 | | ||B||2 /2. 2) . 由于 |x T By| ||B X ||I ||y|| ||B||i ||X ||I ||y|| 从而 | | ||B||i ||x||i ||y|| /((||X |2)2 +(||y||2)2) 易证明 ||x||i ||y|| /((||X ||2)2 +(||y||2) 2) n /2. (显然，不妨假设(||X ||2)2 +(||y||2)2=1, 设HyH =t=cos ()，则y 必为t e 的形式(为什么？) 从 而极值转化为求解如下最大值问题： max ||X ||1,满足约束(||X ||2)2=1 t 2 这样有均值不等式 ||x|h i n ||X ||2= 、、n (1 t 2)1/2, 从而我们需要求解t(1 t 2)1/2的最大值，设t=cos() 可得 t(1 t 2)1/2的最大值为1/2.从而得证。) 因此 11 ||B||1 . n /2. 3) . 由于 b ii =0, i =1,2,…,n, b ij = b ji , n 1 因此 x T By|2=| b ij (X y j X j y i )|2 i 1 j i 2 n (2M)2 |xy j X j Y i | i 1 j i (利用(a 1+a 2+…+a n )2 n((a 1)2+(a 2)2+ …+(a n )2) n (2M)2(n(n 1)/2) | X y j X j yj 2 X T A X y T Ax X T Ay y T Ay T T X X X y T T X y y y clc,clear,close all; A=input('请输入需要进行LU分解的方阵：'); [m n]=size(A); %判断矩阵是否为方阵 if m~=n error('矩阵非方阵，请重新输入！'); end %判断矩阵是否奇异 if det(A)==0 error('矩阵奇异，请重新输入！'); end %判断矩阵顺序主子式是否全不为0 n=size(A,1); flagMat=zeros(n,1); for i=1:n if det(A(1:i,1:i))==0 flagMat(i)=1; end end %以顺序主子式是否含0来决定采用的LU分解方式if any(flagMat)==0 disp('顺序主子式均不为0，采用A=LU'); [L U]=LUFull(A) else disp('顺序主子式存在0，采用PA=LU'); [L U P]=LUPart(A) end %对矩阵A进行LU分解(完全主元法) %L,U矩阵为输出变量；A矩阵为输入变量 U=zeros(size(A));%对U尽可能初始化 U(1,:)=A(1,:); L=eye(size(A));%对L尽可能初始化 L(2:end,1)=A(2:end,1)/A(1,1); n=size(A,1); for i=2:n%Dolittle公式 for j=i:n U(i,j)=A(i,j)-L(i,1:i-1)*U(1:i-1,j); end for k=i:n L(k,i)=(A(k,i)-L(k,1:i-1)*U(1:i-1,i))/U(i,i); end end %对矩阵A进行LU分解(部分主元法) %L,U,P矩阵为输出变量，A矩阵为输入变量 %对输入矩阵A重构 n=size(A,1); B=zeros(n,1); for i=1:n B(i,1)=B(i,1)+i; end A=[A B]; %目标矩阵 tempMat=zeros(size(A)); for j=1:n%第j列 [~,row]=max(abs(A(j:end,j)));%找出第j列最大元素,返回所在行A([row+j-1j],:)=A([j row+j-1],:);%交换最大元素行与第j行 tempMat(j,:)=A(j,:);%将第j行元素复制到tempMat第j行for i=j:n%第i行 if i+1<=n tempMat(i+1,j)=A(i+1,j)/A(j,j); end end %生成新的A(高斯消去) for k=j:n%第k行 if k+1<=n ratio=-A(k+1,j)/A(j,j); A(k+1,j+1:n)=ratio*A(j,j+1:n)+A(k+1,j+1:n); end end end %提取U矩阵 U=triu(tempMat(:,1:n),0);%upper triangle %提取L矩阵 L=tril(tempMat(:,1:n),0);%lower triangle for i=2:n for j=1:(i-1) L(i,j)=L(i,j)./L(j,j); end end 本人本科不是985、211名校，2014年考研，考中国科学院大学计算数学，总成绩386分。【关于数学分析】： 大一大二时，我对老师的数学分析课很感兴趣，所以自认为学得比较认真，理解得也比较深刻。这虽然是我学得最认真的一门课，但也只仅限于完成平时的作业和完成期末的复习，并没有做很多额外的努力，这样四学期下来，对数学分析整体还是有个宏观的把握，比如说起微积分学基本定理，Dini定理，欧拉积分等都会有个大致的印象。 大三时，我旁听了老师给下一级上的数学分析课和数学分析选论，虽然已经是第二次听，但由于以前学得都没有及时复习，老师讲的很多忘记了，很多知识他讲到才想起来。同时，大三这一学年我都在自己钻研一些自己比较感兴趣的数学分析内容，比如反常积分敛散性、级数求和等，因此顺便做了数学分析的课题。但由于大三还是有专业课（实变、抽代等）和其他课程的压力，一直没有系统地开展复习。大三结束，从7月5日起，正式进入考研复习阶段，一开始，我打算做《裴礼文》，计划每天做20页，但是做了1个月不到，发现效率很低，因为对于我来说里面大部分题目都不会，每道题都要花时间。到了后来，每天看10页都已经很不错了，因为之前几天看懂的题目又忘了，又需要不断复习。这时我感到了压力很大，觉得应该以应试为主，《裴礼文》上的题目偏难了，有些大大超过考研难度。 于是，从8月开始，我开始做《题解精粹（钱吉林）》，这本书里面有1000多道各地的考研题，而且对于我来说，有简单题，有少量难题，总的来说比较适合我。所以我每天做40道，这样做了一整个8月，刚好差不多做完。这时，我深深感到之前数学分析课堂上学的真的很不够。从9月份开始，我开始做教材的课后习题，顺便复习一下教材里面的定理。这时做课后题还是需要借助答案，然后发现原来教材后面的题目许多都是考研题，怪不得以前都不会做。这一轮复习到了10月中旬。10月中旬我又把《题解精粹（钱吉林）》在第一遍看时没看透彻或是有跳过的地方补上。 11月，我专门去帝都上了新祥旭的专业课冲刺，主要是老师带着讲解中科院的历年数分真题和成套成套地做与讲模拟题，基本上用一个下午的时间大致做一年的数分真题或模拟卷，再讲解答案解析。到11月中下旬，就把中科院数分的历年真题都研究透彻了（包括2013年的，虽然这年数分题出的相当难）。因为在12月，英语政治复习的压力也很大，就没有时间和精力再待了，总结一点就是感觉收获颇大。这样，数学分析考了128，虽然不高，但也还过得去。 【关于高等代数】： 我在大二学高等代数的时候，因为没有数分的那种喜欢，所以一直使不上劲，上册的内容勉强学懂了，下册的内容在学的时候，到线性变换的时候完全就不懂（一点不谦虚）。虽然期末考得还好，但含水量极大。后来整个大三几乎就没碰过高代。到大三结束7月5号开始做《高等代数考研教案》，每天做10页，觉得这本挺适合我，内容也很全。因为后来到了正定二次型和线性变换那边，本来基础就不好，边看边做，边跳边做，到8月底勉强把它做完，这时终于对整个高代内容有了大概的把握，因为把“矩阵”那章也做了，这章一般上课是不讲的，但考研还挺重要的。 到9月份，开始看课本（真的是看课本，很多以前上课的时候都没看过，例如线性变换、矩阵中的内容）。然后把课后习题做一遍，但补充题基本还是要借助答案。这样做到了10月中旬（国庆放了个假），这时开始又把《考研教案》里曾错过，曾不理解，或曾跳过的内容补上。到了10月，我开始看中科院的考研高代真题，基本上用一个下午的时间大致做一年的高代真题，再看答案解析。到10月中下旬，差不多把高代的真题都研究透彻了（例如2013年的高代还是挺难的）。这样，高等代数考了135，应该算不错。 【关于英语】： 我高中时英语还是挺好的，不过经过大学初的几年，就荒废了，大二上靠着高中时残存中科院矩阵分析与应用大作业
中科院-模式识别考题总结(详细答案)
中科院矩阵分析chapt3
现场质量管理与突破性快速改善
2020-2021年中国科学院大学(中科院)概率论与数理统计考研招生情况、分数线、参考书目及备考经验
中科院矩阵分析课件
中科院矩阵分析与应用大作业
中科院矩阵分析chapt4
中科院-科大真题最完整版+考试攻略
中科院矩阵分析课件.doc
中科院矩阵分析与应用大作业
中科院矩阵分析_第五章
中科院矩阵分析_第二章
中科院矩阵分析_第五章
LU分解(中科院矩阵分析大作业)
中科院计算数学考研必读经验
相关主题
文本预览