python绘制超炫酷动态Julia集示例
作者:微小冷 发布时间:2023-10-04 12:58:49
前言
此Julia非彼Julia,指的是对于某复数 c c c,使得迭代式 f ( z ) = z 2 + c f(z)=z^2+c f(z)=z2+c收敛的复数 z z z的集合。例如,当 c = 0 c=0 c=0时,那么其收敛区间为 z 2 < 1 z^2<1 z2<1的单位圆,对应的 c c c的Julia集便是 cos θ + i sin θ \cos\theta+i\sin\theta cosθ+isinθ。
Mandelbrot集
特别地,当 c = z c=z c=z的初始值时,符合收敛条件的 z z z的便构成大名鼎鼎的Mandelbrot集
在上图中,颜色表示该点的发散速度,可以理解为开始发散时迭代的次数。其生成代码也非常简单,唯一需要注意的是,由于使用了大量的矩阵运算,故使用了cupy
,如果电脑没装cuda
,只需将所有的cp
改为np
即可。
# 这些代码会在后面的程序中反复调用,不再说明
import numpy as np
import time
import matplotlib.pyplot as plt
import cupy as cp
#生成z坐标 x0,y0 为起始点, nx,ny为点数, delta为点距
def genZ(x0, y0, nx, ny, delta):
real, img = cp.indices([nx,ny])*delta
real += x0
img += y0
return real.T+img.T*1j
#获取Julia集,n为迭代次数,m为判定发散点,大于1即可
def getJulia(z,c,n,m=2):
t = time.time()
z,out = z*1, cp.abs(z)
c = cp.zeros_like(z)+c
for i in range(n):
absz = cp.abs(z)
z[absz>m]=0#对开始发散的点置零
c[absz>m]=0
out[absz>m]=i#记录发散点的发散速度
z = z*z + c
print("time:",time.time()-t)
return out
z1 = genZ(-2,-1.5,1000,1000,0.003)
mBrot = getJulia(z1,z1,50)
plt.imshow(mBrot.get(), cmap=plt.cm.jet)
plt.show()
如果对其生成过程感兴趣,那么可以观察一下随着迭代次数的增加,图像的变化情况
代码如下。
from matplotlib import animation
fig = plt.figure()
fig.subplots_adjust(top=1, bottom=0, left=0, right=1)
ax = plt.subplot()
def getJulias(z,c,n,m=2):
z,out = z*1, cp.abs(z)
c = cp.zeros_like(z)+c
J = []
for i in range(n):
z = z*z + c
absz = cp.abs(z)
z[absz>m]=0#对开始发散的点置零
c[absz>m]=0
out[absz>m]=i#记录发散点的发散速度
im = ax.imshow(out.get(),cmap=plt.cm.jet, animated=True)
ax.set_axis_off()
J.append([im])
return J
N = 75 #迭代次数
z1 = genZ(-2,-1.5,1000,1000,0.003)
J = getJulias(z1,z1,N)
ani = animation.ArtistAnimation(fig, J, interval=50, blit=True,repeat_delay=1000)
plt.show()
ani.save('julias.gif',writer='imagemagick')
无限缩放
Mandelbrot集的分形特征意味着我们所生成的图片可以无限放大,但是受到栅格化尺寸的影响,手动的放大并不会更改其真实尺寸,
为了照顾观感,将缩放中心作为图像的中心,所以对genZ
函数进行修改。如果选取(-0.75,-0.2)
作为缩放中心,则其变化如下
代码为
from matplotlib import animation
# 生成z坐标 xy=np.array([xc,yc]) 为起始点,
# nxy=np.array([nx,ny])为点数, delta为点距
def genZbyCenter(xy,nxy,delta):
x0, y0 = xy-np.array(nxy)*delta/2
return genZ(x0,y0,*nxy,delta)
mBrots = []
xy = [-0.75,-0.2]
nxy = [1000,1000]
delta0 = 0.003 #初始宽度
fig = plt.figure()
fig.subplots_adjust(top=1, bottom=0, left=0, right=1)
ax = plt.subplot()
for n in range(50):
z1 = genZbyCenter(xy,nxy,1.1**(-n)*delta0)
out = getJulia(z1,z1,40)
im = ax.imshow(out.get(),cmap=plt.cm.jet, animated=True)
ax.set_axis_off()
mBrots.append([im])
ani = animation.ArtistAnimation(fig, mBrots, interval=50, blit=True)
plt.show()
ani.save('zoom.gif',writer='imagemagick')
Julia集
如果更改c
的值,那么就能得到一个变化着的Julia集,例如,下面选取一条直线
y = x y=x y=x
上面的Julia集,效果如图所示
代码为
z1 = genZ(-2,-1.5,1000,1000,0.003)
fig = plt.figure()
fig.subplots_adjust(top=1, bottom=0, left=0, right=1)
ax = plt.subplot()
mBrots = []
for x in np.arange(0.5,1,0.01):
c = x + x*1j
out = getJulia(z1,c,40)
im = ax.imshow(out.get(),cmap=plt.cm.jet, animated=True)
ax.set_axis_off()
mBrots.append([im])
ani = animation.ArtistAnimation(fig, mBrots, interval=50)
plt.show()
ani.save('julia.gif',writer='imagemagick')
来源:https://blog.csdn.net/m0_37816922/article/details/121880410
猜你喜欢
- 密码算法程序设计实践选的SHA-1。在写的过程中遇到一丢丢关于python移位的问题,记录一下。SHA-1其中第一步需要填充消息。简单阐述一
- Index.asp:程序代码<html><head><meta http-equiv="Conten
- 一、导入所需的库import turtleimport randomfrom math import *二、生成斐波那契数列斐波那契数列是指
- 傅里叶变换是在高数是一个很重要的知识点,今天将结合Python代码实现傅立叶变换。傅立叶变换我们平时是如何去分解一个复杂的问题呢?一个经典的
- 实战场景 本篇博客学习字体反爬,涉及的站点是实习 x,目标站点地址直接百度搜索即可。可以看到右侧源码中出现了很多&ldqu
- 最近项目需要抓包功能,并且抓包后要对数据包进行存库并分析。抓包想使用tcpdump来完成,但是tcpdump抓包之后只能保存为文件,我需要将
- 问题:MySQL 5.0不能使用自动增加字段。具体原因:引起的原因是在MySQL5.0.2版本中的The Server SQL Mode 中
- 1、某汽车网站地址2、使用firefox查看后发现,此网站的信息未使用json数据,而是简单那的html页面而已3、使用pyquery库中的
- Python实现AES算法密码学课程老师留的作业,我觉得用python实现更简单,就用python写了一个加解密的程序。程序分成三个部分,一
- 一个已知管用的方法是,使用session_set_save_handler,接管所有的session管理工作,一般是把session信息存储
- 框架thinkphp 版本:3.2.3 内容:数据库操作1. 多表查找一条数据M('a表')->join("
- asp 在线备份 恢复 sql server 数据库,对于远程没有提供sql server远程连接或打包下载的朋友是个临时解决方法,对于大数
- 1. 直接打印错误try:# your codeexcept KeyboardInterrupt:print("quit"
- python——pip install xxx报错SyntaxError: invalid syntax在安装好python后,进入pyth
- 使用Keras如果要使用大规模数据集对网络进行训练,就没办法先加载进内存再从内存直接传到显存了,除了使用Sequence类以外,还可以使用迭
- matplotlib及相关cmap参数的取值在matplotlib中对于图片的显示有如下方法(这不是重点), 其中有cmap=&ls
- 1. 原地排序:采用sort()方法,按照指定的顺序排列数据后用排序后的数据替换原来的数据(原来的顺序丢失),如:>>>
- 回顾在前面的系列章节中,我们创建了一个数据库并且学着用用户和邮件来填充,但是到现在我们还没能够植入到我们的程序中。 两章之前,我
- 前言由于Django是 * 站,所有每次请求均会去数据进行相应的操作,当程序访问量大时,耗时必然会更加明显,最简单解决方式是使用:缓存,缓存
- 列表是最常用的Python数据类型,它可以作为一个方括号内的逗号分隔值出现。列表的数据项可以是不同的类型,可以是字符串,可以是数字类型,甚至