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
來源：知乎
著作權歸作者所有。商業轉載請聯系作者獲得授權，非商業轉載請注明出處。