ExcelHome技术论坛

 找回密码
 免费注册

QQ登录

只需一步,快速开始

快捷登录

搜索
EH技术汇-专业的职场技能充电站 妙哉!函数段子手趣味讲函数 Excel服务器-会Excel,做管理系统 效率神器,一键搞定繁琐工作
HR薪酬管理数字化实战 Excel 2021函数公式学习大典 Excel数据透视表实战秘技 打造核心竞争力的职场宝典
让更多数据处理,一键完成 数据工作者的案头书 免费直播课集锦 ExcelHome出品 - VBA代码宝免费下载
用ChatGPT与VBA一键搞定Excel WPS表格从入门到精通 Excel VBA经典代码实践指南
楼主: lee1892

[原创] 美丽的数学 之 奇妙的分形

  [复制链接]

TA的精华主题

TA的得分主题

 楼主| 发表于 2013-12-11 11:46 | 显示全部楼层
本帖已被收录到知识树中,索引项:
本帖最后由 lee1892 于 2013-12-14 10:18 编辑

七、一些知识准备

可能会有人去看了顶楼的那个一片漂亮的叶子的帖子,并疑惑于这样的有规律的随机数据是如何产生的,其背后又有着什么样的数学规律。

想要了解个大概的话,可以看看下面的一些内容,或则只是关心好看的效果的话,这一段完全可以跳过~

1、先复习一下复数、复平面

为解决 x^2 = -1,实数显然无解,人们引入了虚数的概念,x = ± 1 * i,这里的 i 表示虚数,于是有了复数的概念 z = x + y*i,其中 x 为实部、y 为虚部。

所有复数的集合就可以构成复平面,与笛卡儿坐标系类似的,y轴表示虚部,x轴表示实部,任何一个复数都可以在复平面上找到一个对应的点。

2、度量空间 Metric Space

度量空间是一个集合,在其中可以定义在这个集合内的元素之间的 “距离”(称为 度量)的概念。最为直观的是三维欧几里得空间,欧几里得度量 定义了两个点之间的距离为其连线的长度。

可以把一个度量空间表达为一个二元组 (M, d),这里 M 是集合,而 d 是在 M 上的 度量,即一个函数
d : M × M → R

3、压缩映射 Contraction Mapping

度量空间 (M,d) 上的压缩映射,是一个从 M 到它本身的函数 f,存在某个实数 0<k<1,使得对于所有 M 内的 xy,都有:
d(f(x),f(y)) ≤ kd(x,y)

满足以上条件的最小的 k 称为 f 的利普希茨常数(Lipschitz Constant)。

一个压缩映射最多有一个不动点。另外,巴拿赫不动点定理(Banach Fixed Point Theorem)说明,非空的完备度量空间上的每一个压缩映射都有唯一的不动点,且对于 M 内的任何 x,迭代函数序列 x, f(x), f(f(x)), f(f(f(x))), ... 收敛于不动点。

4、仿射变换与变换矩阵 Affine Transformation & Transformation Matrix

在几何中,一个向量空间进行一次线形变换并接上一个平移,变换为另一个向量空间。

在线形代数中,线形变换嫩够用矩阵表示。如果 T 是一个把 R1 映射到 R2 的线形变换,且 x 是一个具有 n 个元素的列向量,那么
T(x) = Ax
m×n 的矩阵 A,称为 T 的变换矩阵。

更多变换矩阵的知识,可参考 中文WiKi:http://zh.wikipedia.org/wiki/%E5%8F%98%E6%8D%A2%E7%9F%A9%E9%98%B5

5、复合函数 Composition Function

在数学领域,两个函数的复合函数指将第一个函数作用于参数,然后再将第二个函数作用于所得结果的函数。具体来说,给定两个函数 f: X → Yg: Y → Z,其中 f 的陪域等于 g 的定义域(称为 fg 可复合),则其复合函数记为 g ○ f,以 X 为定义域,Z 为陪域,并将任意 x∈X 映射为 g(f(x))

函数的复合满足结合律:若 fg可复合,gh 可复合,则有:
h ○ (g ○ f) = (h ○ g) ○ f

============================================
顶楼的那个叶子的图案被称为 Barnsley's Fern,使用了 4 个仿射变换,一个变换的矩阵公式为:
9a999429d9c707ff98a6afc2bd1fdda3.png
如果你对矩阵计算熟悉的话就会知道上面的这个公式对应了两个等式:
x = a * x' + b * y' + e
y = c * x' + d * y' + f
上式中,x'、y’ 为上一次迭代计算的结果。

Barnsley's Fern 的参数表如下:
wabcdefp
f10000.16000.01
f20.850.04-0.040.8501.60.85
f30.2-0.260.230.2201.60.07
f4-0.150.280.260.2400.440.07
其中,p 表示的是可能性,即每次迭代计算中,采用 4 种变换中的某一种的可能性。

对应的矩阵变换如下:
97d89acf29651f72c8485dfb56da4bd6.png

489dfa63d783cf2f13a3257c408191c8.png

c8470d468ea7f87a6079f36093003884.png

40bb54270f27210036ac6bb6c7fc2a1a.png


如果愿意的话,可以计算一下上面四个变换的不动点,具体办法可参考 巴拿赫不动点定理

================
以上内容均来自中文或英文 WiKi。

点评

这个完全可以公开,因为不见得就能轻易拿走!!!呵呵  发表于 2014-9-23 20:02

评分

1

查看全部评分

TA的精华主题

TA的得分主题

发表于 2013-12-11 12:48 | 显示全部楼层

数学之美||数学之美||数学之美||数学之美||数学之美

TA的精华主题

TA的得分主题

发表于 2013-12-11 12:50 | 显示全部楼层
[广告] Excel易用宝 - 提升Excel的操作效率 · Excel / WPS表格插件       ★免费下载 ★       ★ 使用帮助

数学之美||数学之美||数学之美||数学之美||数学之美

点评

您这是要表达激动之情嘛{:soso_e128:}  发表于 2013-12-11 13:44

TA的精华主题

TA的得分主题

 楼主| 发表于 2013-12-14 12:26 | 显示全部楼层
本帖最后由 lee1892 于 2013-12-14 13:58 编辑

八、德拉姆曲线 De Rham Curves 0到1的舞蹈

18楼有更新噢~

本贴附件: 奇妙的分形 之 DeRham曲线_By Lee1892.rar (174.98 KB, 下载次数: 368)

考察一个度量空间 (M,d),和它的一个 M 上的压缩映射:
d0: M → M
d1: M → M

由巴纳赫不动点定理,它们存在两个不动点,记为 p0p1。设实数 x ∈ [0,1],并且其二进制扩展表达为:
x = SUM (bk / 2^k) <k=1~∞, bk = 1 或 0>

这时,有:
cx: M → M
其度量函数为:
cx = db1 ○ db2 ○ ... ○ dbk ○ ...

可以发现每一个 cx 都将 d0d1 的普遍吸引域映射到 M 上的一个点,记为 px。由一个实参 x 获得的所有 px 的集合,被称为 德拉姆曲线。

德拉姆曲线连续的条件:
d0(p1)=d1(p0)

嗯,这一段确实太抽象了点儿。接下来的几个实例可以帮助我们稍微形象点儿的理解它。

1、Césaro 曲线
其两个不动点为 p0 = 0p1 = 1,不同的 Césaro 曲线是由一个复数 aa = α + β i)来定义的,|a| < 1  且 |1-a| < 1 (复数的绝对值是其在复平面上代表的点距离原点的距离,即 |z| = sqr(x^2 + y^2) )。

其压缩映射的两个度量函数被定义为复平面上的两个复函数:
d0(z) = az'
d1(z) = a + (1 - a)z'
其仿射变换为:
cesaro curve.PNG

a = 0.3 + 0.3 i
cesaro dwg.PNG
a = 0.5 + 0.5 i 就是 莱维C曲线
levy c dwg.PNG

2、Koch-Peano 曲线
类似的,将上述函数中的变量换为其共轭复数(复数的共轭:z = x + y iz = x - y i 共轭,没有上划线用下划线表示共轭),就变成了 Koch-Peano 曲线了。
d0(z) = az'
d1(z) = a + (1 - a)z'
其仿射变换为:
koch-peano curve.PNG


a = 0.5 + √3 / 6 i 为 Koch 雪花曲线
koch dwg.PNG
a = 0.5 + 0.5 i 为 Peano 曲线 (不好看,不贴了)
a = 0.6 + 0.37 i
koch .6 .37.PNG
a = 0.6 + 0.45 i
koch .6 .45.PNG

3、德拉姆 曲线 的一般化
按上述的仿射变换公式,显然一般化的形式如下:
generalized.PNG
式中包括了全部的向量的线性变化:旋转、偏移、放大缩小等。

由前述的设定,对于 d0 取不动点 p0 = 0 则有 a = d = 0,对于 d1 取不动点 p1 = 1 则有 j = 1 - hm = -l。而由 德拉姆曲线连续的条件 d0(p1)=d1(p0),可得 h = bm = e。将这些变量换以希腊字母,就可以获得 德拉姆曲线的公式:
deRham curve.PNG

deRham dwg.PNG

补充内容 (2014-10-25 00:48):
连续的德拉姆曲线见:http://club.excelhome.net/forum. ... 988&pid=7906838

评分

1

查看全部评分

TA的精华主题

TA的得分主题

发表于 2013-12-14 20:03 | 显示全部楼层

TA的精华主题

TA的得分主题

发表于 2013-12-21 21:48 | 显示全部楼层
[广告] VBA代码宝 - VBA编程加强工具 · VBA代码随查随用  · 内置多项VBA编程加强工具       ★ 免费下载 ★      ★使用手册
网上有本流行书《数学之美》,通读之后,可能与我的行业并无相关,至少我感觉不到美。
但是楼主的这一篇大作,让我由衷的敬仰。
这要是作为大学数学系的入门基础一课,将会是多么好的诱导与兴趣提升。至少实实在在见到了,那些高深公式的实实在在的用处

点评

精彩的建议,教育部应当采纳!!!  发表于 2014-9-23 20:03

评分

1

查看全部评分

TA的精华主题

TA的得分主题

发表于 2014-1-16 11:26 | 显示全部楼层
[广告] VBA代码宝 - VBA编程加强工具 · VBA代码随查随用  · 内置多项VBA编程加强工具       ★ 免费下载 ★      ★使用手册
满满的都是艺术...后悔没学好数学

评分

1

查看全部评分

TA的精华主题

TA的得分主题

发表于 2014-6-10 13:34 | 显示全部楼层

TA的精华主题

TA的得分主题

发表于 2014-6-11 22:08 | 显示全部楼层
[广告] Excel易用宝 - 提升Excel的操作效率 · Excel / WPS表格插件       ★免费下载 ★       ★ 使用帮助
五体投地!...后悔没学好数学,看不懂...

TA的精华主题

TA的得分主题

发表于 2014-6-11 22:34 | 显示全部楼层
膜拜!一直以为数学学的不赖,然而从未将数学与艺术勾连起来,今日是开了眼界了,谢谢楼主。
您需要登录后才可以回帖 登录 | 免费注册

本版积分规则

手机版|关于我们|联系我们|ExcelHome

GMT+8, 2024-11-22 05:32 , Processed in 0.042588 second(s), 8 queries , Gzip On, MemCache On.

Powered by Discuz! X3.4

© 1999-2023 Wooffice Inc.

沪公网安备 31011702000001号 沪ICP备11019229号-2

本论坛言论纯属发表者个人意见,任何违反国家相关法律的言论,本站将协助国家相关部门追究发言者责任!     本站特聘法律顾问:李志群律师

快速回复 返回顶部 返回列表