博客
关于我
matlab vs python: 跑循环的速度对比
阅读量:159 次
发布时间:2019-02-28

本文共 2144 字,大约阅读时间需要 7 分钟。

测试1

matlab代码

N = 20:25;iters = 2.^N;time = zeros(1,length(N));a = 0.111;b = 0.222;for k = 1:length(N)    r = 0;    t1 = clock;        for i = 1:2^N(k)        r = 0.5*a + 0.6*b;    end    t2 = clock;    time(k) = etime(t2,t1);    endplot(iters, time)xlabel('iter')ylabel('time(/s)')

python代码

N = range(20,26)iters = [2**n for n in N]ts = []a, b = 0.111, 0.222for n in N:        t1 = time.time()    for i in range(2**n):        r = 0.5*a + 0.6*b    t2 = time.time()    ts.append(t2-t1)_, ax = plt.subplots()ax.plot(iters, ts)ax.set_xlabel('iter')ax.set_ylabel('time(/s)')

结果对比

将两者数据画到一起,方便对比。
结论:随着循环增多,两者消耗时间都线性增大。对于这个测试案例(两个乘法和一个加法)。python约比matlab慢60倍
在这里插入图片描述

测试2

matlab代码

N = 20:25;iters = 2.^N;time = zeros(1,length(N));a = 0.111;b = 0.222;M = [0.111,0.222;0.111,0.222];for k = 1:length(N)    r = 0;    t1 = clock;        for i = 1:2^N(k)        r = M(1,1)*a + M(1,2)*b;    end    t2 = clock;    time(k) = etime(t2,t1);    endfigure;plot(iters, time)xlabel('iter')ylabel('time(/s)')

python代码

N = range(20,26)iters = [2**n for n in N]ts = []M = np.array([[0.111, 0.222],[0.111, 0.222]])a, b = 0.111, 0.222for n in N:        t1 = time.time()    for i in range(2**n):        r = M[0,0]*a + M[0,1]*b    t2 = time.time()    ts.append(t2-t1)_, ax = plt.subplots()ax.plot(iters, ts)ax.set_xlabel('iter')ax.set_ylabel('time(/s)')

结果对比

将两者数据画到一起,方便对比。
结论:

  • 随着循环增多,两者消耗时间都线性增大。python约比matlab慢110倍
  • 将此测试结果与测试1对比, 可猜想:仅仅是在2*2矩阵中索引一个数,python也要比matlab很多倍,猜想慢110-60=50倍。再通过一个测试3来验证下猜想。
    在这里插入图片描述

测试3

matlab代码

N = 20:25;iters = 2.^N;time = zeros(1,length(N));a = 0.111;b = 0.222;M = [0.111,0.222;0.111,0.222];for k = 1:length(N)    r = 0;    t1 = clock;        for i = 1:2^N(k)        r = M(1,1);    end    t2 = clock;    time(k) = etime(t2,t1);    endfigure;plot(iters, time)xlabel('iter')ylabel('time(/s)')

python代码

N = range(20,26)iters = [2**n for n in N]ts = []M = np.array([[0.111, 0.222],[0.111, 0.222]])a, b = 0.111, 0.222for n in N:        t1 = time.time()    for i in range(2**n):        r = M[0,0]    t2 = time.time()    ts.append(t2-t1)_, ax = plt.subplots()ax.plot(iters, ts)ax.set_xlabel('iter')ax.set_ylabel('time(/s)')

结果对比

猜想正确,仅仅是2*2矩阵索引一个数,python也比matlab慢50倍。
在这里插入图片描述

转载地址:http://wthd.baihongyu.com/

你可能感兴趣的文章
ndk-cmake
查看>>
NdkBootPicker 使用与安装指南
查看>>
ndk特定版本下载
查看>>
NDK编译错误expected specifier-qualifier-list before...
查看>>
Neat Stuff to Do in List Controls Using Custom Draw
查看>>
Necurs僵尸网络攻击美国金融机构 利用Trickbot银行木马窃取账户信息和欺诈
查看>>
Needle in a haystack: efficient storage of billions of photos 【转】
查看>>
NeHe OpenGL教程 07 纹理过滤、应用光照
查看>>
NeHe OpenGL教程 第四十四课:3D光晕
查看>>
Neighbor2Neighbor 开源项目教程
查看>>
neo4j图形数据库Java应用
查看>>
Neo4j图数据库_web页面关闭登录实现免登陆访问_常用的cypher语句_删除_查询_创建关系图谱---Neo4j图数据库工作笔记0013
查看>>
Neo4j图数据库的介绍_图数据库结构_节点_关系_属性_数据---Neo4j图数据库工作笔记0001
查看>>
Neo4j图数据库的数据模型_包括节点_属性_数据_关系---Neo4j图数据库工作笔记0002
查看>>
Neo4j安装部署及使用
查看>>
Neo4j电影关系图Cypher
查看>>
Neo4j的安装与使用
查看>>
Neo4j(1):图数据库Neo4j介绍
查看>>
Neo4j(2):环境搭建
查看>>
Neo4j(3):Neo4j Desktop安装
查看>>