用PolarDB生成“上帝的指纹”

本文涉及的产品
云原生数据库 PolarDB MySQL 版,Serverless 5000PCU 100GB
简介: 分形法则被誉为神性法则, 因为它的公式极其简单, 但是能产生无穷无尽的自相似性. 例如通过分形公式产生的曼德勃罗集, 被成为上帝的指纹.本文介绍了上帝指纹的生成算法, 以及用PolarDB来生成"上帝的指纹".

背景

分形法则被誉为神性法则, 因为它的公式极其简单, 但是能产生无穷无尽的自相似性. 例如通过分形公式产生的曼德勃罗集, 被成为上帝的指纹.

PS: 如何部署PolarDB开源版, 请参考: 《如何用 PolarDB 证明巴菲特的投资理念》

1、分形公式:

z=z^2+c    
z可以是单数、也可以是复数、甚至可以是多维数.

2、复数运算:

(a+bi)(c+di) = (ac-bd)+(ad+bc)i    
(a+bi) + (c+di) = (a+c)+(b+d)i    
(a+bi) - (c+di) = (a-c)+(b-d)i

3、如何产生"上帝的指纹"?

当z的初始值固定时, 取不同的c, z的发散速度是不一样的, 例如某个c, 100次之后z就趋于无穷大, 某个c, 不管迭代多少次, z都在某个范围内跳动.

以一块白布为底图, 把c的坐标映射到这个二维白布中, 用颜色来标记该c值引起的z值发散速度, 越黑表示z值越不会被发散.

我们固定z0=0,那么对于不同的复数c,函数的迭代结果也不同。由于复数c对应平面上的点,因此我们可以用一个平面图形来表示,对于某个复数c,函数f(z)=z^2+cz0=0开始迭代是否会发散到无穷。我们同样用不同颜色来表示不同的发散速度,最后得出的就是Mandelbrot集分形图形:

先来个简单的, 训练一下c和z0固定时, 如何产生z的迭代值:

WITH RECURSIVE t(n, zr, zi, cr, ci) AS (      
    VALUES (1, 0::float8, 0::float8, 0.1::float8, 0.1::float8)      
  UNION ALL      
    SELECT n+1, zr*zr - zi*zi + cr, zr*zi + zi*zr + ci, cr, ci FROM t WHERE n < 100      
)      
SELECT n,zr,zi FROM t;
do language plpgsql $$    
declare    
  zr float8 := 0;    
  zi float8 := 0;    
  cr float8 := 0.1;    
  ci float8 := 0.1;    
  tmpr float8;    
  tmpi float8;    
begin    
  for i in 1..100 loop    
    raise notice '%, %i', zr, zi;    
    tmpr := zr*zr - zi*zi + cr;    
    tmpi := zr*zi + zi*zr + ci;    
    zr := tmpr;    
    zi := tmpi;    
  end loop;    
    raise notice '%, %i', zr, zi;    
end;    
$$;
NOTICE:  0, 0i  
NOTICE:  0.1, 0.1i  
NOTICE:  0.1, 0.12000000000000001i  
NOTICE:  0.0956, 0.12400000000000001i  
NOTICE:  0.09376336, 0.12370880000000001i  
NOTICE:  0.09348770048104961, 0.123198705499136i  
NOTICE:  0.0935620291045716, 0.12303512735871254i  
NOTICE:  0.09361621072599009, 0.12302283233364109i  
NOTICE:  0.09362937763530182, 0.12303386279170858i  
NOTICE:  0.093629128962925, 0.1230391680025096i  
NOTICE:  0.09362777692760627, 0.12304010025679593i  
NOTICE:  0.0936272943412032, 0.1230399421199872i  
...  
基本不发散, 所以这个 0.1,0.1i 的C对应复平面的坐标就用黑色的像素点表示.

把这个函数扩展到整个复数范围。对于复数z0=x+iy,取不同的x值和y值,函数迭代的结果不一样:

对于有些z0,函数值约束在某一范围内;

而对于另一些z0,函数值则发散到无穷。

由于复数对应平面上的点,因此我们可以用一个平面图形来表示,对于哪些z0函数值最终趋于无穷,对于哪些z0函数值最终不会趋于无穷。

我们用深灰色表示不会使函数值趋于无穷的z0;对于其它的z0,我们用不同的颜色来区别不同的发散速度。由于当某个时候|z|>2时,函数值一定发散,因此这里定义发散速度为:使|z|大于2的迭代次数越少,则发散速度越快。这个图形可以编程画出。

f(z)=z^2+(-0.75+0i)时的Julia集:

我们固定z0=0,那么对于不同的复数c,函数的迭代结果也不同。由于复数c对应平面上的点,因此我们可以用一个平面图形来表示,对于某个复数c,函数f(z)=z^2+c从z0=0开始迭代是否会发散到无穷。

我们同样用不同颜色来表示不同的发散速度,最后得出的就是Mandelbrot集分形图形(也就是上帝的指纹):

4、用PolarDB来生成"上帝的指纹":

do language plpgsql $$    
declare    
  zr numeric := 0.0;  -- z0r    
  zi numeric := 0.0;  -- z0i    
  tmpr numeric;    
  tmpi numeric;    
  i int;    
begin    
  <<label_x>>     
  for x in -300..300 loop  -- cr, 表示白布的x轴像素点范围是-300到300    
  <<label_y>>     
  for y in -200..200 loop  -- ci, 表示白布的y轴像素点范围是-200到200    
  <<label_i>>     
  for k in 1..200 loop  -- z的发散速度, i 表示颜色深度; 最多迭代200次, 200次就是黑色, 1次可能就接近白色.      
    tmpr := zr*zr - zi*zi + x::numeric/300.0::numeric;    
    tmpi := zr*zi + zi*zr + y::numeric/200.0::numeric;    
    zr := tmpr;    
    zi := tmpi;     
    i := k;    
    exit label_i when sqrt(zr*zr + zi*zi) > 2;    -- z的迭代次数截止于|z|>2, 因为此时z会无限发散.        
  end loop label_i ;    
  raise notice 'cr:%, ci:%, i:%', x, y, i;    
  zr := 0.0;  -- z0r    
  zi := 0.0;  -- z0i    
  end loop label_y ;    
  end loop label_x ;    
end;    
$$;

cr为x坐标,ci为y坐标,i为颜色深度:

NOTICE:  cr:-300, ci:-200, i:3    
NOTICE:  cr:-300, ci:-199, i:3    
NOTICE:  cr:-300, ci:-198, i:3    
NOTICE:  cr:-300, ci:-197, i:3    
NOTICE:  cr:-300, ci:-196, i:3    
NOTICE:  cr:-300, ci:-195, i:3    
NOTICE:  cr:-300, ci:-194, i:3    
NOTICE:  cr:-300, ci:-193, i:3    
NOTICE:  cr:-300, ci:-192, i:3    
NOTICE:  cr:-300, ci:-191, i:3    
NOTICE:  cr:-300, ci:-190, i:3    
NOTICE:  cr:-300, ci:-189, i:3    
NOTICE:  cr:-300, ci:-188, i:3    
NOTICE:  cr:-300, ci:-187, i:3    
NOTICE:  cr:-300, ci:-186, i:3    
NOTICE:  cr:-300, ci:-185, i:3    
NOTICE:  cr:-300, ci:-184, i:3    
NOTICE:  cr:-300, ci:-183, i:3    
NOTICE:  cr:-300, ci:-182, i:3    
NOTICE:  cr:-300, ci:-181, i:3    
NOTICE:  cr:-300, ci:-180, i:3    
NOTICE:  cr:-300, ci:-179, i:3    
NOTICE:  cr:-300, ci:-178, i:3    
NOTICE:  cr:-300, ci:-177, i:3    
NOTICE:  cr:-300, ci:-176, i:3    
NOTICE:  cr:-300, ci:-175, i:3    
NOTICE:  cr:-300, ci:-174, i:3    
NOTICE:  cr:-300, ci:-173, i:3    
NOTICE:  cr:-300, ci:-172, i:3    
NOTICE:  cr:-300, ci:-171, i:3    
NOTICE:  cr:-300, ci:-170, i:3    
NOTICE:  cr:-300, ci:-169, i:3    
NOTICE:  cr:-300, ci:-168, i:3    
NOTICE:  cr:-300, ci:-167, i:3    
NOTICE:  cr:-300, ci:-166, i:3    
NOTICE:  cr:-300, ci:-165, i:3    
NOTICE:  cr:-300, ci:-164, i:3    
NOTICE:  cr:-300, ci:-163, i:3    
NOTICE:  cr:-300, ci:-162, i:3    
NOTICE:  cr:-300, ci:-161, i:3    
NOTICE:  cr:-300, ci:-160, i:3    
NOTICE:  cr:-300, ci:-159, i:3    
NOTICE:  cr:-300, ci:-158, i:3    
NOTICE:  cr:-300, ci:-157, i:3    
NOTICE:  cr:-300, ci:-156, i:3    
NOTICE:  cr:-300, ci:-155, i:3    
NOTICE:  cr:-300, ci:-154, i:3    
NOTICE:  cr:-300, ci:-153, i:3    
NOTICE:  cr:-300, ci:-152, i:3    
NOTICE:  cr:-300, ci:-151, i:3    
NOTICE:  cr:-300, ci:-150, i:3    
NOTICE:  cr:-300, ci:-149, i:3    
NOTICE:  cr:-300, ci:-148, i:4    
NOTICE:  cr:-300, ci:-147, i:4    
NOTICE:  cr:-300, ci:-146, i:4    
NOTICE:  cr:-300, ci:-145, i:4    
NOTICE:  cr:-300, ci:-144, i:4    
NOTICE:  cr:-300, ci:-143, i:4    
NOTICE:  cr:-300, ci:-142, i:4    
NOTICE:  cr:-300, ci:-141, i:4    
NOTICE:  cr:-300, ci:-140, i:4    
NOTICE:  cr:-300, ci:-139, i:4    
NOTICE:  cr:-300, ci:-138, i:4    
NOTICE:  cr:-300, ci:-137, i:4    
NOTICE:  cr:-300, ci:-136, i:4    
NOTICE:  cr:-300, ci:-135, i:4    
NOTICE:  cr:-300, ci:-134, i:4    
NOTICE:  cr:-300, ci:-133, i:4    
NOTICE:  cr:-300, ci:-132, i:4    
NOTICE:  cr:-300, ci:-131, i:4    
NOTICE:  cr:-300, ci:-130, i:4    
NOTICE:  cr:-300, ci:-129, i:4    
NOTICE:  cr:-300, ci:-128, i:4    
NOTICE:  cr:-300, ci:-127, i:4    
NOTICE:  cr:-300, ci:-126, i:4    
NOTICE:  cr:-300, ci:-125, i:4    
NOTICE:  cr:-300, ci:-124, i:4    
NOTICE:  cr:-300, ci:-123, i:4    
NOTICE:  cr:-300, ci:-122, i:4    
NOTICE:  cr:-300, ci:-121, i:4    
NOTICE:  cr:-300, ci:-120, i:4    
NOTICE:  cr:-300, ci:-119, i:5    
NOTICE:  cr:-300, ci:-118, i:5    
NOTICE:  cr:-300, ci:-117, i:5    
NOTICE:  cr:-300, ci:-116, i:5    
NOTICE:  cr:-300, ci:-115, i:5    
NOTICE:  cr:-300, ci:-114, i:5    
NOTICE:  cr:-300, ci:-113, i:5    
NOTICE:  cr:-300, ci:-112, i:5    
NOTICE:  cr:-300, ci:-111, i:5    
NOTICE:  cr:-300, ci:-110, i:5    
NOTICE:  cr:-300, ci:-109, i:5    
NOTICE:  cr:-300, ci:-108, i:5    
NOTICE:  cr:-300, ci:-107, i:5    
NOTICE:  cr:-300, ci:-106, i:5    
NOTICE:  cr:-300, ci:-105, i:5    
NOTICE:  cr:-300, ci:-104, i:5    
NOTICE:  cr:-300, ci:-103, i:5    
NOTICE:  cr:-300, ci:-102, i:5    
NOTICE:  cr:-300, ci:-101, i:5    
NOTICE:  cr:-300, ci:-100, i:5    
NOTICE:  cr:-300, ci:-99, i:5    
NOTICE:  cr:-300, ci:-98, i:5    
NOTICE:  cr:-300, ci:-97, i:5    
NOTICE:  cr:-300, ci:-96, i:5    
NOTICE:  cr:-300, ci:-95, i:5    
NOTICE:  cr:-300, ci:-94, i:5    
NOTICE:  cr:-300, ci:-93, i:5    
NOTICE:  cr:-300, ci:-92, i:5    
NOTICE:  cr:-300, ci:-91, i:5    
NOTICE:  cr:-300, ci:-90, i:5    
NOTICE:  cr:-300, ci:-89, i:5    
NOTICE:  cr:-300, ci:-88, i:5    
NOTICE:  cr:-300, ci:-87, i:6    
NOTICE:  cr:-300, ci:-86, i:6    
NOTICE:  cr:-300, ci:-85, i:6    
NOTICE:  cr:-300, ci:-84, i:6    
NOTICE:  cr:-300, ci:-83, i:6    
NOTICE:  cr:-300, ci:-82, i:7    
NOTICE:  cr:-300, ci:-81, i:7    
NOTICE:  cr:-300, ci:-80, i:7    
NOTICE:  cr:-300, ci:-79, i:8    
NOTICE:  cr:-300, ci:-78, i:8    
NOTICE:  cr:-300, ci:-77, i:8    
NOTICE:  cr:-300, ci:-76, i:8    
NOTICE:  cr:-300, ci:-75, i:8    
NOTICE:  cr:-300, ci:-74, i:9    
NOTICE:  cr:-300, ci:-73, i:9    
NOTICE:  cr:-300, ci:-72, i:10    
NOTICE:  cr:-300, ci:-71, i:10    
NOTICE:  cr:-300, ci:-70, i:10    
NOTICE:  cr:-300, ci:-69, i:10    
NOTICE:  cr:-300, ci:-68, i:11    
NOTICE:  cr:-300, ci:-67, i:11    
NOTICE:  cr:-300, ci:-66, i:12    
NOTICE:  cr:-300, ci:-65, i:13    
NOTICE:  cr:-300, ci:-64, i:15    
NOTICE:  cr:-300, ci:-63, i:18    
NOTICE:  cr:-300, ci:-62, i:23    
NOTICE:  cr:-300, ci:-61, i:37    
NOTICE:  cr:-300, ci:-60, i:35    
NOTICE:  cr:-300, ci:-59, i:34    
NOTICE:  cr:-300, ci:-58, i:59    
NOTICE:  cr:-300, ci:-57, i:200    
NOTICE:  cr:-300, ci:-56, i:200    
.....

类似的: 固定c (f(z)=z^2+(-0.75+0i)), 取不同的z0值作为起点, 看z的发散速度, 绘制不同z0的二维平面图, 颜色深浅对应z的发散速度, 可以得到julia集合.

参考

http://www.matrix67.com/blog/archives/4570

http://www.matrix67.com/blog/archives/6231

https://www.cnblogs.com/anderslly/archive/2008/10/10/mandelbrot-set-by-fsharp.html

http://www.matrix67.com/blog/archives/292

https://www.eefocus.com/e/500748

分形艺术网: http://www.fxysw.com/

https://zhuanlan.zhihu.com/p/450061289

相关实践学习
使用PolarDB和ECS搭建门户网站
本场景主要介绍基于PolarDB和ECS实现搭建门户网站。
阿里云数据库产品家族及特性
阿里云智能数据库产品团队一直致力于不断健全产品体系,提升产品性能,打磨产品功能,从而帮助客户实现更加极致的弹性能力、具备更强的扩展能力、并利用云设施进一步降低企业成本。以云原生+分布式为核心技术抓手,打造以自研的在线事务型(OLTP)数据库Polar DB和在线分析型(OLAP)数据库Analytic DB为代表的新一代企业级云原生数据库产品体系, 结合NoSQL数据库、数据库生态工具、云原生智能化数据库管控平台,为阿里巴巴经济体以及各个行业的企业客户和开发者提供从公共云到混合云再到私有云的完整解决方案,提供基于云基础设施进行数据从处理、到存储、再到计算与分析的一体化解决方案。本节课带你了解阿里云数据库产品家族及特性。
相关文章
|
2月前
|
存储 缓存 关系型数据库
鱼和熊掌如何兼得?一文解析RDS数据库存储架构升级
阿里云RDS率先推出新型存储类型通用云盘,提供低延迟、低成本、高持久性的用户体验。
鱼和熊掌如何兼得?一文解析RDS数据库存储架构升级
|
7天前
|
SQL 安全 关系型数据库
DMS产品常见问题之DMS的postgis时而掉线如何解决
DMS(数据管理服务,Data Management Service)是阿里云提供的一种数据库管理和维护工具,它支持数据的查询、编辑、分析及安全管控;本汇总集中了DMS产品在实际使用中用户常遇到的问题及其相应的解答,目的是为使用者提供快速参考,帮助他们有效地解决在数据管理过程中所面临的挑战。
|
11天前
|
API 数据库 开发工具
视觉智能平台常见问题之创建人脸数据库失败如何解决
视觉智能平台是利用机器学习和图像处理技术,提供图像识别、视频分析等智能视觉服务的平台;本合集针对该平台在使用中遇到的常见问题进行了收集和解答,以帮助开发者和企业用户在整合和部署视觉智能解决方案时,能够更快地定位问题并找到有效的解决策略。
20 1
|
5月前
|
关系型数据库 定位技术 分布式数据库
沉浸式学习PostgreSQL|PolarDB 18: 通过GIS轨迹相似伴随|时态分析|轨迹驻点识别等技术对拐卖、诱骗场景进行侦查
本文主要教大家怎么用好数据库, 而不是怎么运维管理数据库、怎么开发数据库内核.
839 1
|
1月前
|
SQL 运维 数据库
数据库自治专家DAS在手,安心过春节!
春节期间,企业数据库可能面临的问题高发在稳定性、性能和安全性等方面,比如:会遭遇人手不足,IT支持可能会因运维工程师休假、交通通勤而受到限制,也会存在远端接入不方便、通信不畅通、处理不及时等问题。除此之外,大规模增长的商业活动促销和交易量波动增加,交易系统有可能凌晨依然存在高峰负载压力等问题,这对数据库性能是一场不小的考验,上述种种问题,都需要更精确的、更全方位的守护系统。通过DAS智能助手,能够实现7*24小时的守护,接触DBA焦虑,保障系统平稳有序运行。
332 2
|
3月前
|
SQL 关系型数据库 分布式数据库
西电数据库实验二:XDSQL与PolarDB对比实验
实验二 XDSQL与PolarDB对比实验 姓名:贾瑞  学号:21009201172
|
5月前
|
存储 监控 关系型数据库
PolarDB 常用情景
PolarDB 常用情景
144 0
|
6月前
|
关系型数据库 分布式数据库 数据库
沉浸式学习PostgreSQL|PolarDB 10: 社交、刑侦等业务, 关系图谱搜索
业务场景1 介绍: 社交、刑侦等业务, 关系图谱搜索 - 营销、分销、流量变现、分佣、引爆流行、裂变式传播、家谱、选课、社交、人才库、刑侦、农产品溯源、药品溯源 图式搜索是PolarDB | PostgreSQL在(包括流计算、全文检索、图式搜索、K-V存储、图像搜索、指纹搜索、空间数据、时序数据、推荐等)诸多特性中的一个。 采用CTE语法,可以很方便的实现图式搜索(N度搜索、最短路径、点、边属性等)。 其中图式搜索中的:层级深度,是否循环,路径,都是可表述的。
196 0
沉浸式学习PostgreSQL|PolarDB 10: 社交、刑侦等业务, 关系图谱搜索
|
数据库
阿里云瑶池数据库2022成绩单,有哪些高光时刻?
感谢您的支持与信任!我们将继续韧性向前,携手并进,迎接新的一年。
阿里云瑶池数据库2022成绩单,有哪些高光时刻?
|
运维 安全 数据库
明日14点开播!DAS的数据库自治技术解读
《数据库风向标》是一档聚焦数据库新趋势与新技术的视频栏目,节目每期会请到几位资深技术大咖,与大家共话数据库热点话题。
明日14点开播!DAS的数据库自治技术解读

热门文章

最新文章

相关产品

  • 云原生数据库 PolarDB