PHP项目差分隐私与Laplace

wen PHP项目 2026-07-19 188

本文目录导读：

PHP项目差分隐私与Laplace

差分隐私基础概念
Laplace机制实现
实际应用示例
高级应用：柱状图发布
性能优化与安全考虑

我来详细介绍PHP项目中实现差分隐私和Laplace机制的方案。

差分隐私基础概念

差分隐私的核心目标：在输出统计结果时，保护个体数据是否存在于数据集中不被泄露。

Laplace机制实现

1 基础Laplace分布生成器

<?php
class LaplaceDistribution {
    /**
     * 生成Laplace噪声
     * @param float $location 位置参数μ
     * @param float $scale 尺度参数b
     * @return float
     */
    public static function generateNoise($location = 0.0, $scale = 1.0) {
        $u = mt_rand() / mt_getrandmax() - 0.5;
        return $location - $scale * ($u >= 0 ? log(1 - 2 * abs($u)) : -log(1 - 2 * abs($u)));
    }
    /**
     * 使用Box-Muller方法生成Laplace噪声
     * @param float $scale
     * @return float
     */
    public static function generateBoxMuller($scale = 1.0) {
        $u = mt_rand() / mt_getrandmax();
        $v = mt_rand() / mt_getrandmax();
        return $scale * log($u / $v);
    }
}

2 差分隐私查询器

<?php
class DifferentialPrivacyQuery {
    private $epsilon; // 隐私预算
    private $sensitivity; // 全局敏感度
    public function __construct($epsilon, $sensitivity) {
        $this->epsilon = $epsilon;
        $this->sensitivity = $sensitivity;
    }
    /**
     * 对数值型查询添加Laplace噪声
     * @param callable $queryFunction 原始查询函数
     * @param array $data 数据集
     * @return float
     */
    public function noisyQuery($queryFunction, $data) {
        $trueResult = $queryFunction($data);
        $scale = $this->sensitivity / $this->epsilon;
        return $trueResult + LaplaceDistribution::generateNoise(0, $scale);
    }
}

实际应用示例

1 统计查询保护

<?php
class CensusDataProtector {
    private $epsilon;
    const SENSITIVITY_COUNT = 1.0;
    const SENSITIVITY_SUM = 1.0;
    const SENSITIVITY_AVERAGE = 1.0;
    public function __construct($epsilon = 1.0) {
        $this->epsilon = $epsilon;
    }
    /**
     * 保护后的计数查询
     */
    public function noisyCount($data, $condition = null) {
        $trueCount = $condition ? count(array_filter($data, $condition)) : count($data);
        $scale = self::SENSITIVITY_COUNT / $this->epsilon;
        return max(0, $trueCount + LaplaceDistribution::generateNoise(0, $scale));
    }
    /**
     * 保护后的求和查询
     */
    public function noisySum($data) {
        $trueSum = array_sum($data);
        $scale = self::SENSITIVITY_SUM / $this->epsilon;
        return $trueSum + LaplaceDistribution::generateNoise(0, $scale);
    }
    /**
     * 保护后的平均值查询
     */
    public function noisyAverage($data) {
        $trueAvg = array_sum($data) / count($data);
        $scale = self::SENSITIVITY_AVERAGE / $this->epsilon;
        return $trueAvg + LaplaceDistribution::generateNoise(0, $scale);
    }
}
// 使用示例
$protector = new CensusDataProtector(0.5);
$ages = [25, 30, 35, 40, 45, 50, 55, 60];
$noisyAvg = $protector->noisyAverage($ages);
echo "加噪后的平均年龄: " . round($noisyAvg, 2);

2 数据库查询保护

<?php
class PrivateDatabaseQuery {
    private $pdo;
    private $epsilon;
    public function __construct($pdo, $epsilon = 1.0) {
        $this->pdo = $pdo;
        $this->epsilon = $epsilon;
    }
    /**
     * 保护用户数量查询
     */
    public function getUserCount($conditions = []) {
        $query = "SELECT COUNT(*) as cnt FROM users";
        if (!empty($conditions)) {
            $query .= " WHERE " . implode(" AND ", $conditions);
        }
        $stmt = $this->pdo->query($query);
        $result = $stmt->fetch(PDO::FETCH_ASSOC);
        $trueCount = (int)$result['cnt'];
        // 添加Laplace噪声
        $scale = 1.0 / $this->epsilon;
        $noisyCount = $trueCount + LaplaceDistribution::generateNoise(0, $scale);
        return max(0, round($noisyCount));
    }
    /**
     * 保护平均工资查询
     */
    public function getAverageSalary($departmentId = null) {
        $query = "SELECT AVG(salary) as avg_salary FROM employees";
        if ($departmentId) {
            $query .= " WHERE department_id = :dept_id";
        }
        $stmt = $this->pdo->prepare($query);
        if ($departmentId) {
            $stmt->bindParam(':dept_id', $departmentId, PDO::PARAM_INT);
        }
        $stmt->execute();
        $result = $stmt->fetch(PDO::FETCH_ASSOC);
        $trueAvg = (float)$result['avg_salary'];
        // 假设工资范围在2000-50000之间，敏感度为48000
        $sensitivity = 48000;
        $scale = $sensitivity / $this->epsilon;
        return $trueAvg + LaplaceDistribution::generateNoise(0, $scale);
    }
}

高级应用：柱状图发布

<?php
class PrivateHistogram {
    private $epsilon;
    public function __construct($epsilon) {
        $this->epsilon = $epsilon;
    }
    /**
     * 发布加噪后的柱状图数据
     */
    public function publishHistogram($data, $bins) {
        $histogram = array_fill_keys($bins, 0);
        foreach ($data as $value) {
            foreach ($bins as $bin) {
                if ($value <= $bin) {
                    $histogram[$bin]++;
                    break;
                }
            }
        }
        // 为每个桶添加Laplace噪声
        $scale = 1.0 / $this->epsilon;
        $noisyHistogram = [];
        foreach ($histogram as $bin => $count) {
            $noise = LaplaceDistribution::generateNoise(0, $scale);
            $noisyHistogram[$bin] = max(0, round($count + $noise));
        }
        return $noisyHistogram;
    }
}

性能优化与安全考虑

1 批量处理优化

<?php
class BatchNoiseGenerator {
    private $cache = [];
    public function generateBatchNoise($count, $scale) {
        if (!isset($this->cache[$scale])) {
            $this->cache[$scale] = [];
        }
        if (count($this->cache[$scale]) < $count) {
            $this->fillCache($count, $scale);
        }
        return array_splice($this->cache[$scale], 0, $count);
    }
    private function fillCache($count, $scale) {
        for ($i = 0; $i < $count * 1.5; $i++) {
            $this->cache[$scale][] = LaplaceDistribution::generateNoise(0, $scale);
        }
    }
}

2 安全随机数生成

<?php
class SecureLaplaceDistribution {
    public static function generateSecureNoise($scale = 1.0) {
        // 使用加密安全的随机数生成器
        $bytes = random_bytes(8);
        $u = unpack('d', $bytes)[1] / 2.0;
        return $scale * log($u > 0.5 ? 2 * (1 - $u) : 2 * $u);
    }
}

使用建议：