大量的小气泡。

内容

域

基本域四单元广场,

$ $ \马克斯{(| | x、y | |)} \ le 1 \ z = x + iy $ $

颜色

我可以用HSV colormap。

hsv_bubbles (@ (z) z)

但我更喜欢“周期性parula”parulacolormap与其反向连接,[parula;flipud (parula)]。

泡沫(@ (z) z)

权力

z ^ 2

泡沫(@ (z) z ^ 2)。

z ^ 3

泡沫(@ (z) z ^ 3)。

z ^ 9

泡沫(@ (z) z ^ 9)。

1 / z

泡沫(@ (z) 1. / z)

sqrt (z)

一个复数有两根。一个是正确的半平面。

泡沫(@sqrt)

-√(z)

在左边,另一个是半平面。

泡沫(@ (z) -√(z))

三角函数

罪(z)

泡沫(@sin)

因为(z)

泡沫(@cos)

谭(z)

泡沫(@tan)

床(z)

泡沫(@cot)

指数

exp (z)

这极角是介于1和+ 1弧度。

泡沫(@exp)

exp(π* z) / exp(π)

这些填写整个复杂的飞机。

泡沫(@ (z) exp(π* z) / exp(π))

日志(z)

泡沫(@log)

多项式和理性

z ^ 3 - z

泡沫(@ (z) z。^ 3 z)

5 / (z ^ 5 z / 5)

泡沫(@ (z) .5. / (z。^ 5 z / 5))

必要的奇点。

exp (1 / (8 z ^ 2)

泡沫(@ (z) exp (-1. / (8 * z ^ 2)))

测试

什么是发生在顶部的动画的这篇文章吗?如果你认为你知道,或者即使你只知道答案的一部分,提交评论。我要第一个奖,或最好的解决方案。

代码

也可从https://blogs.mathworks.com/cleve/files/bubbles.m。

类型泡沫

函数泡沫(F) %泡沫。颜色的复值函数F (z), (@sin) % %仓库泡沫泡沫(@ (z) .5. / (z。^ 5 z / 5))如果输入参数个数< 1 F = @ (z) z;端轴(1.5 *[1 1 1 1])轴方形盒cla m = 256;colormap = [parula (m); flipud (parula (m)));n = 25;s = 1:2 / (n - 1): 1;(x, y) = meshgrid(年代);z = x + y * 1;圆= exp((0:32) / 16 *π* 1)/ n;w = F (z); r = abs(w); theta = angle(w)+pi; scale = 20; for k = 1:n for j = 1:n t = m*theta(k,j)/pi/scale; idx = ceil(scale*t+realmin); color = colormap(idx,:); p = w(k,j) + r(k,j)*circle; patch(real(p),imag(p),color) end end titleF = char(F); if (titleF(1) == '@') titleF(1:4) = []; end title(titleF) snapnow end

得到了MATLAB代码

发表与MATLAB®R2019b