matlab函数 bsxfun(高效代码)
来源:互联网 发布:linux 用代码链接网络 编辑:程序博客网 时间:2024/05/17 04:59
bsxfun是一个matlab自版本R2007a来就提供的一个函数,作用是”applies an element-by-element binary operation to arrays a and b, with singleton expansion enabled.”
函数功能:两个数组间元素逐个计算的二值操作
使用方法:C=bsxfun(fun,A,B)
两个数组A合B间元素逐个计算的二值操作,fun是函数句柄或者m文件,也可以为如下内置函数:
@plus 加
@minus 减
@times 数组乘《Simulink与信号处理》
@rdivide 左除
@ldivide 右除
@power 数组幂乘
@max 二值最大值
@min 二值最小值
@rem 余数
@mod 求模
@atan2 四象限反正切
@hypot 平方和的平方根
@eq 等于
@ne 不等于
@lt 小于www.iLoveMatlab.cn
@le 小于或等于
@gt 大于
@ge 大于或等于
@and 逻辑并
@or 逻辑或
@xor 逻辑异或
三种方法对比:
function test_dist(m, n) X = rand(m,n); D = use_bsxfun(X); D = use_repmat(X); D = use_for(X);endfunction D=use_bsxfun(X) v = dot(X,X); D = bsxfun(@plus,v,v')-2*(X'*X);endfunction D=use_repmat(X) v = dot(X,X); D = repmat(v,length(v),1) + repmat(v',1,length(v)) - 2*(X'*X);endfunction D=use_for(X) D = zeros(size(X,2)); for i=1:size(D,1) for j=i+1:size(D,2) D(i,j) = sum((X(:,i)-X(:,j)).^2); end end D = max(D,D');end
应用举例:
假设我们有数据A和B, 每行是一个样本,每列是一个特征。我们要计算高斯核,既:
k(||x-xc||)=exp{- ||x-xc||^2/(2*σ)^2) } 其中xc为核函数中心,σ为函数的宽度参数 , 控制了函数的径向作用范围。
当然可以用双重for实现(如果第一直觉是用三重for的话…)。
K1 = zeros(size(A,1),size(B,1));for i = 1 : size(A,1) for j = 1 : size(B,1) K1(i,j) = exp(-sum((A(i,:)-B(j,:)).^2)/beta); endend
使用2,000×1,000大小的A和B, 运行时间为88秒。
考虑下面向量化后的版本:
sA = (sum(A.^2, 2));sB = (sum(B.^2, 2));K2 = exp(bsxfun(@minus,bsxfun(@minus,2*A*B', sA), sB')/beta);
使用同样数据,运行时间仅0.85秒,加速超过100倍。
如要判断两者结果是不是一样,可以如下
assert(all(all(abs(K1-K2)<1e-12)))
C = bsxfun(fun,A,B) appliesthe element-by-element binary operation specified
- matlab函数 bsxfun(高效代码)
- matlab函数 bsxfun浅谈(转载)
- matlab函数 bsxfun浅谈(转载)
- matlab函数 bsxfun浅谈
- matlab函数 bsxfun浅谈
- matlab函数 bsxfun浅谈
- matlab函数 bsxfun
- matlab函数 bsxfun浅谈
- Matlab bsxfun函数解析
- matlab中 bsxfun函数
- matlab之bsxfun函数
- matlab函数之bsxfun
- Matlab的bsxfun函数
- Matlab bsxfun函数
- bsxfun函数-matlab
- matlab函数之bsxfun
- MATLAB中的bsxfun函数
- matlab学习笔记 bsxfun函数
- TextView中换行符(\\n)显示不正常的问题
- Maven学习记录4——Maven管理多模块项目
- MongoDB 复制集模式Replica Sets
- HDU 5437 Alisha’s Party
- 策略模式
- matlab函数 bsxfun(高效代码)
- MyBatis入门教程(二):实现CRUD
- codeforces #576E Painting Edges 分治+并查集
- 生成带有透明渐变倒影的图片
- LoadRunner的-26612报错解决过程(关联错误)
- ORA-12838的问题
- HDU 1024 Max Sum Plus Plus
- linux中在使用slickedit时怎么修改工程类型,比如将.exe修改成.a库
- kafka集群部署及仍然存在的问题