常微分方程解析理论-正文 复域上的常微分方程理论;应用复变函数论研究微分方程的性状,以及把微分方程的解视为由方程定义的解析函数,并直接从微分方程本身研究解的性质的理论。这是基于A.-L.柯西的基本定理,即在对微分方程作极为广泛的假设下,它的积分是复变数的解析函数。常微分方程解析理论与复变函数理论的发展密切相关。它的先驱性工作是由柯西、(G.F.)B.黎曼、I.L.富克斯、(J.-)H.庞加莱以及P.班勒卫等人所作。
连带勒让德多项式 连带勒让德多项式递推公式推导
解的存在性和惟一性定理 微分方程理论中基本的问题是已给的方程是否有解,早先的数学家们力图通过已知初等函数的有限组合来表示微分方程的解,但在这个观念下大多数微分方程不可积。这实际上是要求方程的大范围通解,是不合适的,因为典型的分析运算与极限过程只要求局部的观点。另一方面,在物理和力学中的问题常是只要求适合某些补充条件的特解。于是柯西提出考虑如下的问题:方程
e=e+(2l+1)(abs(a1).^2+abs(b1).^2);(1)的右端?(z,w)在(z0,w0)点的某个邻域内解析,问是否存在z的解析函数w(z;z0,w0),它在w0点的邻域满足方程(1),并且满足初值条件w(z0;z0,w0)=w0。他证明了在上述假设下,解是存在且惟一。这个定理称为柯西存在性定理。在复域中通常应用幂级数展开式给出惟一的形式解,然后用与某个已知的收敛幂级数相比较的方法(优函数方法)给出形式解的收敛性证明,从而完成存在性和惟一性定理的证明。
奇点 柯西存在性定理所证明的微分方程的解是局部的。即给出了一个解析函数元素,应用外尔斯特拉斯的解析开拓(见常微分方程初值问题)的方法,从z0点的邻域沿一途径Г开拓这个函数元素,如果方程(1)的右端也能沿Γ开拓,则解的开拓元素也满足方程。如果沿着所有可能的途径进行开拓,则得到的所有函数元素构成的在大范围定义了一个单值的或多值的函数。现在重要的问题是在解的整个存在区域上来研究它,而解的存在区域和解的性质是由它的奇点所决定的,这里奇点是指柯西存在性定理不成立的那些点。因此需要研究所考虑的方程的解的奇点的位置和性质。
微分方程的解出现的奇点较解析函数论中的情况要复杂得多。首先当自变量围绕某些点转一圈以后,函数从一个值变为另一个值,称这些点为分支点。代数函数可能具有的奇点称为代数奇点。非代数奇点的分类基于不定区的概念,函数?在z0点的不定区是指以z0为中心的小圆在?映射下的像当圆半径趋于0时的极根。若点z0的不定区由一点组成,则称z0为超越奇点,否则称为本性奇点。富克斯还对微分方程解的奇点提出一种重要的区分,即分为固定奇点和流动奇点。前一种由微分方程本身给出其位置和性质,与方程的个别解无关,也即与通解中所含的任意常数无关。后者则依赖于柯西问题的初始值,也就是依赖于特解的选择,它与任意常数一起变动。例如方程 的解以整数和无穷远点为固定奇点(极点);和 分别有解为 和此时с分别是流动代数分支点,流动对数分支点和流动本性奇点。
班勒卫曾证明如下的定理(称班勒卫定理):若z0是方程(1)的解的奇点,则(z0,w0)不是方程右端?(z,w)的全纯点。
这个定理首次确定解的奇点和方程奇点的关系,同时还说明在方程右端 ?(z, w)的全纯点处除了全纯解之外,不存在非全纯的解。当方程右端是w 的有理函数时,班勒曾卫列举可能出现奇点的种种情况。此外,如果?(z,w)=P(z,w)/Q(z,w),(z0,w0)是P(z,w)和Q(z,w)的全纯点, 但P(z0,w0)=Q(z0,z0)=0,这种不确定的情形下,即使在P(z,w)和Q(z,w)是z 和w 的线性函数的情形,其解在z0点的邻域的性质也相当复杂。
一般地,当对方程的性状加上某些限制以后,也带给解的奇点某些限制,例如线性微分方程的解无流动奇点。1887年班勒卫曾证明,未知函数及其导数代数地出现于方程,而系数是z的解析函数的一阶代数微分方程,它的解无流动超越奇点和流动本性奇点。
反过来,如果对解的奇点作某些限制时,微分方程也要适合某些条件,例如其解无任何奇点的方程必为一个重要的结论是:如果方程(1)的右端是w 的有理函数,其解无流动代数分支点,则方程(1)必化为如下的黎卡提方程
(2) 线性常微分方程 一类很重要的常微分方程,未知函数的阶导数是较低阶导数的线性函数,一般可写成
如果右端恒为零,则称为齐次线性微分方程。如果知道了齐次方程的通解,则能通过参数变动法(或称常数变易法,见初等常微分方程)得到非齐次方程的解。因此线性方程的中心问题是研究齐次方程,而n阶齐次线性方程的通解能由 n个线性独立的特解线性地表示出来。这个基本性质大大简化了对线性方程的研究。此外,在力学和电路理论中有关振动问题常化归为二阶线性方程,纯粹数学中的许多完美思想也是从这类方程的研究中产生,而且常常能展现出n阶线性方程的许多性质。所以大量的工作是关于二阶线性方程的。它的一般形式可写成
(3)已知线性方程的解只有固定奇点,即解w(z)在一点的性质依赖于方程系数 p(z)和 q(z)在该点的性质。许多物理问题引起的微分方程都有奇点,因而对适应这种物理情况的解有较详细的讨论。在奇点领域,方程(3)的解能有如下表示式:设w1(z)和w2(z)是奇点 z0邻域的两个线性独立解,当围绕z0转一周时,它们接受一个线性变换,即 令λ1和λ2是A=的特征根,则当λ1≠λ2时,(3)的解能写为
当λ1=λ2时,则为
下述几个特殊的二阶线性方程在实际应用和理论中都很重要。
富克斯方程 它是奇点全为正则奇点的方程。由于z0为正则奇点的充分必要条件是(z-z0)p(z)和(z-z0)2q(z)在z0点领域全纯,因此富克斯方程可写为
(4)它也是具有正则奇点的的方程,其中p1(z)、q1(z)在αk点全纯;并称
img src="/uploads/90f1ad1c6ce06dde800e740de1ff11df.png" data-original=image/67-7.gif align=abiddle> (5)为在αk点的指标方程,其中,。方程(5)的根称为指标数,记为且有的富克斯关系式这里αn+1=。如果奇点的个数<4且都位于有限平面内,则方程能由奇点的位置和相应的指标数完全确定。特别是当 n=3时即导出超几何方程。对这个方程的研究有着悠久的历史,许多杰出的数学家如L.欧拉、C.F.高斯、E.E.库默尔和黎曼等人都有重要的贡献。这类方程在很多情形中出现,它与共形映射、分方程、连分数和自守函数都有关系;且其理论具有形式上的高度完美性,今设 αk(k=1,2,3)为奇点,()为相应的指标数,则方程可写为
这个形式为黎曼所提出,又称为黎曼方程,它的积分(解)能由黎曼的P函数所表示,通常记为
一个相关的问题是确定一切多值函数,它们仅以给定的αk(k=1,2,3)为奇点,它的奇异性满足一定的要求,在每个奇点附近,此函数有两个独立的值,而任意三个值w1(z)、w2(z)、w3(z)线性相关,这个问题称为黎曼问题。它能化为黎曼方程的积分,一般地可通过超几何函数表示出来,这个问题先后由D.希尔伯特、J.普莱姆利和G.D.伯克霍夫解决和推广。
若富克斯方程的奇点为0、1和,则引入超几何函数中常用的参数之后能导出高斯的标准形式
称为高斯方程或称超几何方程。它的解可表为超几何级数
式中(p)n=p(p+1)(p+2)…(p+n-1)。库默尔于1834年找出24个变换,使得具有三个至多是简单奇点的二阶富克斯方程化为具有不同参数的超几何方程。这24个变换对应着解由超几何级数表示的24个表达式。
勒让德方程 它是形如
的方程。A.-M.勒让德于1785年首先考虑α=n为非负整数的情形。若令t=(1-z)/2,则它能化为以n+1、-n和1为参数的超几何方程,在z=1的全纯解为n阶勒让德多项式
。 贝塞尔方程 它是形如
施瓦兹方程 它是与二阶线性微分方程紧密相关的一类方程, 它由共形地映w上半平面为z平面上圆弧多边形内部的函数所满足,方程为
(6)式中称为施瓦兹导数;α1,α2,…,αn为多边形的角点, P2n-4(w)和2n-4次多项式。方程(6)的解具有一个重要的性质,即当围绕奇点环行一周时,它接受一个分式线性变换 又知二阶线性方程的两个线性独立的解之比亦具有相同的性质,因此方程(6) 的求解问题能化为适当选取的二阶线性方程的求解。设G是一分式线性变换群,?(z)为一单值亚纯函数,如对于任一g∈G有?(g(z))=?(z),则称?(z)是关于群 G的自守函数。自守函数与二阶微分方程有下述的关系:设w=?(z)为自守函数,则z作为w 的函数可用微分方程z〃+uz=0的两个独立解z1(w)和z2(w)之商表示<即的反函数为w=?(z)。
非线性微分方程 由于许多物理系统是非线性的,从而描述它们的微分方程也是非线性的,即未知函数或其导数非线性地出现于方程之中。对于非线性方程一般性质的了解不像线性方程那样完备和深入,而是知道得很少,而且它具有线性方程理论中所未见的新现象。下面只叙述非线性方程理论中的一些事实。
1856年C.A.布里奥和J.-C.布凯考虑如下的方程
(7)式中 F(z,w) 是在某个双圆柱内两个变量的全纯函数。首要的问题是方程(7)是否存在全纯解。他们证明:如果q不是正整数。则(7)在z=0有惟一的全纯解w(z),且w(0)=0。若q=1,p≠0,则不存在全纯解。若p=0,q=1,则有无穷多个全纯解。他们还讨论下面的方程
(8)式中P(x,y)是x和y的常系数多项式,并称(8)为k阶布里奥-布凯方程,或简称BB方程。他们指出,每一椭圆函数满足某个k阶BB方程,并且BB方程具有大范围单值亚纯解的必要条件是代数曲线P(x,y)=0的亏格为0或1。
19世纪末,班勒卫首先讨论了方程式中F(z,w,w┡)是w和w┡的有理函数,系数为z的解析函数。他考虑定出只具有固定分支点和本性奇点的方程。B.O.冈比埃和富克斯对此问题亦作出重要贡献。一般方法是由班勒卫提出,基本技巧是他的α-方法。他们找到了50个不同的类型,但大多数能化为已知的方程,如线性方程或黎卡提方程。只有 6种类型的方程导出新的超越亚纯函数,这些方程是:
< align=center> 等等,并称这些方程为班勒卫方程,它们的解称为班勒卫函数。1913~1914年,P.L.布特鲁对一类二阶方程发展了渐近积分的方法,并指出班勒卫方程的解在某种意义下渐近于外尔斯特拉斯椭圆函数。
常微分方程理论中奈望林纳理论的应用 20世纪20年代芬兰数学家R.奈望林纳创立了亚纯函数值分布理论。不久日本数学家吉田耕作应用此理论于一类非线性常微分方程的研究。50年代H.维蒂希更系统地研究了奈望林纳理论对常微分方程理论的意义,使得这一理论成为研究一类方程解的某些大范围性质(解的增长性,值分布性质,因子分解等)的重要工具。作为柯西存在惟一性定理的直接推论是下述常系数微分方程
(9)的每一非常数亚纯解 w(z)都不取αj(j=1,2,…,n)为值。另方面,根据亚纯函数皮卡定理,任一非常数亚纯函数能取所有的复值为值,至多除去两个例外。因此,如果方程(9)具有非常数亚纯解,则必有方程(9)的右端对w的次数≤2。对此,在1913年J.马尔姆奎斯特得到了重要的推广,他证明了下述的马尔姆奎斯特定理:设方程(1)的右端是z和w的有理函数,如果方程存在全平面单值超越亚纯解,则(1)必为黎卡提方程。1933~1934年吉田耕作应用奈望林纳理论给出这个定理一个漂亮的证明,并且大大推进了结果。由于微分方程的解更多出现为有限多值的解析解,即代数体函数解,他还考虑了方程
(10)的代数体解存在的必要条件,其中P(z,w)和Q(z,w)分别是w 的p次和q次多项式,系数是z的有理函数。他证明:若方程(10)存在v值超越代数体解,则必有p≤2nv和q≤2n(v-1)。特别地,当 n=v=1时即是马尔姆奎斯特定理。
等等。
此外,对于代数微分方程亦有相应的结果,数学工作者对相当广泛的高阶代数微分方程存在“较快”增长的代数体函数解的必要条件亦得到形式的马尔姆奎斯特型定理。近年来奈望林纳理论还被用来研究常微分方程复振荡理论、解的增长性估计和解的因子分解等。
这个可以用递归函数,计算n阶勒让德多项式的值的C函数为
式中ck(k=0,1,2)是常数,uk(z)(k=1,2,3)是在z0点邻域的洛朗级数。这个表示式的作用在于将解的单值解析部分和多值解析部分明显地表示出来。另一方面在大多数物理问题中,奇异性比较“弱”,出现较弱奇异性的点称为正则奇点,其定义如下:若在z0点,uk(z)(k=1,2,3)只有极点,则称z0为正则的;若uk(z)中至少有一个以z0为本性奇点,则称z0是非正则的。legendre(int
n,
x)
{if(n
==
0)
1;
==
x;
((2
n
-legendre(n
-1,
x)
legendre(n
-2,
x))
/(二)岩(矿)石标本的磁参数测定n;
重力的三个公式分别为:G=mg,g=G/m,m=Gif(n == 0)/g。式中g=9.8N/Kg,它表示的物理意义是:质量是1Kg的物体受到的重力为9.8N。
由于地球的吸引而使物体受到的力,叫做重力。方向总是竖直向下,不一定是指向地心的(只有在赤道和两极指向地心)。地面上同一点处物体受到重力的大小跟物体的质量m成正比,同样,当m一定时。
else if(n>1) a=((2n-1)xp(n-1,x)-(n-1)p(n-2,x))/n;//少了个,改为else if(n>1) a=((2n-1)xp(n-1,x)-(n-1)p(n-2,x))/n
(4)按一定比例进行重复测量,计算视磁化率和剩余磁化强度的相对误。17物体所受重力的大小与重力加速度g成正比,用关系式G=mg表示。通常在地球表面附近,g值约为9.8N/kg,表示质量是1kg的物体受到的重力是9.8N。9.8N是一个平均值,在赤道上g小,g=9.79N/kg,在两极上g,g=9.83N/kg。N是力的单位,字母表示为N,1N大约是拿起两个鸡蛋的力。89年法国大革命后,于1790年宣布要对当时相当混乱的度量衡制度进行改革。科学院组成了一个由拉格朗日为首的委员会。委员会建议以从赤道到北极的一千万分之一为长度基本单位――米,这个方案于1791年被法国国民议会通过。于是就要着手实际测量从赤道到北极的长度。勒让德参加了(2)在标本中心建立起空间直角坐标系,如图2-4-2所示。测量,并且是经度局的一名成员。1813年拉格朗日逝世,勒让德接替他成为经度局的主席。他在数学上的贡献,勒让德多项式就是在计算地球形状时的一项创造。
(一)CZM-2质子磁力仪磁测资料的整理内容及方法
x1.检查与验算
printf("请输入n,x的值:");检查记录中各项填写是否完整、清楚,并且验算平均读数。
2.计算读数
将各测点的平均读数Ti减去早基点
(0点)的平均读数T0计算读数ΔTi,即ΔTi=Ti-T0。3.日变改正
从日变曲线图上查出日变影响值Δ'1。
具体方法为:在时间轴(横轴)上定出早基点的读数时刻t0,过t0点作一条平行于纵轴的直线,过该直线与日变曲线的交点作一条平行于横轴的直线如图2-4-1所示的虚线所示,此虚线即为该台仪器的日变曲线的零值线,就可以查出各个测点的日变影响值Δ1i',而日变改正值Δi=-Δ1i'。对于MP-4质子磁力仪可直接将日变观测仪与测量仪对接,自动进行日变改正,也可用上述方法改正。日变改正的精度用仪器噪声水平衡量。
图2-4-1 日变曲线
4.正常场改正及高度改正
应用地球物理教学实习指导
式中:gmn和hmn为高斯系数(可查表求得);Pmn(cosθ)为n次m阶缔和勒让德多项式;λ为各节点的经度值;θ为各节点的余纬度值。
应用地球物理教学实习指导
由地磁场正常梯度公式,即
应用地球物理教学实习指导
可计算高度改正值δT,式中地球半径r=6371200m,Δh是测点与总基点的高,当高相30m、地磁场T0=5000nT时,垂向变化为-0.72nT,因此要进行高度改正。高度改正以总基点高程算起,一般T0=5000nT时,每42m高改正1nT,若高相不大的平原地带这项改正可以忽略。
综上所述,终求取ΔT异常的计算步骤是:
(1)由GPS测量给出的测量点的大地经纬度坐标λ和θ,代入式(2-4-9)求得X,Y,Z。
(3)将各测点经过各项改正的总磁场强度T减去T0值求得ΔT。
5.综合影响改正值的计算方法
质子磁力仪不受温度变化的影响,无零点漂移。但两次对基点消除日变Δ1影响后,仍有值,此值Δ2称综合影响。即
应用地球物理教学实习指导
此Δ2'应按时间线性分配到各个测点上,其改正值为
应用地球物理教学实习指导
式中:Δt为两次基点读数的时间;ti,t0分别为测点和早基的读数时间。
但本实习对总的Δ2有要求:
若 倍质量检查均方时,则本次测量作废;
应用地球物理教学实习指导
倍质量检查均方时 倍质量检查均方时,按上述方法在测点上作线性分配;
应用地球物理教学实习指导
倍质量检查均方时,可忽略不记而不作改正。
6.总基点改正
7.正常梯度改正
由于本实习测区范围小(南北方向仅0.5km),因此,不进行地磁场的正常梯度改正。
综上所述,各测点异常值的算式应是
应用地球物理教学实习指导
全部数据的计算,必须百分之百的进行检查,并且计算者、检查者必须签名以示负责。
记录和计算均要求内容完善、记录真实、字迹清楚、格式统一、不能撕页,原始记录不能涂改,也不能用橡皮擦改,如有记录有错只能用细线划去在旁边另写。
应指出,教学实习中测量地磁场总强度的总均方误设计为±5nT,如经质量检查不合格则应查明原因,部分或全部返工。
磁参数系指磁化率κ和剩余磁化强度Mr,κ是纯量,Mr是向量,因此为测量出标本磁性的大小和方向,必须采集定向标本,即采集前在露头上标出磁北方向和铅直向下的方向,如图2-4-2所示。用MP-4型质子磁力仪(高斯位置)测磁参数的方法如下。
1.准备工作
图2-4-2 采集定向标本示意图 图2-4-3 特制的测量磁性的仪器架
2.测定步骤
(1)读取T0(未置标本的读数);
(2)将标本置于标本托板上,设Z轴朝斜下方时的读数为T1,朝斜上方时读数为T2对X轴和Y轴也可以读取T3,T4和T5,T6;探头中心位于标本磁距X,Y,Z分量的高斯位置;
(3)读取T0'(取下标本);
(4)量取标本中心到探头中心的距离R(以m为单位);
(5)量取标本的体积(以m3为单位)。
3.计算磁参数的公式
计算标本视磁化率的公式为
应用地球物理教学实习指导
式中:R以m为单位;V以m3为单位;T0,T1,T2,T3,T4,T5,T6以nT为单位。
计算剩余磁化强度大小的公式为
应用地球物理教学实习指导
计算剩余磁化强度偏角φ和倾角θ的公式分别为
应用地球物理教学实习指导
4.技术及精度要求
(1) 应分别大于等于T0;
(2)标本中心与探头中心距离R应小于等于0.45m;
(3)|T0-T'0|≤2nT;
clc
elsealpha=5.2182;
e=0;
f=0;
syms x f1 f2
if l==1
f1=sin(x)./x-cos(x);
f1dot=diff(f1,x,l);
f2=(sin(x)+icos(x))./x-cos(x)+isin(x);
f2dot=diff(f2,x,l);
f1=(pix/2).besselj(l+0.5,x);
f1dot=(pix/2).besselj(l-0.5,x)-l(pix/2).besselj(l+0.5,x)./x;
f2=(pix/2).(besselj(l+0.5,x)-ibessely(l+0.5,x));
f2dot=(pix/2).(besselj(l-0.5,x)-ibessely(l-0.5,x))-l(pix/2).(besselj(l+0.5,x)-ibessely(l+0.5,x))./x;
end
g1=inline(vectorize(f1), x);
g1dot=inline(vectorize(f1dot), x);
g2=inline(vectorize(f2), x);
g2dot=inl}else if (n==1)ine(vectorize(f2dot), x);
a1=(g1dot(nalpha).g1(alpha)-ng1(nalpha).g1dot(alpha))./(g1dot(nalpha).g2(alpha)-ng1(nalpha).g2dot(alpha));
b1=(ng1dot(nalpha).g1(alpha)-g1(nalpha).g1dot(alpha))./(ng1dot(nalpha).g2(alpha)-g1(nalpha).g2dot(alpha));
f=f+(2l+1)real(a1+b1);
end
kext=2/alpha.^2.f
首先,由于已经给定了次数而且很低,就不需要积分。Pl(n=1.55;x)中x次幂就是x^l,P2(x)=0.5(3x^2-1),P1(x)=x,P0(x)=1.于是f(x)=4/3P2(x)+3P1(x)+18/3P0(x)
是2return(c)
if(n递归公式
return1 (n=0)
pn(x)=x (n=1)
((2n-1)xpn-1(x)-(n-1)pn-2(x))/n (n>1)
例如:
#include
{int f;
if(n<0)
{f = -1;
printf("error, n should be larger than 0");
}else if(n==0)
{f = 1;
{f = x;
}else if (n>1)
{f=((2n-1)xp(n-1, x)-(n-1)p(n-2,x))/n;
}return f;
}void main ()
{int n,x;
scanf("%d%d",&n,&x);
printf("结果为:%d
getchar();
}扩展资料:
由于递归引起一系列的函数调用,并且可能会有一系列的重复计算,递归算法的执行效率相对较低。当某个递归算法能较方便地转换成递推算法时,通常按递推算法编写程序。
例如上例计算斐波那契数列的第n项的函数fib(n)应采用递推算法,即从斐波那契数列的前两项出发,逐次由前两项计算出下一项,直至计算出要求的第n项。
参考资料来源:
是2return(c)
语法不熟悉造成的问题,造成stack overfelselow....
注意是:
if(xxx)
else if(xxx)
才
和上面的不一样 在执行-(n完if(n==1) return x;以后,因为return这个关键字,所以已经返回了,后面的所有语句都不再执行...而if(n==1) p=x;却不一样,这句执行完后还会接下去执行else后面的语句...所以你那样写,程序并不会按照你想的那样"跑"...
把代码改成这样,应该就OK,我测试时OK的...
double 基点联测时所求出的各基点相对于总基点的磁场值,即为总基点的改正值,如基点A比总基点高50nT,则总基点改正值Δ3=50nT;如基点A比总基点低30nT,则总基点改正值Δ3=-30nT。f(double x, int n)
{double p;
p = 1;
else if(n == 1)
p = x;
p = ((2n - 1)xf(x, n - 1) - (n - 1)f(x, n - 2))/n;
return p;
版权声明:本文内容由互联。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌抄袭侵权/违法违规的内容, 请发 a13828211729@163.com 邮箱删除。