本例改编自 Box 和 Draper (1987),使用 Stochastic Optimization.jmp。一个化学反应将化学物质“A”转换为化学物质“B”。化学物质“B”的生成量为反应时间和反应温度的函数。时间越长,温度越高,生成的“B”就越多。但是,时间过长,温度过高还会导致一些“B”转换为另一化学物质“C”。多长的反应时间和多高的反应温度将使生成的“B”最多、“A”和“C”最少呢?反应应在高温下快速进行还是在低温下慢速进行呢?
化学反应
目标是使生成的化学物质“B”最多。一个方法是进行实验并拟合作为时间和温度函数的反应产量(化学物质“B”的量)的响应曲面模型。但是,由于已知的化学反应模型,基于阿伦尼乌斯定律,可以直接计算反应产量。产量列包含产量的公式。该公式是反应时间(小时)和反应速度 k1k2 的函数。反应速度是反应温度(开氏度)和已知的物理常数 θ1θ2θ3θ4 的函数。因此,产量反应时间反应温度的函数。
1.
选择帮助 > 样本数据库,然后打开 Stochastic Optimization.jmp
2.
选择图形 > 刻画器
3.
产量分配给 Y,预测公式
4.
点击展开中间公式。点击确定
5.
从“预测刻画器”红色小三角菜单,选择优化和意愿 > 最大化意愿
刻画器将最大化产量,并将图形设置为反应时间反应温度的最优值。
产量最大值
最大产量约为 0.62,此时反应时间为 0.115 小时,反应温度为 540 开氏度,即在高温下快速反应。(由于优化过程的起始值是随机的,您的结果可能略有不同。)
在生产环境中,有时无法精确控制过程输入。若输入(反应时间反应温度)具有随机变异,产量会有什么变化呢?此外,若产量具有规格限,规格外需要被丢弃的批次所占百分比是多少?“模拟器”还可以帮助我们调查在指定反应时间反应温度变异的条件下产量的变异和缺陷率情况。
8.
“模拟器”初始设置中所示,填充因子参数使得温度服从标准差为 1 的正态加权分布,时间服从标准差为 0.03 的正态加权分布。均值参数默认为当前因子值。
“模拟器”初始设置
产量具有下规格限 0.55,将它设置为列属性并在“模拟器”初始设置上显示为红线。为输入因子设置随机变异后,就可以运行模拟来研究产量的变异和缺陷率。
10.
点击模拟按钮。
模拟结果
模拟结果中所示,预测的产量为 0.62,但是若因子具有给定的变异,平均产量为 0.60,标准差为 0.03。
对于反应温度反应时间的其他设置而言缺陷率如何呢?假定您将反应温度改为 535,然后将反应时间设置为使产量最大化的值。
11.
最大产量记住的设置
13.
现在,将反应温度的均值改为 535。
14.
反应时间图中的红色虚线移至最大化产量的值(大约 0.16)。将反应时间的均值改为 0.16。
15.
点击模拟
温度为 535 时的缺陷率
温度为 535 时的缺陷率中所示,缺陷率降低到约为 1.8%,这比 6.0% 低多了。因此,您看到的是使产量最大化的该固定(无变异性)设置与存在因子变异时缺陷率最小化的设置不同。
通过运行模拟实验,您可以找到使缺陷率最小化的反应温度反应时间设置。为此,您模拟反应温度反应时间设计的每个点的缺陷率,然后拟合缺陷率的预测模型并使它最小化。
结果显示在“高斯过程”模型拟合的结果中。由于模拟中随机抽签,您的结果可能略有不同。“高斯过程”平台自动打开“预测刻画器”。
“高斯过程”模型拟合的结果
19.
要查找使缺陷率最小化的温度时间设置,请从“预测刻画器”红色小三角菜单中选择优化和意愿 > 最大化意愿
最小缺陷率的设置
使缺陷率最小化的设置约是反应温度 = 526,反应时间 = 0.3。
20.
点击传输回因子设置按钮。
这将温度时间设置为使在原始“刻画器”报表窗口中缺陷率最小化的哪些设置。
22.
最小缺陷设置
24.
指定新设置后,点击模拟按钮以估计新设置下的缺陷率。
更低的缺陷率
在新设置下,缺陷率为 0.05%,比使产量最大化的设置下的 6.0% 低很多。这几乎降低了大约 120 倍。我们回想一下,第一种设置下的平均产量为 0.60,而新平均产量为 0.59。在缺陷率降低 120 倍的情况下,平均产量减少 0.01 是完全可以接受的。
25.
点击记住的设置单选按钮可以查看每种设置的刻画器。
设置比较
化学家现在知道了将哪个设置用于质量过程。若因子没有变异,则使产量最大化的设置为高温下快速反应。但是,若该过程输入有类似于我们模拟的变异,使产量最大化的设置会导致高的缺陷率。因此,要在有因子变异的情况下使缺陷率最小化,设置应为低温下慢速反应。