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

你可能感兴趣的文章
nacos源码 nacos注册中心1.4.x 源码 nacos源码如何下载 nacos 客户端源码下载地址 nacos discovery下载地址(一)
查看>>
nacos源码 nacos注册中心1.4.x 源码 spring cloud alibaba 的discovery做了什么 nacos客户端是如何启动的(二)
查看>>
nacos源码 nacos注册中心1.4.x 源码 如何注册服务 发送请求,nacos clinet客户端心跳 nacos 注册中心客户端如何发送的心跳 (三)
查看>>
Nacos源码分析:心跳机制、健康检查、服务发现、AP集群
查看>>
nacos看这一篇文章就够了
查看>>
Nacos简介、下载与配置持久化到Mysql
查看>>
Nacos简介和控制台服务安装
查看>>
Nacos管理界面详细介绍
查看>>
Nacos编译报错NacosException: endpoint is blank
查看>>
nacos自动刷新配置
查看>>
nacos运行报错问题之一
查看>>
Nacos部署中的一些常见问题汇总
查看>>
NACOS部署,微服务框架之NACOS-单机、集群方式部署
查看>>
Nacos配置Mysql数据库
查看>>
Nacos配置中心中配置文件的创建、微服务读取nacos配置中心
查看>>
Nacos配置中心集群原理及源码分析
查看>>
nacos配置在代码中如何引用
查看>>