拉丁超立方抽样 Latin hypercube sampling，java 代码

以下是使用Java实现拉丁超立方抽样的示例代码：

import java.util.*;

public class LatinHypercubeSampling {

    public static void main(String[] args) {
        int n = 10; // 样本数量
        int m = 5; // 变量数量
        double[][] samples = generateLatinHypercubeSamples(n, m);
        for (int i = 0; i < n; i++) {
            System.out.println(Arrays.toString(samples[i]));
        }
    }

    public static double[][] generateLatinHypercubeSamples(int n, int m) {
        double[][] samples = new double[n][m];
        Random rand = new Random();
        double[] intervals = new double[n + 1];
        double delta = 1.0 / n;
        for (int i = 0; i <= n; i++) {
            intervals[i] = i * delta;
        }
        for (int j = 0; j < m; j++) {
            for (int i = 0; i < n; i++) {
                double randValue = intervals[i] + rand.nextDouble() * delta;
                samples[i][j] = randValue;
            }
            shuffle(samples, j);
        }
        return samples;
    }

    public static void shuffle(double[][] samples, int j) {
        Random rand = new Random();
        for (int i = samples.length - 1; i > 0; i--) {
            int index = rand.nextInt(i + 1);
            double temp = samples[index][j];
            samples[index][j] = samples[i][j];
            samples[i][j] = temp;
        }
    }
}

说明：

- `generateLatinHypercubeSamples`方法用于生成拉丁超立方抽样，其中`n`为样本数量，`m`为变量数量。
- `shuffle`方法用于对每个变量的样本进行随机打乱，避免样本之间存在相关性。
- 示例代码中仅输出了样本结果，可以根据需要进行进一步处理。