基于GMRES（m）法的双连通区域数值保角变换的计算法

吕毅斌，赖富明，王樱子，武德安

(1.昆明理工大学理学院，云南昆明650500）

(2.昆明理工大学计算中心，云南昆明650500）

(3.电子科技大学数学科学学院，四川成都611731）

基于GMRES（m）法的双连通区域数值保角变换的计算法

吕毅斌1，赖富明1，王樱子2，武德安3

(1.昆明理工大学理学院，云南昆明650500）

(2.昆明理工大学计算中心，云南昆明650500）

(3.电子科技大学数学科学学院，四川成都611731）

本文研究了基于模拟电荷法的双连通区域的数值保角变换问题.利用限制Krylov子空间最大维数的算法-GMRES（m）算法，求解基于模拟电荷法的双连通区域数值保角变换中的约束方程，获得了模拟电荷和变换半径，构造了近似保角变换函数.数值实验表明了本文算法的有效性.

模拟电荷法；双连通区域；Krylov子空间；GMRES（m）法

1　引言

保角变换是复变函数的一个基本问题，广泛应用于物理学与工学［1-3］.保角变换的主要求解方法：解析法和数值计算法.解析法指出了变换函数的存在，但是只能给出一些特殊区域的变换函数表达式.基于实际问题的复杂性，必须通过数值计算方法求变换函数.目前的保角变换分为两类：一类是单连通区域的保角变换，另一类是多连通区域的保角变换.

1969年，德国的Steinbigler提出了在计算空气中回转对称电极周围电场时可以用设置在电极内部的若干个虚设的模拟电荷来计算电极表面上电荷分布的电场的计算方法，形成了模拟电荷法（Charge Simulation Method）的基本思想［4］.自20世纪80年代起，天野要等对模拟电荷法、数值保角变换、数值保角变换的模拟电荷法以及数值保角变换中模拟电荷点的计算方法做了许多研究工作，提出了基于模拟电荷的保角变换计算法（天野法）［5-8］.在复杂的数值保角变换的计算中，相比传统的数值计算法，基于模拟电荷法的数值保角变换计算法具有计算精度高、误差评价容易、计算时间短等优势.

广义极小残量法（Generalized Minimum Residual Methods）［9-11］是在Krylov子空间进行迭代的计算方法，但是GMRES法的储存量和正交化工作量会随着迭代次数的增加而增加，因此，一般需要利用限制Krylov子空间最大维数的算法-GMRES（m）法求解计算问题.本文利用GMRES（m）求解出了双连通区域模拟电荷法中的约束方程，得到了模拟电荷和变换半径，从而得到近似保角变换函数，最后利用数值实验验证了算法的有效性和结果的正确性.

2　基于模拟电荷法的数值保角变换计算法

本节主要阐述基于模拟电荷法的双连通区域保角变换的数值计算法［12-13］.对于z平面上的一个有限双连通区域D，它由两条封闭的Jordan曲线C1和C2所围成.通过保角变换把它变成了w平面的一个圆环µ＜|ω|＜1，C1和C2分别是外部边界和内部边界.

图1：基于模拟电荷法的双连通区域数值保角变换

在不失一般性的情况下，假定f（0）=0，f（z）满足正规化条件f（∞）=∞，f＇（∞）＞0时，可以表示成

g（z）是Dirichlet型势场问题

的解，h（z）是g（z）的共轭调和函数［14］.由模拟电荷法可知g（z）可以用C1和C2所围成的区域外部配置N个电荷点作为极的对数势场的一次结合

高度近似.因此g（z）在区域D内的共轭调和函数h（z）可以用

近似.另外，µ由M近似.ζi（i=1，2，···，N）为电荷点，分布在给定区域的外部，即N/2个点分布在边界C1的外部，另外N/2个点分布在边界C2的内部.根据（2.2）式，双连通区域的数值保角变换是单连通区域的内部保角变换与外部保角变换的组合，因此边界条件可以看成内部保角变换边界条件与外部保角变换边界条件的组合.未知电荷Qi可以用边界上选择的约束点zj在满足Dirichlet边界条件时进行求解，即

又由于g（∞）=0，h（∞）=0，根据（2.3）和（2.4）式可推导出

上述（2.5）式中zj（j=1，2，···，N/2）取在边界C1上.类似的，（2.6）式中zj（j=N/2+ 1，N/2+2，···，N）取在边界C2上.

联立（2.5）-（2.7）式，化为N+1维约束方程

其中aij=log|zj-ζi|.最后，利用zi，ζi，Qi，M计算双连通保角变换.

3　基于GMRES（m）法的模拟电荷的求解

将约束方程（2.8）式写成标准线性方程组Ax=b的形式，其中

约束方程的系数矩阵A非对称的，广义极小残量法（GMRES）是当前求解大型非对称线性方程非常有效的算法之一.它是在Krylov子空间

的基础上发展起来的算法，通过寻找近似解xm∈x0+Km（A，r0），来逼近约束方程（2.8）的准确解.但是GMRES法在求解约束方程时，随着电荷点数的增加，计算量和储存量都会随着迭代次数的增加而大大增加.为了避免这种情况，每m步必须重新启动，这就是GMRES（m）法［15］.GMRES（m）算法是GMRES算法的进一步改进，根据文献［16］，可以得到求解模拟电荷的GMRES（m）法，其算法步骤如下：

Algorithm 1 GMRES（m）Algorithm

上述算法中，x0是给定的初始近似解，取为零向量，r0=b-Ax0为初始向量，ε是给定精度.

4　数值实验

在MATLAB 7.0环境下，以椭圆为边界的双连通区域C1：x2/a21+y2/b21=1，C2：x2/a22+y2/b22=1为例，C1和C2为边界.利用模拟电荷法对双连通区域的保角变换进行数值实验，保角变换误差由

确定，检验GMRES（m）法的有效性，具体步骤如下：

步骤1设定电荷点ζi（i=1，2，···，N），约束点zi（i=1，2，···，N）以及各个参数；

步骤2根据边界条件确定关于模拟电荷Qi（i=1，2，···，N）和变换半径M的约束方程；

步骤3通过GMRES（m）法求解约束方程，得到模拟电荷Qi（i=1，2，···，N）和变换半径M，从而构造近似保角变换函数f（z）.

例1 a1=7，b1=5，a2=5，b2=1时.分别用Method 1和Method 2表示天野法和基于GMRES（m）法的双连通区域保角变换的计算法.图2给出的是这两种方法数值保角变换误差.由图2可知电荷点数增加时，误差会随着电荷点数的增加而减小.当电荷点数多于90时，Method 2的保角变换误差要小于Method 1.电荷点数N=200时，Method 2的保角变换误差为4.1729×10-6，而Method 1的保角变换误差为8.7557×10-6，说明本文所采用的算法可以得到更高精度的模拟电荷Qi（i=1，2，···，N）和变换半径M.图3给出的模拟电荷点数N为200时，模拟电荷点的分布情况.

图2：数值保角变换误差（例1）

图3：模拟电荷点分布（例1）

为进一步验证本文算法的有效性，取以

为边界的内部区域的等高线.图4中，粗实线表示的边界，细线表示的等高线.电荷点数取200，利用GMRES（m）法求解约束方程，从而构造变换函数.图5为区域及其等高线的映射结果，可以看出，对于C1和C2的所围成的区域内的任意部分，经过保角变换，对应的仍然是变换后所围成区域的内部.且边界经过保角变换对应的同心圆也是变换后图像的边界，从而验证了基于GMRES（m）法的双连通区域保角变换的计算法的有效性.

图4：边界及区域等高线（例1）

图5：边界及区域等高线保角变换（例1）

例2 a1=8，b1=4，a2=4，b2=2时.类似于图2和图3，图6和图7分别给出了两种方法的数值保角变换误差和模拟电荷点的分布.从图6中可以看出，电荷点数越多，两种方法的保角变换误差越小.但是，当电荷点数大于180时，Method 2的保角变换误差明显要小于Method 1.电荷点数为300，Method 2的保角变换误差已经达到3.7814×10-13.图8中细线表示以

为边界的内部区域的等高线，粗线表示边界.从图9可以看出，Method 2可以得到了很好的保角变换结果.

图6：数值保角变换误差（例2）

图7：模拟电荷点分布（例2）

图8：边界及区域等高线（例2）

图9：边界及区域等高线保角变换（例2）

5　结束语

本文利用GMRES（m）法求解出了双连通区域保角变换模拟电荷法中的约束方程，进而提出了基于GMRES（m）法的双连通区域数值保角变换的计算法.并通过数值实验验证了新算法的有效性，且用等高线模拟了保角变换的计算结果.本算法同样可以运用解决其它的多连通区域数值保角变换问题.

［1］曹伟杰.保形变换理论及其应用［M］.上海：上海科学技术文献出版社，1988.

［2］朱满座．数值保角变换及其在电磁理论中的应用［D］.西安：西安电子科技大学，2008.

［3］林为干.微波理论与技术［M］.北京：科学出版社，1979.

［4］Singer H.A charge simulation method for calculation of high voltage field［J］.IEEE Trans.Power App.Sys.，1974，9（3）：1660-1668.

［5］Amano K.Numerical conformal mappings of exterior domains based on the charge simulation method［J］.Trans.Inform.Proc.Soc.Japan，1998，29（1）：62-72（in Japanese）.

［6］Amano K.Numerical conformal mappings of interior domains based on the charge simulation method ［J］.Trans.Inform.Proc.Soc.Japan，1988，29（7）：697-699（in Japanese）.

［7］Amano K.A bidirectional method for numerical conformal mapping based on the charge simulation method［J］.Trans.Inform.Proc.Soc.Japan，1991，28（4）：473-482（in Japanese）.

［8］Amano K.Numerical conformal mapping of doubly-connected domains based on the charge simulation method［J］.Trans.Inform.Proc.Soc.Japan，1988，29（7）：914-919（in Japanese）.

［9］Golub G H，Van Loan C F.Matrix computation（4th ed.）［M］.北京：人民邮电出版社，2014.

［10］全忠，向淑晃.基于GMRES的多项式预处理广义极小残差法［J］.计算数学，2006（4）：365-376.

［11］潘春平，王红玉，赵伟良.一种求解鞍点问题的广义对称超松弛迭代法［J］.数学杂志，2011，31（3）：569-574.

［12］Amano K.A charge simulation method for numerical conformal mapping of interior，exterior and doubly-connected domains［J］.Comput.Appl.Math.，1994，53（30）：353-370.

［13］Lu Y，Wu D，Wang Y，et al.The accuracy improvement of numerical conformal mapping using the modified Gram-Schmidt method［C］.19th Intern.Confer.Indus.Engin.Engin.Manag.，Berlin Heidelberg：Springer，2013：555-56.

［14］戴济能，杜金元.单准周期的Riemann边值问题［J］.数学杂志，2006，26（5）：579-584.

［15］Saad Y.A Flexible inner-outer preconditioned GMRES algorithm［J］.SIAM J.Sci.Comp.，1993，14（2）：461-469.

［16］Habu M，Nodera F.GMRES（m）algorithm with changing the restart cycle adaptively［C］.Proc. Alg.2000 Confer.Sci.Comp.，2000：254-263.

THE GMRES（M）METHOD FOR NUMERICAL CONFORMAL MAPPING OF DOUBLY-CONNECTED DOMAIN

LU Yi-bin1，LAI Fu-ming1，WANG Ying-zi2，WU De-an3
（1.School of Science，Kunming University of Science and Technology，Kunming 650500，China）
（2.Computer Center，Kunming University of Science and Technology，Kunming 650500，China）
（3.School of Mathematical Sciences，University of Electronic Science and Technology，Chengdu 611731，China）

In this paper，we study the GMRES（m）method for numerical conformal mapping based on charge simulation method of doubly-connected domain.In this method，using the GMRES（m）method which limited the number of Krylov subspace's maximum dimension，the linear equation of charge simulation method is solved，and the approximate conformal mapping function is constructed by using the charges and conformal mapping radius.Numerical results show that the proposed method is effective.

charge simulation method；doubly-connected domains；Krylov subspace；GMRES（m）method

MR（2010）主题分类号：65E05；30C30O241.85

A

0255-7797（2016）05-1028-07

2015-07-17接收日期：2016-01-19

国家自然科学基金资助（11461037）.

吕毅斌（1972-），男，黑龙江牡丹江，副教授，主要研究方向：科学计算和图像处理.

王樱子.

2010 MR Subject Classification：65E05；30C30