package com.moral.util;

public class EmissionDataUtil {

    //计算污染源强度（参数：c：设备点测到是数据，x：源点到【设备点与风向垂线交点】的距离，y：设备点到【设备点与风向垂线交点】的距离，u：风速）
    public static Double  getPollutionSourceIntensity(Double c, Double x, Double y, Double u){
        double e = 2.718281828;
        double r1 = 0.281846;
        double r2 = 0.127190;
        double x1 = 0.914370;
        double x2 = 0.964435;
        double π = Math.PI;
        double diffusionCoefficient_y = r1*Math.pow(x,x1);//扩散系数y
        double diffusionCoefficient_z = r2*Math.pow(x,x2);//扩散系数z
        double molecule = c*2*π*u*diffusionCoefficient_y*diffusionCoefficient_z;
        double power = -1/2*(Math.pow(y,2)/Math.pow(diffusionCoefficient_y,2));
        double denominator = Math.pow(e,power);
        double pollutionSourceIntensity = molecule/denominator;
        return pollutionSourceIntensity;
    }

    public static void main(String[] args) {
        double result = EmissionDataUtil.getPollutionSourceIntensity(100.0,10.0,5.0,3.0);
        System.out.println(result);
    }

}