Please wait a minute...
浙江大学学报（理学版）
浙江大学学报（理学版）  2017, Vol. 44 Issue (2): 161-167    DOI: 10.3785/j.issn.1008-9497.2017.02.007
数学与计算机科学     
单位圆到任意多边形区域的Schwarz Christoffel变换数值解法
崔建斌1, 姬安召2, 王玉风2, 于江涛2, 周华龙2
1. 陇东学院 数学与统计学院, 甘肃 庆阳 745000;
2. 陇东学院 能源工程学院, 甘肃 庆阳 745000
Numerical solution method for Schwarz Christoffel transformation from unit circle to arbitrary polygon area
CUI Jianbin1, JI Anzhao2, WANG Yufeng2, YU Jiangtao2, ZHOU Hualong2
1. Mathematics and Statistics Institute, Longdong University, Qingyang 745000, Gusu Province, China;
2. Energy Engineering Institute, Longdong University, Qingyang 745000, Gusu Province, China
 全文: PDF(892 KB)   HTML  
摘要： Schwarz Christoffel变换技术在处理某些工程问题时具有重要作用.从黎曼存在定理出发，建立了单位圆到任意多边形区域的映射函数Schwarz Christoffel变换模型，采用Levenberg-Marquardt算法求解含约束条件的非线性映射函数Schwarz Christoffel变换模型参数系统.针对映射函数中出现的奇异积分问题，对映射函数进行2次参数变换，将其化为高斯雅克比型积分，以积分路径中的奇异点为界，缩短积分路径，对子路径采用修正高斯积分方法进行计算.通过指数变换、连乘变换和累加变换，使任意初值问题均可进行迭代计算并满足初值的约束条件.提出以边长绝对误差和顶点绝对误差为迭代计算的收敛条件，并保证了映射函数的精度.给出了11顶点多边形区域映射函数的求解算例，4种方案的计算结果表明，Schwarz Christoffel变换数值解法操作简单、精度高、收敛快.
关键词： 修正高斯雅克比型积分;  单位圆;  初值变换;  Levenberg-Marquardt算法;  Schwarz Christoffel变换    
Abstract: Schwarz Christoffel transformation technique has an important role in dealing with engineering problem. Based on the Riemann's existence theorem, a transform model of Schwarz Christoffel from unit circle to arbitrary polygon area was established. Levenberg-Marquardt algorithm was used to solve parameters of the nonlinear mapping function of Schwarz Christoffel with constraint conditions. As to the mapping function in the singular integral problem, the mapping function was twice transformed to Gauss Jacobi integral, then the singular point in the integral path was taken as the boundary; The length of integral path was narrowed, and the modified Gaussian integral was used to calculate the sub path. By the exponential transform, multiplicative transform and accumulation transform process, the arbitrary initial value could be iteratively calculated and the results satisfied the constraints of the initial value. The convergence conditions of the iterative copulation based on the absolute errors of length and absolute error of vertex were put forward to ensure the precision of the mapping function. The solution of the mapping function with 11 vertex polygon region was calculated by four calculation schemes. Results of the four schemes showed that the numerical solution of Schwarz Christoffel transform is simple with high precision and fast convergence.
Key words: modified Gauss Jacobi type integral    unit circle    initial value transformation    Levenberg-Marquardt algorithm    Schwarz Christoffel transformation
收稿日期: 2016-03-01 出版日期: 2017-07-08
CLC:  O241  
基金资助: 甘肃省科技计划资助项目（1606RJZM092，1606RJYM259，1506RJYM324）.
作者简介: 崔建斌(1972-),ORCID:http://orcid.org/0000-0002-6693-3415,男,硕士,副教授,主要从事数值分析与数据挖掘研究,E-mail:cuijb0658@163.com.
服务　　
把本文推荐给朋友
加入引用管理器
E-mail Alert
RSS
作者相关文章　　
崔建斌
姬安召
王玉风
于江涛
周华龙

引用本文:

崔建斌, 姬安召, 王玉风, 于江涛, 周华龙. 单位圆到任意多边形区域的Schwarz Christoffel变换数值解法[J]. 浙江大学学报（理学版）, 2017, 44(2): 161-167.

CUI Jianbin, JI Anzhao, WANG Yufeng, YU Jiangtao, ZHOU Hualong. Numerical solution method for Schwarz Christoffel transformation from unit circle to arbitrary polygon area. Journal of ZheJIang University(Science Edition), 2017, 44(2): 161-167.

链接本文:

https://www.zjujournals.com/sci/CN/10.3785/j.issn.1008-9497.2017.02.007    或    https://www.zjujournals.com/sci/CN/Y2017/V44/I2/161

[1] 王刚,许汉珍,顾王明,等.数值许瓦尔兹-克力斯托夫变换与数值高斯-雅可比型积分[J].海军工程学院学报,1994,1(2):25-34. WANG G,XU H Z,GU W M,et al. Numerical schwarz-christoffel transformation and numerical gauss-jacobi quadrature[J]. Journal of Naval Academy of Engineering,1994,1(2):25-34.
[2] 王刚,陆小刚,顾王明.槽形内域中的数值许瓦尔兹-克力斯托夫保角变换[J].海军工程学院学报,1995,1(4):16-23. WANG G,LU X G,GU W M. Numerical schwarz-christoffel conformal mapping in channel region[J].Journal of Naval Academy of Engineering,1995,1(4):16-23.
[3] 祝江鸿.隧洞围岩应力复变函数分析法中的解析函数求解[J].应用数学和力学,2013,34(4):345-354. ZHU J H. Analytic functions in stress analysis of the surrounding rock for caverns with the complex variable theory[J]. Applied Mathematics and Mechanics,2013,34(4):345-354.
[4] 祝江鸿,杨建辉,施高萍,等.单位圆外域到任意开挖断面隧洞外域共形映射的计算方法[J].岩土力学,2014,35(1):175-183. ZHU J H,YANG J H,SHI G P,et al. Calculating method for conformal mapping from exterior of unit circle to exterior of cavern with arbitrary excavation cross-section[J]. Rock and Soil Mechanics,2014,35(1):175-183.
[5] 皇甫鹏鹏,伍法权,郭松峰.基于边界点搜索的洞室外域映射函数求解法[J].岩石力学,2011,32(5):1418-1424. HUANGFU P P,WU F Q,GUO S F. A new method for calculating mapping function of external area of cavern with arbitrary shape based on searching points on boundary[J]. Rock and Soil Mechanics,2011,32(5):1418-1424.
[6] 朱大勇,钱七虎,周早生.复杂形状洞室映射函数的新解法[J].岩石力学与工程学报,1999,18(3):279-282. ZHU D Y,QIAN Q H,ZHOU Z S. New method for calculating mapping function of opening with complex shape[J]. Chinese Journal of Rock Mechanics and Engineering,1999,18(3):279-282.
[7] 王润富.一种保角映射法及其微机实现[J].河海大学学报,1991,19(1):86-89. WANG R F. A method of conformal mapping and its computer implementation[J]. Journal of Hohai University,1991,19(1):86-89.
[8] HU C H. A software package for computing schwarz-christoffel conformal transformation for doubly connected polygonal regions[J]. ACM Transactions on Mathematical Software,1998,24(3):317-333.
[9] COSTAMAGNA E. Numerical inversion of the schwarz-christoffel conformal transformation:Strip-line case studies[J]. Microwave and Optical Technology Letters,2001,28(3):179-183.
[10] 崔建斌,姬安召,鲁洪江,等.Schwarz Christoffel变换数值解法[J].山东大学学报:理学版,2016,51(4):104-111. CUI J B,JI A Z,LU H J,et al. Numerical solution of Schwarz Christoffel transform[J]. Journal of Shandong University:Natural Science,2015,50(12):104-111.
[11] COSTAMAGNA E. A new approach to standard schwarz-christoffel formula calculations[J]. Microwave and Optical Technology Letters,2002,32(3):196-199.
[12] DRISCOLL T A. Improvements to the schwarz-christoffel toolbox for matlab[J]. ACM Transactions on Mathematical Software,2005,31(2):239-251.
[13] HOUGH D M. Asymptotic gauss jacobi quadrature error estimation for schwarz-christoffel integrals[J]. Journal of Approximation Theory,2007,146(2):157-173.
[14] NATARAJAN S, BORDAS S P A, MAHAPATRA R D. Numerical integration over arbitrary polygonal domains based on schwarz-christoffel conformal mapping[J]. International Journal for Numerical Methods in Engineering,2009,80(1):103-134.
[15] CROWDY D. The schwarz-christoffel mapping to bounded multiply connected polygonal domains[J]. Proceedings of the Royal Society,2005,146(2061):2653-2678.
[16] 刘浩.大规模非线性方程组和无约束优化方法研究[D].南京:南京航空航天大学,2008. LIU H. Research on Methods for Large-scale Nonlinear Equations and Unconstrained Optimization[D]. Nanjing:Nanjing University of Aeronautics and Astronautics,2008.
[1] 房亮, 刘三阳. 一类非线性矩阵方程的正定解[J]. 浙江大学学报（理学版）, 2019, 46(1): 1-8.
[2] 初鲁, 鲍亮, 董贝贝. 求解非埃尔米特正定方程组的广义LHSS迭代法[J]. 浙江大学学报（理学版）, 2018, 45(6): 694-697,706.
[3] 张理涛, 张国辉, 赵莹超. 基于块分裂求解线性互补问题的新模系同步多分裂方法[J]. 浙江大学学报（理学版）, 2018, 45(6): 685-693.
[4] 杨延涛. Lp空间中Lipschitz强单调算子方程解的迭代算法[J]. 浙江大学学报（理学版）, 2018, 45(4): 405-408.
[5] 张勇, 杨雪玲, 舒永录. 一类大气混沌模型的动力学分析及数值仿真[J]. 浙江大学学报（理学版）, 2018, 45(1): 18-22.
[6] 李光耀, 杨连喜, 徐晨东. 一种基于离散插值的多项式曲线逼近有理曲线的方法[J]. 浙江大学学报（理学版）, 2017, 44(6): 705-710.
[7] 张勇, 胡永才, 舒永录. 新三维非线性混沌系统的动力学特性分析[J]. 浙江大学学报（理学版）, 2016, 43(4): 401-405.
[8] 王学彬. 混合边界条件下广义二维多项时间分数阶扩散方程的解析解[J]. 浙江大学学报（理学版）, 2016, 43(4): 406-410.
[9] 李军成, 谢炜. 自动满足C2连续的带参数五次Hermite插值样条[J]. 浙江大学学报（理学版）, 2016, 43(2): 175-180.
[10] 梁仙红 . 求解 Banach空间中的非线性方程的 修正的 Chebyshev迭代方法 [J]. 浙江大学学报（理学版）, 2000, 27(1): 8-19.
[11] 吴　彪 . 冷却模型中固化层的数值计算 [J]. 浙江大学学报（理学版）, 2000, 27(1): 42-46.
[12] 刘才明 . 三棱镜折射率测量结果的不确定度评定 [J]. 浙江大学学报（理学版）, 2000, 27(1): 53-58.
[13] 孙方裕,叶　绿. 一个并行迭代的加速法 [J]. 浙江大学学报（理学版）, 2000, 27(1): 4-.
[14] 阮火军 . 递归 IFS维数公式的推广 [J]. 浙江大学学报（理学版）, 2000, 27(3): 243-246.
[15] 胡桂武，刘晓斌. 非线性不可微方程的迭代解法[J]. 浙江大学学报（理学版）, 1999, 26(1): 1-6.