奇摄动Volterra型积分微分方程非线性边值问题.pdf

桥梁非线性动力响应方法及volterra级数非线性方法探讨作者：冉晴等来源：《建筑科技与经济》2014年第01期摘要：强震下桥梁的破坏和倒塌，造成的人员伤亡和经济社会效益是不容忽视的；另一方面，随着桥梁结构体系越来越复杂和延性抗震设计理念的转变，都使得桥梁非线性抗震验算越来越受到工程师的重视。

本文回顾了桥梁非线性动力学问题的发展，总结了目前桥梁非线性抗震验算方法的优、缺点，着眼于非线性系统理论及最新成果，通过volterra级数在非线性系统理论中的成果，探讨了volterra级数在桥梁工程的运用的可行性及优越性，并提出了volterra 级数运用于桥梁结构的关键性问题。

关键词：桥梁；非线性动力学；volterra级数；非线性频率响应函数1.引言我国广阔、复杂的地貌，造就了多山多河的地形，为了方便人们的出行，桥梁作为可以跨谷跨河的工程结构物得到了广阔的发展。

其中桥梁动力学问题分析作为桥梁分析中重要的一环，关系到桥梁功能的正常使用和生命财产安全，被研究学者和设计工程师重点关注着。

在2008年我国汶川8.0级地震，有24条高速公路、6140座桥梁受损，导致了69225人遇难，4600多万人受灾[1]，造成难以估计的损失。

只进行桥梁线性动力分析已经不满足需求，因此破坏性地震下桥梁的非线性动力验算愈来愈得到工程人员的重视。

另一方面，现代桥梁体系延性设计的理念的转变、减隔震支座的采用和大量大跨桥梁复杂新体系的出现，都决定着桥梁非线性动力响应分析的必要性。

2.桥梁非线性动力响应发展经过了各界研究学者的1个多世纪的探索，从1900年日本提出的采用静力等效力来模拟地震力，到现代能够比较完整的考虑地震三大主要效应（峰值、频率、持时）的时域分析方法和频域分析方法，对结构地震动力的线性响应考虑也愈趋细致。

非线性动力响应不再满足叠加原理，表现为非常复杂的力学行为，如分叉、混沌现象。

这些理论和计算远远不理想，本文并不探讨。

第３６卷第４期２０２２年７月兰州文理学院学报(自然科学版)J o u r n a l o fL a n z h o uU n i v e r s i t y ofA r t s a n dS c i e n c e (N a t u r a l S c i e n c e s )V o l ．３６N o ．４J u l ．２０２２收稿日期:２０２２Ｇ０３Ｇ１３基金项目:国家自然科学基金(１１８６１０６８);新疆维吾尔自治区自然科学基金杰出青年基金项目(２０２２D ０１E １３)作者简介:罗紫洋(１９９６Ｇ),男,新疆昌吉人,在读硕士,研究方向:微分方程理论及数值模拟．E Ｇm a i l :４６６９４９８４１＠q q ．c o m．∗通讯作者:张新东(１９８１Ｇ),男,江苏徐州人,教授,硕士生导师,研究方向:微分方程理论及数值模拟．E Ｇm a i l :L i Ｇa o yu a n １１２６＠１６３．c o m ㊀㊀文章编号:２０９５Ｇ６９９１(２０２２)０４Ｇ００１０Ｇ０５V o l t e r r a 型微分Ｇ积分方程的数值格式构造及理论分析罗紫洋,安文静,张新东(新疆师范大学数学科学学院,新疆乌鲁木齐８３００１７)摘要:本文研究了一类V o l t e r r a 型微分Ｇ积分方程的数值格式构造及其理论分析．格式构造方面,利用有限差分方法进行时间和空间离散,对于积分项采用复合梯形求积公式进行处理．最后,给出了数值格式的稳定性分析和误差估计,其误差的收敛阶为Ο(τ＋h ４),其中τ为时间步长,h 为空间步长．关键词:V o l t e r r a 型微分Ｇ积分方程;有限差分;复合梯形求积公式;稳定性;误差分析中图分类号:O ２４１．８２㊀㊀㊀文献标志码:AN u m e r i c a l S c h e m eC o n s t r u c t i o na n dT h e o r e t i c a lA n a l y s i s f o r I n t e g r a l ＧD i f f e r e n t i a l E q u a t i o n s o fV o l t e r r aT y pe L U OZ i Ｇy a n g ,A N W e n Ｇj i n g ,Z HA N G X i n Ｇd o n g(S c h o o l o fM a t h e m a t i c a l S c i e n c e s ,X i n j i a n g N o r m a lU n i v e r s i t y ,U r u m qi ８３００１７,C h i n a )A b s t r a c t :I n t h i s p a p e r ,an u m e r i c a l s c h e m ec o n s t r u c t i o na n dt h e o r e t i c a l a n a l y s i sw a sm a i n l yd e v e l o p e d f o rs o l v i n g t h e V o l t e r r ai n t e g r a l Ｇd i f f e r e n t i a le q u a t i o n s ．I nt e r m so fs c h e m ec o n Ｇs t r u c t i o n ,t h e f i n i t ed i f f e r e n c ea p p r o x i m a t i o ni s s u e df o rt i m ea n ds pa c ed e r i v a t i v e ,a n dt h e c o m p o u n d t r a p e z o i d a l q u a d r a t u r e f o r m u l a i s s u e d f o r i n t e g r a l t e r m．F i n a l l y,t h eu n c o n d i t i o n a l s t a b i l i t y a n d c o n v e r g e n c ew e r e g i v e n ．I na d d i t i o n ,t h ee r r o r a n a l y s i sw a s c a r r i e do u t ,w h i c h s h o w e d t h a t t h e c o n v e r g e n c e o r d e rw a s O (ι＋h ４),w h e r e ιw a s t h e t i m e s t e p ,a n d h w a s t h e s p a c e s t e p．K e y wo r d s :V o l t e r r ai n t e g r a l Ｇd i f f e r e n t i a le q u a t i o n ;f i n i t ed i f f e r e n c e ;c o m p o u n dt r a p e z o i d a l q u a d r a t u r e f o r m u l a ;s t a b i l i t y ;e r r o r a n a l ys i s 0㊀引言V o l t e r r a 型微分Ｇ积分方程(V I D E s )源于１９世纪末２０世纪初,是在竞争或人口增长模型中所引入的一类具有积分内核的方程．V I D E s 中的积分项具有记忆性质[１],这一性质是其在物理领域有着广泛应用的重要原因．例如:粘弹性方程[２]㊁具有记忆性的热传导方程[３Ｇ４],以及核反应堆中的热交换过程等．由于求解此类方程的解析解较为困难,因此研究此类方程的数值解具有重要的理论价值和实际意义．处理微分Ｇ积分方程的主要方法为差分法．１８９７年,V o l t e r r a [５]在其著作中对该类方程的数值解法已有所研究．１９７４年,B r u n n e r [６]采用隐式R u n g e ＧK u t t a 方法求解V o l t e r r a 型微分Ｇ积分方程,并分析了此类方程数值解的稳定性．１９９２年,陈传淼等[７]利用内积近似一类微分Ｇ积分方程的积分项,并对空间项采用有限元方法,得到了相应的误差估计．１９９３年,汤涛[８]运用梯形求积技巧构造了一类微分Ｇ积分方程中积分项的数值格式,使其时间误差收敛阶可以达到Ο(τ３/２),并给出具体的理论分析．２００５年,S h a h m o r a d [９]基于T a u方法分析了线性F r e d h o l m ＧV o l t e r r a 型微分Ｇ积分方程的有效误差．２００６年,A m i r a l i ye v 等[１０]在均匀网格中构建V o l t e r r a 型微分Ｇ积分方程的数值格式,并证明该格式的时间收敛阶为一阶均匀收敛．２００８年,徐大[１１]在希尔伯特空间中分析验证了有限差分方法求解线性V o l t e r r a 型方程的稳定性;同年,T a r i 和S h a h m o r a d [１２]基于勒让德多项式研究二维线性F r e d h o l m 积分方程,并给出了数值格式的误差边界．２０１０年,W a z w a z [１３]提出结合L a p l a c e 变换和A d o m i a n 分解法,构造非线性V o l t e r r a 型微分Ｇ积分方程的数值格式,并通过实验证明该方法的有效性;与此同时,Z a r e b Ｇn i a [１４]利用S i n c 方法求解V o l t e r r a 型微分Ｇ积分方程,并验证该方法所具有的优良性,不仅收敛速度较快,且不存在当使用其他数值方法时常见的不稳定问题．２０２１年,C i m e n 和C a k i r [１５]基于正交规则和积分恒等式构造了一类F r e d h o l m 型微分Ｇ积分方程一种新的差分格式,并分析证明该格式的稳定性和收敛性．２０２２年,S a n t r a 和M o h a Ｇpa t r a [１６]利用复合梯形公式逼近V o l t e r r a 型微分Ｇ积分方程中的积分项,使其误差时间收敛阶为Ο(τ)．本文考虑如下V o l t e r r a 型微分Ｇ积分方程(V I D E ):∂u (x ,t )∂t －∂２u (x ,t )∂x ２＝ʏt０K (x ,t －s )u (x ,s )d s ＋f (x ,t ),(x ,t )ɪΩ;u (x ,０)＝φ(x ),∀x ɪ[０,L ];u (０,t )＝g １(t ),u (L ,t )＝g ２(t ),∀t ɪ(０,T ],ìîíïïïïïïïïïï(１)其中:Ω＝[０,L ]ˑ(０,T ],空间L ＞０,且时间T ＞０;积分核K (x ,t －s )在Ω上光滑且有界;f ,φ,g １和g ２为给定的光滑函数．1㊀紧致差分格式推导为方便起见,本文中不同地方的常数C 可以代表不同的数值．给定两个正整数M 和N ,定义h＝L/M 为空间变量x 的步长,τ＝T /N 为时间变量t 的步长,所以x j ＝j h ,j ＝０,１, ,M ,t n ＝n τ,n ＝０,１, ,N ．用u n j 和U nj 分别表示函数u 在点(x j ,t n )处的精确解和数值解．为了方便表达及书写,给出如下算子定义,δ２xU j ＝U j －１－２U j ＋U j ＋１h ２,H U j ＝１１２(U j －１＋１０U j ＋U j ＋１)＝(１＋h ２１２δ２x )U j ,㊀㊀j ＝１,２, ,M －１;U j ,j ＝０,M ．ìîíïïïï(２)由复合梯形求积公式可得ʏbag (t )d t ＝ðn －１i ＝０ʏt i ＋１t ig (t )d t ＝τ２ðn －１i ＝０[g (t i )＋g (t i ＋１)]＋R ,(３)其中:R 为整个区间的截断误差．由T n ＝τ２ðn －１i ＝０[g (t i )＋g (t i ＋１)]和中值定理可得R ＝ʏbag (t )d t －T n ＝－b －a １２τ２g ᵡ(ξ),ξɪ[a ,b ]．(４)引理１[１７]㊀假设s (x )ɪC ６[a ,b ],则有１１２(s x x (x j ＋１)＋１０s x x (x j )＋s x x (x j －１))－１h２(s (x j ＋１)－２s (x j )＋s (x j －１))＝h ４２４０s (６)(ξj ),其中:ξj ɪ(x j －１,x j ＋１),j ＝１,２, ,M －１．接下来对方程(１)进行逐项分析．对于方程(１)的时间项来说,利用T a y l o r 展开,可得∂u (x j ,t n )∂t n＝u (x j ,t n ＋１)－u (x j ,t n )τ＋O (τ)．(５)对于方程(１)的空间项来说,H ∂２u (x j ,t n ＋１)∂x ２j＝δ２x u (x j ,t n ＋１)＋O (h ４)．(６)对于方程(１)的积分项而言,由式(３)和(４)可得１１第４期罗紫洋等:V o l t e r r a 型微分Ｇ积分方程的数值格式构造及理论分析ʏt n ０K (x j ,t n －s )u (x j ,s )d s ＝ðn －１i ＝０ʏt i ＋１t iK (x j,t n－s )u (x j,s )d s ＝τ２ðn －１i ＝０[K (x j ,t n －t i )u (x j ,t i )＋K (x j ,t n －t i ＋１)u (x j ,t i ＋１)]＋O (τ２)．(７)如果对方程(１)作用H 算子,可得H ∂u (x j ,t n )∂t n－H ∂２u (x j ,t n ＋１)∂x ２j ＝Hʏt n ０K (x j ,t n －s )u (x j ,s )d s ＋H f (x j ,t n ),进而由式(５)㊁(６)和式(７)可得H u n ＋１j－τδ２x u n ＋１j ＝H u nj＋H τ２２ðn －１i ＝０(K n －i －１j u i ＋１j ＋K n －i j u i j )éëêêùûúú＋H τf n j ＋R nj ,(８)其中:f nj ＝f (x j ,t n ),R n j ɤC (τ＋h ４)．省去式(８)中的截断误差,当j ＝１,２, ,M－１且n ＝１,２, ,N 时得到方程(１)的离散格式如下:H U n ＋１j －τδ２x U n ＋１j ＝H U nj ＋㊀H τ２２ðn －１i ＝０(K n －i －１j U i ＋１j ＋K n －i j U i j )éëêêùûúú＋H τf n j ,U ０j ＝φ(x ),U n ＋１０＝g １(t n ＋１),U n ＋１M ＝g ２(t n ＋１)．ìîíïïïïïï(９)2㊀稳定性分析本节将利用能量不等式的方法讨论数值格式的稳定性,并证明数值格式为无条件稳定．首先给出一些符号说明和引理．定义V h 为空间网格函数,V h ＝{V |V ＝(V ０,V １, ,V M ),V ０＝g １(t n ＋１),V M ＝g ２(t n ＋１)}．以下内积和范数在证明中将用到,(U ,V )＝h ðM －１j ＝１U jV j , U ２＝(U ,U )．引理１[１８]㊀假设U ɪV h ,则有－(δ２xU n H U n )ȡ２３δx U n ２．定理１㊀紧致差分格式(９)是无条件稳定的,并且满足如下不等式:H U n ＋１ ɤC H U ０ ＋m a x ０ɤp ɤnH f p (),其中:C 为常数．证明㊀设m a x (x ,t )ɪΩK (x ,t ) ＝K x t ．由式(９)可得H U n ＋１－τδ２xU n ＋１＝H U n ＋H τ２２ðn －１i ＝０(K n －i －１U i ＋１＋K n －i U i )éëêêùûúú＋H τf n ,(１０)对式(１０)两边同乘H U n ＋１,并在Ω上积分可得(H U n ＋１,H U n ＋１)－τ(δ２xU n ＋１,H U n ＋１)＝(H U n ,H U n ＋１)＋τ２２ðn －１i ＝０(H (K n －i －１U i ＋１),H U n ＋１)＋τ２２ðn －１i ＝０(H (K n －i U i ),H U n ＋１)＋τ(H f n ,H U n ＋１)．(１１)由引理１可知－(δ２xU n H U n )ȡ０．由式(１１)可得如下不等式,(H U n ＋１,H U n ＋１)ɤ(H U n ,H U n ＋１)＋τ２２ðn －１i ＝０(H (K n －i －１U i ＋１),H U n ＋１)＋τ２２ðn －１i ＝０(H (K n －i U i ,H U n ＋１))＋τ(H f n ,H U n ＋１)．(１２)接下来将用归纳法证明．由式(１２)可知,当n ＝０时,(H U １,H U １)ɤ(H U ０,H U １)＋τ(H f ０,H U １)．再由柯西Ｇ施瓦茨不等式可得H U １ ɤ H U ０ ＋τ H f ０ ɤC H U ０ ＋ H f ０ ()．假设对任意k ,k ɤn 时结论成立,即 H U k ɤC H U ０ ＋m a x ０ɤp ɤnH f p ()．(１３)由式(１２),当k ＝n ＋１时,有(H U n ＋１,H U n ＋１)ɤ(H U n ,H U n ＋１)＋τ２２ðn －１i ＝０(H (K n －i －１U i ＋１),H U n ＋１)＋τ２２ðn －１i ＝０(H (K n －i U i ),H U n ＋１)＋τ(H f n ,H U n ＋１)．由上述不等式,结合式(１３)可得 H U n ＋１ ɤ２１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀兰州文理学院学报(自然科学版)㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀第３６卷H U n＋τ２K x t２ðn －１i ＝０H U i ＋１ ＋τ２K x t２ðn －１i ＝０H U i ＋τ H f n ɤC H U ０＋m a x ０ɤp ɤn －１H f p()＋τK x tT C H U ０＋m a x ０ɤp ɤn －１H f p ()＋τ H f nɤ(C ＋τK x tT C ) H U ０ ＋(C ＋τK x tT C ＋τ)m a x ０ɤp ɤnH f p ɤC H U ０＋m a x ０ɤp ɤnH f p(),其中,C 为常数．定理１证明完毕．3㊀误差估计本节将讨论数值格式的收敛性,并推导出数值格式的误差收敛阶为O (τ＋h ４)．记e nj ＝u n j －U n j ,j ＝０,１, ,M 且n ＝０,１, ,N ．定理２㊀假设u nj 为方程(１)的精确解,U n j 为离散格式(９)的数值解,则如下不等式成立, H e n ＋１jɤC (τ＋h ４),其中,C 为常数．证明㊀设m a x (x ,t )ɪΩK (x ,t ) ＝K x t ．用式(８)减去式(９),再由e nj 的定义可得H e n ＋１j －τδ２x e n ＋１j＝H e n j＋H [τ２２ðn －１i ＝０(K n －i －１j e i ＋１j ＋K n －i j e i j )]＋R n j ．(１４)将式(１４)两边同乘H e n ＋１j,并在Ω上积分可得(H e n ＋１j ,H e n ＋１j )－τ(δ２x e n ＋１j ,H e n ＋１j)＝(H e n j ,H e n ＋１j)＋τ２２ðn －１i ＝０(H (K n －i －１j e i ＋１j ),H e n ＋１j )＋τ２２ðn －１i ＝０(H (K n －i j e i j ),H e n ＋１j )＋(R n j ,H e n ＋１j)．(１５)由引理１和式(１５),有(H e n ＋１j ,H e n ＋１j )ɤ(H e n j ,H e n ＋１j)＋τ２２ðn －１i ＝０(H (K n －i －１j e i ＋１j ),H e n ＋１j )＋τ２２ðn －１i ＝０(H (K n －i j e i j ),H e n ＋１j )＋(R n j ,H e n ＋１j)．(１６)类似定理１的证明,将采用归纳法．由式(１６)可知,当n ＝０时,可得(H e １j ,H e １j )ɤ(H e ０j ,H e １j )＋(R ０j ,H e １j ),再由柯西Ｇ施瓦茨不等式及R ０j ɤC (τ＋h ４)可得 H e １j ɤ H e ０j ＋R ０j ɤC (τ＋h ４)．假设对任意k ,k ɤn 时,结论成立,即 H e kj ɤC (τ＋h ４)．(１７)由式(１６),当k ＝n ＋１时,有(H e n ＋１j ,H e n ＋１j )ɤ(H e n j ,H e n ＋１j)＋τ２２ðn －１i ＝０(H (K n －i －１j e i ＋１j ),H e n ＋１j )＋τ２２ðn －１i ＝０(H (K n －i j e i j ),H e n ＋１j )＋(R n j ,H e n ＋１j)．由上述不等式,结合式(１７)可得 H e n ＋１jɤ H e nj＋τ２K x t ２ðn －１i ＝０H e i ＋１j ＋τ２K x t ２ðn －１i ＝０H e ij ＋ R n j ɤC (τ＋h ４)＋τK x tT (C (τ＋h ４))＋C (τ＋h ４)ɤm a x {C ,τK x tT C ,C }(τ＋h ４)ɤC (τ＋h ４),其中,C 为常数．定理２证明完毕．4㊀结论本文通过有限差分方法构造了一类V o l t e r r a型微分Ｇ积分方程的数值格式,分析了格式的稳定性和误差估计,得到其误差收敛阶为O (τ＋h ４)．相较于解决该问题原有的数值格式,在空间收敛阶上有进一步的提升．在今后的研究工作中,将以构建收敛速度更快,误差更低的数值格式为目标,并将此方法应用到分数阶微分Ｇ积分方程的数值求解．参考文献:[１]MO L F E R S D O R F L V．O ni d e n t i f i c a t i o no fm e m o r yk e r n e l s i n l i n e r t h e o r y of h e a t c o n d u c t i o n [J ]．M a t h e Ｇm a t i c a lM e t h o d s i n A p pl i e dS c i e n c e s ,１９９４,１７:９１９Ｇ９３２．[２]R C N A R D Y M．M a t h m e a t i c a l a n a l ys i so fv i s c o e l a s t i c f l o w s [J ]．A n n u a lR e v i e w o fF l u i d M e c h a n i c ,１９８９,２１:２１Ｇ３６．[３]G U R T I N M E ,P I P K I N AC ．A g e n e r a l t h e o r y of h e a t c o n d u c t i o nw i t h f i n i t ew a v e s pe e d [J ]．A r c h i v ef o rR a Ｇt i o n a lM e c h a n i c s a n dA n a l y s i s ,１９６８,３１:１１３Ｇ１２６．[４]M I L L E R R K．A ni n t eg r o Ｇd i f f e r e n t i a le qu a t i o nf o r g i r dh e a tc o n d u c t o r s w i t h m e m o r y [J ]．J o u r n a lo f ３１第４期罗紫洋等:V o l t e r r a 型微分Ｇ积分方程的数值格式构造及理论分析M a t h e m a t i c a l A n a l y s i sa n d A p p l i c a t i o n s,１９７８,６６:３１３Ｇ３３２．[５]V O L T E R R A V．S o p r aa l c u n e q u e s t i o n i d i i n v e r s i o n e d i i n t e g r a l i d e f i n i t i[J]．A n n a l i d iM a t e m a t i c aP u r ae d A p p l i c a t a,１８９７,２５:１３９Ｇ１７８．[６]B R U N N E R H．I m p l i c i t R u n g eＧK u t t am e t h o d s o f o p t iＧm a lo r d e rf o r V o l t e r r ai n t e g r oＧd i f f e r e n t i a le q u a t i o n s [J]．M a t h e m a t i c s o fC o m p u t a t i o n,１９８４,４２:９５Ｇ１０９．[７]C H E NC M,T HOM E E V,WA H L B I NLB．F i n t e e lＧe m e n t a p p r o x i m a t i o n o f a p a r a b o l i c i n t e g r oＧd i f f e r e n t iＧa l e q u a t i o nw i t haw e a k l y s i n g u l a rk e r n e l[J]．M a t h eＧm a t i c a l o fC o m p u t a t i o n,１９９２,５８:５８７Ｇ６０２．[８]T A N G T．Af i n i t ed i f f e r e n c es h e m ef o r p a r t i a l i n t eＧg r oＧd i f f e r e n t i a l e q u a t i o n sw i t haw e a k l y s i n g u l a rk e rＧn e l[J]．A p p l i e d N u m b e r i c a l M a t h e m a t i c s,１９９３,１１:３０９Ｇ３１９．[９]S HA HMO R A DS．N u m e r i c a l s o l u t i o no f t h e g e n e r a l f o r ml i n e a rF r e d h o l mＧV o l t e r r a i n t e g r oＧd i f f e r e n t i a l eＧq u a t i o n sb y t h eT a um e t h o dw i t ha ne r r o r e s t i m a t i o n [J]．A p p l i e d M a t h e m a t i c sa n d C o m p u t a t i o n,２００５,１６７:１４１８Ｇ１４２９．[１０]AM I R A L I Y E V G M,S E V G I N S．U n i f o r m d i f f e rＧe n c em e t h o df o r s i ng u l a r l yp e r t u r b e dV o l t e r r a i n t eＧg r oＧd i f f e r e n t i a l e q u a t i o n s[J]．A p p l i e d M a t h e m a t i c sa n dC o m p u t a t i o n,２００６,１７９:７３１Ｇ７４１．[１１]X U D．S t a b i l i t y o f t h ed i f f e r e n c et y p e m e t h o d sf o r l i n e a rV o l t e r r ae q u a t i o ni n H i l b e r ts p a c e s[J]．N uＧm e r i s c h eM a t h e m a t i k,２００８,１０９:５７１Ｇ５９５．[１２]T A R IA,S HA HMO R A DS．Ac o m p u t a t i o n a lm e t hＧo d f o r s o l v i n g t w oＧd i m e n s i o n a l l i n e a rF r e d h o l mi nＧt e g r a le q u a t i o n so fs e c o n d k i n d[J]．T h e A n z i a mJ o u r n a l,２００８,４９:５４３Ｇ５４９．[１３]WA Z WA Z A M．T h ec o m b i n e dL a p l a c et r a n s f o r mＧA d o m i a nd e c o m p o s i t i o n m e t h o df o rh a n d l i n g n o nＧl i n e a r V o l t e r r ai n t e g r oＧd i f f e r e n t i a le q u a t i o n s[J]．A p p l i e d M a t h e m a t i c sa n dC o m p u t a t i o n,２０１０,２１６:１３０４Ｇ１３０９．[１４]Z A R E B N I A M．S i n c n u m e r i c a l s o l u t i o n f o r t h eV o l tＧe r r ai n t e g r oＧd if f e r e n t i a le q u a t i o n[J]．C o mm u n i c aＧt i o n s i n N o n l i n e a rS c i e n c ea n d N u m e r i c a lS i m u l aＧt i o n,２０１０,１５:７００Ｇ７０６．[１５]C I M E N E,C A K I R M．Au n i f o r m n u m e r i c a lm e t h o df o rs o l v i ng s i n g u l a r l yp e r t u r b e d F r e dh o l mi n t e g r oＧd i f fe r e n t i a l p r o b l e m[J]．C o m p u t a t i o n a l a n dA p p l i e dM a t h e m a t i c s,２０２１,４２:１Ｇ１４．[１６]S A N T R AS,MOHA P A T R AJ．An o v e l f i n i t e d i f f e rＧe n c e t e c h n i q u ew i t h e r r o r e s t i m a t ef o r t i m e f r a c t i o nＧa l p a r t i a li n t e g r oＧd i f f e r e n t i a le q u a t i o n o f V o l t e r r at y p e[J]．J o u r n a l o f C o m p u t a t i o n a la n d A p p l i e dM a t h e m a t i c s,２０２２,４００:１１３７４６．[１７]L U O M,X U D,L I L M．A c o m p a c td i f f e r e n c e s c h e m ef o r a p a r t i a li n t e g r oＧd i f f e r e n t i a le q u a t i o nw i t ha w e a k l y s i n g u l a rk e r n e l[J]．A p p l i e d M a t h eＧm a t i c a lM o d e l l i n g,２０１５,３９:９４７Ｇ９５４．[１８]A K B A R M．C o m p a c t f i n i t e d i f f e r e n c e s c h e m e f o r t h e s o l u t i o n f o r a t i m e f a r a c t i o n a l p a r t i a l i n t e g r oＧd i f f e rＧe n t i a l e q u a t i o n w i t ha w e a k l y s i n g u l a rk e r n e l[J]．M a t h e m a t i c a l M e t h o d si n t h e A p p l i e d S c i e n c e s,２０１７,４０:７６２７Ｇ７６３９．[责任编辑:赵慧霞]４１㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀兰州文理学院学报(自然科学版)㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀第３６卷。

奇摄动Volterra型积分微分方程非线性边值问题
吴钦宽;林平健;孙福树;尤兴华
【期刊名称】《自然科学进展》
【年(卷),期】2005(015)003
【摘要】利用微分不等式理论,研究了一类Volterra型积分微分方程非线性边值问题.在适当条件下,构造出问题的上下解,得出解的存在性和渐近估计.
【总页数】5页(P375-379)
【作者】吴钦宽;林平健;孙福树;尤兴华
【作者单位】南京工程学院基础部,南京,210013;南京工程学院基础部,南
京,210013;南京工程学院基础部,南京,210013;南京工程学院基础部,南京,210013【正文语种】中文
【中图分类】O175
【相关文献】
1.某一类型积分微分方程非线性边值问题的奇摄动 [J], 金丽;王国灿
2.奇摄动Volterra型积分微分方程的非线性边值问题 [J], 金丽;王国灿
3.二阶Volterra型积分微分方程非线性边值问题的奇摄动 [J], 王国灿;金丽
4.二阶Volterra 型积分微分方程奇摄动非线性边值问题解的惟一性 [J], 王国灿;丁传华
5.奇摄动Volterra型积分微分方程非线性边值问题的渐近估计 [J], 王国灿;金丽因版权原因，仅展示原文概要，查看原文内容请购买。

第四章Volterra模型在本章中,我们首先回顾最常用的非线性建模方法.第 4.1.2节讨论多项式Volterra模型建模方法的特性,第4.2节详细解说记录下来的完全电非线性和电热非线性,以及建模非线性级数展开中所用的项.4.3节描述了如何应用Volterra分析方法计算一个普通发射放大器的失真,并且作为第一个研究案例,第 4.4介绍了在一个BJT CE放大器中,影响IM3失真的所有项的分析.同样的分析(对一个MESFET放大器)在第4.5节中介绍.4.1非线性建模为了能够分析功放的非线性行为,我们需要为实际的非线性电路,无源匹配以及偏置元件建立精确的模型.但此模型的获取又有一定难度.众所周知,N阶失真的数量与I-V和Q-V波形的N阶导数成比例(参照[1-3]).因此,为了达到精确的失真模拟,对有源元件的I-V和 Q-V波必须建模,因此不仅是直流值,高阶导数同样是正确且连续的.(为了便于参考,在早期的仿真模型中,第一阶导数可以不连续).此外,电容是很容易建模的,因此电荷没有存储,这将导致完全容性结点的非物理整流和自偏压.因此,特别是如果电容值同时取决于两个终端电压,将电容模拟等效于电荷平衡是十分重要的[5].由于无源元件在高频具有分布特性,很难在射频频率上对其建模.有损耗的传输线难以在时域上进行建模,一般而言,无源元件的建模在频域上更加精确.尽管如此,有些无源元件的频域仿真模型在高次谐波时也可能是不准确的,例如传输线宽上的阶跃变化,在电路分析中,可以采用测试电路的标准终端阻抗值来进行分析.简而言之,为了得到精确的失真模拟,从模拟模型中需要:1. 对N阶失真模拟来说, I-V 和Q-V波形的N阶导数必须足够精确.2. 结点阻抗的频率响应对于最高相关谐波,必须是正确的.同时,在基带频率,偏置阻抗和热阻抗的正确模型也是十分需要的.3. 如果可以得到主要失真源的组成信息,将十分有用.4.1.1非线性仿真模型通常来说,功放和发射机设计师使用两种非线性模型:一种是用于系统仿真的行为黑盒子模型,另一种是用于电路仿真的器件模型.根据建模方法的不同,可以对这两类模型进行更深层分割:可以是解析的,基于一些预先确定的和物理学的用参数表示的模型函数，或者完全根据实验的,将测量数据列表并以内插值替换的,或者用简单的曲线或物理意义不清楚的多项式表示的模型.在表4.1中有所介绍.表4.1 功率放大器的非线性模型行为基带模型广泛用于模拟和优化整个发射机和收发器,并且增加了新的功能,例如模拟记忆效应,在[6]中介绍.然而,行为模型描述的或者是一个已经存在的放大器,或者行为模型源自到目前为止不存在的放大器规范说明,但是在设计一个新的功率放大器中,行为模型的使用遭到限制.这里简明地介绍了最常用的行为模型的特性,仅用于参考.简单的静态的AM-AM 和 AM-PM波形不能够模拟记忆效应,但是基于调制频率的AM-AM和AM-PM波形的模型已经被开发出来.如图4.1所示,在Saleh模型中, AM-AM和AM-PM非线性模块的输入和输出端都增加了线性滤波器.在Blum 和Jeruchim(在[7]中描述)模型中,用快速傅立叶算法及足够的抽样来找到用于修改AM-AM表的瞬时调制频率.有一种Volterra型的行为模型被称作Volterra输入输出图(VIOMAP).它是普通S参数的非线性的概念性扩展,包括谐波响应,并且被成功应用于单音负载下拉仿真中[9,10].图4.1 (a)功率放大器的AM-AM和AM-PM波形(b)基于滤波器和无记忆非线性的由频率决定的非线性模型.器件模型描述了半导体设备的动作,以及无源和分布式元件的合适模型,可以建立并优化功放的模型.早期的半导体模型是解析的,所采用的等式首先来源于半导体物理学,然后将其简化以减少仿真时间.这些基于等式的模型的一个基本问题是,所选用的函数和控制参数固定了I-V 和Q-V特性的可能形状,并且可能没有足够的自由度来模拟例如I C-V CE曲率.例如,在基本的Gummel- Poon (GP) BJT SPICE模型中,集电极电流的简化形式如下:其中, 基本的指数仅可被三个控制参数修改:IS依比例决定电流, VAF (所谓的早期电压[12], 如图4.2所示)构成输出电导的一个极其简化的模型, IKF (所谓的拐点电流)降低高电流时的增益[13, 14].这个简单的等式涵盖了整个I-V平面,同时固定了导数dn I C/d V n,因而固定了非线性行为. SPICE GP模型可以适当地用于模拟基带信号[15],但是特别是对于过分简单化和固有性线的输出阻抗模型, 不能用SPICE GP模型进行精确的失真仿真,这将在本书后面说明.更好的物理模型已经被开发出来,比如BJT的Mextram和VBIC ,以及用于LDMOS 的摩托罗拉MET模型,这都是久经考验的模型.与早期的SPICE模型其比,这些模型的性能大大地提高了.后者对于找到正确的直流偏置十分重要,这是因为自我加热使得I-V波形产生一个大的差值.假如热模型有足够的时间常量来模拟缓慢加热包(主要影响直流偏置和芯片表面的微秒范围热记忆),后面的模型也可用于模拟热记忆效应.图4.2 在BJT中使用早期电压VAF模拟输出阻抗模型的额外自由度增加了其复杂性以及控制参数的数量.在一个极限中,MOS BSIM模型有数十个参数来单独控制比例特性.因此,模型的复杂性趋于失控,并且其配置愈加复杂,对错误愈加敏感.另一种设备建模的方法是放弃等式,而采用列表的测试数据或者完全根据实际以验的函数来代替.现在,任何形式的I-V和Q-V特性的模型可以被建立,这是通过Root模型得到的方法,称作”设备最了解”模型[5].在内插列表数据时,存在一些技术问题,这是因为内插的多项式容易使数据点之间产生振荡,因此派生出高阶非物理波动.然而,由于预定函数不需要压力,列表模型使用灵活.Volterra模型是一种经验模型,它不依赖于半导体物理学.其非线性描述为多项式,系数可以通过对I-V和Q-V函数微分得到,也可通过将多项式直接填入测量数据表里得到.在此我们使用的是后者,在接下来的章节中,我们将会更深层次地对Volterra模型的特性进行研究.4.1.2 Volterra模型的特性多项式模型并不自动地对模拟快速响应,相反,它可能严重地会聚在高于原始设置范围的信号电平上.然而,多项式模型允许使用高效的Volterra分析程序.然而,采用Volterra模拟方法的主要动机并不是看中了其速度优势,而是它能提供一个极好的分析工具来进行分析.主要的失真机制可以用与在普通交流噪声分析中采用的相同的方法来进行分析,由于非线性分析,多重的混合机制同样可被识别,例如可以帮助谐波终端阻抗的设计等.因此, Volterra分析是少数可以帮助理解记忆效应和帮助设计优化的方法之一.尽管如此,仍需承认多项式模型存在一些缺点.首先,多项式模型要遭受在适宜的带宽范围外,其响应接近无穷大这一事实.传统的非线性建模函数与此正好相反,它在整个偏置范围内平滑,有限的表现是我们设计的特性,因为它帮助使信号收敛,并且信号摆动不必要进行推理的了解.因此, Volterra分析并不是一个非常普通的工具.由于速度原因, Volterra分析被用于快速失真分析和模拟器中的低噪放型小信号电路(例如Voltaire XL [22]和SPICE的早期版本)中,或者甚至做为独立的模拟器使用[23].然而,为了功放能被成功地进行分析,可靠的早期信息仍是必需的.第二,预先得需要实际的大信号直流偏置电压.大信号动作经常会引起直流工作点的移位,它同时影响增益和非线性的数量.此信号感应引起的直流移位会减缓谐波平衡模拟中收敛的速度,并且在非反复Volterra计算程序中,只能对它进行估计,而非完全地预测.为了克服这点,我们需要检查直流移位是不是很大,或者在实际的大信号工作点使得多项式模型合适.第三,在多项式函数的装配中,需要知道输入和输出电压摆动的范围.多项式模型的实际功率是除开失真成分的其它部分.适宜的范围越大,低阶多项式的精确度越小.因此,沿着最大信号振幅安装是合适的,为确保多项式模型的精度,不能太超过此范围,并且由于多项式响应可能在合适的范围外是完全非物理的,也不能超过一个较小的范围.在这种情况下,就需要对输入输出轨道有较好地评估.总之, Volterra分板并不是一种简单使用的独立的仿真方法,但当它用在与其它仿真方法(例如谐波平衡)并联使用时,此方法提供了更多调试功能.在本书中,研究的案例仅限于单级晶体管放大器,并且对Volterra分析半解析地计算(象征性地来源于每一个失真源到所有结点电压的转移函数). 象征性的分析并不是必需的,它限制了对CE或CS放大器固定的结构以及双音测试信号的分析.而Volterra分析却能达到,通过在(用标准修改结分析矩阵以及非线性电流源表示的)任何电路上运用通常非线性交流分析方法.4.2非线性I-V和Q-V特性大多数的晶体管模型是以Π模型或T模型为基础.这里使用的是Π模型,本节中描述了用BJT,异质结BJT(HBT)和场效应管(FET)的Π模型表示的典型传导(I-V)和电容(Q-V)的非线性特性.在这里将BJT作为一个案例,但同样的模型也可用于FET晶体管,只是多项式系数设置不同.前面已经讲过, Volterra模型是以I-V和Q-V曲线的多项式建模为基础.测量这些曲线也许会有些困难,详见第五章.这样的电荷不能直接进行测量,我们必须依靠交流测量的电容以及对所得电容值得到的电荷等式求积分.用类似的方法,I-V曲线可以通过由S参数测量得到的和值进行大部分重造,但是实际的I-V曲线是一个更安全的出发点.这里介绍的模型是电热模型,这意味着其结温是一自由变量.然而,直流温度上升包含在偏置点中,并且只考虑由动态自我加热引起的温度变化.由于功耗是电压和电流的产物,我们认为结温中的交流成分已经是一个二阶现象.因此,一个三阶的模型仅仅包括温度的一次方,这意味着电容元件的温度依靠性被认为是线性的.4.2.1 特性在大多数被报导的BJT/HBT Volterra级数分析中,集电极电流被认为只是基极电压的函数[25-27],此考虑方法抓住了主要指数的输入输出非线性但是认为输出电导是常数.在MESFET Volterra级数分析中,漏电压的效应通常用的一个多项式来实现, (参照[28]),但即使如此,也难以捕捉所有的非线性特性.在(4.1)式中,等式是,和结温T的一个三维简单函数，就像.通过扩展大信号I-V函数到一个三输入的泰勒级数(在直流工作点,和周围),很容易得到一个多项式模型.因此,交流电流的电热三阶级数展开可以写作:其中, ,,并且K ixxx是元素xxx的i阶非线性系数.( 可以用来标志).由等式可见,第一行只受影响,第二行只受影响(例如非线性输出电导).尽管如此,第三行又列出了和的向量积.最后,第四行列出了与两个终端电压混合在一起的温度变化.图4.3中证明了电非线性的影响,其中,绘制了在三个基极电压处的集电极电流,它是三个不同基极电压处集电极电压的函数.如果除外的所有系数都是零,我们可以得到如图4.3(a)所示的三条等间距的水平线.由于线精确地保持水平,输出电导为零,且集电极电压不影响电流数量.此外,由于线之间等间距,跨导是线性的.然而,如果或偏离了零,在I-V平面的线距离将会变得不等,这表明跨导是非线性的.图4.3 证明集电级电流非线性.垂直的轴是集电极电流,水平轴是电压.(a)线性响应(b)非零 (c)非零 (d) 非零 (e) 非零 (f)非零图4.3(b)证明了的影响,仅仅存在和.与图4.3(a)相比,图4.3(b)中的线有一个非零的斜率,它与成比例且不依赖与.图4.3(b)仍是完全线性的, 图4.3(c)证实了输出电导的非线性,其电流的斜率随而变化.在这种情况下,仅有一个非零值,和可被用于模拟输出电导的曲线效应，例如饱和和击穿.图4.3(d-f)分别图解说明了,和的截项,它模拟了基极和集电极非线性物性的交互作用.为了帮助比较,图4.3(d-f)中的稍细线是临摹图(c)的( 和都有非零值).(对应于项)在图4.3(d)中是非零的,这是由于其线的斜率变化不仅仅受图4.3(c)中集电极电压的作用,也受基极电压的作用.这对于图4.2所示的模拟早期效应是十分必需的.相似的推理也可应用于和,如图4.3(e, f)所示,确定输出电导的形状,分别是和的函数.图4.4对不同建模方法的I-V曲线做了更多的比较.如果集电极电流被模拟为基极电压和线性的一维函数，就产生了一条直的I-V线,如图4.4中细虚线所示.用SPICE Gummel-Poon模型模拟的I-V曲线同样也是直的细的线,但是其斜率和输出电导随集电极电流变化,如图4.2所示.实际上,由于准饱和和截止影响,在大信号或半大信号的情况下,BJT的I-V曲线决不是直线.曲率可以通过使用和的一维多项式模拟,如图4.4中粗实线所示.然而,初步的现象(例如早期效应)在没有引起斜率决定于的值截项时不能被模拟.如图中粗虚线所示,并对应于(4.3)的完全级数展开.饱和和截止的开关同样取决于基极电压,这使得截项的使用强制地避免了I-V平面角落处的重大误差,如图4.4所示.图4.4 三个Volterra模型和Gummel-Poon SPICE模型的I-V特性最后讨论了集电极电流的电热效应,如(4.3)中最后三项所示.在图4.5(a, b)中,描述了一个二次项,它模拟了电流中由温度决定的移位.值得注意的是,不能由小信号参数和得到,相反的,需要实际的电流测量.是一个包含温度和集电极电压影响的三次项,如图4.5(c)所示.它在本质上模拟了温度对输出电导的依靠性.同样的,如图4.5(d)所示,模拟了温度和基极电压的结合效应.由于曲线的斜率反应了跨导，可被认为是受温度影响的跨导的一个变化.图4.5 电热非线性系数的影响. 非零在(a) 轴(b)轴.(c)非零影响(d)非零影响特性是FET型晶体管中唯一重要的传导性非线性.在BJT中,存在另外的两个非线性:由指数引起的非线性和非线性.电导的影响通常更重大,并且也容易模拟.理论上, 等式可以粗略地用电流增益β除，但可用一些方式简化.由于基极电流并不是强烈地依赖于集电极电压,我们可以用一个仅由和决定的两维模型:这里,系数与之前有着相似的含义.线性项通过进行模拟,和模拟其指数曲率.此外, 模拟了由自我加热产生的移位,可以看成线性项的温度依靠.本征基极电阻较难模拟.它是内部基极点和外在基极点之间的串联电阻,它的值取决于基极区的电流拥挤,也取决于内部值.因此,它可以模拟为一个由电阻器电压(),内部基极电压和结温控制的三维电导.所有的K项(k=1,2,...)都是零,但是电流拥挤效应是用和间的截项模拟的,如(4.5)所示.不论如何,通常较小,且在下面例子中被模拟为一个线性电导.4.2.3电容模型如前面解释,将电容建模成多项式电荷,然后将其关于时间进行区分以得到位移电流.电荷可能并且经常是由多于一个的端口电压所控制,这使得我们必须使用一个类似于(4.3)的多维多项式.电荷同样可以模拟为一个电容,在这种情况下,电荷不出现在控制结点之间,而出现在一些其它的结点间.在下面的例子中,只假定了一个控制电压,式(4.6)描述了基极到发射极的电荷,它是基极到发射极的电压以及温度的函数.从这个等式中,可以轻松地得到对应测量电荷C pi和非线性电流源.只需将电荷等式(4.6)分别关于和时间进行区分即可.在(4.8)式中, ω仅是失真音调的频率;因此,电容并不会引起直流失真电流但是在谐波频率处失真最严重.等式(4.7)指出由温度决定的电荷项K2CPIT不能来源于电荷测量;尽管如此,一个时变的结温可能引起一个与它成比例的电流.另外, (4.6)式中的第一项C pi描述了小信号电容, K2CPI和K3CPI定义了它的有关电的非线性. K3CPIT描述了是控制电压和结温函数的电荷,由于C = d Q/d v, K3CPIT的作用可以看成是电荷值的温度决定性.如(4.6)所示,一个线性的C-V趋势K2CPI引起了二次电荷非线性.同样的,与v2 (K3CPI)成比例的电荷引起立方的非线性.不同类型的电容有着不同的特性,如图4.6(a)所示.如 (4.27)所示,基极-射极电容C pi是指数的,因此它是高度的非线性的.BJT和FET中有偏的P-N结或肖特基结仅是稍微的非线性,它们可以通过增加反偏压进行更深一层的线性化. MOSFET型晶体管有着特殊的栅电容,例如, C GS在门限电压的周围dip.如果MOSFET动作接近于关断,此dip会引起大量的二阶非线性.图4.6 (a)归一化的电容(b) 和的电荷4.3共射极BJT/HBT放大器模型现在我们使用直接的方法来计算一个共射极BJT/HBT放大器的IM3成分,使用第2.5.2节所列出的步骤.这样进行分析:首先为电路建立一个模型,通过一个线性的交流分析建立基本的幅度.然后,通过使用第4.3.2节所示的步骤计算二阶电流电压以及三阶电流电压.4.3.1线性分析图4.7所示的是一个共射极BJT放大器,它包括输入阻抗Z IN(混合匹配网络和偏置电路),基极-射极电导g pi和电容C pi,反馈电容CBC,输出电容CCE,输出电导跨导,负载阻抗Z L以及发射极阻抗Z E.输入和负载阻抗不仅包含匹配网络的阻抗,也包括偏置网络和包寄生的阻抗,Z IN由前级的输出阻抗和本征基极电阻r bb组成,如图4.7所示.图4.7 共射极BJT放大器线性化一阶电路为减少等式的数量,输入电压源可以用诺顿等效源来代替且用下列简化符号通过替换图4.7中所有阻抗并在结点1至3应用基尔霍夫电流定律,可以得到如下矩阵等式.通过使用Cramer法则,可以得到结果以及对的响应.因此,在基极,发射极和集电极的电压如下所示:导纳矩阵的行列式写作:基极到发射极,集电极到发射极的电压分别是和.最后,由于和经常用于计算失真,例如由gm 成分产生的失真,很容易得到它们的比值:线性化分析的目的是为了获得所有非线性元件的基本电压幅度,这样我们就可以继续计算这些元件内部产生的非线性电流.在此之前,我们需要对信号摆动进行一些观察.BJT的指数响应是极端的非线性,并且在没有过度失真的情况下,不能承受高于10到30 mV的信号幅度.这听起来不像一个功率放大器,但是两件事情恰好帮助改变此情况.首先,器件并不是完全指数的,但是当驱动到高注入时,BJT线性化可用(4.1)中的参数IKF来模拟.第二,放大器有一些反馈机制来减小BE结中的信号电平.串联发射极阻抗引起一个线性化的串联反馈,CBC引起一个并联反馈.CBC的作用十分重大,这是由于强烈的电容性反馈降低了基极阻抗,因此也减小了BE电压摆动和从驱动激励产生的失真数量.4.3.2非线性分析在本节中展示了一个CE BJT放大器的非线性模型,并由它得到IM3失真.该电路有三个两输入和一个三输入的I-V和Q-V非线性,通过鉴定7个一阶系数,二阶系数和三阶系数模拟得到,其中18个是电系数,其它9个与动态温度变化相关.最后,此分析将介绍IM3音调,它是一个矢量和,由以下组成:7个由立方电非线性引起的项,21个由级联二次方非线性(修正包络信息向上转换到IM3中)引起的项,21个二次谐波向下混合到IM3的项,最后,5个立方的和24个级联的二阶电热项.这些看起来也许很多,但它描述了产生失真的不同机制的真实幅度;它同样清楚地证明了,只要立方项是解析的,大量的信息就会失去.若电路较大,分析的阶数越高,则需要对数据进行压缩,但是原理还是一样的:我们希望知道,通过从直流或谐波波段混合失真并使总失真最小化(或最小化其记忆效应),能产生多少IM3总数;我们希望知道在这些谐波波段，由失真电压引起了什么样的非线性和实际阻抗.由于项数众多,我们不单个地对这些项进行讨论.完整的分析见附件C,接下来用一些例子讲述计算步骤.4.3.2.1二阶失真电流图4.8所示的电路用于解决二阶响应,图4.7中的线性输入电压在此被短路, 添加二阶失真电流源与所有非线性电路元件并行.以T结束的电流是电热电流,将在第4.3.2.5节中讨论.如同以前, ZIN, ZE和ZL将包寄生,偏置阻抗和匹配阻抗结合在一起.为了计算自我加热效应,瞬时功耗用计算,第3.4节所示的热阻抗用于计算频率ω2–ω1处的瞬时温度波动.可以对不同的电路元件使用不同的温度,但是它们在物理上靠近基极区,此处使用共同的温度.然而对于大型设备,将此设备分成较小的并行设备是十分有利的,这样可以看出不同的温度变化.非线性被模拟为,和温度组成的三维函数,它包括和非线性以及所有高达三阶的交叉项. 和是基极到射极电压和温度的函数,非线性由集电极到基极电压和温度控制.图4.8 含电流源电路的二阶响应表示我们通过计算二阶失真电流开始进行分析.举一个例子,在ω2–ω1处由引起的二阶包络电流是:使用表2.5举另一个例子,由非线性引起的二阶包络电流是:它结合了二阶输入非线性,输出非线性和输入输出交叉项的作用,见I-V模型(4.3).从表2.5中可以看出,相量频率和可能的常量项值的选择取决于音频:例如,的乘积在2ω1处产生一个音调.上述的相量音调被选中,所以能在包络频率ω2–ω1产生失真.计算基频ω1和ω2的相量和用(4.12)-(4.14)式.4.3.2.2跨导倒数转移函数和二阶电压接下来,我们需要把不同结点的失真电流转换为失真电压.此处,我们选用了一个象征性的分析,因此很容易得到从结点X和Y到结点Z之间的非线性电流源的转移函数。

一类第二种非线性Volterra 积分程积分数值解法1前言微分程和积分程都是描述物理问题的重要数学工具，各有优点.相对于某种情况来说，对于某种物理数学问题，积分程对于问题的解决比微分程更加有优势，使对问题的研究更加趋于简单化，在数学上，利用积分形式讨论存在性、唯一性往往比较便，结果也比较完美,所以研究积分程便得越来越有用,日益受到重视. 积分程的发展，始终是与数学物理问题的研究息息相关.一般认为，从积分发展的源头可以追溯到国外的数学家克莱茵的著作《古今数学思想》，该书是被认为第一个清醒的认为应用积分程求解的是Abel.Abel 分别于1833年和1826年发表了两篇有关积分程的文章，但其正式的名称却是由数学家du Bois-Raymond 首次提出的，把该问题的研究正式命名为积分程。

所以最早研究积分程的是Abel,他在1823年从力学问题时首先引出了积分程，并用两种法求出了它的解,第一的积分程便是以Abel 命名的程.该程的形式为：⎰=-baax f dt t x t )()()(ϕ,该程称为广义Abel 程,式中a 的值在(0,1)之间.当a=21时，该式子便成为)()(x f dt tx x x a =-⎰ϕ.在此之前，Laplace 于1782年所提出的求Laplace 反变换问题，当时这个问题就要求解一个积分程.但是Fourier 其实已经求出了一类积分程的反变换，这就说明在早些时候积分程就已经在专业性很针对的情况下得到了研究，实际上也说明了Fourier 在研究反变换问题是就相当于解出了一类积分程.积分程的形成基础是有两位数学家Fredholm 和Volterra 奠定的，积分程主要是研究两类相关的程，由于这两位数学家的突出贡献，所以这两个程被命名为Fredholm程和Volterra程。

后来又有德国数学家D.Hilbert进行了重要的研究，并作出了突出的贡献，由于D.Hilbert领头科学家的研究，所以掀起了一阵研究积分程的热潮，并出现了很多重要的成果，后来该理论又推广到非线性部分。

伴有边界摄动非线性积分微分方程系统的奇摄动
吴钦宽
【期刊名称】《吉林大学学报（理学版）》
【年(卷),期】2009(47)5
【摘要】研究伴有边界摄动的二阶非线性积分微分方程组的奇摄动问题.在适当的条件下,利用对角化技巧证明了解的存在性,并构造了解的渐近展开式,给出了余项的一致有效估计.
【总页数】6页(P881-886)
【作者】吴钦宽
【作者单位】南京工程学院,数学研究所,南京,211167
【正文语种】中文
【中图分类】O175.1
【相关文献】
1.伴有边界摄动非线性系统初值问题的奇摄动 [J], 黄蔚章;陈育森
2.伴有边界摄动的积分微分方程组的奇摄动 [J], 林宗池
3.伴有边界摄动的Volterra型积分微分方程组的奇摄动 [J], 吴钦宽
4.伴有边界摄动的积分微分方程组的奇摄动 [J], 林宗池
5.伴有边界摄动二阶非线性系统的奇摄动 [J], 陈育森
因版权原因，仅展示原文概要，查看原文内容请购买。

第一类Volterra积分方程论文：第一类Volterra积分方程数值方法的研究【中文摘要】第一类Volterra积分方程是很重要的一类积分方程，它是在二十世纪发展并成熟起来的。

物理，力学等领域中的许多实际问题都可以通过转化为第一类Volterra积分方程来求解。

当核函数是连续或具有弱奇性时，通常精确解很难给出。

因此，Tolterra 积分方程的数值解法占有了很重要的地位，通过研究它们有很多有益的分析结果得以实现。

本文正是考虑在数据没有扰动的情况下第一类Volterra积分方程的数值解法。

本文结构如下：第一章主要介绍第一类Volterra积分方程的历史背景，国内外研究现状以及发展趋势。

第二章是一些求解第一类Volterra积分方程的预备理论，包括不适定问题，本文所需要使用的正则化方法：Tikhonov，正则化方法，全变差正则化方法等知识。

第三章研究在数据没有扰动的情况下，求解第一类、Volterra积分方程。

主要利用配置点方法，包括方法的格式构造以及收敛性分析。

第四章数值实验，主要利用Tikhonov正则化方法及全变差正则化方法，正则化参数选取方法为L-曲线法。

【英文摘要】The first-kind Volterra integral equations are a very important kind of integral equa-tions. It has been developed and matured since the twentieth century. Many practicalproblems about physics and mechanics can be solved by changing into the first-kindVolterra integral equations. Whenthe kernel function is continuous or weakly singu-lar, the exact solution is always di?cult to work out. Therefore, the numerical methodsof the first-kind Volterra integral equations play a very important role in mathematics. Byresearching the first-kind Volterra integral equations, there are many wonderful analysis .This article considers the numerical methods of the first-kind Volterra integral equationswhen the data are undisturbed and disturbed.This structure is as follows:In chapterⅠ, we introduce the background , the domestic and foreign researchingsituation and the developping tendency of the first-kind Volterra integral equations. Thisarticle lists some classical methods of solving thefirst-kind Volterra integral equations.In chapterⅡ,we show some preparatory theory of solving the first-kind Volterraintegral equations, ill-posed problems, the regularization methods using in this article:Tikhonov regularization method,total variation regularization method .In chapterⅢ, we research the numerical methods of the first-kind Volterra integralequations when the data is undisturbed, The format structure and convergence analysisare also introduced in this article.In chapterⅣ, we give a numerical experiment based on Tikhonov regularizationmethodand total variation regularization method.The regularization parameter method isthe L-curve method.【关键词】第一类Volterra积分方程不适定问题离散的正则化方法 Tikhonorv正则化方法全变差正则化方法【英文关键词】The first-kind Volterra integral equation Ill-posed problem Discrete regular-ization method Tikhonov regularization method Total variation regularization method【备注】索购全文在线加好友：1.3.9.9.3.8848同时提供论文写作一对一指导和论文发表委托服务【目录】第一类Volterra积分方程数值方法的研究中文摘要3-4Abstract4第1章绪论7-12 1.1 第一类Volterra积分方程的发展历史7-8 1.2 第一类Volterra积分方程的基本理论8-12第2章预备理论12-29 2.1 不适定问题的简介12 2.2 不适定定问题的正则化方法12-29 2.2.1 Titkhonov正则化方法14-16 2.2.2 离散的正则化方法16-19 2.2.3 全变差正则化方法19-26 2.2.4 正则化参数的选取方法26-29第3章数据没有扰动情况下第一类Volterra方程的数值解法29-46 3.1 线性第一类Volterra积分方程的离散正则化方法29-35 3.1.1 方法的格式构造29-30 3.1.2 方法的理论分析30-35 3.2 非线性第一类Volterra积分方程的离散正则化方法35-45 3.2.1 方法的格式构造35-40 3.2.2 方法的理论分析40-45 3.3 本章小结45-46第4章数值实验46-53 4.1 Tikhonov正则化方法的数值实验47 4.2 全变差正则化方法的数值实验47-52 4.3 本章小结52-53结论53-54参考文献54-59致谢59-60攻读学位期间发表的学术论文60。

Science &Technology Vision 科技视界0引言Volterr 型积分不等式在微分方程,积分方程的定性研究中有着非常重要的作用。

2004年Pachpatte [1]研究了一类线性Volterra -Fredholm 型积分不等式:u (t )≤k +h (t )h (α)∫a (t,s )[f (s )u (s )+sh (α)∫c (s,σ)u (σ)dσ]ds+h (β)h (α)∫b (t,s )u(s )ds (1)解的估计。

2008年Ma [2]研究了一类非线性Volterra-Fredholm 型积分不等式:u (t )≤k +α(t )α(t 0)∫σ1(s )[f (s )w (u (s ))+sα(t 0)∫σ2(τ)w (u (τ))dτ]ds+α(T )α(t 0)∫σ1(s )[f (s )w (u (s ))+s α(t 0)∫σ2(τ)w (u (τ))dτ]ds ,t ∈[t 0,T ](2)解的估计。

本文研究了一类幂形式的非线性Volterra-Fredholm 型积分不等式:u (t )≤k +α(t )α(t 0)∫σ1(s )[f (s )u 12(s )+s α(t 0)∫σ2(τ)u 12(τ)dτ]ds+α(T )α(t 0)∫σ1(s )[f (s )u 12(s )+sα(t 0)∫σ2(τ)u 12(τ)dτ]ds ,t ∈[t 0,T ](3)解的估计。

1主要结果本文中R +=[0,+∞],I ∈[t 0,T ];C i (M ,S )为定义在(M ,S )上的i 次连续可微的函数集,其中i =1,2,…;令C 0(M ,S )=C (M ,S )。

定理令u (t ),f (t ),σ1(t ),σ2(t )∈C (I,R +),α∈C 1(I ,I ),α(t )是定义在[t 0,T ]上的连续单调不减函数且α(t )≤t 。