Matlab代码优化之道
一、 遵守Performance Acceleration的规则
关于什么是“Performance Acceleration”请参阅matlab的帮助文件。
1、只有使用以下数据类型,matlab才会对其加速:
logical,char,int8,uint8,int16,uint16,int32,uint32,double而语句中如果使用了非以上的数据类型则不会加速,如numeric,cell,structure,single,function handle,java classes,user classes,int64,uint64
2、matlab不会对超过三维的数组进行加速。
3、当使用for循环时,只有遵守以下规则才会被加速:
a、for循环的范围只用标量值来表示;
b、for循环内部的每一条语句都要满足上面的两条规则,即只使用支持加速的数据类型,只使用三维以下的数组;
c、循环内只调用了内建函数(build-in function)。
4、当使用if、elseif、while和switch时,其条件测试语句中只使用了标量值时,将加速运行。
5、不要在一行中写入多条操作,这样会减慢运行速度。即不要有这样的语句:
x = a.name; for k=1:10000, sin(A(k)), end;
6、当某条操作改变了原来变量的数据类型或形状(大小,维数)时将会减慢运行速度。
7、应该这样使用复常量x = 7 + 2i,而不应该这样使用:x = 7 + 2*i,后者会降低运行速度。
二、 遵守三条规则
1、尽量避免使用循环
改进这样的状况有两种方法:
a、尽量用向量化的运算来代替循环操作。
如
i=0;
for t = 0:.01:10
i = i+1;
y(i) = sin(t);
end
替换为:
t = 0:.01:10;
y = sin(t);
速度将会大大加快。最常用的使用vectorizing技术的函数有:All、diff、ipermute、permute、reshapen squeeze、any、find、logical、prod、shiftdim、sub2ind、cumsum、ind2sub、ndgrid、repmat、sort、sum 等。
请注意matlab文档中还有这样一句补充:
“Before taking the time to vectorize your code, read the section on Performance Acceleration.You may be able to speed up your program by just as much using theMATLAB JITAccelerator instead of vectorizing.”。
b、在必须使用多重循环时下,如果两个循环执行的次数不同,则在循环的外环执行循环次数少的,内环执行循环次数多的。这样可以显著提高速度。
2、预分配矩阵空间
a、预分配矩阵空间,即事先确定变量的大小,维数。这一类的函数有zeros、ones、cell、struct、repmat等。
b、当要预分配一个非double型变量时使用repmat函数以加速,如将以下代码:
A = int8(zeros(100)); 换成:A = repmat(int8(0), 100, 100);
c、当需要扩充一个变量的大小、维数时使用repmat函数。
3、优先使用matlab内建函数
a、优先使用matlab内建函数,将耗时的循环编写进MEX-File中以获得加速。者利用matlab提供的工具将程序转化为C语言、Fortran语言。
b、使用Functions而不是Scripts 。
4、改用更有效的算法
文章最后发布于: 2015-05-08 20:56:00, 来自CSDN的迁移
参考文章
相关文章
- matlab读取多幅图片,并对读取的图片降采样和双三次插值
- Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动
- 风电最大化消纳的热电联产机组联合优化控制(Matlab代码实现)
- 基于贝叶斯推理估计稳态 (ST) 和非稳态 (NS) LPIII 模型分布拟合到峰值放电(Matlab代码实现)
- 使用带外空间信息选择毫米波波束(Matlab代码实现)
- 【图像处理】从点云数据中提取边界(识别和追踪)(Matlab代码实现)
- 【MATLAB】关系运算与逻辑运算
- 【图像处理】基于matlab蚁群聚类图像边缘检测
- 【图像处理】基于MATLAB的小波变换DWT图像检索(实现)
- 基于LS最小二乘法的OFDM信道估计误码率matlab仿真
- 基于ACGWO混沌灰狼优化算法的MATLAB对比仿真,对比标准的GWO
- 基于PSO粒子群优化的SVM(PSO-SVM)的短期电力负荷预测matlab仿真
- 基于kmeans算法的数据聚类matlab仿真
- Matlab Newton‘s method
- 【故障诊断】用于轴承故障诊断的候选故障频率优化克改进包络频谱研究(Matlab代码实现)