【matlab】求1/(1-z)^2的级数展开式

最新推荐文章于 2023-01-09 13:46:52 发布

404Gloria

最新推荐文章于 2023-01-09 13:46:52 发布

阅读量1.3k

点赞数

分类专栏： matlab 文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/S_Glora/article/details/128586119

版权

matlab 专栏收录该内容

1 篇文章 0 订阅

订阅专栏

文章详细介绍了如何对1/(1-z)^2进行泰勒展开和洛朗展开，利用已知的1/(1-z)展开式，推导出1/(1-z)^2的级数表示，并通过MATLAB代码展示了在不同收敛域内的图形表示。

摘要由CSDN通过智能技术生成

求 $1/(1-z)^2$ 的级数展开式。

已知 $1/(1-z) = 1+z+z^2+z^3+...+z^n,|z|<1$
则 $1/(1-z)^2=1+2z+3z^2+...+nz^{n-1},|z|<1$

令 $t = 1/ z$ ，可得洛朗展开
$t^2/(1-t)^2=t^2*[1/(1-t)^2]=t^2+2t^3+3t^4+...+nt^{n+1}$
即为
$1/z^2+2/z^3+3/z^4+...+n/z^{n+1}$

以下为matlab代码：

m=30;
r=(0:2*m)'/m;
theta=pi*(-m:m)/m;
z=r*exp(i*theta);
z(find(z==1))=NaN;
figure(1)
cplxmap(z,1./(1-z).^2);
2
title('1/(1-z)^2');
colorbar('vert');
%泰勒展开如下：
z1=z;
z1(abs(z1)>=1)=NaN; %选定收敛域
f1=1;u1=1;
for k=1:100
 u1=u1.*z1;
 f1=f1+u1.*(k+1);
end
figure(2)
subplot(1,2,1);cplxmap(z1,f1);%Taylor
%洛朗展开如下：
z2=z;
z2(abs(z2)<=1)=NaN; %选定收敛域
f2=1./z2./z2;u2=1./z2./z2;
for k=1:100
 u2=u2./z2;
 f2=f2+u2.*(k+1);
end
figure(2)
subplot(1,2,2);cplxmap(z2,f2);%Laure