毕业设计(论文)--异形脊波导传输特性分析.doc

题 目 异形脊波导传输特性分析
学生姓名 学号
所在学院 物理与电信工程学院
专业班级 电子1204班
指导教师
完成地点 陕西理工学院

2016年 6月 5日

陕西理工学院毕业设计

异形脊波导传输特性分析

作者：
（陕西理工学院 ,物理与电信工程学院, 电子信息工程系2012级4班，陕西 汉中 723000）
指导老师：

[摘要]在现代微波工程中，为了满足微波传输系统性能的某些需求，需要不断探索和研究具有特殊截面形状的各类新型波导，多脊波导的理论研究对解决实际问题有很好的参考价值。由于其能实现微波小型化及具有良好的宽带特性等优点，具有很好的工程应用价值。有限差分法突破了解析法只能分析结构简单的波导器件，并且应用其处理脊形波导结构具有不需要复杂的物理模型等许多优点。本文从泰勒级数出发推导出用有限差分法计算脊波导的截止频率，以单脊波导为例，编程计算及验证了单脊波导内TM波的截止频率，并分析计算了新型脊波导-梯形脊波导的TM波的截止频率，所得结果为多脊波导的设计与制作提供理论参考。
[关键词]有限差分；单脊波导；梯形脊波导；截止频率；TM波

The analysis of transmission characteristics about the shaped trapezium ridge waveguide

Author： LuYang
（Grade 11, Class 1, Major electronics and information engineering, School of Physics and Electronic Information Engineering, Shaanxi University of Technology, Hanzhong 723000, Shaanxi）
Tutor： Chunjiang Shuai

Abstract:In modern microwave engineering, in order to meet some of the requirements of microwave transmission system performance, we need to continue to explore and research new types of waveguide having a special cross-sectional shape, more ridge waveguide theory to solve practical problems for a good reference value. Due to its size and to achieve microwave has good broadband characteristics, etc., and has good engineering application value. Finite difference method to break the analytical method can only analyze a simple structure waveguide devices, and applications which deal with ridge waveguide structure has many advantages not require complex physical models. From the Taylor series starting deduced ridge waveguide by finite difference method cutoff frequency to a single ridge waveguide, for example, calculation and verification program within a single ridge waveguide cutoff frequency field intensity distribution and TM waves, and analysis intensity distribution shaped trapezium ridge waveguide TM wave and the cut-off frequency, the results provide a theoretical reference for the design and production of multi-ridge waveguide.
Key words: Finite-Difference;single ridge waveguide; trapezium ridge waveguide;Cutoff frequency;TM wave

目录
1 引言 1
1.1 脊波导的发展 1
1.2 课题意义 1
1.3 数值法简介 2
2 有限差分基本原理 3
2.1 基本差分公式 3
2.2差分方程的求解过程 4
2.3 波导中的电磁场方程 6
2.4 亥姆霍兹方程的差分表达式 8
3 FORTRAN语言 12
3.1 简介 12
3.2 基本程序结构 12
4 脊波导TM波传输特性分析 13
4.1 单脊波导传输特性分析 13
4.2 梯形脊波导传输特性分析 16
5 总结 23
致谢 24
参考文献 25
附录 A 外文文献 26
附录 B 中文翻译 30
附录 C 程序设计 34

1 引言
电磁波在有界空间的传播，就被称为导波系统中的电磁波。所谓导波系统是指引导电磁波沿一定方向传播的装置，它不仅是能量和信息的载体和传递工具，而且是构成各种高频、微波元件和仪器的基础。脊波导是其中的一种常见的导波结构，它是矩形波导的一种变形，又称凸缘波导[1]。
近年来，微波因其本身的特点，在信息获取、传输中所占地位愈显重要，在这样的大环境下，各种高质量的微波器件（包括用脊波导制作的各种器件）会有很大的实际需求，在雷达及卫星通讯设备中，最常用的是矩形波导和圆波导，在现代微波工程中矩形脊波导具有较低的截止频率,较宽的工作带宽,低特性阻抗等优点,使得矩形脊波导在微波和毫米波器件中被广泛应用[2]，如：宽带脊波导滤波器、宽带定向耦合器、双工器、变频器、移相器、混频器、低噪声放大器、脊波导缝隙天线阵等。
尽管脊波导的应用越来越广泛,然而有关它的理论分析计算的资料却较少，因此如何计算脊形波导中传输的电磁波的截止频率是一个值得研究讨论的问题。经过对论文的研究发现，对非规则脊波导研究的文献很少，在本文中将采用有限差分法这一经典的数值方法结合计算软件Visual fortran来分析梯形脊波导的传输特性。
1.1 脊波导的发展
早在1947年，Cohn Seymour. B在论文“Properties of ridged waveguide” 中对脊波导特性进行了初步研究，利用横向谐振技术给出了求解矩形脊波导的截止波长和特性阻抗的方程和特性图，从此脊波导便一直受到广泛的关注。
1955年，Hopfer 在论文“The Design of Ridged waveguides”中给出了一些不同外观比例的脊波导各个参数的计算方法，还给出了一些设计曲线。到了60年代，Pyle用准静态法得到了任意外形比例的脊波导的特征值，其结果比Hopfer的更精确。以上文章对脊波导的研究属于开创性的，但都只能解决脊波导的TEn0模的特征值问题，都对其它高次模的解无能为力，特别是当脊很薄和脊沟很窄时，其精度受很大限制。70年代初Konishi等把用脊波导制作的12GHz低噪声变频器等器件成功地用于卫星通信，并对矩形脊波导的不连续性作了分析计算，使得矩形脊波导在微波和毫米波器件中被广泛应用。80年代，Monsour用模式分析技术得到了双脊波导的TE模的近似闭式表达式，其中用到了横向谐振的概念。90年代，J.Helszajn用有限元法分析了脊波导的特性阻抗，给出了精确的曲线图，并提出了梯形脊波导，并分析了其阻抗特性。2000年Javis给出了任意高度的双脊波导设计曲线，他们所用的方法都属于与上述文章类似的等效电路法[3]。
国内在脊波导的研究领域起步较晚，但发展迅速。1997年黄彩华在文献“矩形变形脊波导主模截止波长和特性计算”中，用等效电路法得到了矩形变形脊波导的主模和特性阻抗。1999年李永明等在论文“有限元外推法在波导本征值问题中的应用”中用有限元法对矩形双脊波导进行了研究，得到了主模截止波长相关数据[4]。2004年任列辉等在“用电磁场算子理论分析脊波导的传输特性”中，用电磁场算子理论求解脊波导的本征值，在此基础上讨论了脊波导的的传输特性。王萍在“脊波导各种参数的计算”一文中，用等效电路法对脊波导进行了分析，得到了截止波长，特性阻抗，功率容量等数据。2006年薛德强在论文“下限单脊矩形波导传输特性研究”中，用有限差分法对下限单脊矩形波导进行了研究。2007年陈小强教授等在论文“两种新型双脊波导传输特性的研究”中提出了三角形和倒梯形两种新型对称双脊波导，并计算了TE模式的传输特性[4]。
尽管波导历史悠久,但时至今日仍有各种新的导波结构不断涌现。如今种类繁多的波导广泛应用于现代通讯,雷达及微波,毫米波,光电系统中。其种类包括了传输线类波导;金属类空心波导(如圆波导,矩形波导,椭圆波导,三角形波导及近年来新兴的五边形波导等);微带线类波导(如屏蔽微带线,带状线等);介质类波导(如阶跃型光纤,渐变折射率型光纤,椭圆光纤,蝴蝶结光纤,光纤光栅等) [5]。各种新材料波导更是层出不穷,等离子体金属波导,等离子体介质波导,非线性波导,光子晶体光纤等成为波导研究的重要领域。纵观今日的各类波导,工作频率已从射频贯穿到光频,波导横截面尺寸从米量级直到亚微米级。

1.2 课题意义
在现代微波工程中脊波导具有较低的截止频率，较宽的工作带宽，低特性阻抗等优点，使得脊波导被广泛应用。它常用作宽频带传输系统及宽频带测试系统，比如宽带脊波导滤波器、宽带定向耦合器、双工器等等；也常作为微波管的输出波导，比如：变频器、移相器、脊波导缝隙天线阵等等。为了满足微波传输系统性能的某些需求，我们还需要不断探索和研究具有特殊截面形状的各种新型波导。最近几十年来各种结构形状的脊波导应运而生，它们可以看作是矩形波导的变形，本文所研究的梯形脊波导就是其中一种情形。
长期以来分析波导问题的解析方法和数值方法有很多，有限差分法和有限元法对复杂问题具有较强的适应性，近年来,为了解决大规模电磁场数值计算问题，基于有限差分和有限元的区域分解法已受到人们的广泛重视。
本文就是采用有限差分法对梯形脊波导的传输特性进行了分析，得到了它们的截止频率的一些变化规律，给脊波导在今后工程应用中提供了一定的数据参考价值。
1.3 数值法的简介
所使用的解析法和数值法都是基于实际问题的解析模型进行求解，其过程是从基本的麦克斯韦方程组或波动方程出发，加上特定的边界条件和本构参数，推导出微分方程或积分方程的体系，然后尽可能地采用解析手段具体问题具体分析，最后则可以编制出一个紧凑而且高效的计算程序[6]。
数值法是指直接将待求解的数学方程进行离散化处理，将无限维的连续性问题化为有限维的离散问题，将解析方程的求解问题化为代数方程的计算问题的一类方法。用高性能的计算机，可以直接以数值的、程序的形式代替解析形式来描述电磁场问题。在纯数值法中，通常以差分代替微分、用有限求和代替积分，将问题化为求解差分方程或代数方程问题。其优点是原则上它可适用于任何复杂的边界问题的求解，既可用于求解第一、二类可解的问题，又可用于求解因边界复杂由它们无法解决的问题，原则上可以求得所需要的任意精度，仅受限于计算机容量、速度和舍入误差、但随着大容量、高速度巨型电子计算机的发展和数值计算方法本身的发展，可以预料，许多过去只能靠实验来分析和设计的电磁场边值问题都可以靠数值法来实现，因而数值法有着广阔的应用前景[7]。
数值法，包括有限差分法、矩量法、边界元法及有限元法等。矩量法及边界元法的特点是计算精度高，但对于复杂边界形状不易导出特定的积分，不能模拟任意的媒质分布参数问题。有限元法适应性强，可以用于各种复杂的媒质分布等问题。有限差分法由于其直接从麦克斯韦方程组出发，不需要任何导出方程，这就避免了使用更多的数学工具，使其成为所有电磁场计算方法中最简单的一种。还有一种数值方法就是时域有限差分法,它是在有限差分法的基础上提出的，也是当今重要的电磁场数值计算方法之一，但由于其通常对计算机的内存配置和CPU的运算能力要求较高，其应用受到一定限制。
2 有限差分基本原理
2.1 基本差分公式
所我们要求的电位函数，它在区域内满足下面的拉普拉斯方程

(2-1)
在边界上，它服从以下条件：
(2-2)
式中为边界点的函数。这类问题一般称为第一类边值问题或称狄里赫利问题。
为了用差分方法求解电位分布，先在平面分别作两族平行于轴和轴的直线,线间的距离为，于是各直线的和坐标分别为：

式中为正整数，取值1、2、……。这样区域就被许多边长为的正方形所覆盖，在图2-1中示出了这种情况。各正方形的顶点被称为网格的节点，从图可以看到，各节点所处位置有所不同。一些节点（例如节点）恰落在边界上，我们把它叫做边界节点。有些节点到边界的距离不足（例如节点）,这些节点叫做不规则节点。但是大部分节点到边界的距离大于，例如图上的点，它们属于规则节点。差分法就是求这些离散节点处的近似值。以后将用符号来表示节点，而将节点处的值用代表[8]。
图2-1 将场域划分成网格

i=1,2,3…
j=1,2,3…
a
b
0
D
（x0，y0）

0
1
2
3
4
h1
h2
h3
h4
图2-2 不等距网格

当用差商来替代偏导数，可以从(2-1)式导出它的差分表达式。在直角坐标系如果采用图2-2所示的不等距网格时，求出的差分表达式如下
(2-3)
对于常用的等距网格情况，。根据式（2-3）就得到
(2-4)
在轴对称情况，拉普拉斯方程就成为下面的形式
(2-5)
它的差分表达式对不等矩网格为
图2-3 轴对称情况的不等距网格
h1
h3
h2
r0
轴线
r
Z
１
２
３
４

(2-6)
式中各量见图2-3。
在等距网格情况，。如果令轴线处，而点落行并且，则。根据（2-6）式就得出
(2-7)
假使点落在的轴上，可以求出差分格式为
(2-8)

2.2差分方程的求解过程
0
100
100
100
0
0
0
0
0
0
0
0
下面用一个计算金属槽内电位分布的例子来说明差分方程的求解过程。设金属槽上盖板电位等于100，侧壁与底壁的电位为0。将场域用正方形网格划分，要求的是场域内部九个节点的电位值。对网格节点1，根据式（2-4）的等距差分公式可以得到

也可以扩展成以下的形式

图2-4 金属槽结构
对其它标号为2~9的各个节点也可写出类似的扩展方程，将它们组合起来，就得到下面的依矩阵形式表示的线性代数方程组

图1－7 差分法计算框图
这是一个包含九个变量的方程组，可采用线性代数中塞德尔迭代法求解。但是采用这种方法，必须先给各个变量赋以初值，这样就可从第一至第九个方程依次求出的改善值。不断重复种计算，直至相邻两次算出的各个值收敛到我们需要的精度为止。
在实际的计算过程中，当然用不着先列出上面所述的方程组，它仅在网格节点很少时候才有可能。注意到前面列出的各个差分公式中，网格中某节点的电位只和它四邻节点电位有关这一特点，例如图1-6中节点5的电位仅和它四周相邻点2、4、6、8的电位有关。就可以从左至右，从上至下一行一行地扫描逐一地求出各网格节点的电位[9]。
为判别计算的收敛程度，一般采用以下两种方法：
（一） 对所有网格点，找出相继的两次迭代的最大相对误差
(2-9)
式中表示第次迭代。
（二） 找出相继两次迭代形成的各网格节点的平均相对误差
(2-10)
对同一问题，由第二种方法求出的收敛误差比第一种方法小1~2个数量级，所以如果用第一种方法是收敛精度误差控制在，则用第二种方法时就应控制在[6]。
根据以上所说，可以写出用差分方法计算静电电位问题的框图如图2-5。

启动
给定边界位置信息

给定边界电位和给其它节点赋初值

图2－5 差分法计算框图

是
打印计算结果

计算收敛到指定精度乎？

全区扫描结束乎？

k＝k +1

迭代次数k＝0

是
否
否
是
是
否
否
j＝j +1

i＝i +1

用不规则点差分公式计算

用规则点的五点差分公式计算

第j行扫描结束乎？

j＝1

i＝1

规则点乎？

2.3 波导中的电磁场方程
当波在任意截面的波导中传播时，可以将其中传播的电磁波分成TE波和TM波[10]。现在将分别把各个场分量关系总括如下：
对TM波，如令，则其它各个场分量为：

， (2-11)
， (2-12)
以上诸式中的横向电磁场分量可以归并成以下简洁的形式：
(2-13)
(2-14)
对TE波，可令，则其它各个场分量为

， (2-15)
， (2-16)
各个横向分量也可合并如下：
(2-17)
(2-18)
在上面各式中，为相位常数，它等于
(2-19)
被称为截止波数，他由下式决定
(2-20)
式中，——波导的截止频率；
——光速。
2.4 亥姆霍兹方程的差分表达式
无论是TE波还是TM波，上面定义的函数均满足亥姆霍兹方程
(2-21)
但是对于不同波型有不同的边界条件。对TM波，；对TE 波,，代表波导截面的边界线。
用求拉普拉斯方程的差分方程一样的方法，不难求出在不等距网格情况下，亥姆霍兹方程的差分表达式为
(2-22)
式中的、、、定义见图2-6。
对等距网格情况，，就得到
(2-23)
当用式（2-21）或（2-22）求解亥姆霍兹方程时，和解拉普拉斯方程有所不同。前者需要给定值，但这个值在计算前是不知道的。为消除这个困难，可采用下面的措施，就是在计算开始前，除了给内部网格节点赋初值外，也给一个初值，有了这个初值，就可进行值迭代计算而将各节点值改善到一个较好的精度，如果松弛系数选得适当，只要几次迭代就行了。然后用各节点的这些值用下面讲的方法求出改进的值。再利用这个值去求，这种过程循环重复，直到和各节点值均收敛到应用的精度为止[11]。
为了从各网格节点值求，可采用以下方法。用乘式（2-20）后并将整个等式对波导横截面进行面积分，得到：

或者写成

式中为对横向坐标的拉普拉斯算子。由上式解出
(2-24)
利用这个公式求，直观上也可看出比直接用式（2-21）来求更合理，因为它意味着在整个场域内的平均，这就可能使计算形成的误差部分地被抵消。
（2-23）式可用下面的差分近似式来表示
(2-25)
上式中的代表对所有网格节点求和，代表第（i，j）节点所占有的面积[15]。对于区域内的规则点，例如图2-7中的g点，。但对于面积上的W点，；

1
2
3
4

0

图2-6 不等距网格

图2-7 求边界节点的

对于边界上的角点z，，这些关系由图上的阴影区容易看出。在TM波情况，由于边界上，所以沿边界的乘积，这是在式（2-24）中所有其他内部点均等于而可约去，式（2-24）就简化成
(2-26)
在以后对波导计算时，编制的计算机程序所采用的方法同于计算电位的超松弛法，但稍作改变如下：
将式（2-22）改写成下面的形式
(2-27)
式中叫做余量，表示0点的第次迭代求出的值，各点1，2，…等位置见图2-2。再来求0点的新值，它应使得式（2-26）的余量为零，即
(2-28)
为了加速收敛，乃令第次迭代求出的值为，它等于
(2-29)
式中为松弛系数。用式（2-26）减去式（2-25）可得

或者
(2-30)
将式（2-27）代入式（2-28）就有
(2-31)
计算过程是在给定了后，先由（2-29）式计算各节点余量，再根据式（2-30）求出改善后的，用这种方法算出场区各点的值。重复这样的计算，经对指定的迭代次数后，再由式（2-24）或（2-25）求出较精确的值。利用这个新的值，重复上述过程直至值收敛为止[16]。
在求出后，就可以算出截至频率。因为

由上面的公式可解出
(2-32)
根据以上所述，整个用差分法求解波导问题的计算框图示如图2-8。

给边界点赋值和给内部节点赋初值

输入的初始估计值

进行指定次数的迭代以改善场量u的值

计算新的值

是
值是否收敛

给边界点和内部点赋值

否

总迭代次数大于指定次数乎？

否
给边界点赋值

是

图2-8求解波导问题的计算框图

3 FORTRAN语言
FORTRAN语言，是第一个面向过程的高级程序设计语言，主要用于科学计算，也可用于数据处理和仿真。FORTRAN是英文FORmula TRANslator的缩写，原意是公式翻译。FORTRAN语言可使程序员用一种非常接近于常用数学表达式和英语自然语言的方式编制计算机程序。自1956年开始使用以来，一直在国际上广泛流行，是使用最广泛的程序设计语言之一。
3.1 简介
FORTRAN语言问世以来，根据需要几经发展，先后推出了不同的版本，其中最流行的是1958年出现的FORTRANⅡ和1962年出现的FORTRANⅣ。1966年美国标准化协会(ANSI)公布了两个美国标准文本：
·标准FORTRAN（X3.9-1966），大致相当于FORTRANⅣ。
·标准基本FORTRAN（X3.10-1966），大致相当FORTRANⅡ。
1972年国际标准化组织(ISO)接受了美国标准，在稍加修改后公布了ISO FORTRAN标准，即《程序设计语言FORTRAN ISO 1539-1972》，它分为三级，即：
· 完全的(一级)FORTRAN，相当于FORTRANⅣ。
· 中间的(二级)FORTRAN，介于FORTRANⅡ和FORTRANⅣ之间。
· 基本的(三级)FORTRAN，相当于FORTRANⅡ。
FORTRANⅣ(即FORTRAN 66)流行了十几年，几乎统治了所有的数值计算领域，许多应用程序和程序库都是用FORTRANⅣ语言编写的。美国标准化协会(ANSN)在1976年对ANSI FORTRAN(X3.9-1966)进行了修订，预定在1977年通过，为了区别于FORTRAN 66，新标准定名为FORTRAN 77。实际上到1978年4月才由ANSI正式公布作为新的美国国家标准。即FORTRAN(X3.9-1978)。1980年，FORTRAN 77被接受为国际标准，即《程序设计语言FORTRAN ISO 1539-1980》，该标准分为全集和子集。
中国制订的FORTRAN标准，基本上采用了国际标准，于1983年5月公布执行，标准号为GB3057-82。
FORTRAN 77标准完成后，新版本的修订工作也在同一时间开始进行。这个版本进行了15年，最后在1992年正式由国际标准组织ISO公布，它就是FORTRAN 90。FORTRAN 90对以往的FORTRAN语言标准作了大量的改动，使之成为一种功能强大、具有现代语言特征的计算机语言。其主要特色是加入了面向对象的概念及工具、提供了指针、加强了数组的功能、改良了旧式FORTRAN语法中的编写“版面”格式。
FORTRAN 95标准在1997年同样由ISO公布，它可以视为是FORTRAN 90的修正版，主要加强了FORTRAN在并行运算方面的支持。同时一些公司纷纷推出Visual Fortran，这为工程技术界进行科学计算和编写面向对象的工程实用软件的用户提供了极大的方便。熟悉VB或VC的读者可以很容易地掌握Visual Fortran的使用，进一步开发出自己专业领域的Windows下的界面友好的工程应用软件[13]。
3.2 基本程序结构
FORTRAN程序由一个或几个相对独立的程序段组成，其中必须有一个主程序段(PROGRAM),还可以有（也可以没有）子程序段（subroutine）、函数段(FUNCTION)或数据段(BLOCK DATA)。程序段由语句组成。语句分为可执行语句和非执行语句两类[14]。可执行语句是在程序执行时能导致系统硬件作出一个实际操作的语句，包括赋值语句、控制语句、输入输出语句等。非执行语句是在程序编译时为编译程序提供有关信息的语句，包括说明语句、格式语句、数值语句、参数语句、函数定义语句和程序段语句等。语句是由常数、变量、运算符和专用定义符等按事先规定的格式书写。语句标号则提供了引用单个语句的方法。

4 脊波导TM波传输特性分析
4.1 对称单脊波导传输特性分析
对于脊形波导中传播的TM波，；在边界应满足的条件。此外，就基波而言，它的分量称于截面中心线，所以只要研究截面的一半区域就可以了，此时中心处的边界条件为，为垂直中心线的外法线方向[15]。由于沿金属边界，故而计算时可以应用式（2-25）。
根据图2-10框图编制的程序列在下面(见附录C)，其中ALFA变量代表松弛系数，SQKHOLD代表上一计算得到的值，SQKHNEW代表新算出的值，为计算内部节点时常用的一个函数子程序。[16]如图4-1为对称单脊波导的立体模型，图4-2所示为对称单脊波导的横截面图。

图4-1 对称单脊波导立体模型 图4-2 对称单脊波导截面图

b
c
d
e
HH
a
对于对称单脊波导，由于其对称性，我们仅仅需要研究其一半电位分布来求取截止频率。[17，18]待计算的对称单脊波导的一半截面示如图4-3，采用正方形网格对波导进行划分。HH为波导侧壁高度。

图4-3 对称单脊波导的网格划分图

4.2对称单脊波导的TM波的截止频率分析

由于工程上的所需，改变脊的高度与长度，波导的TM波的截止频率会有不同的变化[19]。如图4-3所示，取=5.0mm，=2.5mm,=2.5mm,=2.5mm,=1.25mm,=1.25mm。划分网络=10；=5；=5-10；=5；=5-10。
计算所得截止频率数与表4-1理论近似，计算结果正确，程序无误。
表4-1计算值与参考值比较f（GHz）

本文
3.2226
参考值
3.2223
(1) 取=5不变时， 改变/的比值所得到的TM波的截止频率如表4-2所示，据此绘制出如图4-4所示的对称单脊波导TM波截止频率随脊的高度比变化曲线。
表4-2对称单脊波导的TM波随高度比变化的截止频率（GHz）
b/HH
0.1
0.2
0.3
0.4
0.5
0.6
0.7
f
20.184
22.389
25.307
28.990
32.226
33.487
33.806

图4-4 对称单脊波导TM波截止频率随脊的高度比变化曲线图
由图4-5可看出，随脊高比（）的增大，TM波的截止频率呈上升趋势，且比例为0.5时开始趋势较为平缓。
(2) 当=5不变时，改变/的比值得到的TM波的截止频率如表4-3所示，据此绘制出如图4-5所示的对称单脊波导TM波截止频率随脊的宽度比变化曲线。
表4-3对称单脊波导的TM波随宽度比变化的截止频率（GHz）
d/c
0.1
0.2
0.3
0.4
0.5
0.6
0.7
f
24.430
26.331
28.847
32.226
35.946
37.275
37.460

图4-5对称单脊波导TM波截止频率随脊的宽度比变化曲线
由图4-7可以看出，随脊的宽度比（）的增大，TM波的截止频率呈上升趋势，比列为0.5时趋势比较平缓。

4.2.1梯形脊波导传输特性分析

图4-6梯形脊波导的实际模型 图4-7 梯形脊波导的剖面图

4.2.2梯形脊波导的TM波的截止频率分析
由于工程上的所需，改变脊的高度比，波导的TM波的截止频率会有不同的变化。如图4-8所示=10.0mm，=5.0mm,=5.0mm,=28.0mm,=5.0mm,=8.0mm，θ=45。。划分网络=20；=10；=10-20；=10；=8。
保持d不变，改变b/hh,得到的TM波的截止频率如表4-4所示，据此绘制出如图4-8所示的梯形脊波导TM波截止频率随脊的角度变化曲线。
表4-4对称梯形脊波导的TM波随高度比变化的截止频率（GHZ）

d/c
0.1
0.2
0.3
0.4
0.5
0.6
0.7
f
5.94
3.37
1.98
0.947
0.808
0.741
0.695

图4-8梯形脊波导TM波截止频率随脊的高度比变化曲线
由图4-8可以看出，随脊高比（）的增大，TM波的截止频率呈下降趋势，且比例为0.5时开始趋势较为平缓。
保持b不变，改变d/c得到的TM波的截止频率如表4-5所示，据此绘制出如图4-9所示的梯形脊波导TM波截止频率随脊的角度变化曲线。

b/HH
0.1
0.2
0.3
0.4
0.5
0.6
0.7
f
18.56
15.24
12.93
9.47
8.36
8.29
8.24
表4-5对称梯形脊波导的TM波随高度比变化的截止频率（GHz）

图4-9梯形脊波导TM波截止频率随脊的宽度比变化曲线
由图4-9可以看出，随脊高比（）的增大，TM波的截止频率呈下降趋势，且比例为0.5时开始趋势较为平缓。

5 总结
本文主要运用有限差分法（FDM）分析了脊波导的TM波的传输特性，用FORTRAN程序语言编制系列程序，对对称单脊与梯形脊波导中电磁场进行模拟，从而得到在这些波导中传输的电磁波的截止特性。主要做了以下工作：
（1）对有限差分法的基本概念及原理进行了综合性论述，通过求解波导截止频率给出其相关公式推导；
（2）用泰勒级数法推出直角坐标系下拉普拉斯方程的差分格式，并从直角坐标系中TE波电磁场纵横关系出发，以对称单脊波导为例，采用FDM算法，用FORTRAN程序语言编制系列程序分析计算其截止频率，并对梯形脊波导进行模拟，从而得到在梯形波导中传输的电磁波的截止特性。
（3）通过对各种脊形波导的结构变化分析，绘制出其TM截止频率随高度比与长度比变化曲线，
得到其传输特性。

致谢

本论文的写作从开题、需求调研、搜集资料、分析设计到最后成文，历时数月。其间，得到了老师和同学及朋友的各种帮助。在此，我表达衷心地感谢。感谢导师帅春江一直以来的精心指导,从立题到完成学位论文,每一步都凝聚着他的心血，老师以其严谨求实的治学态度、高度的敬业精神、孜孜以求的工作作风和大胆创新的进取精神对我产生了重要影响。他渊博的知识、开阔的视野和敏锐的思维给了我深深的启迪，将我带入一个崭新的领域,导师对学术的严谨态度和执著追求的精神,将永远激励我努力奋斗。导师宽以待人和豁达的生活作风,将使我终生受益。深深地感谢帅老师为我创造的一切学习的条件。
感谢物电学院领导的支持和鼓励,感谢电子信息工程实验室提供的条件,感谢电子信息工程系的同学的关心和支持。