用Maxima画出一些有趣的图
Maxima可以画出Chaos、Duffing 、Fern、Lorenz、Rossler 、Portraits 、Mandelbrot、Staircase、Triangles等有趣的图...
Chaos
orbits(x^2+a, 0, 100, 400, [a,-1,-1.53], [x,-1.6,-0.8],
[nticks, 400], [style,dots], WEB_IMAGE);
Duffing √
a:1.0;
b:0.3;
c:0;
f:10.0;
eqs:[y,-(a*x^3+c*x+b*y)+f*cos(t)];
t_range:[t, 0, 40, 0.01];
inits:[0.0, 0.1];
sol:rk(eqs,[x,y],inits, t_range)$
plot2d([discrete, makelist([sol[i][2], sol[i][3]], i, 1, length(sol))], WEB_IMAGE);
Fern
a1: matrix([0.85,0.04],[-0.04,0.85])$
a2: matrix([0.2,-0.26],[0.23,0.22])$
a3: matrix([-0.15,0.28],[0.26,0.24])$
a4: matrix([0,0],[0,0.16])$
p1: [0,1.6]$
p2: [0,1.6]$
p3: [0,0.44]$
p4: [0,0]$
w: [85,92,99,100]$
ifs(w, [a1,a2,a3,a4], [p1,p2,p3,p4], [5,0], 50000, [style,dots], WEB_IMAGE);
Lorenz √
d:10;
r:28;
b:8/3;
eqs:[d*(y-x),r*x-y-x*z, x*y-b*z];
t_range:[t, 0, 40, 0.01];
inits:[0.0, 0.1,0.0];
sol:rk(eqs,[x,y,z],inits, t_range)$
plot2d([discrete, makelist([sol[i][2], sol[i][3]], i, 1, length(sol))], WEB_IMAGE);
plot2d([discrete, makelist([sol[i][3], sol[i][4]], i, 1, length(sol))], WEB_IMAGE);
plot2d([discrete, makelist([sol[i][2], sol[i][4]], i, 1, length(sol))], WEB_IMAGE);
Rossler √
a:0.2;
b:0.2;
c:5.7;
eqs:[-y-z, x+a*y, b+x*z-c*z];
t_range:[t, 0, 200, 0.01];
inits:[0.0, 0.01, 0.01];
sol:rk(eqs,[x,y,z],inits, t_range)$
plot2d([discrete, makelist([sol[i][2], sol[i][3]], i, 1, length(sol))], WEB_IMAGE);
plot2d([discrete, makelist([sol[i][2], sol[i][4]], i, 1, length(sol))], WEB_IMAGE);
plot2d([discrete, makelist([sol[i][3], sol[i][4]], i, 1, length(sol))], WEB_IMAGE);
Portraits √
numPortraits:8;
eqs:[sin(x)+y, sin(y)-x];
t_range:[t, 0, 10, 0.01];
for k:1 thru numPortraits do (
inits:[0.0, k*0.1],
sol[k]:rk(eqs,[x,y],inits, t_range)
)$
sol_list:makelist(
[discrete, makelist([sol[k][i][2], sol[k][i][3]], i, 1, length(sol[k]))],
k, 1, numPortraits)$
plot2d(sol_list, WEB_IMAGE)$
simple Mandelbrot
mandel(reMin, reMax, imMin, imMax, resolution, iterations):=block(
[reZ, imZ, absZ, i, pixels, reTemp],
pixels:[[0,0]],
for x:reMin step resolution thru reMax do [
for y:imMin step resolution thru imMax do [
reZ:0,
imZ:0,
absZ:0,
for i:1 step 1 thru iterations while absZ < 4 do [
reTemp:reZ^2 - imZ^2 + x,
imZ:2*reZ*imZ + y,
reZ:reTemp,
absZ:reZ^2+imZ^2
],
if absZ < 4 then pixels:append(pixels,[[x,y]])
]
],
return(pixels)
)$
plot2d([discrete, mandel(-2,0.5,-1,1, 0.02, 16)],[style,
[points,1,1,1]], WEB_IMAGE);
Staircase
load("dynamics")$
staircase(cos(y^2), 1, 11, [y, 0, 1.2], WEB_IMAGE);
Triangles
chaosgame([[0, 0], [1, 0], [0.5, sqrt(3)/2]], [0.1, 0.1], 1/2,
30000, [style, dots], WEB_IMAGE);
Maple中可以试着画出:
Duffing
Lorenz
Rossler
Portraits
用Maxima画出一些有趣的图相关推荐
- df python 增加数据_美国确诊超100万!教你用Python画出全球疫情动态图(附数据下载)...
CDA数据分析师 出品 [导语]:今天我们教你用Python绘制全球疫情动态图,技术部分请看第二部分. 获取数据:扫描下方公众号回复关键字"全球疫情" 目前,全球新冠疫情还十 ...
- 生信漫谈如何画出美美的SeqLogo图
画出美美的SeqLogo图,下面介绍三种不同的方法,适用不同的人群的需求. 1.直接上在线绘制SeqLogo图的网站WebLogo 3: https://weblogo.threeplusone.co ...
- matlab绘制两个函数,用matlab程序 在同一图中绘制两个函数,这两个函数分别是: , ,,matlab怎样同时画出两个函数图...
导航:网站首页 > 用matlab程序 在同一图中绘制两个函数,这两个函数分别是: , ,,matlab怎样同时画出两个函数图 用matlab程序 在同一图中绘制两个函数,这两个函数分别是: , ...
- java矢量图_你见过js画出这样的矢量图以及动画特效吗?
js 画出几十种矢量图 raphaeljs 超炫丽的动画效果 3.动态曲线图 ---- js 画图 raphaeljs demo http://raphaeljs.com/analytics.ht ...
- 根据伪代码画出流程图和盒图以及根据流程图判断是否为结构化流程图,并且为其设计一个等价结构化程序。
一.首先附上作业图: 二.解决实例 1.画出程序流程图和盒图. 流程图 ...
- 美国确诊超100万!教你用Python画出全球疫情动态图
目前,全球新冠疫情还十分严峻.截止4月29日,疫情突破多个标志性的统计节点.其中,全球累计确诊近296万例,累计死亡超20万例:美国累计确诊超100万例,西班牙和意大利累计确诊均超过20万例. 01 ...
- 推荐:一个画出高大上的机器学习图的神器(论文必备)
本人在写论文的时候,很多图片是用matplotlib和seaborn画的,但是,我还有一个神器,Scikit-plot,通过这个神器,画出了更加高大上的机器学习图,本文对Scikit-plot做下简单 ...
- 同一界面画出多个QQ图(R语言)
题目:一名研究者用光子吸收法测量了妇女骨骼中无机物含量,对三根骨头主侧和非主侧记录了测量值,数据框"T1bones.txt"中的第2至第7列记录了相应数据.对各个变量做qq图,在同 ...
- Python零基础学习 | 分享一个使用Python画出太极阴阳八卦图的代码
大家好,我是王某人. 一.写在前面 前几天在Python技术交流群有个小伙伴分享一个使用Python画出太极阴阳八卦图的代码,这里拿出来给大家分享下,一起学习下. 不过这里他的代码是有点问题的,不过不 ...
- 不需要懂得编程,但却可以使用ggplot2画出论文级别的图?
你有没有遇到过这样的烦恼,你需要画一些论文级别的图,并且你知道R中的ggplot2是一个很好的选择,可以画出符合你要求的图.但是由于你不熟悉ggplot2的使用,你需要上网倒弄一番,了解与你图相关的代 ...
最新文章
- HBase概念学习(八)开发一个类twitter系统之表设计
- [全]php-redis函数使用
- 数据存储方式_视频监控系统的数据存储方式的概念及应用
- java flatmap_Java 8 Steam API map和flatMap方法使用详解
- 杂七杂八DataTable获取Table某一列成List转换Table成List
- securecrt上传下载文件命令
- 【自我救赎--牛客网Top101 4天刷题计划】 第一天 热身运动
- 使用HTML+CSS设计个人简历
- 最全最新cpu显卡天梯图_2019.12月CPU和显卡性能天梯图
- 服务器cpu哪个型号可以超频,CPU超频史上的14个神器:你用过几个?
- React lazyLoad懒加载
- PNP问题学习笔记1
- asp.net paypal信用卡支付功能
- 介绍一下国家葡萄产业体系,列出全世界最重要的葡萄育种单位。
- Firefox浏览器基本使用
- springcloud + oauth2
- 青海湖 游走遇见一路风景
- 嵌入式远程机房数采监控报警系统
- 收藏这些网站,做设计和PPT不再为找素材发愁
- POJ - 3179 Corral the Cows(二分,离散化,前缀和)
热门文章
- c 语言 登陆窗口界面,c/c++语言实现登陆界面
- 2022-02-09 Android setLayoutParams动态修改RelativeLayout控件layout_width、layout_marginTop等。实际应用实例:控件随手指移动。
- 老鸟必备 | 如何画出优秀的架构图
- 国外广告联盟前期需要准备的事情
- 苹果换原装电池_苹果手机换电池客户必看!苹果原装电池科普鉴别!
- Android电子白板
- 30天自制操作系统——第三天
- 线性代数的那些事(二)行列式与逆
- CSS布局:多种方案实现固定页脚(sticky footer)
- 利用批处理文件快速设置IP地址