kaiyu
2020-12-07 cf42a18aba4da77141dd0ea65918444636f899e1
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
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);
    }
 
}