python实现差分隐私Laplace机制详解
作者:kyton123 时间:2022-01-08 16:37:41
Laplace分布定义:
下面先给出Laplace分布实现代码:
import matplotlib.pyplot as plt
import numpy as np
def laplace_function(x,beta):
result = (1/(2*beta)) * np.e**(-1*(np.abs(x)/beta))
return result
#在-5到5之间等间隔的取10000个数
x = np.linspace(-5,5,10000)
y1 = [laplace_function(x_,0.5) for x_ in x]
y2 = [laplace_function(x_,1) for x_ in x]
y3 = [laplace_function(x_,2) for x_ in x]
plt.plot(x,y1,color='r',label='beta:0.5')
plt.plot(x,y2,color='g',label='beta:1')
plt.plot(x,y3,color='b',label='beta:2')
plt.title("Laplace distribution")
plt.legend()
plt.show()
效果图如下:
接下来给出Laplace机制实现:
Laplace机制,即在操作函数结果中加入服从Laplace分布的噪声。
Laplace概率密度函数Lap(x|b)=1/2b exp(-|x|/b)正比于exp(-|x|/b)。
import numpy as np
def noisyCount(sensitivety,epsilon):
beta = sensitivety/epsilon
u1 = np.random.random()
u2 = np.random.random()
if u1 <= 0.5:
n_value = -beta*np.log(1.-u2)
else:
n_value = beta*np.log(u2)
print(n_value)
return n_value
def laplace_mech(data,sensitivety,epsilon):
for i in range(len(data)):
data[i] += noisyCount(sensitivety,epsilon)
return data
if __name__ =='__main__':
x = [1.,1.,0.]
sensitivety = 1
epsilon = 1
data = laplace_mech(x,sensitivety,epsilon)
for j in data:
print(j)
来源:https://blog.csdn.net/qq_38242289/article/details/80798952
标签:python,差分隐私,Laplace机制
![](/images/zang.png)
![](/images/jiucuo.png)
猜你喜欢
转换字符串单词的第一个字母为大写
2007-10-18 10:50:00
asp如何做一个全面的服务器探测器?
2010-07-12 18:55:00
Django之创建引擎索引报错及解决详解
2023-06-11 10:55:15
![](https://img.aspxhome.com/file/2023/5/59125_0s.png)
Python+Qt身体特征识别人数统计源码窗体程序(使用步骤)
2021-06-03 10:40:54
![](https://img.aspxhome.com/file/2023/5/66745_0s.jpg)
Python读取YAML文件过程详解
2021-04-21 07:12:25
python使用wxPython打开并播放wav文件的方法
2023-08-23 10:30:58
70+漂亮且极具亲和力的导航菜单设计推荐[译]
2009-02-23 13:15:00
![](https://img.aspxhome.com/file/UploadPic/up/2009022313285860.gif)
以PHP代码为实例详解RabbitMQ消息队列中间件的6种模式
2023-06-08 12:15:53
![](https://img.aspxhome.com/file/2023/2/55492_0s.png)
Js实现仿msn的右下角popup提示窗口
2007-12-27 20:30:00
![](https://img.aspxhome.com/file/UploadPic/200712/27/20071227212229674s.gif)
XML:OpenSearch 浏览器指定搜索应用
2010-05-04 19:37:00
![](https://img.aspxhome.com/file/UploadPic/20105/4/2008101720221064077801-78s.png)
Python实现用手机监控远程控制电脑的方法
2021-06-22 07:57:49
asp定时生成静态HTML的代码
2010-07-02 12:29:00
Python Unittest自动化单元测试框架详解
2023-09-13 23:13:56
![](https://img.aspxhome.com/file/2023/2/67822_0s.png)
ORACLE 回收站当前状态查询整理
2023-07-14 09:23:01
使用go gin来操作cookie的讲解
2023-09-12 14:21:49
那些被我遗忘掉的XHTML标签们
2008-06-07 14:27:00
php广告加载类用法实例
2023-11-14 14:56:53
python SQLAlchemy的Mapping与Declarative详解
2022-12-04 02:48:37
论Asp与XML的关系
2008-03-05 12:01:00
Python内置数据结构列表与元组示例详解
2021-08-17 21:28:14
![](https://img.aspxhome.com/file/2023/1/78801_0s.png)