matlab共形映射程序,Matlab在共形映射中的应用

共形映射是复变函数中重要的概念之一。共形映射的方法,解决了动力学,弹性理论,静电场与磁场等方面的许多实际问题。由于该理论独特的思维方法、抽象的理论基础,使得学生学习这方面知识时存在着一定的困难,借助于MATLAB这一强大的计算和画图软件,可以使得共形映射理论的若干问题,计算得以简化,抽象的理论可视化。1Matlab在分式线性映射中的运用在最简单的共形映射中,需要用到分式线性函数。其定义是:在扩充的z平面上任意给定三个不同的点z1,z2,z3,在扩充的w平面也任意给定三个不同的点w1,w2,w3,则存在唯一的分式线性映射,把z1,z2,z3分别映射成w1,w2,w3。313121323121::zzzzzzzzwwwwwwww=(1)此即所求分式线性映射。以下分两种情况运用Matlab求解分式线性映射。(1)如果z1,z2,z3或w1,w2,w3中有一个为无穷远点,不妨设w3=,其他各点均为有限点,则显然有31312121:11:zzzzzzzzwwww=例1:求把z1=2,z2=2i,z3=1分别映为w1=3;w2=1;w3=inf的分式线性映射。程序及运行结果如下:symsz1z2z3w1w2w3z1=2;z2=2*i;z3=1;w1=3;w2=1;w3=inf;A=(z3-z1)/(z3-z2)A=-0.2000-0.4000ia=Aa=-0.2000-0.4000ia=A*w1-w2a=-1.6000-1.2000ib=z1*w2-A*w1*z2b=-0.4000+1.2000ic=A-1c=-1.2000-0.4000id=z1-A*z2d=1.2000+0.4000isymszw=(a*z+b)/(c*z+d)w=(z*((6*i)/5+8/5)+2/5-(6*i)/5)/(z*((2*i)/5+6/5)-6/5-(2*i)/5)即所求分式线性映射为w=(z*((6*i)/5+8/5)+2/5-(6*i)/5)/(z*((2*i)/5+6/5)-6/5-(2*i)/5)(2)若z1,z2,z3或w1,w2,w3不存在无穷远点,按照一般的Matlab求方程根的命令即可。例2:求将点z1=2,z2=i,z3=-2分别映射为w1=-1,w2=i,w3=1的分式线性映射。解:根据(1)式及z1=2,z2=i,z3=-2,w1=-1,w2=i,w3=1Matlab命令为:s1=sym(['(w+1)/(w-i)/2*(1-i)-(z-2)/(z-i)/-4*(-2-i)']);solve(s1,'w')ans=-(3*z^2*i+20*z-12*i)/(9*z^2+4)即所求分式线性映射为w=-(3*z^2*i+20*z-12*i)/(9*z^2+4)2Matlab在共形映射图形中的应用借助Matlab我们可以通过画出曲线的图形或是曲线在某共形映射下的像,来理解不同映射的特点。例3:作出圆周z=r在映射zwfzz1()=1=+下的像并作出w的实部与虚部的等值线。解:首先绘制映射w的图形如图1所(a)w的实部的等值线(b)w的虚部的等值线图21()w=fz的实部与虚部的等值线-2.5-2-1.5-1-0.500.511.522.5-1.5-1-0.500.511.5图1映射1()w=fz的图形示,绘制映射曲线的M文件shiyan.m为:Functionshiyan(r)r=r;t=0:0.01*pi:2*pi;z=r*exp(i*t);w=z+1./z;plot(w);title‘(w=z+1/z’);axisequal;运行shiyan(2)的结