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

登录注册

应《网络安全法》要求，自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]

回复此楼

» 收录本帖的淘帖专辑推荐

Funny

杂趣

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

求救！comsol如何把运行到一半的程序结果导出来？已经有3人回复
怎么在按MFC的打开之后打开一个图片让图片显示在MFC窗体里？已经有10人回复
52单片机的程序问题，程序刚开始可以执行，可是执行一段时间后就不灵敏了！求高人指点已经有3人回复
用MPICH2运行并行程序，代码中的printf语句只有在所有进程都结束后才会打印出来么？已经有5人回复
vc6.0编译32位软件在win7的64位电脑中运行的问题已经有3人回复
投稿要求图片为黑白的，如何把matlab运行的程序改成黑白的？已经有8人回复
MS一提交任务，就出现“由于另一个程序正在运行中,此操作无法完成...” 已经有18人回复
用matlab运行代码出不来图像！！已经有3人回复
关于matlab GUI打包生产exe文件，并使其在没有安装matlab的计算机上运行的问题。。已经有4人回复
【求助】关于vb调用dll，做成exe后在别人电脑上无法运行（在线等）已经有14人回复
MS能带计算出现“另一程序正在运行”怎么回事？？--求助啊，弄了一下午都没解决啊。。已经有16人回复
四核电脑最多能同时运行几个C++程序？已经有3人回复
配置一台1万元以内的计算机跑linux系统，运行C语言分子动力学模拟已经有9人回复
同样的MATLAB程序，为何在别人的电脑上运行不了已经有21人回复
【求助】如何将源代码在linux下运行？已经有8人回复
【求助】此程序没有错误，怎么运行不出来？已经有6人回复
【求助】vc6.0运行一个exe文件对图像进行canny边缘检测时就会出现以下错误，请教如何� 已经有5人回复
【求助】此程序也是没有问题了，但是依然运行不出结果已经有7人回复
【求助】此段C语言代码运行不了已经有9人回复
【讨论】如何过写代码运行心理上这一关已经有4人回复

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

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

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

易炀

金虫 (职业作家)

应助: 8 (幼儿园)
金币: 293.4
散金: 1154
红花: 26
沙发: 4
帖子: 3134
在线: 636小时
虫号: 3309299
注册: 2014-07-06
性别: GG
专业: 应用高分子化学与物理

matlab行不

[ 发自小木虫客户端 ]

回复此楼

援雅琴以变调兮，奏愁思之不可长

5楼2014-09-17 10:56:59

已阅回复此楼关注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 (幼儿园)
金币: 4676.2
红花: 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 个回答