自由曲面坐标系stokes公式证明和计算实例.doc

自由曲面坐标系Stokes公式证明和计算实例 杨科 四川 成都 610017 028E-mail: more2010e@ 摘 要： 用符号逻辑的方式,证明在无穷多个任意参数曲面坐标系[包括正交曲线坐标系和非正交曲线坐标系] Stokes公式的存在, 将空间环路积分和曲面积分推广到无穷多个任意参数曲面坐标系,实现任意曲面积分, 实现向量场(静电场、磁场、流体场、引力场等)在任意自由空间(闭合路径、曲面)的积分计算. 关键词: 向量场 自由参数曲面坐标系 Stokes公式 证明 计算实例 中图分类号：O17/O412.3 1.1自由曲面坐标系Stokes公式证明: Stokes公式 设光滑或分片光滑的有向曲面S的正向边界L+为光滑或分段光滑的闭合曲线.如果函数P(x,y,z),Q(x,y,z),R(x,y,z) [构成向量场A] 在有向曲面S上有一阶连续偏导数,则 (1) 其中rotA为向量场A的旋度,n为有向曲面S的单位外法向量 证明: 定义曲面S的参数表达式:[a sin(u)cos(v),b sin(u) sin(v),c cos(u)] (2) 其中a,b,c为任意非零常数或一阶连续可导函数(以u,v为自变量的三角函数)表达式;设定参数u,v的变化范围[0,/n - ],[0,2],其中n为任意常数,并且n 1;为任意常数或连续函数(以v为自变量的三角函数)表达式,使曲面S非闭合. [12] 定义边界曲线L的参数表达式:[ cos(v), sin(v),] (3) 其中,,为依存于a,b,c的常数(其中a0,b0)或一阶连续可导函数(以v为自变量的三角函数)表达式; 因参数v的变化范围为[0,2],边界曲线L闭合. [12] 向量场A在边界曲线L的环路积分: (4) 定义并计算偏导数矩阵M,获取曲面S的切平面法向量: M = = (5) 从(5)式分别提取i,j,k项系数,获得曲面S的切平面法向量: [,,] (6) 计算向量场A的旋度rotA,并将其从直角坐标形式(7)转变为参数曲面坐标形式(8): (7) = (8) 旋度rotA(8)与曲面S的切平面法向量(6)的空间点积对变量u,v的积分: (9) = 即(4)式=(9)式: 亦可表述为 (1), 证毕 1.2球面坐标系Stokes公式证明: Stokes公式 设光滑或分片光滑的有向曲面S的正向边界L+为光滑或分段光滑的闭合曲线.如果函数P(x,y,z),Q(x,y,z),R(x,y,z) [构成向量场A] 在有向曲面S上有一阶连续偏导数,则 (1) 其中rotA为向量场A的旋度,n为有向曲面S的单位外法向量 证明: 定义球面S的参数表达式:[sin(u) cos(v), sin(u) sin(v), cos(u)] (2) 设定参数u,v的变化范围[0,/n - ],[0,2], 其中n为任意常数,并且n1; 为任意常数或连续函数(以v为自变量的三角函数)表达式,使球面S非闭合. 将变量u的边界值带入曲面CS的参数表达式,获得边界曲线CL参数表达式: [] (3) 因参数v的变化范围为[0,2],边界曲线L闭合. 向量场A在边界曲线L的环路积分: (4) 定义并计算偏导数矩阵M,获取球面S的切平面法向量: (5) M = = 从(5)式分别提取i,j,k项系数,获得球面S的切平面法向量: [,,] (6) 计算向量场A的旋度rotA,并将其从直角坐标形式(7)转变为球面坐标形式(8): (7) (8) 旋度rotA(8)与球面S的切平面法向量(6)的空间点积对变量u,v的积分: (9) = 即(4)式=(9)式: 亦可表述为 (1), 证毕 1.3自由曲面坐标系Stokes公式证明[Maple 11计算机代数系统格式]: 符号表达系统: 向量场V,向量场V的旋度cV1,cV2, 任意曲面CS[非闭合],曲面CS的闭合曲线边界CL,曲面CS的切平面法向量[A,B,C] restart; # 系统复位 with(linalg): # 加载线性代数符号分析库 CS:=[a*sin(u)*cos(v),b*sin(u)*sin(v),c*cos(u)]; # 定义任意参数曲面CS,其中a,b,c为任意常数或一阶连续可导函数(以u,v为自变量的三角函数)表达式,并且a0,b0,c0. [12] rgu:=[0,Pi/n-theta]; # 参照于任意闭合参数曲面,参数u的变