下载文档原格式
工程认证背景下基于OBE 理念的工程数学教学改革与探索郭凯红 刘婷婷*(辽宁大学信息学院 辽宁沈阳 110036)摘要： 针对工程数学课程的性质与特点，以探究式教学为主线，同时融合混合式、案例驱动式等教学方法开展教学。
利用先进教育理论联系工程实践，探索工程认证背景下基于成果导向教育（Outcomes Based Educa‐tion，OBE ）的工程数学课程创新性教学模式，以培养学生的独立思考能力和创新性思维能力为出发点，研究和探讨了工程数学课程教学改革的关键问题与实施方案。
关键词： 成果导向教育 教学模式 工程数学 工程认证中图分类号： G642.0;TB11-4文献标识码： A文章编号： 1672-3791(2023)10-0192-04Reform and Exploration of Engineering Mathematics TeachingBased on the Theory of OBE under the Background ofEngineering AccreditationGUO Kaihong LIU Tingting*(College of Information, Liaoning University, Shenyang, Liaoning Province, 110036 China)Abstract: According to the properties and characteristics of Engineering Mathematics, teaching is carried out with inquiry-based teaching as a main clue, along with the blended, case-driven and other teaching methods. Using ad‐vanced educational theories and engineering practice, this paper explores the innovative teaching model of Engi‐neering Mathematics based on OBE in the context of engineering accreditation, and studies and discusses the key problems and implementation plans of the teaching reform of Engineering Mathematics with the cultivation of stu‐dents' independent thinking ability and innovative thinking ability as the starting point.Key Words: OBE; Teaching model; Engineering mathematics; Engineering accreditation工程数学课程（包括线性代数、概率论与数理统计等）作为工科专业最重要的基础课程之一，根本任务是在教学中培养学生严谨的逻辑思维能力、独立思考能力与创新能力。
*1.【圆函数】e (φ)=cos φi +sin φj .*2.a.弧长的微分ds =以点M 为界,当ds 位于s 增大一方时取正号;反之取负号.b.矢性函数的微分的模,等于(其矢端曲线的)弧微分的绝对值.矢性函数(其矢端曲线的)弧长s 的导数d r /ds 在几何上为一切单位矢量,恒指向s 增大的一方.+3.证明||.ds d d r t dt=证,d dx dy dz dtdt dtr i j k dt =++d dt r =由于ds 与dt 有相同的符号,故有.ds d dt dt r ===由此可知:矢端曲线的切向单位矢量.d d ds d d dt dt dt dtd r s r r r ==*4.【二重矢积】公式：a ×(b ×c )=(a ·c )b -(a ·b )c .+5.矢性函数A (t)的模不变的充要条件是.d d A A t•=0证假定|A |=常数,则有A 2=|A |2=常数.两端对t 求导[左端用导数公式],就得到.d d A A t •=0反之,若有.d d A A t •=0则有,d dt A =20从而有A 2=|A |2=常数.所有有|A |=常数.定常矢量A (t)与其导矢相互垂直.*6.''.A B A dt t B B A d ×=×+×∫∫''.A B A dt t B B A d •=•−•∫∫+7.一质点沿曲线r =rcos φi +rsin φj 运动,其中r,φ均为时间t 的函数.求速度v 在矢径方向及其垂直方向上的投影v r 和v φ.解将r 写成r =r e (φ),则有()().d dr d r dt dt v d r e e t ϕϕϕ==+1由此可知：,.r dr d v v r dt dtϕϕ==[使用圆函数e (φ)，则e (φ)及e 1(φ)之方向即为矢径方向及与之垂直的方向.]*8.【矢量线】A =A x i +A y j +A z k 为单值、连续且有一阶连续导数。
第二十二章 曲面积分4 场论初步一、场的概念概念：若对全空间或其中某一区域V 中每一点M ，都有一个数量(或向量)与之对应，则称V 上给定了一个数量场(或向量场).温度场和密度场都是数量场. 若数量函数u(x,y,z)的偏导数不同时为0, 则满足方程u(x,y,z)=c(常数)的所有点通常是一个曲面.曲面上函数u 都取同一个值时，称为等值面，如温度场中的等温面.重力场和速度场都是向量场. 设向量函数A(x,y,z)在三坐标轴上投影分别为：P(x,y,z), Q(x,y,z), R(x,y,z), 则A(x,y,z)=(P(x,y,z), Q(x,y,z), R(x,y,z)), 其中P , Q, R 为定义区域上的数量函数，且有连续偏导数.设向量场中的曲线L 上每点M 处的切线方向都与向量函数A 在该点的方向一致，即P dx =Q dy =Rdz, 则称曲线L 为向量场A 的向量场线. 如, 电力线、磁力线等都是向量场线.二、梯度场概念：梯度是由数量函数u(x,y,z)定义的向量函数grad u=⎪⎪⎭⎫⎝⎛∂∂∂∂∂∂z u y u x u ,,, 且grad u 的方向是使lu∂∂达到最大值的方向, 其大小为u 在这个方向上的方向导数. 所以可定义数量场u 在点M 处的梯度grad u 为在M 处最大的方向导数的方向，及大小为在M 处最大方向导数值的向量. 因为方向导数的定义与坐标系的选取无关，所以梯度定义也与坐标系选取无关. 由梯度给出的向量场，称为梯度场. 又数量场u(x,y,z)的等值面u(x,y,z)=c 的法线方向为⎪⎪⎭⎫⎝⎛∂∂∂∂∂∂z u y u x u ,,, 所以 grad u 的方向与等值面正交, 即等值面法线方向. 引进符号向量： ▽=⎪⎪⎭⎫⎝⎛∂∂∂∂∂∂z y x ,,. 将之视为运算符号时, grad u=▽u.基本性质：若u,v 是数量函数, 则 1、▽(u+v)=▽u+▽v ；2、▽(uv)=u(▽v)+(▽u)v. 特别地▽u 2=2u(▽u)；3、若r=(x,y,z), φ=φ(x,y,z), 则d φ=dr ▽φ；4、若f=f(u), u=u(x,y,z), 则▽f=f ’(u)▽u ；5、若f=f(u 1,u 2,…,u n ), u i =u i (x,y,z) (i=1,2,…,n), 则▽f=i ni iu u f∑=∇∂∂1. 证：1、▽(u+v)=⎪⎪⎭⎫ ⎝⎛∂+∂∂+∂∂+∂z v u y v u x v u )(,)(,)(=⎪⎪⎭⎫⎝⎛∂∂+∂∂∂∂+∂∂∂∂+∂∂z v z u y v y u x v x u ,, =⎪⎪⎭⎫⎝⎛∂∂∂∂∂∂z u y u x u ,,+⎪⎪⎭⎫⎝⎛∂∂∂∂∂∂z v y v x v ,,=▽u+▽v. 2、▽(uv)=⎪⎪⎭⎫ ⎝⎛∂∂∂∂∂∂z uv y uv x uv )(,)(,)(=⎪⎪⎭⎫⎝⎛∂∂+∂∂∂∂+∂∂∂∂+∂∂z v u v z u y v u v y u x v u v x u ,,=⎪⎪⎭⎫ ⎝⎛∂∂∂∂∂∂z v u y v u x v u,,+⎪⎪⎭⎫ ⎝⎛∂∂∂∂∂∂v z u v y u v x u ,,=u ⎪⎪⎭⎫ ⎝⎛∂∂∂∂∂∂z v y v x v ,,+⎪⎪⎭⎫⎝⎛∂∂∂∂∂∂z u y u x u ,,v=u(▽v)+(▽u)v. 当u=v 时，有▽u 2=▽(uv)=u(▽v)+(▽u)v =2u(▽u).3、∵dr=dx+dy+dz, ▽φ=⎪⎪⎭⎫⎝⎛∂∂∂∂∂∂z y x ϕϕϕ,,, ∴dr ▽φ=(dx+dy+dz)⎪⎪⎭⎫ ⎝⎛∂∂∂∂∂∂z y x ϕϕϕ,,=dz z dy y dx x ∂∂+∂∂+∂∂ϕϕϕ=d φ. 4、∵▽f=⎪⎪⎭⎫ ⎝⎛∂∂∂∂∂∂z f y f x f ,,=⎪⎪⎭⎫⎝⎛∂∂∂∂∂∂z u du df y u du df x u du df ,,, 又▽u=⎪⎪⎭⎫⎝⎛∂∂∂∂∂∂z u y u x u ,,, f ’(u)=du df, ∴f ’(u)▽u=⎪⎪⎭⎫ ⎝⎛∂∂∂∂∂∂z u y u x u du df ,,=⎪⎪⎭⎫ ⎝⎛∂∂∂∂∂∂z u du df y u du df x u du df ,,=▽f. 5、▽f =⎪⎪⎭⎫ ⎝⎛∂∂∂∂∂∂z f y f x f ,,=⎪⎪⎭⎫⎝⎛⎪⎪⎭⎫ ⎝⎛∂∂∂∂⎪⎪⎭⎫ ⎝⎛∂∂∂∂⎪⎪⎭⎫ ⎝⎛∂∂∂∂∑∑∑===n i i i n i i i n i i i z u u f y u u f x u u f 111,,=∑=⎪⎪⎭⎫ ⎝⎛∂∂∂∂∂∂∂∂∂∂∂∂ni i i i i i i z u u f y u u f x u u f 1,,=∑=⎪⎪⎭⎫ ⎝⎛∂∂∂∂∂∂∂∂n i i i i iz u y u x u u f1,,=i n i iu u f∑=∇∂∂1.例1：设质量为m 的质点位于原点, 质量为1的质点位于M(x,y,z), 记OM=r=222z y x ++, 求rm的梯度. 解：rm∇=⎪⎭⎫ ⎝⎛-r z r y r x r m ,,2.注：若以r 0表示OM 上的单位向量，则有r m∇=02r rm -, 表示两质点间引力方向朝着原点, 大小是与质量的乘积成正比, 与两点间的距离的平方成反比. 这说明引力场是数量函数r m 的梯度场. 所以称rm为引力势.三、散度场概念：设A(x,y,z)=(P(x,y,z), Q(x,y,z), R(x,y,z))为空间区域V 上的向量函数, 对V 上每一点(x,y,z), 定义数量函数D(x,y,z)=zRy Q x P ∂∂+∂∂+∂∂, 则 称D 为向量函数A 在(x,y,z)处的散度，记作D(x,y,z)=div A(x,y,z).设n 0=(cos α, cos β, cos γ)为曲面的单位法向量, 则=n 0dS 就称为曲面的面积元素向量. 于是得高斯公式的向量形式：⎰⎰⎰VdivAdV =⎰⎰⋅SdS A .在V 中任取一点M 0, 对⎰⎰⎰VdivAdV 应用中值定理,得⎰⎰⎰VdivAdV =div A(M*)·△V=⎰⎰⋅SdS A , 其中M*为V 中某一点，于是有div A(M*)=VdSA S∆⋅⎰⎰. 令V 收缩到点M 0(记为V →M 0) 则M*→M 0, 因此div A(M 0)=VdSA SM V ∆⋅⎰⎰→0lim.因⎰⎰⋅SdS A 和△V 都与坐标系选取无关，所以散度与坐标系选取无关.由向量场A 的散度div A 构成的数量场，称为散度场.其物理意义：div A(M 0)是流量对体积V 的变化率，并称它为A 在点M 0的流量密度.若div A(M 0)>0, 说明在每一单位时间内有一定数量的流体流出这一点，则称这一点为源.反之，若div A(M 0)<0, 说明流体在这一点被吸收，则称这点为汇. 若向量场A 中每一点皆有div A=0, 则称A 为无源场.向量场A 的散度的向量形式为：div A=▽·A.基本性质：1、若u,v 是向量函数, 则▽·(u+v)=▽·u+▽·v ； 2、若φ是数量函数, F 是向量函数, 则▽·(φF)=φ▽·F+F ·▽φ；3、若φ=φ(x,y,z)是一数量函数, 则▽·▽φ=222222zy x ∂∂+∂∂+∂∂ϕϕϕ.证：1、记u(P 1(x,y,z),Q 1(x,y,z),R 1(x,y,z)), v(P 2(x,y,z),Q 2(x,y,z),R 2(x,y,z)), 则▽·(u+v)=zR R y Q Q x P P ∂+∂+∂+∂+∂+∂)()()(212121 =⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂+⎪⎪⎭⎫⎝⎛∂∂+∂∂+∂∂z R y Q x P z R y Q x P 222111=▽·u+▽·v. 2、▽·(φF)=z R y Q x P ∂∂+∂∂+∂∂)()()(ϕϕϕ=zR z R y Q y Q x P x P ∂∂+∂∂+∂∂+∂∂+∂∂+∂∂ϕϕϕϕϕϕ =φ⎪⎪⎭⎫⎝⎛∂∂+∂∂+∂∂z R y Q x P +(P ,Q,R)⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂z y x ϕϕϕ=φ▽·F+F ·▽φ. 3、∵▽φ=⎪⎪⎭⎫⎝⎛∂∂∂∂∂∂z y x ϕϕϕ,,, ∴▽·▽φ=⎪⎭⎫ ⎝⎛∂∂∂∂+⎪⎪⎭⎫ ⎝⎛∂∂∂∂+⎪⎭⎫ ⎝⎛∂∂∂∂z z y y x x ϕϕϕ=222222zy x ∂∂+∂∂+∂∂ϕϕϕ.注：算符▽的内积▽·▽常记作△=▽·▽=222222zy x ∂∂+∂∂+∂∂，称为拉普拉斯算符, 于是有▽·▽φ=△φ.例2：求例1中引力场F=⎪⎭⎫⎝⎛-r z r y r x r m,,2所产生的散度场.解：∵r 2=x 2+y 2+z 2, ∴F=3222)(z y x m ++-(x,y,z),▽·F=-m ⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛∂∂+⎪⎭⎫ ⎝⎛∂∂+⎪⎭⎫ ⎝⎛∂∂333r z z r y y r x x =0.注：由例2知，引力场内每一点处的散度都为0(除原点处外).四、旋度场概念：设A(x,y,z)=(P(x,y,z),Q(x,y,z),R(x,y,z))为空间区域V 上的向量函数, 对V 上每一点(x,y,z), 定义向量函数F(x,y,z)=⎪⎪⎭⎫⎝⎛∂∂-∂∂∂∂-∂∂∂∂-∂∂y P x Q x R z P z Q y R ,,, 称之为向量函数A 在(x,y,z)处的旋度, 记作rot A.设(cos α,cos β,cos γ)是曲线L 的正向上的单位切线向量t 0的方向余弦, 向量ds =(cos α,cos β,cos γ)ds= t 0dl 称为弧长元素向量. 于是有 斯托克斯公式的向量形式：⎰⎰SdS rotA ·=⎰Lds A ·.向量函数A 的旋度rot A 所定义的向量场，称为旋度场.在流量问题中，称⎰L A ·为沿闭曲线L 的环流量. 表示流速为A 的不可压缩流体在单位时间内沿曲线L 的流体总量，反映了流体沿L 时的旋转强弱程度. 当rot A=0时，沿任意封闭曲线的环流量为0，即流体流动时不成旋涡，这时称向量场A 为无旋场.注：旋度与坐标系的选择无关. 在场V 中任意取一点M 0,通过M 0作平面π垂直于曲面S 的法向量n 0, 且在π上围绕M 0作任一封闭曲线L, 记L 所围区域为D ，则有⎰⎰SrotA ·=⎰⎰DdS n rotA 0·=⎰LA ·. 又由中值定理有 ⎰⎰DdS n rotA 0·=(rotA ·n 0)M*μ(D)=⎰LA ·, 其中 μ(D)为区域D 的面积, M*为D 中的某一点. ∴(rotA ·n 0)M*=)(·D A Lμ⎰.当D 收缩到点M 0(记作D →M 0)时, 有M*→M 0, 即有 (rotA ·n 0)0M =)(·limD A LMD μ⎰→ .左边为rot A 在法线方向上的投影，即为旋度的另一种定义形式. 右边的极限与坐标系的选取无关，所以rot A 与坐标系选取无关.物理意义：⎰⎰DdS n rotA 0·=(rotA ·n 0)M*μ(D)=⎰LA ·, 表明向量场在曲面边界线上的切线投影对弧长的曲线积分等于向量场旋度的法线投影在曲面上对面积的曲面积分. 即流体的速度场的旋度的法线投影在曲面上对面积的曲面积分等于流体在曲面边界上的环流量.刚体旋转问题：设一刚体以角速度ω绕某轴旋转，则角速度向量ω方向沿着旋转轴，其指向与旋转方向的关系符合右手法则，即右手拇指指向角速度ω的方向，其它四指指向旋转方向. 若取定旋转轴上一点O 作为原点，则刚体上任一点P 的线速度v 可表示为v=ω×r, 其中r=OP 是P 的径向量. 设P 的坐标为(x,y,z),便有r=(x,y,z),设ω(ωx ,ωy ,ωz ), ∴v=(ωy z-ωz y,ωz x-ωx z,ωx y-ωy x), ∴rot v=(2ωx ,2ωy ,2ωz )=2ω或ω=21rot v.即线速度向量v 的旋度除去21, 就是旋转的角速度向量ω. 也即 v 的旋度与角速度向量ω成正比.基本性质：rot A=▽×A. 1、若u,v 是向量函数, 则 (1)▽×(u+v)=▽×u+▽×v ；(2)▽(u ·v)=u ×(▽×v)+v ×(▽×u)+(u ·▽)v+(v ·▽)u ； (3)▽·(u ×v)=v ·(▽×u)-u ·(▽×v)；(4)▽×(u ×v)=(v ·▽)u-(u ·▽)v+(▽·v)u-(▽·u)v.2、若φ是数量函数, A 是向量函数, 则▽×(φA)=φ(▽×A)+▽φ×A.3、若φ是数量函数, A 是向量函数, 则 (1)▽·(▽×A)=0, ▽×▽φ=0,(2)▽×(▽×A)=▽(▽·A)-▽2A =▽(▽·A)-△A.证：1、记u(P 1(x,y,z),Q 1(x,y,z),R 1(x,y,z)), v(P 2(x,y,z),Q 2(x,y,z),R 2(x,y,z)),则(1)▽×(u+v)=⎪⎪⎭⎫⎝⎛∂+∂-∂+∂∂+∂-∂+∂∂+∂-∂+∂yP P xQ Q xR R zP P zQ Q yR R )()(,)()(,)()(212121212121=⎪⎪⎭⎫ ⎝⎛∂∂-∂∂∂∂-∂∂∂∂-∂∂y P xQ xR zP zQ yR 111111,,+⎪⎪⎭⎫⎝⎛∂∂-∂∂∂∂-∂∂∂∂-∂∂y P x Q x R z P z Q yR 222222,,=▽×u+▽×v. (2)∵▽(u ·v)=▽(P 1P 2+Q 1Q 2+R 1R 2)=⎪⎪⎭⎫⎝⎛∂++∂∂++∂∂++∂z R R Q Q P P y R R Q Q P P x R R Q Q P P )(,)(,)(212121212121212121 = ⎝⎛∂∂+∂∂+∂∂+∂∂+∂∂+∂∂,122112211221x RR x R R x Q Q x Q Q x P P x P P,122112211221y RR y R R y Q Q y Q Q y P P y P P ∂∂+∂∂+∂∂+∂∂+∂∂+∂∂⎪⎭⎫∂∂+∂∂+∂∂+∂∂+∂∂+∂∂z R R z R R z Q Q z Q Q z P P z P P 122112211221.又u ×(▽×v)=u ×⎪⎪⎭⎫⎝⎛∂∂-∂∂∂∂-∂∂∂∂-∂∂y P x Q x R z P z Q yR 222222,, = ⎝⎛∂∂+∂∂-∂∂-∂∂,21212121xRR z P R y P Q xQ Q ⎪⎪⎭⎫∂∂+∂∂-∂∂-∂∂∂∂+∂∂-∂∂-∂∂z Q Q y R Q x R P z P P x R P z P P y P R x Q R 2121212121212121,. v ×(▽×u)= ⎝⎛∂∂+∂∂-∂∂-∂∂,12121212xR R zP R yP Q xQ Q ⎪⎪⎭⎫∂∂+∂∂-∂∂-∂∂∂∂+∂∂-∂∂-∂∂z Q Q y R Q x R P z P P x R P z P P y P R x Q R 1212121212121212,. (u ·▽)v=⎪⎪⎭⎫⎝⎛∂∂+∂∂+∂∂z R y Q x P 111v =⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂∂∂+∂∂+∂∂∂∂+∂∂+∂∂z R R y R Q x R P z Q R y Q Q x Q P z P R y P Q x P P 212121212121212121,,(v ·▽)u=⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂∂∂+∂∂+∂∂∂∂+∂∂+∂∂z R R yR Q xR P zQ R yQ Q xQ P zP R yP Q xP P 121212121212121212,,; ∴▽(u ·v)=u ×(▽×v)+v ×(▽×u)+(u ·▽)v+(v ·▽)u. (3)∵▽·(u ×v)=▽·(Q 1R 2-R 1Q 2,R 1P 2-P 1R 2,P 1Q 2-Q 1P 2) =zP Q Q P y R P P R xQ R R Q ∂-∂+∂-∂+∂-∂)()()(212121212121=y P R y R P y R P y P R x R Q x Q R x Q R x R Q ∂∂-∂∂-∂∂+∂∂+∂∂-∂∂-∂∂+∂∂1221122112211221zQP z P Q z P Q z Q P ∂∂-∂∂-∂∂+∂∂+12211221.又v ·(▽×u)=v ·⎪⎪⎭⎫⎝⎛∂∂-∂∂∂∂-∂∂∂∂-∂∂y P xQ xR zP zQ yR 111111,,=yP R xQ R xR Q zP Q zQ P yR P ∂∂-∂∂+∂∂-∂∂+∂∂-∂∂121212121212;u ·(▽×v)=yPR x Q R x R Q z P Q z Q P yR P ∂∂-∂∂+∂∂-∂∂+∂∂-∂∂212121212121;∴▽·(u ×v)=v ·(▽×u)-u ·(▽×v).(4)∵▽×(u ×v)=▽×(Q 1R 2-R 1Q 2,R 1P 2-P 1R 2,P 1Q 2-Q 1P 2)=⎪⎪⎭⎫⎝⎛∂-∂-∂-∂∂-∂-∂-∂∂-∂-∂-∂y Q R R Q x R P P R x P Q Q P z Q R R Q z R P P R y P Q Q P )()(,)()(,)()(212121212121212121212121= ⎝⎛∂∂+∂∂+∂∂-∂∂-∂∂-∂∂-∂∂+∂∂,1221122112211221zP R zR P zR P zP R yQ P yP Q yP Q yQ P,1221122112211221x QP x P Q x P Q x Q P z R Q z Q R z Q R z R Q ∂∂+∂∂+∂∂-∂∂-∂∂-∂∂-∂∂+∂∂⎪⎪⎭⎫∂∂+∂∂+∂∂-∂∂-∂∂-∂∂-∂∂+∂∂y R Q y Q R y Q R y R Q x P R x R P x R P x P R 1221122112211221; 又(v ·▽)u=⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂∂∂+∂∂+∂∂∂∂+∂∂+∂∂z R R yR Q xR P zQ R yQ Q xQ P zP R yP Q xP P 121212121212121212,,; (u ·▽)v=⎪⎪⎭⎫⎝⎛∂∂+∂∂+∂∂∂∂+∂∂+∂∂∂∂+∂∂+∂∂z R R y R Q x R P z Q R y Q Q x Q P z P R y P Q xP P 212121212121212121,,;(▽·v)u=⎪⎪⎭⎫⎝⎛∂∂+∂∂+∂∂z R y Q xP 222u =⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂∂∂+∂∂+∂∂∂∂+∂∂+∂∂z R R y Q R x P R z R Q y Q Q x P Q z R P y Q P xP P 212121212121212121,,; (▽·u)v=⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂∂∂+∂∂+∂∂∂∂+∂∂+∂∂z R R yQ R xP R zR Q yQ Q xP Q zR P yQ P xP P 121212121212121212,,; ∴▽×(u ×v)=(v ·▽)u-(u ·▽)v+(▽·v)u-(▽·u)v. 2、记φ=φ(x,y,z), A=A(P(x,y,z),Q(x,y,z),R(x,y,z)), 则▽×(φA)=⎪⎪⎭⎫ ⎝⎛∂∂-∂∂∂∂-∂∂∂∂-∂∂y P xQ xR zP zQ yR )()(,)()(,)()(ϕϕϕϕϕϕ=⎪⎪⎭⎫ ⎝⎛∂∂-∂∂-∂∂+∂∂∂∂-∂∂-∂∂+∂∂∂∂-∂∂-∂∂+∂∂P yyP Q xxQ R xxR P zzP Q zzQ R yyR ϕϕϕϕϕϕϕϕϕϕϕϕ,,=φ⎪⎪⎭⎫ ⎝⎛∂∂-∂∂∂∂-∂∂∂∂-∂∂y P xQ x R zP z Q yR ,,+⎪⎪⎭⎫ ⎝⎛∂∂-∂∂∂∂-∂∂∂∂-∂∂P yQ xR xP zQ zR yϕϕϕϕϕϕ,,=φ(▽×A)+▽φ×A.3、记φ=φ(x,y,z), A=A(P(x,y,z),Q(x,y,z),R(x,y,z)), 则(1)▽·(▽×A)=▽·⎪⎪⎭⎫ ⎝⎛∂∂-∂∂∂∂-∂∂∂∂-∂∂y P xQ x R zP z Q yR ,,=⎪⎪⎭⎫⎝⎛∂∂-∂∂∂∂+⎪⎭⎫ ⎝⎛∂∂-∂∂∂∂+⎪⎪⎭⎫ ⎝⎛∂∂-∂∂∂∂y P x Q z x R z P y z Q y R x=⎪⎪⎭⎫ ⎝⎛∂∂∂∂-⎪⎭⎫ ⎝⎛∂∂∂∂+⎪⎭⎫ ⎝⎛∂∂∂∂-⎪⎭⎫ ⎝⎛∂∂∂∂+⎪⎭⎫ ⎝⎛∂∂∂∂-⎪⎪⎭⎫ ⎝⎛∂∂∂∂y P z x Q z x R y z P y z Q x y R x =⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛∂∂∂∂-⎪⎭⎫ ⎝⎛∂∂∂∂+⎥⎦⎤⎢⎣⎡⎪⎪⎭⎫ ⎝⎛∂∂∂∂-⎪⎭⎫ ⎝⎛∂∂∂∂+⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛∂∂∂∂-⎪⎪⎭⎫ ⎝⎛∂∂∂∂z Q x x Q z y P z z P y x R y y R x =0. ▽×▽φ=▽×⎪⎪⎭⎫ ⎝⎛∂∂∂∂∂∂z y x ϕϕϕ,,=⎪⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛∂∂∂∂-⎪⎪⎭⎫ ⎝⎛∂∂∂∂⎪⎭⎫ ⎝⎛∂∂∂∂-⎪⎭⎫ ⎝⎛∂∂∂∂⎪⎪⎭⎫ ⎝⎛∂∂∂∂-⎪⎭⎫ ⎝⎛∂∂∂∂x y y x z x x z y z z y ϕϕϕϕϕϕ,,=0. (2)▽×(▽×A)=▽×⎪⎪⎭⎫ ⎝⎛∂∂-∂∂∂∂-∂∂∂∂-∂∂y P xQ x R zP z Q yR ,,= ⎪⎪⎭⎫⎝⎛⎪⎪⎭⎫ ⎝⎛∂∂-∂∂∂∂-⎪⎭⎫ ⎝⎛∂∂-∂∂∂∂⎪⎪⎭⎫ ⎝⎛∂∂-∂∂∂∂-⎪⎪⎭⎫ ⎝⎛∂∂-∂∂∂∂⎪⎭⎫ ⎝⎛∂∂-∂∂∂∂-⎪⎪⎭⎫ ⎝⎛∂∂-∂∂∂∂z Q y R y x R z P x y P x Q x z Q y R z x R z P z y P x Q y ,, =⎪⎪⎭⎫ ⎝⎛∂∂∂+∂∂-∂∂-∂∂∂∂∂∂+∂∂-∂∂-∂∂∂∂∂∂+∂∂-∂∂-∂∂∂z y Q y R x R z x P y x P x Q z Q y z R x z R z P y P x y Q 222222222222222222,,; 又▽(▽·A)=▽⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂z R yQ xP=⎪⎪⎭⎫⎝⎛⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂∂∂⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂∂∂⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂∂∂z R y Q x P z z R y Q x P y z R y Q x P x ,,, =⎪⎪⎭⎫⎝⎛∂∂+∂∂∂+∂∂∂∂∂∂+∂+∂∂∂∂∂∂+∂∂∂+∂∂222222222222,,z R y z Q x z P z y R y Q x y P x z R y x Q x P ; ▽2A=△A=⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂∂+∂∂+∂∂∂+∂∂+∂∂222222222222222222,,z R y R x R z Q y Q x Q z P y P x P ;∴▽×(▽×A)=▽(▽·A)-▽2A =▽(▽·A)-△A.五、管量场与有势场概念：对无源场A ，即div A=0，由高斯公式知，此时沿任何闭曲面的曲面积分都为0，这样的向量场称为管量场. 因为 在向量场A 中作一向量管，即由向量线围成的管状曲面， 用断面S 1, S 2截它，以S 3表示所截出的管的表面，即得到 由S 1, S 2, S 3围成的封闭曲面S ，于是有⎰⎰⋅SdS A =⎰⎰⋅外侧1S dS A +⎰⎰⋅外侧2S dS A +⎰⎰⋅外侧3S dS A =0. 又由向量线与曲面S 3的法线正交知,⎰⎰⋅外侧3S dS A =0.∴⎰⎰⋅外侧1S dS A +⎰⎰⋅外侧2S dS A =0, 即⎰⎰⋅内侧1S dS A +⎰⎰⋅外侧2S dS A . 等式说明，流体通过向量管的任意断面流量相同,∴称场A 为管量场. 如例2，由梯度rm ∇所成的引力场F 是管量场.概念：对无旋场A ，即rot A=0，由斯托克斯公式知，这时在空间单连通区域内沿任何封闭曲线的曲线积分都等于0，该向量场称为有势场. 因为当rot A=0时，由定理22.7推得此时空间曲线积分与路线无关, 且有u(x,y,z), 使得du=Pdx+Qdy+Rdz, 即grad u=(P ,Q,R), u 称为势函数. 所以，若向量场A 的旋度为0，则必存在某势函数u ，使得grad u=A. 这也是一个向量场是某个数量场的梯度场的充要条件. 例1中引力势u=r m 就是势函数. ∴▽u=F=-⎪⎭⎫⎝⎛r z r y r x r m ,,2. 又▽×▽u ≡0, ∴▽×F=0, 它也是引力场F 是有势场的充要条件.若向量场A 既是管量场，又是有势场，则称其为调和场.例2中的引力场F 就是调和场. 若A 是一个调和场，则必有 ▽·A=0, ▽u=A. 显然▽·▽u=▽2u=△u=0, 即必有势函数u 满足222222z uy u x u ∂∂+∂∂+∂∂=0, 这时称函数u 为调和函数. 习题1、若r=222z y x ++, 计算▽r, ▽r 2, ▽r1, ▽f(r), ▽r n (n ≥3). 解：∵x r ∂∂=r x , y r ∂∂=r y , z r ∂∂=r z, ∴▽r=⎪⎭⎫ ⎝⎛r z r y r x ,,=r1(x,y,z); 记u=r 2=x 2+y 2+z 2, ∵x u ∂∂=2x, y u ∂∂=2y, zu ∂∂=2z, ∴▽r 2=▽u=2(x,y,z);记v=r1, ∵x v ∂∂=-3r x , y v ∂∂=-3r y , z v∂∂=-3rz , ∴▽r 1=▽v=31r -(x,y,z);∵x f ∂∂=f ’(r)r x , y f ∂∂=f ’(r)ry , z f∂∂=f ’(r)r z , ∴▽f(r)=f ’(r)r 1(x,y,z); ∴▽r n =nr n-1⎪⎭⎫ ⎝⎛r z r y r x ,,=nr n-2(x,y,z), (n ≥3).2、求u=x 2+2y 2+3z 2+2xy-4x+2y-4z 在O(0,0,0), A(1,1,1), B(-1,-1,-1)处的梯度，并求梯度为0的点. 解：∵x u ∂∂=2x+2y-4, y u ∂∂=4y+2x+2, zu∂∂=6z-4,∴在O(0,0,0), grad u=(-4,2,-4); 在A(1,1,1), grad u=(0,8,2); 在B(-1,-1,-1), grad u=(-8,-4,-10);又由2x+2y-4=0, 4y+2x+2=0, 6z-4=0, 解得x=5, y=-3, z=32, ∴在(5,-3,32), |grad u|=0.3、证明梯度的基本性质1~5. 证：见梯度的基本性质.4、计算下列向量场A 的散度与旋度：(1)A=(y 2+z 2,z 2+x 2,x 2+y 2)；(2)A=(x 2yz,xy 2z,xyz 2)；(3)A=⎪⎪⎭⎫⎝⎛++xy z zx y yz x . 解：(1)∵P=y 2+z 2, Q=z 2+x 2, R=x 2+y 2； ∴div A=x ∂∂(y 2+z 2)+y ∂∂(z 2+x 2)+z ∂∂(x 2+y 2)=0；又y ∂∂(x 2+y 2)-z ∂∂(z 2+x 2)=2y-2z; z ∂∂(y 2+z 2)-x∂∂(x 2+y 2)=2z-2x; x∂∂(z 2+x 2)-y ∂∂(y 2+z 2)=2x-2y. ∴rot A=2(y-z,z-x,x-y).(2)∵P=x 2yz, Q=xy 2z, R=xyz 2； ∴div A=x ∂∂(x 2yz)+y ∂∂(xy 2z)+z∂∂(xyz 2)=6xyz ；又y ∂∂(xyz 2)-z ∂∂(xy 2z)=x(z 2-y 2); z ∂∂(x 2yz)-x∂∂(xyz 2)=y(x 2-z 2); x∂∂(xy 2z)-y ∂∂(x 2yz)=z(y 2-x 2). ∴rot A=(x(z 2-y 2),y(x 2-z 2),z(y 2-x 2)).(3)A=⎪⎪⎭⎫ ⎝⎛++xy z zx y yz x . ∵P=yz x , Q=zxy, R=xy z ；∴div A=⎪⎪⎭⎫ ⎝⎛∂∂yz x x +⎪⎭⎫ ⎝⎛∂∂zx y y +⎪⎪⎭⎫ ⎝⎛∂∂xy z z =xyzx yz 111++； 又⎪⎪⎭⎫ ⎝⎛∂∂xy z y -⎪⎭⎫ ⎝⎛∂∂zx y z =22xy z xz y -; ⎪⎪⎭⎫ ⎝⎛∂∂yz x z -⎪⎪⎭⎫ ⎝⎛∂∂xy z x =22yz x y x z-; ⎪⎭⎫ ⎝⎛∂∂zx y x -⎪⎪⎭⎫ ⎝⎛∂∂yz x y =z x y z y x 22-. ∴rot A=⎪⎪⎭⎫⎝⎛---x y y x z x x z y z z y xyz 222222,,1.5、证明散度的基本性质1~3. 证：见散度的基本性质.6、证明旋度的基本性质1~3. 证：见旋度的基本性质.7、证明：场A=(yz(2x+y+z),zx(x+2y+z),xy(x+y+2z))是有势场并求其势函数.证：P=yz(2x+y+z), Q=zx(x+2y+z), R=xy(x+y+2z),y ∂∂[xy(x+y+2z)]-z∂∂[zx(x+2y+z)]=x 2+2xy+2xz-x 2-2xy-2xz=0; z ∂∂[yz(2x+y+z)]-x∂∂[xy(x+y+2z)]=2xy+y 2+2yz-2xy-y 2-2yz=0; x∂∂[zx(x+2y+z)]-y ∂∂[yz(2x+y+z)]=2xz+2yz+z 2-2xz-2yz-z 2=0.∴对空间任一点(x,y,z)都有rot A=(0,0,0)=0i+0j+0k=0, ∴A 是有势场. 由d[xyz(x+y+z)]=yz(2x+y+z)dx+xz(x+2y+z)dy+xy(x+y+2z)dz 知， 其势函数为u(x,y,z)=xyz(x+y+z)+C.8、若流体流速A=(x 2,y 2,z 2), 求单位时间内穿过81球面x 2+y 2+z 2=1, x>0,y>0,z>0的流量.解：设S 为所给81球面，S 1, S 2, S 3分别是S 在三个坐标面上的投影, 则 所求流量为：⎰⎰⋅SdS n A 0+⎰⎰⋅11S dS n A +⎰⎰⋅22S dS n A +⎰⎰⋅33S dS n A =⎰⎰⎰⎪⎭⎫ ⎝⎛球体81V divAdV=⎰⎰⎰++Vdxdydz z y x )(2=⎰⎰⎰++103202sin )cos sin sin cos (sin 2dr r d d ϕϕθϕθϕϕθππ=⎰⎥⎦⎤⎢⎣⎡++2021)sin (cos 421πθθθπd =83π.注：其中n 0, n 1, n 2, n 3分别是S, S 1, S 2, S 3的单位法矢，显然有A|n i (i=1,2,3)，∴A ·n i =0，从而⎰⎰⋅iS i dS n A =0 (i=1,2,3), 于是所求流量为：⎰⎰⋅SdS n A 0=83π.9、设流速A=(-y,x,c) (c 为常数)，求环流量： (1)沿圆周x 2+y 2 =1, z=0；(2)沿圆周(x-2)2+y 2 =1, z=0.解：(1)圆周x 2+y 2 =1, z=0的向径r 适合方程r=costi+sintj+0k(0≤t ≤2π). ∵A ·dr=(-sinti+costj+ck)·(-sinti+costj+0k)dt=dt, ∴所环流量为⎰⋅c dr A =⎰π20dt =2π.(2)圆周(x-2)2+y 2 =1, z=0的向径r=(2+cost)i+sintj+0k (0≤t ≤2π); ∵A ·dr=[-sinti+(2+cost)j+ck]·(-sinti+costj+0k)dt=(2cost+1)dt, ∴所环流量为⎰⋅c dr A =⎰+π20)1cos 2(dt t =2π.。
