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

以下是一个简单的 Java 代码示例，用于生成拉丁超立方抽样：

import java.util.Random;

public class LatinHypercubeSampling {

    public static void main(String[] args) {
        int n = 10; // 样本点数
        int k = 3; // 维度

        double[][] lhs = generateLatinHypercube(n, k);
        for (int i = 0; i < n; i++) {
            for (int j = 0; j < k; j++) {
                System.out.print(lhs[i][j] + " ");
            }
            System.out.println();
        }
    }

    public static double[][] generateLatinHypercube(int n, int k) {
        double[][] lhs = new double[n][k];
        double[] fractions = new double[n];

        Random rand = new Random();

        // 初始化分数数组（从0到1之间的小数）
        for (int i = 0; i < n; i++) {
            fractions[i] = (i + rand.nextDouble()) / n;
        }

        // 对于每个维度，将样本点分配到分数区间中
        for (int j = 0; j < k; j++) {
            double[] sequence = new double[n];
            for (int i = 0; i < n; i++) {
                sequence[i] = fractions[i];
            }
            shuffle(sequence);

            for (int i = 0; i < n; i++) {
                lhs[i][j] = sequence[i];
            }
        }

        return lhs;
    }

    // Fisher-Yates 随机置换算法
    public static void shuffle(double[] array) {
        Random rand = new Random();
        for (int i = array.length - 1; i > 0; i--) {
            int j = rand.nextInt(i + 1);
            double temp = array[i];
            array[i] = array[j];
            array[j] = temp;
        }
    }

}

这个代码示例生成了一个大小为10x3的拉丁超立方抽样，并将其打印到控制台上。请注意，这只是一个非常简单的示例，实际使用中可能需要进行更多的参数设置和优化。