% 此函数用于生成拉丁超立方体或 sobol % 基于用户提供的约束和所需大小的准随机集% 用于敏感性分析或 GA 总体初始化% 注意：需要统计工具箱%% 示例用法： % 在这个例子中，4 个设计变量是通过一个矩阵隐式定义的% 表示上限和下限（分别为第一行和第二行） % % 任何需要整数约束的变量都通过% 标志向量，其索引对应于它们各自的索引% 变量。 True 表示变量受整数约束。 % 在下面的例子中，只有第三个变量有整数约束% 应用% % 调用这个函数，在这个例子中将创建一个 200 x N 的变量在应用的界限和约束内设置的 % sobol。 % % myDOEbnds = [-1.75, 32.4, -1, -234.2; 3.95、112.56、9、156.37]； % myDOEintcon = [0 0 1 0]; % myDOEset = createdoe(200,myDO

sobol+matlab+代码使用 Monte Carlo 抽样的 Sobol 敏感性分析 此代码演示了使用蒙特卡罗采样的 Sobol 灵敏度分析方法。 采样方法是MATLAB中的Sobol序列。 测试函数是gmath函数 有关详细说明，请参阅以下参考资料。 参考： [1] Sobol, IM “非线性数学模型的全球敏感性指数及其蒙特卡罗估计。” 模拟中的数学和计算机 55(1), 2001: 271-280。 [2] IM Sobol、S. Tarantola、D. Gatelli、SS Kucherenko、W. Mauntz，在修复全局敏感性分析中的非必要因素时估计近似误差，可靠性工程与系统安全，92(7) 2007：957-960。

sobol 算子程序.txt

内容中包含3个例子说明在sobol采样后基于方差的全局灵敏度分析，不明白的地方可以参照说明书中的指导

Sobol序列代码

Sobol敏感性分析-Sobol.rar 采用网上的一篇博客编写了Sobol敏感性分析程序，博客地址为 blog.csdn.net/xiaosebi1111/article/details/46517409

Sobol全局敏感性分析Matlab代码，输出一阶敏感度Sol_1及总敏感度Sol_t。 对简单函数来说，自己仿照构造一个目标函数Sobol_obj即可； 如果分析对象是Matlab外部的模型，其实就不需要Sobol_obj了，代码中的kp就是模型参数值，output、c_out_1和c_out_2就是kp代入模型算出来的模型输出，用实际数据替换一下就行了。 本程序参考文献：Bilal, N. (2014). Implementation of Sobol's Method of Global Sensitivity Analysis to a Compressor Simulation Model. International Compressor Engineering Conference. Paper 2385.

Sobol随机序列是一种低差异序列。分布均匀性好，耗时少。