matlab软件imag函数_「复变函数与积分变换」基本计算代码

9c89d72e6efb02890d0d7eaf30a01ca4.png
使用了Matlab代码,化简平时遇到的计算问题,也可以用于验算结果
来自211工科专业2学分复变函数与积分变换课程

求复角主值

sym(angle(待求复数))%公式
sym(angle(1+sqrt(3)*i))%举例

代入化简

equation?tex=z%3D%5Cfrac%7B1%7D%7Bz%7D 代入关于z的函数f(z)中并化解,用于公式法
计算无穷远点的留数
syms f z 
f=1/(z*(z+1));%某一函数f(z)
f=subs(f,z,1/z);%可以把1/z换成任意关系式
f=simplify(f/(z^2));

一维积分计算

简单的积分计算,matlab中不加;即可显示符号型变量,使用double()函数可以转化为数值型显示。

syms x z(x) %定义关于x的函数z
u=@(x)3*(i*x^2-x^2);%某函数
jf1=int(u,x,[0,i]);%这里可以是定积分、不定积分、变上限积分,只需要在[]内更改
jf1%符号型显示
double(jf1)%数值型显示

级数求和

用于计算级数是否收敛、条件收敛、绝对收敛。对目标级数s(n)判断是否收敛直接使用symsum()级数求和函数,如果需要判断是否绝对收敛时使用abs()函数取模即可。

syms s n
s=@(n)(8^n)/factorial(n);%factorial()阶乘函数
s=@(n)abs(((-1)^n)/n);
qh=symsum(s,n,1,inf)
double(qh)

二次积分

连续计算两次积分的基础计算,详细见「概率论」部分

低薪小明:「概率论」基础计算代码​zhuanlan.zhihu.com
378a54c69f2400050827d7b9e6a4d776.png
syms x y z(x,y) 
u=@(x,y)6-x-y;
jf1=int(u,y,[2,3-x]); 
jf2=int(jf1,x,[0,1]);
%jf=int(z,x,[0,inf]);
jf1
jf2
%double(jf)

一阶、二阶微分计算

经典题型:由u(x,y)构造解析函数f(z),首先对u相应的求导,然后由于解析函数的性质对v积分,最终化简得到f(z).

clc
syms x y u 
u=@(x,y)x^2-y^2+x*y;
ux=diff(u,x)%u对x求导
uy=diff(u,y)%u对y求导
vx=diff(u,x,2)
vy=diff(u,y,2)
vjfy=int(ux,y)
vx=diff(vjfy,x)

普通求导并代入计算

结合高阶求导公式计算结果,先求导,再将相应的z0代入求值,非常简便。

syms f z
f=1/sin(z);
df=diff(f,z)
subs(df,z,0)

普通极限

显而易见,这用来求一个函数的极限

  1. 用于判断极限是否存在
  2. 判断是否收敛
  3. 判断奇点类型
  4. 求留数
  5. 等等
syms f n
f=@(n)(1+1/n)^n;
limit(f,n,inf)%这里也可以计算在某一点的极限,不一定是无穷∞

三角极限

用于含三角函数的留数计算,不常用

syms f z 
for k=0:1:10f=@(z)(z-k*pi)/sin(z);limit(f,z,k*pi)
end

对圆积分

在圆上的积分,用于第三章、第五章的积分计算,不过使用柯西积分公式、高阶求导公式或留数公式等等也能算出。

syms seita
u=@(seita)exp(seita.*j);%圆的指数是方程
f=@(seita)exp(-1*j*seita);%被积函数
I=int(f*diff(u,seita),seita,0,pi);

保形映射公式计算

针对第六章知道几对映射点求保形映射后的点,只是简单的化简计算,在代码中输入初始点就行。

clc
syms z w(z);
w=(z+2)*i/(z-2);
w1=subs(w,z,0)
scatter(real(w1),imag(w1));
hold on
w2=subs(w,z,1)
scatter(real(w2),imag(w2));
hold on
w3=subs(w,z,-1)
scatter(real(w3),imag(w3));
hold on
w4=subs(w,z,-i)
scatter(real(w4),imag(w4));
hold on
%zuo=expand((w+i)/((w+3*i))*4)
%you=expand(z/(2*(z-1)))

保形映射画图

针对第六章对复平面某一区域求保形映射后的区域,在第一行的x和y分别输入21个表示该区域边界的离散点,并在循环中输入映射函数,即可得到映射后图形。

不过这里有一个限制条件:初始区域只能是方形区域,对圆形区域或者三角区域无法处理,欢迎大家提出宝贵的改进意见。

x=-10:1:10;y=linspace(-pi,0,21);w=zeros(1,441);
for j=1:20for k=1:21w((j-1)*21+k)=exp(x(j)+y(k)*i);end
end
realw=real(w);
imagw=imag(w);
scatter(realw,imagw);

ec79e6d3b5877fbb221b019ed48bdfec.png

4bff80d9a17a1d47be36911f78a159c9.png

7dff5d56c7d453fdd79be5a1871c0a98.png

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/388893.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

数据科学 python_为什么需要以数据科学家的身份学习Python的7大理由

数据科学 pythonAs a new Data Scientist, you know that your path begins with programming languages you need to learn. Among all languages that you can select from Python is the most popular language for all Data Scientists. In this article, I will cover 7 r…

[luoguP4142]洞穴遇险

https://www.zybuluo.com/ysner/note/1240792 题面 戳我 解析 这种用来拼接的奇形怪状的东西,要不就是轮廓线\(DP\),要不就是网络流。 为了表示奇数点(即\((xy)\%21\))的危险值,把该点拆为两个点,连一条边长…

飞信虚拟机

做完了一个图片处理软件,突然想到上次上网看到C#程序脱离.NET FRAMEWORK运行的文章,于是决定自己动手试一下。 之前看到有用别的方法来实现的,但我还是选择了现在比较流行的软件飞信中带的VMDotNet,也就是所谓的.NET FRAMEWORK虚拟机吧。相信有很多人也已…

django的contenttype表

https://blog.csdn.net/aaronthon/article/details/81714496 这篇文章已经非常详细了,供自己以后忘了...回看...... 总结: 当一张表和多个表FK关联,并且多个FK中只能选择其中一个或其中n个时,可以利用contenttype,固定用三个字段…

视频播放问题和提高性能方案

1.Five symptoms of poor video performance 1.1 视频加载缓慢 ​Perceived Wait Time Time to first frame (TTFF): ​ 播放开始所需的adaptive bitrate(ABR)流媒体段的数量。(我们稍后将对此进行更详细的讨论。) ​ 视频请求发送到视频加载之间的时间(即接收到足够的数据…

rabbitmq 不同的消费者消费同一个队列_RabbitMQ 消费端限流、TTL、死信队列

消费端限流1. 为什么要对消费端限流假设一个场景,首先,我们 Rabbitmq 服务器积压了有上万条未处理的消息,我们随便打开一个消费者客户端,会出现这样情况: 巨量的消息瞬间全部推送过来,但是我们单个客户端无法同时处理这…

动量策略 python_在Python中使用动量通道进行交易

动量策略 pythonMost traders use Bollinger Bands. However, price is not normally distributed. That’s why only 42% of prices will close within one standard deviation. Please go ahead and read this article. However, I have some good news.大多数交易者使用布林…

css3 变换、过渡效果、动画

1 CSS3 选择器 1.1 基本选择器 1.2 层级 空格 > .itemli ~ .item~p 1.3 属性选择器 [attr] [attrvalue] [attr^value] [attr$value] [attr*value] [][][] 1.4 伪类选择器 :link :visited :hover :active :focus :first-child .list li:first-child :last-chi…

webservice 启用代理服务器

您会发现你写完了一个webservice在调用的时候发现怎也没办法调用,一个简单的webservice怎么不能使用,一肚子的怨恨,哈哈您可能没有为webservice设置代理。 下面就给您写个调用的用例和大家分享下。其实很简单,但是你没有想到的时…

mysql常用的存储引擎_Mysql存储引擎

什么是存储引擎?关系数据库表是用于存储和组织信息的数据结构,可以将表理解为由行和列组成的表格,类似于Excel的电子表格的形式。有的表简单,有的表复杂,有的表根本不用来存储任何长期的数据,有的表读取时非…

android studio设计模式和文本模式切换

转载于:https://www.cnblogs.com/judes/p/9437104.html

高斯模糊为什么叫高斯滤波_为什么高斯是所有发行之王?

高斯模糊为什么叫高斯滤波高斯分布及其主要特征: (Gaussian Distribution and its key characteristics:) Gaussian distribution is a continuous probability distribution with symmetrical sides around its center. 高斯分布是连续概率分布,其中心周…

C# webbrowser 代理

百度,google加自己理解后,将所得方法总结一下: 方法1:修改注册表Software//Microsoft//Windows//CurrentVersion//Internet Settings下 ProxyEnable和ProxyServer。这种方法适用于局域网用户,拨号用户无效。 1p…

C MySQL读写分离连接串_Mysql读写分离

一 什么是读写分离MySQL Proxy最强大的一项功能是实现“读写分离(Read/Write Splitting)”。基本的原理是让主数据库处理事务性查询,而从数据库处理SELECT查询。数据库复制被用来把事务性查询导致的变更同步到集群中的从数据库。当然,主服务器也可以提供…

golang 编写的在线redis 内存分析工具 rma4go

redis 内存分析工具 rma4go redis是一个很有名的内存型数据库,这里不做详细介绍。而rma4go (redis memory analyzer for golang) 是一个redis的内存分析工具,这个工具的主要作用是针对运行时期的redis进行内存的分析,统计redis中key的分布情…

从Jupyter Notebook到脚本

16 Aug: My second article: From Scripts To Prediction API8月16日:我的第二篇文章: 从脚本到预测API As advanced beginners, we know quite a lot: EDA, ML concepts, model architectures etc…… We can write a big Jupyter Notebook, click “Re…

【EasyNetQ】- 使用Future Publish调度事件

许多业务流程要求在将来某个日期安排事件。例如,在与客户进行初次销售联系后,我们可能希望在将来的某个时间安排跟进电话。EasyNetQ可以通过其Future Publish功能帮助您实现此功能。例如,这里我们使用FuturePublish扩展方法来安排未来一个月的…

Java这些多线程基础知识你会吗?

0、并发和并行、进程核线程、多进程和多线程的区别: (这里的时间和时刻上的概念同物理上的一样) 并发:在一段时间内多个任务同时执行,或者说是在一段很短的时间内可以执行多条程序指令,微观上看起来好像是可…

MySQL set names 命令_mysql set names 命令和 mysql 字符编码问题

先看下面的执行结果:(rootlocalhost)[(none)]mysql>show variables like character%;---------------------------------------------------------------------------------------| Variable_name | Value |---------------------------------------------------…

设置Proxy Server和SQL Server实现数据库安全

首先,我们需要了解一下SQL Server在WinSock上定义协议的步骤: 1. 在”启动”菜单上,指向”程序/Microsoft Proxy Server”,然后点击”Microsoft Management Console”。 2. 展开”Internet Information Service”,再展开运行Proxy…