zoukankan      html  css  js  c++  java
  • 曼德勃罗(Mandelbrot)集合与其编程实现

    一、从科赫雪花谈起

    设想一个边长为1的等边三角形(例如以下图所看到的)。取每边中间的三分之中的一个,接上去一个形状全然类似的但边长为其三分之中的一个的三角形,结果是一个六角形。如今取六角形的每个边做相同的变换,即在中间三分之中的一个接上更小的三角形,以此反复,直至无穷。

    外界的变得原来越细微曲折,形状接近理想化的雪花。

    这个图形的名字叫柯赫雪花。

    1904年,瑞典数学家柯赫首次描写叙述了一种后来以其名字命名的曲线——柯赫曲线,而柯赫雪花就是由柯赫曲线所围成的一个封闭的图形。

    柯赫雪花有非常多奇异的性质。

    比如。柯赫雪花的周长趋于无穷。而柯赫雪花的面积却是有限的!此外,柯赫雪花也体现了一种(递归的)自类似性,即当把图形的局部放大。它会呈现出一种与总体(或总体的部分)之间的惊人形似性。


    二、分形

    虽然类似于柯赫雪花这种自类似图形早已经出现。可是真正把这些问题系统地发展成一门学问,则是在二十世纪下半页。数学家本华·曼德勃罗(Benoit B. Mandelbrot)把这些部分与总体以某种方式类似的形体称为分形(fractal)。1975年,他创立了分形几何学(Fractal Geometry)。

    在此基础上。形成了研究分形性质及其应用的科学,称为分形理论(Fractal Theory)。因此,曼德勃罗也被称为“分形学之父”。如今。分形理论已经发展成为一门十分风靡和活跃的新理论、新学科,特别是把分形理论和混沌理论结合之后,更是衍生出一大片广阔的研究天地。分形的世界与我们寻常所研究的几何学中非常多直观的常识之间具有巨大的冲突。比如,寻常我们说一个几何图形的维度,那么这个维度一般都应该是整数。比如平面上一个矩形就是二维的,空间中一个球体就是三维的。可是在分析几何中。几何图形的维度都是不是整数而是分数。这个维度又称为豪斯多夫维或豪斯多夫-贝塞科维奇维(Hausdorff-Becikovich Dimesion),它是由数学家豪斯多夫于1918年引入的。

    通过豪斯多夫维能够给一个随意复杂的点集合比方分形赋予一个维度。


    三、曼德勃罗集合(Mandelbrot Set)

    曼德勃罗集合(Mandelbrot Set)或曼德勃罗复数集合。是一种在复平面上组成分形的点的集合,因由曼德勃罗提出而得名。曼德博集合能够使复二次多项式 进行迭代来获得

    当中,c是一个复參数。对于每个c。从 z = 0 開始对fc(z)进行迭代。序列 的值或者延伸到无限大,或者仅仅停留在有限半径的圆盘内(这与不同的參数c有关)。曼德布洛特集合就是使以上序列不延伸至无限大的全部c点的集合。

    最后。我们给出一个利用C语言生成Mandelbrot集合并绘制图形的程序(该程序来自文献【1】):

    #include <stdio.h>
    #include <stdlib.h>
    #include <complex.h>
    
    #define width_size      800
    #define height_size     600
    #define Maxval          255
    
    static const float orig_x = width_size * 2/3;
    static const float orig_y = height_size * 1/2;
    
    typedef struct _pixel {
        unsigned char r;
        unsigned char g;
        unsigned char b;
    } pixel;
    
    static const pixel dim_gray = { 105, 105, 105 };
    
    static unsigned char iteration(int x, int y)
    {
        const int limit = Maxval + 1;
        int i;
        complex c = ((x - orig_x) / (width_size / 3)) +
        ((orig_y - y) / (height_size / 2)) * I;
        complex z = 0;
        
        for (i = 0; i < limit; i++) {
            /* basic formula */
            z = z * z + c;
            if (creal(z) > 2 || cimag(z) > 2)
                break;
        }
        return (unsigned char) (i == limit ? 0 : i);
    }
    
    
    int main()
    {
        FILE *f = fopen("mandelbrot.ppm", "w+");
        
        /* PPM header */
        fprintf(f,
                "P6
    "      /* PPM magic number */
                "#Mandelbrot Set
    "
                "%d "       /* width, in ASCII decimal */
                "%d
    "      /* height, in ASCII decimal */
                "%d
    ",     /* maximum color value, in ASCII decimal */
                width_size, height_size, Maxval);
        
        /* Write every pixel generated by Mandelbrot Set */
        for (int i = 0; i < height_size; i++) {
            for (int j = 0; j < width_size; j++) {
                unsigned char iter = iteration(j, i);
                if (iter) {
                    pixel p = {
                        .r = iter,
                        .g = (float) abs(j - orig_x) / width_size * Maxval,
                        .b = (float) abs(i - orig_y) / height_size * Maxval };
                    fwrite(&p, sizeof(pixel), 1, f);
                } else {
                    fwrite(&dim_gray, sizeof(pixel), 1, f);
                }
            }
        }
        
        fclose(f);
        return 0;
    }

    上述程序所生成的图像结果例如以下图所看到的,须要补充说明的是:该图像文件格式为ppm,在Windows下你能够使用Photoshop来查看这种类型的图像文件,在OS X系统下你能够使用免费的GIMP软件来查看它。



    參考文献:

    【1】http://blog.linux.org.tw/~jserv/archives/2011/09/_mandelbrot_set.html


  • 相关阅读:
    Dynamics CRM修改密码界面
    Dynamics CRM新加了组织后提示数据加密错误的解决方法
    ADFS修改默认访问端口
    Dynamics CRM与ADFS安装到同一台服务器后ADFS服务与Dynamics CRM沙盒服务冲突提示808端口占用问题
    Dynamics CRM各个版本的元数据浏览解决方案
    Dynamics CRM9.0安装CRM的时候提示Microsoft.Crm.Setup.Server.RegisterSandboxServiceAction操作失败
    IIS误删了默认网站,恢复方法
    grafana配置邮件告警
    Git连接Gitlab并提交内容
    docker自定义bridge
  • 原文地址:https://www.cnblogs.com/llguanli/p/7307589.html
Copyright © 2011-2022 走看看