CN107122572A - 一种模拟行波管注波互作用的欧拉方法 - Google Patents
一种模拟行波管注波互作用的欧拉方法 Download PDFInfo
- Publication number
- CN107122572A CN107122572A CN201710433995.4A CN201710433995A CN107122572A CN 107122572 A CN107122572 A CN 107122572A CN 201710433995 A CN201710433995 A CN 201710433995A CN 107122572 A CN107122572 A CN 107122572A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mfrac
- mtr
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000003993 interaction Effects 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000006870 function Effects 0.000 claims abstract description 15
- 230000008569 process Effects 0.000 claims abstract description 15
- 238000002347 injection Methods 0.000 claims abstract description 13
- 239000007924 injection Substances 0.000 claims abstract description 13
- 238000005290 field theory Methods 0.000 claims abstract description 6
- 238000004088 simulation Methods 0.000 claims abstract description 6
- 238000010894 electron beam technology Methods 0.000 claims description 7
- 239000012530 fluid Substances 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 claims description 2
- 238000005859 coupling reaction Methods 0.000 claims description 2
- 238000009826 distribution Methods 0.000 claims description 2
- 239000002245 particle Substances 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 230000009467 reduction Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 6
- 238000003672 processing method Methods 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000010363 phase shift Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 229940028444 muse Drugs 0.000 description 1
- GMVPRGQOIOIIMI-DWKJAMRDSA-N prostaglandin E1 Chemical compound CCCCC[C@H](O)\C=C\[C@H]1[C@H](O)CC(=O)[C@@H]1CCCCCCC(O)=O GMVPRGQOIOIIMI-DWKJAMRDSA-N 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明属于行波管模拟技术,公开了一种模拟行波管注波互作用的欧拉方法。针对现有行波管欧拉非线性理论计算精度低的问题,本发明的方法在基于拉格朗日场论模型的基础上,通过对电子相位采用一种新的处理方式:对电子相位进行傅里叶一阶展开,然后利用贝塞尔母函数关系式,建立了基于欧拉坐标系下的非线性注波互作用理论模型。本发明的欧拉方法具有计算速度快、消耗内存小的特点,能够精确模拟螺旋线行波管中的中度注波互作用过程,为行波管非线性非线性的奠定基础。
Description
技术领域
本发明属于行波管模拟技术领域,具体涉及一种模拟行波管注波互作用的欧拉方法。
背景技术
行波管是使用最广泛的真空电子器件之一,广泛应用于卫星通讯、雷达、电子对抗等领域。随着我国航天事业的飞速发展,对空间行波管放大器的研制和生产提出了非常迫切的需求。然而空间行波管放大器中的各项非线性特性(比如:群时延、相位失真、谐波、三阶互调等)会直接影响卫星和航天器的整体性能。因此非常有必要建立行波管非线性理论来研究这些非线性特性的产生机理和抑制方法。
目前,基于欧拉坐标系下的非线性理论是模拟螺旋线行波管内注波互作用(电子注与高频电磁场的相互作用)过程的重要方法,它具有计算速度快、消耗内存小的特点。然而现有的基于欧拉坐标系下的非线性理论是采用等效线路的方法建立起来的,在模拟行波管内注波互作用的过程中,存在着精度低的问题,因而无法精确模拟螺旋线行波管中的中度注波互作用过程。
发明内容
本发明的目的是为了解决现有欧拉非线性理论计算精度低的问题,以拉格朗日场论理论模型为基础,提出了一种模拟行波管注波互作用的欧拉方法。
本发明的技术方案是:一种模拟行波管注波互作用的欧拉方法,包括以下步骤:
A.在基于拉格朗日坐标系下的场论模型中,将离散的电子注处理为连续分布的流体,对电子相位进行傅里叶一阶展开,然后将展开式代入基于拉格朗日坐标系下的一维运动方程和一维场方程,可以得到基于欧拉坐标系下的运动方程和场方程;
B.利用贝塞尔母函数关系式对步骤A得到的方程组中的积分进行处理;
C.利用步骤B得到的微分方程组(非线性注波互作用理论模型),设置互作用初始条件,求解建立的微分方程组,直到注波互作用结束,即可完成一次注波互作用过程的模拟。
本发明的有益效果:以拉格朗日场论理论模型为基础,本发明的方法对电子相位采用一种新处理方式:对电子相位进行傅里叶一阶展开。将朗格朗日体系的理论转化为欧拉体系理论,构建了一种基于欧拉坐标系下的非线性注波互作用理论模型。本发明的欧拉方法具有计算速度快、消耗内存小的特点,能够精确模拟螺旋线行波管中的中度注波互作用过程,为行波管非线性理论的研究奠定了基础。
附图说明
图1是本发明的模拟行波管注波互作用的欧拉方法流程示意图。
图2是无翼片品型夹持杆高频结构图。
图3是本发明与传统欧拉非线性理论模型以及BWIS-1D代码的功率对比图。
图4是本发明与传统欧拉非线性理论模型以及BWIS-1D代码的增益对比图。
图5是本发明与传统欧拉非线性理论模型以及BWIS-1D代码的相移对比图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。
本发明的行波管注波互作用的模拟方法的流程示意图如图1所示,具体包括如下步骤:
A.基于拉格朗日坐标系的场论模型,将离散的电子注处理为连续分布的流体,对电子相位采用一种新的处理方式:对电子相位进行傅里叶一阶展开。然后将展开式代入基于拉格朗日坐标系下的一维运动方程和一维场方程,可以得到基于欧拉坐标系下的运动方程和场方程;
联立一维拉格朗日坐标下的运动方程和电子相位方程,化简后可以得到:
公式(1.1)为联立后的运动方程。其中,ψn(z)表示第n个粒子的相位,m0和q0分别表示电子的质量和电荷,η(=q0/m0)表示荷质比,ω表示角频率,γ表示相对论因子,υ0表示电子的初速度,F(z)表示归一化场幅值,I表示电子注电流,b表示电子注半径,Kc表示线路的耦合阻抗,z表示轴向位置,Vpc表示归一化相速,上标*表示变量取共轭,IA(=m0c0 3/q0 2),c0表示光速,R'n表示等离子体频率降低因子,kz表示冷腔传播常数,I1(·)表示一阶第一类的变态贝塞尔函数。
一维拉格朗日坐标下的场方程为:
其中:α表示线路衰减,
下面对电子相位采用一种新的处理方式:对电子相位进行傅里叶一阶展开。
其中,初始电子相位分布为ψn|z=0=φ,A0(z)和A1(z)分别表是电子相位的直流分量和一阶分量。
将电子相位的一阶展开式代入运动方程和场方程,化简后可以得到基于欧拉坐标系下运动方程和场方程,分别表示为如下:
运动方程:
其中,分别对K2(z)、和F(z)取共轭便得到和F*(z)。
场方程:
B.利用贝塞尔母函数关系式对方程组中的积分进行处理;
利用贝塞尔母函数关系式处理方程中的积分,对方程进行化简。贝塞尔母函数关系式如下:
其中:Jn(·)表示n阶贝塞尔函数。
利用贝塞尔母函数关系式,方程中的积分可以表示成如下形式:
其中:|·|是绝对值符号,表示取变量的模,ζ1是A1(z)的幅角。
将(1.8)-(1.11)代入方程组(1.4)-(1.6)后并化简,最终构建了一种基于欧拉坐标系下的非线性注波互作用理论模型:
C.利用步骤B得到的微分方程组(非线性注波互作用理论模型),设置互作用初始条件,求解微分方程组,直到注波互作用结束,即可完成一次注波互作用过程的模拟;
这里,可以利用龙格库塔法逐步求解建立的微分方程组。为求解建立的微分方程组,必须设置初始条件。各微分方程对应的初始条件设置如下:
对高频场方程,高频场幅值的初始条件为:
公式(1.11)中,pin和θin分别为输入信号的功率和相位。
对于运动方程,由于在入口处,电子注还未进行速度调至和密度调制,因此相位的零阶、一阶分量及其导数的初始条件均为0:
BWIS-1D代码是由实验验证过的拉格朗日理论模型,因此这里认为BWIS-1D代码是更加准确的模型,因此可以通过对比本发明的以及传统的欧拉非线性模型与BWIS-1D代码的差异来验证本发明的准确性。
本发明采用图2的高频结构,下面将对本发明提出的欧拉非线性理论模型的准确性进行对比验证。图3-图5分别对比了本发明的欧拉非线性理论模型、一维的基于拉格朗日体系的注波互作用仿真软件(BWIS-1D)以及传统的欧拉非线性理论模型(MUSE模型),得到了它们的功率、增益以及相移的对比曲线图。从对比图中可以看出:在中度互作用区(轴向位置0-8cm),本发明的欧拉非线性模型和BWIS-1D代码的功率曲线、增益曲线几乎重合。在中度互作用区(轴向位置0-8cm),本发明的欧拉非线性模型和BWIS-1D代码的相移曲线几乎重合。而在整个互作用区域内,传统的欧拉非线性模型都与前两个模型在增益、相移和功率上存在着很大的差异。可以看出,本发明的欧拉非线性模型不但比传统的欧拉非线性模型更加准确,而且还能精确的分析中度注波互作用过程。
本发明通过对电子相位采用一种新的处理方式:对电子相位进行傅里叶一阶展开,然后利用贝塞尔母函数关系式,建立了基于欧拉坐标系下的非线性注波互作用理论模型。本发明的欧拉方法具有计算速度快、消耗内存小的特点,能够精确模拟螺旋线行波管中的中度注波互作用过程,为行波管非线性非线性的奠定基础。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (4)
1.一种模拟行波管注波互作用的欧拉方法,包括以下步骤:
A.在基于拉格朗日坐标系下的场论模型中,将离散的电子注处理为连续分布的流体,对电子相位进行傅里叶一阶展开,然后将展开式代入基于拉格朗日坐标系下的一维运动方程和一维场方程,可以得到基于欧拉坐标系下的运动方程和场方程;
B.利用贝塞尔母函数关系式对步骤A得到的方程组中的积分进行处理;
C.利用步骤B得到的微分方程组(非线性注波互作用理论模型),设置互作用初始条件,求解建立的微分方程组,直到注波互作用结束,即可完成一次注波互作用过程的模拟。
2.根据权利要求1所述的模拟行波管注波互作用的欧拉方法,其特征在于,步骤A的具体过程如下:
联立一维拉格朗日坐标下的运动方程和电子相位方程,化简后可以得到:
<mrow>
<mfrac>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<msub>
<mi>&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>dz</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&eta;</mi>
<mi>&omega;</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;v</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
<mi>Re</mi>
<mo>&lsqb;</mo>
<mn>2</mn>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>K</mi>
<mn>2</mn>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>-</mo>
<mi>i</mi>
<msub>
<mover>
<mi>Q</mi>
<mo>^</mo>
</mover>
<mi>q</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munder>
<mo>&Sigma;</mo>
<mi>n</mi>
</munder>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.1</mn>
<mo>)</mo>
</mrow>
</mrow>
公式(1.1)为联立后的运动方程。其中,ψn(z)表示第n个粒子的相位,m0和q0分别表示电子的质量和电荷,η(=q0/m0)表示荷质比,ω表示角频率,γ表示相对论因子,υ0表示电子的初速度,F(z)表示归一化场幅值,I表示电子注电流,b表示电子注半径,Kc表示线路的耦合阻抗,z表示轴向位置,Vpc表示归一化相速,上标*表示变量取共轭,θ(=ωz(1/Vpcc0-1/υ0)),IA(=m0c0 3/q0 2),c0表示光速,R'n表示等离子体频率降低因子,kz表示冷腔传播常数,I1(·)表示一阶第一类的变态贝塞尔函数;
一维拉格朗日坐标下的场方程为:
<mrow>
<mo>(</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mi>&alpha;</mi>
<mo>)</mo>
<mi>F</mi>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>IK</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munder>
<mo>&Sigma;</mo>
<mi>n</mi>
</munder>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:α表示线路衰减,
对电子相位进行傅里叶一阶展开:
<mrow>
<msub>
<mi>&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&phi;</mi>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mi>i&phi;</mi>
</msup>
<mo>+</mo>
<msubsup>
<mi>A</mi>
<mn>1</mn>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>i&phi;</mi>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,初始电子相位分布为ψn|z=0=φ,A0(z)和A1(z)分别表是电子相位的直流分量和一阶分量;
将电子相位的一阶展开式代入运动方程和场方程,化简后可以得到基于欧拉坐标系下运动方程和场方程,分别表示为如下:
运动方程:
<mrow>
<mfrac>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>dz</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>&pi;</mi>
</mfrac>
<mfrac>
<mrow>
<mi>&eta;</mi>
<mi>&omega;</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;v</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
<mi>Re</mi>
<mo>&lsqb;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>K</mi>
<mn>2</mn>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mi>d</mi>
<mi>&phi;</mi>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.4</mn>
<mo>)</mo>
</mrow>
</mrow>
1
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>dz</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&eta;</mi>
<mi>&omega;</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&gamma;v</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
<mo>&lsqb;</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>K</mi>
<mn>2</mn>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>i</mi>
<mi>&phi;</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>&phi;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>8</mn>
<msup>
<mi>&pi;</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<msub>
<mover>
<mi>Q</mi>
<mo>^</mo>
</mover>
<mi>q</mi>
</msub>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>i</mi>
<mi>&phi;</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>&phi;</mi>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mi>d</mi>
<mi>&phi;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<msup>
<mi>F</mi>
<mo>*</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>i</mi>
<mi>&phi;</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>&phi;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>8</mn>
<msup>
<mi>&pi;</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<msubsup>
<mover>
<mi>Q</mi>
<mo>^</mo>
</mover>
<mi>q</mi>
<mo>*</mo>
</msubsup>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>i</mi>
<mi>&phi;</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>&phi;</mi>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mi>d</mi>
<mi>&phi;</mi>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,分别对K2(z)、和F(z)取共轭便得到和F*(z)。
场方程:
<mrow>
<mo>(</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mi>&alpha;</mi>
<mo>)</mo>
<mi>F</mi>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<msub>
<mi>IK</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mi>d</mi>
<mi>&phi;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.6</mn>
<mo>)</mo>
</mrow>
</mrow>
3.根据权利要求2所述的模拟行波管注波互作用的欧拉方法,其特征在于,步骤B的具体过程如下:
利用贝塞尔母函数关系式处理方程中的积分,对方程进行化简。贝塞尔母函数关系式如下:
<mrow>
<msup>
<mi>e</mi>
<mrow>
<mi>i</mi>
<mi>x</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
</mrow>
</msup>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mi>&infin;</mi>
</munderover>
<msub>
<mi>J</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>i</mi>
<mi>n</mi>
</msup>
<msup>
<mi>e</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:Jn(·)表示n阶贝塞尔函数。
利用贝塞尔母函数关系式,方程中的积分可以表示成如下形式:
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mi>d</mi>
<mi>&phi;</mi>
<mo>=</mo>
<msub>
<mi>iJ</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mn>2</mn>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>i&zeta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>iA</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mi>d</mi>
<mi>&phi;</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>iJ</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mn>2</mn>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>i&zeta;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>iA</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.9</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>i</mi>
<mi>&phi;</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>&phi;</mi>
<mo>=</mo>
<msub>
<mi>J</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mn>2</mn>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>iA</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.10</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>i&psi;</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>i</mi>
<mi>&phi;</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>&phi;</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>J</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mn>2</mn>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>i&zeta;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>iA</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.11</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:|·|是绝对值符号,ζ1是A1(z)的幅角。
将(1.8)-(1.11)代入方程组(1.4)-(1.6)后并化简,最终构建了一种基于欧拉坐标系下的非线性注波互作用理论模型:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>dz</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&eta;</mi>
<mi>&omega;</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&gamma;v</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>J</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mo>|</mo>
<mn>2</mn>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>iA</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mfrac>
<mrow>
<mi>i</mi>
<mi>&eta;</mi>
<mi>&omega;</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&gamma;v</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
<msub>
<mover>
<mi>Q</mi>
<mo>^</mo>
</mover>
<mi>q</mi>
</msub>
<msub>
<mi>J</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mo>|</mo>
<mn>2</mn>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>J</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mo>|</mo>
<mn>2</mn>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mi>i</mi>
<mi>&phi;</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mi>&eta;</mi>
<mi>&omega;</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&gamma;v</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
<msup>
<mi>F</mi>
<mo>*</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>K</mi>
<mn>2</mn>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>(</mo>
<msub>
<mi>J</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mo>|</mo>
<mn>2</mn>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>i&zeta;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>iA</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mfrac>
<mrow>
<mi>i</mi>
<mi>&eta;</mi>
<mi>&omega;</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&gamma;v</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
<msubsup>
<mover>
<mi>Q</mi>
<mo>^</mo>
</mover>
<mi>q</mi>
<mo>*</mo>
</msubsup>
<msub>
<mi>J</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mo>|</mo>
<mn>2</mn>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>J</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mo>|</mo>
<mn>2</mn>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>i&zeta;</mi>
<mn>1</mn>
</msub>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.12</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>z</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>iIK</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>J</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mn>2</mn>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>i</mi>
<mo>&lsqb;</mo>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>&zeta;</mi>
<mn>1</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
</msup>
<mo>-</mo>
<mi>&alpha;</mi>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.13</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>dz</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>i</mi>
<mi>&eta;</mi>
<mi>&omega;</mi>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&gamma;v</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mfrac>
<mo>|</mo>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<msub>
<mi>J</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mn>2</mn>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mo>&lsqb;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>i&zeta;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mi>i</mi>
<mi>&theta;</mi>
<mo>+</mo>
<msub>
<mi>iA</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>F</mi>
<mo>*</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>i&zeta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mi>i</mi>
<mi>&theta;</mi>
<mo>-</mo>
<msub>
<mi>iA</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.14</mn>
<mo>)</mo>
</mrow>
</mrow>
4.根据权利要求3所述的模拟行波管注波互作用的欧拉方法,其特征在于,步骤C所述的设置互作用初始条件具体如下:
对高频场方程,高频场幅值的初始条件为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>Re</mi>
<mi> </mi>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>=</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<msub>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</msqrt>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>Im</mi>
<mi> </mi>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>=</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<msub>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</msqrt>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.15</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,pin和θin分别为输入信号的功率和相位。
对于运动方程,由于在入口处,电子注还未进行速度调至和密度调制,因此相位的零阶、一阶分量及其导数的初始条件均为0:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>=</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>A</mi>
<mn>0</mn>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>=</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.16</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>=</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>A</mi>
<mn>1</mn>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>=</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1.17</mn>
<mo>)</mo>
</mrow>
</mrow>
3
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710433995.4A CN107122572A (zh) | 2017-06-09 | 2017-06-09 | 一种模拟行波管注波互作用的欧拉方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710433995.4A CN107122572A (zh) | 2017-06-09 | 2017-06-09 | 一种模拟行波管注波互作用的欧拉方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107122572A true CN107122572A (zh) | 2017-09-01 |
Family
ID=59729165
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710433995.4A Pending CN107122572A (zh) | 2017-06-09 | 2017-06-09 | 一种模拟行波管注波互作用的欧拉方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107122572A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108268708A (zh) * | 2017-12-29 | 2018-07-10 | 东南大学 | 获知热状态下真空电子器件参数变化状况的方法及系统 |
CN112417744A (zh) * | 2020-10-27 | 2021-02-26 | 电子科技大学 | 一种基于稳态频域模型的速调管非线性注波互作用模拟方法 |
CN117217133A (zh) * | 2023-11-09 | 2023-12-12 | 北京大学 | 改善激波稳定性的数值模拟方法、装置、设备和存储介质 |
-
2017
- 2017-06-09 CN CN201710433995.4A patent/CN107122572A/zh active Pending
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108268708A (zh) * | 2017-12-29 | 2018-07-10 | 东南大学 | 获知热状态下真空电子器件参数变化状况的方法及系统 |
CN108268708B (zh) * | 2017-12-29 | 2021-04-06 | 东南大学 | 获知热状态下真空电子器件参数变化状况的方法及系统 |
CN112417744A (zh) * | 2020-10-27 | 2021-02-26 | 电子科技大学 | 一种基于稳态频域模型的速调管非线性注波互作用模拟方法 |
CN117217133A (zh) * | 2023-11-09 | 2023-12-12 | 北京大学 | 改善激波稳定性的数值模拟方法、装置、设备和存储介质 |
CN117217133B (zh) * | 2023-11-09 | 2024-02-23 | 北京大学 | 改善激波稳定性的数值模拟方法、装置、设备和存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104850697B (zh) | 基于ansys和adams的大型天线动力学建模方法 | |
CN107122572A (zh) | 一种模拟行波管注波互作用的欧拉方法 | |
CN103199525A (zh) | 一种基于等效节点注入电流变换的配电网潮流计算方法 | |
CN101944145B (zh) | 可去除伪直流模式的微波管高频电路有限元仿真方法 | |
CN103763230A (zh) | 一种改进的自适应盲源分离方法 | |
CN104035750A (zh) | 一种基于fpga的实时模板卷积实现方法 | |
CN106650218A (zh) | 基于ceemd算法和希尔伯特变换的谐波分析方法 | |
Guerin | Farfield radiation of induct-cutoff pressure waves | |
CN102298658B (zh) | 一种行波管注波互作用的模拟方法 | |
CN105760597A (zh) | 基于DG算法的二维色散介质Crank-Nicolson完全匹配层实现算法 | |
CN107944113A (zh) | 一种计算三维高速平动目标电磁散射场的方法 | |
CN106777552A (zh) | 一种能量耦合器的热匹配优化设计方法 | |
CN104933261A (zh) | 一种高效序列拉丁超立方试验设计方法 | |
CN104102471B (zh) | 一种应用fpga定点化技术实现指数cordic算法收敛域扩张的方法 | |
CN106156475A (zh) | 电大尺寸目标的瞬态电磁特性快速提取方法 | |
CN107609233A (zh) | 一种行波管通用注波互作用模型的不连续场匹配方法 | |
CN105303022B (zh) | 快速获取目标电磁散射特性的高斯波束方法 | |
CN107229050A (zh) | 一种基于极坐标格式的雷达成像优化方法 | |
Chen et al. | On the internal resonance modes of time domain surface integral equations for acoustic transmission problems | |
Kim et al. | Two dimensional Riemann problems for the nonlinear wave system: Rarefaction wave interactions | |
CN107391825B (zh) | 一种模拟行波管线性注波互作用的欧拉方法 | |
CN108054070B (zh) | 一种耦合腔行波管群时延的抑制方法 | |
CN109033686B (zh) | 一种行波管返波振荡模拟方法 | |
Xu et al. | Windows simulator: An advanced three-dimensional finite-element S-parameter-simulation tool for microwave tubes | |
CN103763045A (zh) | 一种适用于三维水声信道仿真的高斯束方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170901 |
|
RJ01 | Rejection of invention patent application after publication |