呂毅斌，賴富明，王櫻子，武德安

(1.昆明理工大學理學院，云南 昆明 650500； 2.昆明理工大學計算中心，云南 昆明 650500； 3.電子科技大學數(shù)學科學學院，四川 成都 611731)

基于改進譜修正迭代法的數(shù)值保角變換計算法

呂毅斌1，賴富明1，王櫻子2，武德安3

(1.昆明理工大學理學院，云南 昆明 650500； 2.昆明理工大學計算中心，云南 昆明 650500； 3.電子科技大學數(shù)學科學學院，四川 成都 611731)

通過改進的譜修正迭代法，計算基于模擬電荷法的雙連通區(qū)域數(shù)值保角變換的約束方程，得到了新的電荷點和變換半徑，構(gòu)造了近似保角變換函數(shù).進而提出了新算法，并通過數(shù)值實驗檢驗了新算法的有效性.

數(shù)值保角變換；雙連通區(qū)域；模擬電荷法；譜修正迭代法

0 引言

保角變換的理論來源于復(fù)變函數(shù)論，其在電磁理論、電場與熱場理論、流體力學等許多領(lǐng)域均有廣泛應(yīng)用.[1-3]保角變換的主要求解方法有解析法和數(shù)值法.解析法的基礎(chǔ)是Riemann存在唯一性定理，它指出任意一個邊界多于一點的單連通區(qū)域總存在一個解析函數(shù)將其映射為單位圓盤.類似地，對于雙連通區(qū)域也存在一個解析函數(shù)將其映射成同心圓.Riemann存在唯一性定理證明了變換函數(shù)的存在性，但卻不能給出變換函數(shù)的具體表達式.相比于解析法，保角變換的數(shù)值法在解決實際工程問題的過程中具有獨特的優(yōu)勢，因此吸引許多學者對其做了大量研究.目前，保角變換的數(shù)值法主要有積分方程式法[4-6]、正交多項式法[7]和模擬電荷法[8-10].基于模擬電荷法的數(shù)值保角變換計算法是由日本數(shù)學家Amano提出來的，它具有計算簡單、不用數(shù)值積分、計算精度高等優(yōu)點.然而，Amano并沒有給出求解基于模擬電荷法的雙連通區(qū)域保角變換約束方程的具體算法.

本文針對以橢圓為邊界的雙連通區(qū)域，給出了模擬電荷點和約束點，進而構(gòu)造了約束方程.對約束方程進行分析，發(fā)現(xiàn)其是病態(tài)非奇異的.針對這一特性，采用改進的譜修正迭代法[11-12]對其進行求解，計算得到新的電荷點和變換半徑，從而構(gòu)造了高精度的近似保角變換函數(shù).相比于條件預(yù)優(yōu)法、投影法、剛性常微分方程法等求解方法，改進的譜修正迭代法能很好地改善方程組的病態(tài)性，且計算結(jié)果也是無偏的，適合于計算良態(tài)、病態(tài)和秩虧方程.最后，在改進的譜修正迭代法的基礎(chǔ)上，通過等高線模擬了雙連通區(qū)域數(shù)值保角變換的計算結(jié)果.

1 基于模擬電荷法的雙連通區(qū)域數(shù)值保角變換計算法[13-17]

設(shè)D為Z平面上由曲線C1和C2圍成的有限雙連通區(qū)域，通過基于模擬電荷法的保角變換將之映射為W平面的一個圓環(huán)μ<|w|<1，C1和C2分別是外部和內(nèi)部邊界(見圖1).

(+表示模擬電荷點，·表示約束點)

不失一般性，假定f(0)=0.當f(z)滿足正規(guī)化條件f(∞)=∞，f′(∞)>0時，f(z)可以表示成

(1)

(2)

h(z)和g(z)為共軛調(diào)和函數(shù)，根據(jù)模擬電荷法，h(z)和g(z)可以分別由

(3)

(4)

(5)

(6)

(7)

聯(lián)立(5)—(7)式，得到以Qi(i=1，2，…，N)和logM為未知量的線性方程

(8)

其中系數(shù)矩陣中的一般項aij=log|zj-ζi|.

最后，利用zi，ζi，Qi，M計算雙連通區(qū)域的數(shù)值保角變換.

2 改進的譜修正迭代法在雙連通區(qū)域數(shù)值保角變換上的應(yīng)用

方程(8)稱為基于模擬電荷法的雙連通區(qū)域數(shù)值保角變換的約束方程，簡稱為約束方程.將其寫成標準線性方程的形式

Ax=b，

(9)

其中A∈R(N+1)×(N+1)，x∈RN+1，b∈RN+1，N表示模擬電荷點數(shù).由于cond(A)=‖A‖‖A‖-1較大且det(A)≠0，因此方程(9)是病態(tài)非奇異的，b的微小擾動會引起解x的很大變化，從而不能得到高精度的未知電荷Q1，Q2，…，QN和保角變換半徑M.

改進的譜修正迭代法[18]能夠很好地改善上述方程的病態(tài)性.對方程(9)進行預(yù)處理，兩邊同時乘以AT有

ATAx=ATb，

(10)

其中ATA是對稱正定的.顯然x=(ATA)-1ATb，但是cond(ATA)很大，因此對ATA求逆會不穩(wěn)定.

為克服這一缺點，假定E是(N+1)×(N+1)的單位矩陣，在(10)式的兩邊同時加上μx，得到

(ATA+μE)x=ATb+μx，

(11)

其中μ為修正因子(0<μ<1).進一步得到改進的譜修正迭代法的迭代公式

xk=(ATA+μE)-1(ATb+μxk-1).

(12)

另一方面，迭代(12)式可以改寫成

xk=μ(ATA+μE)-1xk-1+(ATA+μE)-1ATb.

(13)

設(shè)x*為約束方程的精確解，那么x*=μ(ATA+μE)-1x*+(ATA+μE)-1ATb.根據(jù)(13)式，

(14)

表1 改進的譜修正迭代算法

表1中，MaxIter和SolRD分別表示最大迭代次數(shù)和相對誤差.ATA+μE為對稱正定矩陣，所以可以按LDLT分解，其中L為單位下三角矩陣，D為對角矩陣.LDLT分解大約需要n3/6次乘除法，大約是直接LU分解法計算量的一半.相比于LTL分解，LDLT分解不需要開方計算，使得結(jié)果精度更高.[19]

綜合前面的分析，基于改進的譜修正迭代法的雙連通數(shù)值保角變換計算法分為以下幾步：

步驟2 構(gòu)造約束方程，利用改進的譜修正迭代法計算約束方程(見表1)，得到電荷Q1，Q2，…，QN和保角變換半徑M.

3 數(shù)值實驗

圖2 雙連通區(qū)域模擬電荷點分布

圖3 誤差曲線

根據(jù)前面求出的電荷和保角變換半徑，可構(gòu)造近似保角變換函數(shù)f(z).取區(qū)域D的等高線，如圖4，其中粗實線代表邊界C1和C2，細線代表等高線.圖5表示由近似保角變換函數(shù)映射后得到的圖像，從中可以看出，邊界C1和C2變換后的圖像為同心圓，等高線與邊界組成的內(nèi)部區(qū)域仍然對應(yīng)同心圓圍成的內(nèi)部區(qū)域.進一步證明了基于改進譜修正迭代法的雙連通區(qū)域數(shù)值保角變換計算法的有效性.

圖4 區(qū)域及其等高線

圖5 區(qū)域及其等高線的保角變換

4 結(jié)論

本文提出了一種求解雙連通區(qū)域數(shù)值保角變換約束方程的新算法，并通過數(shù)值實驗驗證了新算法的可行性.利用新方法得到的數(shù)值保角變換的精度比較高，且通過等高線模擬了雙連通區(qū)域數(shù)值保角變換的計算結(jié)果.這種方法同樣可考慮用于多連通區(qū)域的數(shù)值保角變換.

[1] 林為干.微波理論與技術(shù)[M].北京：科學出版社，1979：64-120.

[2] 曹偉杰.保形變換理論及其應(yīng)用[M].上海：上?？茖W技術(shù)文獻出版社，1988：193-231.

[3] 朱滿座.數(shù)值保角變換及其在電磁理論中的應(yīng)用[D].西安：西安電子科技大學，2008.

[4] SYMM G T.An integral equation method in conformal mapping[J].Numer Math，1966，9(3)：250-258.

[5] SANGAWIA A W K，MURIDA A H M，NASSERD M M S.Linear integral equations for conformal mapping of bounded multiply connected regions onto a disk with circular slits[J].Appl Math Comput，2011，218(5)：2055-2068.

[6] SYMM G T.Conformal mapping of doubly-connected domains[J].Numer Math，1969，13(5)：448-457.

[7] KOKKINOS C A，PAPAMICHAEL N，SIDERIDIS A B.An orthonormalization method for the approximate conformal mapping of multiply-connected domains [J].IMAJ Numer Anal，1990，9：343-359.

[8] AMANO K.Numerical conformal mappings of exterior domains based on the charge simulation method [J].Trans Inform Process Soc Japan，1998，29(1)：62-72.

[9] AMANO K.Numerical conformal mappings of interior domains based on the charge simulation method [J].Trans Inform Process Soc Japan，1988，29(7)：697-699.

[10] AMANO K.Numerical conformal mapping of doubly-connected domains based on the charge simulation method [J].Trans Inform Process Soc Japan，1988，29(7)：914-919.

[11] 王新洲，劉丁酉.最小二乘估計中法方程的迭代解法[J].湖北民族學院學報(自然科學版)，2002，20(3)：1-4.

[12] 劉斌，龔健雅，江萬壽，等.基于嶺參數(shù)的譜修正迭代法及其在有理多項式參數(shù)求解中的應(yīng)用[J].武漢大學學報(信息科學版)，2012，37(4)：399-402.

[13] 鄭志強.單位圓到任意曲線保角變換的近似計算方法[J].應(yīng)用數(shù)學和力學，1992(5)：449-457.

[14] AMANO K.A charge simulation method for the numerical conformal mapping of interior，exterior and doubly-connected domains [J].Journal of Computational and Applied Mathematics，1994，53：354-361.

[15] LU Y，WU D，WANG Y，et al.The accuracy improvement of numerical conformal mapping using the modified Gram-Schmidt method [M]//QI E R，SHEN J，DOU R L.The 19th International Conference on Industrial Engineering and Engineering Management.Berlin：Springer，2013：555-563.

[16] 姚國梅，呂毅斌，王櫻子.數(shù)值保角變換的新算法[J].價值工程，2014，31：308-310.

[17] 代榮恒，呂毅斌，王櫻子.基于LMS法的數(shù)值保角變換計算法[J].曲阜師范大學學報(自然科學版)，2015，41(5)：23-27.

[18] 潘朝毅.譜修正法的收斂分析及其改進[J].四川教育學院學報，2009，25(5)：112-113.

[19] GOLUGB G H，VAN LOAN C F.Matrix computations[M].Baltimore：The Johns Hopkins University Press，2012：153-165.

(責任編輯：李亞軍)

Theimprovementofiterationmethodbycorrectingcharacteristicvaluesfornumericalconformalmapping

LU Yi-bin1，LAI Fu-ming1，WANG Ying-zi2，WU De-an3

(1.Faculty 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 Science，University of Electronic Science and Technology，Chengdu 611731，China)

A new method for numerical conformal mapping of doubly-connected domain is considered.This method calculates the constraint equations of numerical conformal mapping based on charge simulation method of doubly-connected domain using the improvement of iteration method by correcting characteristic values，and the approximate conformal mapping function is constructed using new charge points and conformal mapping radius.The efficiency of the proposed method is illustrated by some numerical results.

numerical conformal mapping；doubly-connected domain；charge simulation method；the iteration method by correcting characteristic values

1000-1832(2017)03-0042-05

10.16163/j.cnki.22-1123/n.2017.03.010

2015-12-02

國家自然科學基金資助項目(11461037).

呂毅斌(1972—)，男，副教授，主要從事科學計算與圖像處理研究；通信作者：王櫻子(1972—)，女，碩士，講師，主要從事科學計算與數(shù)學應(yīng)用軟件設(shè)計研究.

O 241 [學科代碼] 110·6199

A