Matlab的差分隐私laplace噪声实现

typedef double db;
int sgn( db x) {return x < - EPS ? - 1 : x > EPS ;}

db Uniform(db x=.0){
	return x+(db)(1.0+rand())/(RAND_MAX+2.0);
}

db Lap( db miu = .0, db lambda = 1.0) {
	db U = Uniform(- 0.5) ;
	return miu - lambda * sgn(U) * log( fabs(1.0 - 2* fabs(U)) ) ;
}

作者：可以不用真名吗
链接：https://www.zhihu.com/question/52487363/answer/140322524
来源：知乎
著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。