【SEU程序设计课笔记】 17 - 2020/12/10 - Write C++ programs - 分形绘图

Teddy van Jerry

于 2020-12-10 15:29:03 发布

阅读量300

点赞数

分类专栏：【SEU程序设计课（C/C++）】文章标签： c++ 分形绘图 EasyX

本文链接：https://blog.csdn.net/weixin_50012998/article/details/110947338

版权

【SEU程序设计课（C/C++）】专栏收录该内容

39 篇文章 2 订阅

订阅专栏

This is Fractal Designer 0.1.

题目

绘制一个分形图形，以之前的复数类为基础内容。

我之前的复数（Complex）类详见我的博客【C++ 程序】 TVJ Complex Calculator (v 2.1) 复数计算器中 Complex.h 和 Complex.cpp。
这个复数类功能比较强大，定义了多种运算，可以直接使用，包括+（加法），^（乘方），.modulus()（取模）。
最新版的复数类（Complex）见我的博客【C++ 程序】复数类（Complex）。

平面上的每一个点，都对应了一个复数坐标，下面判断该点是否发散：
原来的点设为 $z_0$ ，下面对于每一个 $z$ ,进行如下操作100次：
$z_{n+1}=(z_n)^2+z_0\text{,}$
如果 $|z_{100}|\leqslant10$ ，则认为不发散，否则认为发散。

界面范围： $x$ : $- 2.3$ ~ $0.9$ ； $y$ : $- 1.2$ ~ $1.2$ 。

初级版

代码

#include <iostream>
#include <graphics.h> // the EasyX graphics library
#include "Complex.h"
using namespace std;

IMAGE blue, black;

int main(int argc, char** argv)
{
	initgraph(801, 601);
	for (double i = 800; i >= 0; i--)
	{
		for (double j = 0; j <= 600; j++)
		{
			Complex this_point((i - 575) / 250, (j - 300) / 250);
			Complex z0 = this_point;
			bool test = true;
			for (int k = 0; k != 100; k++)
			{
				if (this_point.modulus() > 10)
				{
					test = false;
					break;
				}
				this_point = (this_point ^ 2) + z0; // the bracket can not be omitted due to the priority
			}
			if (!test)
			{
				loadimage(&black, _T("Black.jpg"), 1, 1);
				putimage(i, j, &black);
			}
			else
			{
				loadimage(&blue, _T("Blue.jpg"), 1, 1);
				putimage(i, j, &blue);
			}
		}
	}
	system("pause");
	return 0;
}

// This is Fractal Designer 0.1 by Teddy van Jerry.