运行这段代码到出图片,你的电脑用了多久? - 程序语言 - 第3页小木虫论坛-学术科研互动平台

ส็็็็็็็็็็็็็็็็็็็็

21楼2014-09-17 22:19:49

easezh

木虫 (著名写手)

应助: 24 (小学生)
金币: 10265.8
红花: 12
帖子: 1542
在线: 438.9小时
虫号: 1197825
注册: 2011-01-29
专业: 模式识别

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

25.9740456 s
------
i7-3770 3.4G
32G RAM
64 WIN7
运行时，操作系统下还有若干其他程序在运行。
因此我觉得可以更快，但这个似乎已经够快了。

Untitled-1.png

22楼2014-09-18 01:10:09

mshwangg

至尊木虫 (正式写手)

程序强帖: 5
应助: 206 (大学生)
金币: 10702.8
散金: 100
红花: 19
帖子: 597
在线: 195.4小时
虫号: 576702
注册: 2008-06-21
专业: 物理学I

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

21楼: Originally posted by cooooldog at 2014-09-17 22:19:49
这个配置,只有这个性能...

真这样啊？那看来需要弄个像样点的台式机跑计算了
本本还是做做doc，ppt啥的好了

23楼2014-09-18 09:10:09

龙骧楼主

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 538.8
帖子: 36
在线: 6.3小时
虫号: 1645510
注册: 2012-02-26
性别: GG
专业: 化工冶金

我的46秒多

24楼2014-09-30 18:25:26

戴钢盔的猪头

木虫 (知名作家)

应助: 200 (高中生)
金币: 5083.9
散金: 1310
红花: 18
帖子: 6001
在线: 319.9小时
虫号: 2854040
注册: 2013-12-07
性别: GG
专业: 波谱分析与成像分析

引用回帖:

2楼: Originally posted by mshwangg at 2014-09-15 23:08:18
牛人啊，树叶挺漂亮
还可以画别的吗？分享下

教材上有的

[ 发自手机版 http://muchong.com/3g ]

25楼2014-09-30 23:39:23

fzhqhd

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 2401.3
帖子: 98
在线: 50.9小时
虫号: 1151273
注册: 2010-11-20
性别: GG

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

上传图片:

Mathematica.jpg

26楼2014-10-01 17:49:35

chyanog

金虫 (小有名气)

程序强帖: 1
应助: 12 (小学生)
金币: 808.1
红花: 5
帖子: 85
在线: 52.8小时
虫号: 1540955
注册: 2011-12-17
性别: GG
专业: 计算数学与科学工程计算

引用回帖:

25楼: Originally posted by 戴钢盔的猪头 at 2014-09-30 23:39:23
教材上有的
...

是什么教材啊？

27楼2014-10-06 22:58:37

chyanog

金虫 (小有名气)

程序强帖: 1
应助: 12 (小学生)
金币: 808.1
红花: 5
帖子: 85
在线: 52.8小时
虫号: 1540955
注册: 2011-12-17
性别: GG
专业: 计算数学与科学工程计算

★
小木虫: 金币+0.5, 给个红包，谢谢回帖
jjdg: 程序强帖+1, 强帖 2014-10-18 00:23:51

说实话，代码写得有点糟，不符合Mathematica style，而且速度太慢。
同样的算法一样的效果其实只需不到半秒甚至不到0.1s

CODE:

Clear["`*"];

t0 = AbsoluteTime[];

ifs[prob_, A_, init_, maxIter_] := 

  FoldList[#2.{#[[1]], #[[2]], 1} &, init, 

   RandomChoice[prob -> A, maxIter]];

init = {0., 0.};

prob = N@{85, 6, 8, 1};

A = N@{{{0.83, 0.03, 0}, {-0.03, 0.86, 1.5}}, {{0.2, -0.25, 0}, {0.21,

       0.23, 1.5}}, {{-0.15, 0.27, 0}, {0.25, 0.26, 0.45}}, {{0, 0,  0}, {0, 0.17, 0}}};

iter = 10^5;

pts = ifs[prob, A, init, iter];

Graphics[{PointSize[0.002], Green, Point[pts]}]

AbsoluteTime[] - t0

赞一下(4人)

28楼2014-10-06 23:21:13

Mr__Right

专家顾问 (著名写手)

专家经验: +31
应助: 317 (大学生)
金币: 14456.3
散金: 500
红花: 54
帖子: 2716
在线: 950.6小时
虫号: 1972612
注册: 2012-09-04
性别: GG
专业: 应用数学方法
管辖: 程序语言

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

28楼: Originally posted by chyanog at 2014-10-06 23:21:13
说实话，代码写得有点糟，不符合Mathematica style，而且速度太慢。
同样的算法一样的效果其实只需不到半秒甚至不到0.1s

Clear;
t0 = AbsoluteTime[];
ifs :=
FoldList[#2.{#[], #;

init = {0., 0.};
...

你mathematica怎么学的? 这么厉害

文章乃身外之物,要多考虑编辑、审稿人和读者的感受。

29楼2014-10-11 15:30:09