版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

本帖产生 1 个程序强帖，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

应助: 237 (大学生)
金币: 6101.9
散金: 1114
红花: 39
帖子: 1380
在线: 553.8小时
虫号: 506699
注册: 2008-02-18
专业: 模式识别

[交流] 运行这段代码到出图片,你的电脑用了多久? 已有14人参与

测试一下Mathematica的效率:

CODE:

ClearAll["Global`*"];

t0 = AbsoluteTime[];

list = {{0, 0}};

last = {{0}, {0}};

GetNextPoint[pt_] := Module[{r,

   t1 = {{0.83, 0.03}, {-0.03, 0.86}}, p1 = {0, 1.5},

   t2 = {{0.2, -0.25}, {0.21, 0.23}}, p2 = {0, 1.5},

   t3 = {{-0.15, 0.27}, {0.25, 0.26}}, p3 = {0, 0.45},

   t4 = {{0, 0}, {0, 0.17}}, p4 = {0, 0}},

  r = Random[];

  If[r <= 0.85, t1.pt + p1, 

   If[r <= 0.91, t2.pt + p2, If[r <= 0.99, t3.pt + p3, t4.pt + p4]]]]

For[i = 0, i < 100000, i++, last = GetNextPoint[last];

 list = Append[list, First[Transpose[last]]];]

h = ListPlot[list, PlotStyle -> {PointSize[0.002], Red}, 

   Axes -> False];

range = First /@ Differences /@ (PlotRange /. Options[h]);

Show[h, AspectRatio -> Last@range/First@range]

Print[AbsoluteTime[] - t0]

回复此楼

» 猜你喜欢

售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有6人回复
面上没中，邀请各位路过的虫友分析一下分数已经有7人回复
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生已经有3人回复
找导师已经有8人回复
国社科又开始会评了，不知道这次命运如何已经有14人回复
麻烦专家们看看评委们的意见（F口面上）已经有13人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有14人回复
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够已经有17人回复
学科评审组评审是指会评吗？已经有5人回复
小白求助投论文要求的highlights应该如何写已经有4人回复

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

1楼 2014-09-15 21:03:46

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 34 个回答

mshwangg

至尊木虫 (正式写手)

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

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

牛人啊，树叶挺漂亮
还可以画别的吗？分享下

赞一下(2人)

回复此楼

2楼2014-09-15 23:08:18

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

walk1997

金虫 (著名写手)

应助: 1 (幼儿园)
金币: 4677.7
红花: 22
帖子: 1066
在线: 798.1小时
虫号: 416039
注册: 2007-06-29
性别: GG
专业: 粒子物理学和场论

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

21.0682050
mshwangg 的机器是AMD的cpu？

赞一下

回复此楼

3楼2014-09-17 08:43:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zzs510

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 639
散金: 6
帖子: 21
在线: 41.6小时
虫号: 1284673
注册: 2011-05-03
性别: GG
专业: 无机材料化学

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

ClearAll["Global`*"];
t0 = AbsoluteTime[];
list = {{0, 0}};
last = {{0}, {0}};
GetNextPoint[pt_] := Module[{r,
t1 = {{0.83, 0.03}, {-0.03, 0.86}}, p1 = {0, 1.5},
t2 = {{0.2, -0.25}, {0.21, 0.23}}, p2 = {0, 1.5},
t3 = {{-0.15, 0.27}, {0.25, 0.26}}, p3 = {0, 0.45},
t4 = {{0, 0}, {0, 0.17}}, p4 = {0, 0}},
r = Random[];
If[r <= 0.85, t1.pt + p1,
If[r <= 0.91, t2.pt + p2, If[r <= 0.99, t3.pt + p3, t4.pt + p4]]]]
For[i = 0, i < 100000, i++, last = GetNextPoint[last];
list = Append[list, First[Transpose[last]]];]
h = ListPlot[list, PlotStyle -> {PointSize[0.002], Red},
Axes -> False];
range = First /@ Differences /@ (PlotRange /. Options[h]);
Show[h, AspectRatio -> Last@range/First@range]
Print[AbsoluteTime[] - t0]

赞一下

回复此楼

4楼2014-09-17 10:31:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 34 个回答