在|z|<1的环域展开
在命令行窗口输入
![](https://img-blog.csdnimg.cn/img_convert/fe1c34a8a1079a6747e000ad3fdaecd7.png)
查看taylor函数文档:https://www.mathworks.com/help/releases/R2021b/symbolic/sym.taylor.html
点击enter键,得到结果:
![](https://img-blog.csdnimg.cn/img_convert/7c70f4d6747d09e1072bfe7beedc0224.png)
整理得:
![](https://img-blog.csdnimg.cn/img_convert/ca7c33a034ce1278cd30be483ee2d908.png)
就是函数 1/(1-z)^2 (|z|<1)的级数展开式
在1<|z|<∞的环域展开
![](https://img-blog.csdnimg.cn/img_convert/c8288b00db3c495f69e0da3aa0803683.png)
命令行窗口输入:
![](https://img-blog.csdnimg.cn/img_convert/2ddc068c3093f7f4d31f393d01ea0302.png)
得到结果:
![](https://img-blog.csdnimg.cn/img_convert/508d868708b42979d8ce63a439a086c9.png)
整理得:
![](https://img-blog.csdnimg.cn/img_convert/da70b2be36298dab195ce2eca0eaef85.png)
就是函数 1/(1-z)^2 (1<|z|<∞)的级数展开式。
绘图
根据上述得到的展开式编写M程序。
在M程序文件(“.m”后缀):
%1/(1-z)^2幂级数展开及绘图表示
m=30;%与图形的平滑度(精度)相关
r=(0:2*m)'/m; %模
theta=pi*(-m:m)/m; %幅角
z=r*exp(1i*theta); %Euler欧拉表示
z(z==1)=NaN;%挖去奇点
f=1./(1-z).^2;
figure(1)%打开图形窗口
cplxmap(z,f); %原函数
colorbar('vert');%此 MATLAB 函数 在当前坐标区或图的右侧显示一个垂直颜色栏
z1=z;
z1(abs(z1)>=1)=NaN; %选定收敛域,|z|<1
f1=1;%通项和
u1=1;%通项次幂
for k=1:100%for循环,求101个项数的和
u1=u1.*z1;%累乘得到下一个高次幂(不含系数)
f1=f1+(k+1).*u1;%累加求和,逼近原函数
end
figure(2)%打开图形窗口
subplot(1,2,1);%在各个分块位置创建坐标区
cplxmap(z1,f1);%Taylor,绘制一个复变量的函数
z2=z;
z2(abs(z2)<=1)=NaN; %选定收敛域,|z|>1
f2=1./z2.^2;%通项和
u2=1./z2.^2;%通项次幂
for k=1:100%for循环,求101个项数的和
u2=u2./z2;%累除得到下一个高次幂(不含系数)
f2=f2+(k+1).*u2;%累加求和,逼近原函数
end
figure(2)%打开图形窗口
subplot(1,2,2);%在各个分块位置创建坐标区
cplxmap(z2,f2);%Lauren,洛朗展开
运行得到输出图形:
![](https://img-blog.csdnimg.cn/img_convert/4aec4377bea6b66d1b04c6e09528e37b.png)
原函数的图像
![](https://img-blog.csdnimg.cn/img_convert/b66290d7de3adc0a9d74f2a88c09d4db.png)
展开式图像:在|z|<1的环域展开(左),在1<|z|<∞的环域展开(右,部分)
![](https://img-blog.csdnimg.cn/img_convert/b5bfe35e9a02f531828952cde990bf8e.png)
在工作区查看原函数和展开式的部分值并对比,分别点击f,f1,f2
![](https://img-blog.csdnimg.cn/img_convert/782faedeef6f163d6f8dc13946168aea.png)
f= 1/(1-z)^2
![](https://img-blog.csdnimg.cn/img_convert/63a28cc2b66b680cd6c316a05c36fdbf.png)
在|z|<1的环域展开,与上图一致
![](https://img-blog.csdnimg.cn/img_convert/5b8488b1a1609f6098bd5652cfff8e14.png)
f= 1/(1-z)^2
![](https://img-blog.csdnimg.cn/img_convert/1b771bd23dbbca362bec413d10ce8f55.png)
在1<|z|<∞的环域展开,与上图一致
结束语
本文是作者第一篇文章,末流985大二学生的一次作业分享,希望得到大家的认同,如有错误或更好的解决办法,欢迎在评论区交流Thanks♪(・ω・)ノ。
我曾到CSDN上找相关文章完成作业,但是没找到自己想要的,完成以后想着帮到其他有需要的人。
可能有些步骤细节没有展示,有不懂的可以私信。