博客
关于我
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/

你可能感兴趣的文章
Neo4j图数据库的介绍_图数据库结构_节点_关系_属性_数据---Neo4j图数据库工作笔记0001
查看>>
Neo4j安装部署及使用
查看>>
Neo4j电影关系图Cypher
查看>>
Neo4j的安装与使用
查看>>
Neo4j(1):图数据库Neo4j介绍
查看>>
Neo4j(2):环境搭建
查看>>
Neo4j(4):Neo4j - CQL使用
查看>>
NEO改进协议提案1(NEP-1)
查看>>
Neo私链
查看>>
NervanaGPU 项目使用教程
查看>>
Nerves 项目教程
查看>>
nessus快速安装使用指南(非常详细)零基础入门到精通,收藏这一篇就够了
查看>>
Nessus漏洞扫描教程之配置Nessus
查看>>
Nest.js 6.0.0 正式版发布,基于 TypeScript 的 Node.js 框架
查看>>
nested exception is org.apache.ibatis.builder.BuilderException: Error parsing Mapper XML.
查看>>
nestesd exception is java .lang.NoSuchMethodError:com.goolge.common.collect
查看>>
nestJS学习
查看>>
net core 环境部署的坑
查看>>
NET Framework安装失败的麻烦
查看>>
Net 应用程序如何在32位操作系统下申请超过2G的内存
查看>>