matlab vs python: 跑循环的速度对比
本文共 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/
你可能感兴趣的文章
MySQL Join算法与调优白皮书（二）
查看>>
Mysql order by与limit混用陷阱
查看>>
Mysql order by与limit混用陷阱
查看>>
mysql order by多个字段排序
查看>>
MySQL Order By实现原理分析和Filesort优化
查看>>
mysql problems
查看>>
mysql replace first,MySQL中处理各种重复的一些方法
查看>>
MySQL replace函数替换字符串语句的用法（mysql字符串替换）
查看>>
mysql replace用法
查看>>
Mysql Row_Format 参数讲解
查看>>
mysql select， from ，join ，on ，where groupby,having ,order by limit的执行顺序和书写顺序
查看>>
MySQL Server 5.5安装记录
查看>>
mysql server has gone away
查看>>
mysql slave 停了_slave 停止。求解决方法
查看>>
MySQL SQL 优化指南：主键、ORDER BY、GROUP BY 和 UPDATE 优化详解
查看>>
MYSQL sql语句针对数据记录时间范围查询的效率对比
查看>>
mysql sum 没返回,如果没有找到任何值，我如何在MySQL中获得SUM函数以返回'0'？
查看>>
mysql Timestamp时间隔了8小时
查看>>
Mysql tinyint(1)与tinyint(4)的区别
查看>>
mysql union orderby 无效
查看>>