分型图z=z*z+c的画法
python代码:
import numpy as np
from PIL import Image
from numba import jit
MAXITERS = 200
RADIUS = 100@jit
def color(z, i):v = np.log2(i+1-np.log2(np.log2(abs(z))))/5if v < 1.0:return v**4, v**2.5, velse:v = max(0, 2-v)return v, v**1.5, v**3@jit
def iterate(c):z = 0jfor i in range(MAXITERS):if z.real*z.real +z.imag*z.imag > RADIUS:return color(z, i)z = z*z + creturn 0, 0, 0def main(xmin, xmax, ymin, ymax, width, height):x = np.linspace(xmin, xmax, width)y = np.linspace(ymax, ymin, height)z = x[None, :] +y[:, None]*1jred, green, blue = np.asarray(np.frompyfunc(iterate, 1, 3)(z)).astype(np.float)img = np.dstack((red, green, blue))Image.fromarray(np.uint8(img*255)).save('mandelbrot.png')main(-2.1, 0.8, -1.16, 1.16, 1200, 960)
效果图:
分型图z=z*z+c的画法相关推荐
- Scratch编程与数学之神奇的分型图-勾股树!
注:这个程序脚本需要用到积木模块:动作.画笔.控制.事件.自制积木(适合有一定编程经验,且学过勾股定理,想编程进阶的孩子) 勾股树,又可以称为毕达哥拉斯树.它是古希腊数学家毕达哥拉斯根据勾股定理画出的 ...
- hdu4282 x^z+y^z+x*y*z=k 解的个数
题意: x^z + y^z + x*y*z = k; (x < y ,z > 1),给你一个k问有多少组解. 思路: 暴力枚举z,y,然后二分查找x.注意一点最好用 ...
- 将字符串中的字母全部替换成字母的下一个字母,要是最后一位是z或Z则替换为a或A
第九题 功能描述:将字符串中的字母全部替换成字母的下一个字母, 要是最后一位是z或Z则替换为a或A. 输入:aBxyZ 输出:bCyzA /*** 功能描述:将字符串中的字母全部替换成字母的下一个字母 ...
- k线顶分型 python_顶底分型-(K线分类及顶底分型的一种数学原理 源码 贴图)...
好股票软件下载网(www.goodgupiao.com)提示:您正在下载的是:顶底分型-(K线分类及顶底分型的一种数学原理 源码 贴图) 参考缠论,研究了很多天终于将顶底分型进行了具体的数学量化,涵盖 ...
- C语言 一条语句max = x>y? (x>z?x:z) : (y>z?y:z)实现求 x, y, z 最大值
从键盘输入三个整数,按由小到大的顺序输出. #include "stdio.h" main() {int i,j,k,max;scanf("%d%d%d",&a ...
- ML:通过数据预处理(分布图/箱型图/模型寻找异常值/热图/散点图/回归关系/修正分布正态化/QQ分位图/构造交叉特征/平均数编码)利用十种算法模型调优实现工业蒸汽量回归预测(交叉训练/模型融合)之详
ML之LightGBM:通过数据预处理(分布图/箱型图/模型寻找异常值/热图/散点图/回归关系/修正分布正态化/QQ分位图/构造交叉特征/平均数编码)利用十种算法模型调优实现工业蒸汽量回归预测(交叉训 ...
- PHP正则表达式之\A、\Z、 \z、^、$
文章目录 \A.\Z. \z 断言与 ^.$ 之间的区别 单行模式和多行模式的差异 换行符匹配的差异 示例代码及解读 文中提到的模式修饰符 D (PCRE_DOLLAR_ENDONLY) m (PCR ...
- 主图顶底分型自动提示
V00:= HIGH<REF(HIGH,1) AND LOW<REF(LOW,1);{//高点低于昨日,低点低于昨日} V01:= HIGH<REF(HIGH,1) AND LOW& ...
- matplotlib柱状图、面积图、直方图、散点图、极坐标图、箱型图
一.柱状图 1.通过obj.plot() 柱状图用bar表示,可通过obj.plot(kind='bar')或者obj.plot.bar()生成:在柱状图中添加参数stacked=True,会形成堆叠 ...
最新文章
- sublime怎么编译java_文本编辑利器:Sublime Text3 编译运行Java程序
- MFC ComboBox
- 日报 18/07/22 您的设计模式!终于有时间开写~ ~
- cocos2d-x初探学习笔记(7)--CCProgressTimer
- RFileWriteStream 写入汉字
- 【项目管理】变更管理与过程改进
- Problem F: 时间类的构造和输出
- 无人机集群定位与编队方式研究总结(最详细的鸟群模型分析)
- 让两个Div并排显示
- Python数据分析——matplotlib折线图和散点图
- 腾讯 2016届实习生招聘笔试
- 加拿大首个以女性名字命名的工程学院:商界领袖Gina Cody向蒙特利尔康考迪亚大学捐赠1500万加元,创造历史
- mysql for mac 重置密码_Mac版,mysql重置密码
- 关于H5中的Canvas API的探索
- 洪水填充算法_洪水填充算法分析
- Kotlin协变和逆变
- 你可曾了解过,数据分析师与商业数据分析师?
- H.264视频中SPS/PPS缺失时实现补帧
- 政务使用区块链技术,网络安全风险不容小觑
- 学UI设计有必要考Adobe国际认证证书吗?
热门文章
- 电子科技大学计算机调剂要求,2016电子科技大学考研调剂信息发布(3.16更新)
- 路由器上网行为管理功能浅谈 -- 路由端
- openstack octavia部署 ussuri U版
- android拍照模糊,解决Android拍照并显示在ImageView中变模糊
- 2019最赞的大屏可视化,一定是这款双十一的酷屏
- 开源示波器ADALM2000介绍、原理图学习
- 2019/11/11 还是没能脱单噻
- 小程序毕业设计 汽车维修保养预约小程序毕设作品开题报告参考
- event_dlist
- .NET-1.理论知识(历史介绍和了解)