CaiXiang/SPC
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

新增 MR MS 控制图

Browse Source
This commit is contained in:
caixiang
2022-07-29 10:28:14 +08:00
parent fde6a8b684
commit e7a542d386
14 changed files with 628 additions and 36 deletions
Download Patch File Download Diff File Expand all files Collapse all files

90
ym-process-inspection/src/main/java/com/cnbm/processInspection/controller/ProcessInspectionController.java
View File

@@ -1,34 +1,20 @@
package com.cnbm.processInspection.controller;
import com.cnbm.basic.dto.FactoryDTO;
import com.cnbm.basic.dto.UnitDTO;
import com.cnbm.basic.entity.ProductFeatures;
import com.cnbm.common.spc.math.StandardDiviation;
import com.cnbm.common.spc.util.DataUtils;
import com.cnbm.common.utils.Result;
import com.cnbm.common.vo.R;
import com.cnbm.influx.config.InfluxClient;
import com.cnbm.influx.constant.Constant;
import com.cnbm.influx.param.PageInfo;
import com.cnbm.influx.param.QueryDataParam;
import com.cnbm.influx.param.Range;
import com.cnbm.influx.param.Tag;
import com.cnbm.influx.template.Event;
import com.cnbm.processInspection.dto.GraphArg;
import com.cnbm.processInspection.dto.InterpretationListArg;
import com.cnbm.processInspection.dto.XbarRGraphData;
import com.cnbm.processInspection.dto.XbarSGraphData;
import com.cnbm.processInspection.graphAnalyzed.MeanStandardDeviationGraph;
import com.cnbm.qualityPlanning.common.StatisticalControlledTest;
import com.cnbm.qualityPlanning.entity.ControlLimit;
import com.cnbm.qualityPlanning.entity.Point;
import com.cnbm.qualityPlanning.entity.SpecificationLimit;
import com.cnbm.qualityPlanning.entity.XbarSPoint;
import com.influxdb.client.domain.WritePrecision;
import com.influxdb.query.FluxRecord;
import com.influxdb.query.FluxTable;
import com.cnbm.processInspection.graphAnalyzed.mr.MeanRGraph;
import com.cnbm.processInspection.graphAnalyzed.ms.MeanStandardDeviationGraph;
import org.slf4j.Logger;
import org.slf4j.LoggerFactory;
import org.springframework.beans.factory.annotation.Autowired;
import org.springframework.web.bind.annotation.PostMapping;
import org.springframework.web.bind.annotation.RequestBody;
import org.springframework.web.bind.annotation.RequestMapping;
@@ -37,9 +23,7 @@ import org.springframework.web.bind.annotation.RestController;
import java.math.BigDecimal;
import java.time.Instant;
import java.util.ArrayList;
import java.util.Date;
import java.util.List;
import java.util.Random;
@RestController
@RequestMapping("/processInspection")
@@ -58,6 +42,12 @@ public class ProcessInspectionController {
productFeatures.setName("LTWeight");
productFeatures.setSampleSize(BigDecimal.TEN);
MeanStandardDeviationGraph meanStandardDeviationGraph = new MeanStandardDeviationGraph(productFeatures);
//判读方案
List<InterpretationListArg> interpretationScheme = new ArrayList<>();
interpretationScheme.add(new InterpretationListArg(1,null,null));
interpretationScheme.add(new InterpretationListArg(5,3,2));
meanStandardDeviationGraph.isNeedInterpretation(interpretationScheme);
QueryDataParam queryDataParam = new QueryDataParam();
queryDataParam.setMeasurement("Weight");
@@ -76,6 +66,40 @@ public class ProcessInspectionController {
return R.ok("成功",xbarSGraphData);
}
@PostMapping("/XbarRGraphTest")
public R<XbarRGraphData> XbarRGraphTest() throws Exception {
ProductFeatures productFeatures = new ProductFeatures();
productFeatures.setSl(new Float(5));
productFeatures.setUsl(new Float(10));
productFeatures.setLsl(new Float(1));
productFeatures.setName("LTWeight");
productFeatures.setSampleSize(BigDecimal.TEN);
MeanRGraph meanRGraph = new MeanRGraph(productFeatures);
//判读方案
List<InterpretationListArg> interpretationScheme = new ArrayList<>();
interpretationScheme.add(new InterpretationListArg(1,null,null));
interpretationScheme.add(new InterpretationListArg(5,3,2));
meanRGraph.isNeedInterpretation(interpretationScheme);
QueryDataParam queryDataParam = new QueryDataParam();
queryDataParam.setMeasurement("Weight");
queryDataParam.setRange(new Range(DataUtils.getBeforeDate(10).toInstant(), Instant.now()));
meanRGraph.initialDate(queryDataParam);
XbarRGraphData xbarSGraphData = new XbarRGraphData(
meanRGraph.getMrGraphEntity(),
meanRGraph.getXbarCL(),
meanRGraph.getRCL(),
meanRGraph.getSpecificationLimit(),
meanRGraph.getProcessCapacity(),
new StandardDiviation(meanRGraph.getXigma(),meanRGraph.getTotalXigma())
);
return R.ok("成功",xbarSGraphData);
}
//筛选条件顺序：
// ① 先工厂（产品表） => ② 检验类型（产品表）=> ③ 产品名（产品表）=> ④ 工序（产品表）=> ⑤ 特性（检验参数name）（ product_features 表） ；；
@@ -109,4 +133,32 @@ public class ProcessInspectionController {
return R.ok("成功",xbarSGraphData);
}
@PostMapping("/XbarRGraph")
public R<XbarRGraphData> XbarRGraph(@RequestBody GraphArg graphArg) throws Exception {
//new 对象
ProductFeatures productFeatures = graphArg.getProductFeatures();
MeanRGraph meanRGraph = new MeanRGraph(productFeatures);
//如果要检验，，那么set 判读方案
if(graphArg.getInterpretationScheme()!=null){
meanRGraph.isNeedInterpretation(graphArg.getInterpretationScheme());
}
QueryDataParam queryDataParam = new QueryDataParam();
queryDataParam.setMeasurement(Constant.measurement);
queryDataParam.setRange(new Range(graphArg.getBegin().toInstant(), graphArg.getEnd().toInstant()));
meanRGraph.initialDate(queryDataParam);
XbarRGraphData xbarSGraphData = new XbarRGraphData(
meanRGraph.getMrGraphEntity(),
meanRGraph.getXbarCL(),
meanRGraph.getRCL(),
meanRGraph.getSpecificationLimit(),
meanRGraph.getProcessCapacity(),
new StandardDiviation(meanRGraph.getXigma(),meanRGraph.getTotalXigma())
);
return R.ok("成功",xbarSGraphData);
}
}

9
ym-process-inspection/src/main/java/com/cnbm/processInspection/dto/InterpretationListArg.java
View File

@@ -15,4 +15,13 @@ public class InterpretationListArg {
private Integer number;
private Integer arg1;
private Integer arg2;
public InterpretationListArg() {
}
public InterpretationListArg(Integer number, Integer arg1, Integer arg2) {
this.number = number;
this.arg1 = arg1;
this.arg2 = arg2;
}
}

45
ym-process-inspection/src/main/java/com/cnbm/processInspection/dto/XbarRGraphData.java Normal file
View File

@@ -0,0 +1,45 @@
package com.cnbm.processInspection.dto;
import com.cnbm.common.spc.math.StandardDiviation;
import com.cnbm.processInspection.graphAnalyzed.mr.MRGraphEntity;
import com.cnbm.qualityPlanning.entity.ControlLimit;
import com.cnbm.qualityPlanning.entity.ProcessCapability;
import com.cnbm.qualityPlanning.entity.SpecificationLimit;
import io.swagger.annotations.ApiModel;
import io.swagger.annotations.ApiModelProperty;
import lombok.Data;
/**
* @Desc: ""
* @Author: caixiang
* @DATE: 2022/7/22 14:18
*/
@Data
@ApiModel(value = "均值极差控制图 结果类")
public class XbarRGraphData {
@ApiModelProperty(value = "控制图list数据")
private MRGraphEntity mrGraphEntity;
@ApiModelProperty(value = "xbar控制图 控制限")
private ControlLimit XBarCL;
@ApiModelProperty(value = "R控制图 控制限")
private ControlLimit RCL;
@ApiModelProperty(value = "工艺规格限")
private SpecificationLimit SL;
@ApiModelProperty(value = "工序能力")
private ProcessCapability processCapability;
@ApiModelProperty(value = "标准差/总体标准差")
private StandardDiviation standardDiviation;
public XbarRGraphData(MRGraphEntity mrGraphEntity, ControlLimit xBarCL, ControlLimit rCL, SpecificationLimit sl, ProcessCapability processCapability, StandardDiviation standardDiviation){
this.mrGraphEntity = mrGraphEntity;
this.XBarCL = xBarCL;
this.RCL = rCL;
this.SL = sl;
this.processCapability = processCapability;
this.standardDiviation = standardDiviation;
}
}

2
ym-process-inspection/src/main/java/com/cnbm/processInspection/dto/XbarSGraphData.java
View File

@@ -1,7 +1,7 @@
package com.cnbm.processInspection.dto;
import com.cnbm.common.spc.math.StandardDiviation;
import com.cnbm.processInspection.graphAnalyzed.MSDGraphEntity;
import com.cnbm.processInspection.graphAnalyzed.ms.MSDGraphEntity;
import com.cnbm.qualityPlanning.entity.ControlLimit;
import com.cnbm.qualityPlanning.entity.ProcessCapability;
import com.cnbm.qualityPlanning.entity.SpecificationLimit;

50
ym-process-inspection/src/main/java/com/cnbm/processInspection/graphAnalyzed/mr/MRGraphEntity.java Normal file
View File

@@ -0,0 +1,50 @@
package com.cnbm.processInspection.graphAnalyzed.mr;
import com.cnbm.qualityPlanning.entity.XbarRPoint;
import com.cnbm.qualityPlanning.entity.XbarSPoint;
import java.util.List;
/**
* @Desc: "xbar-R 结果类"
* @Author: caixiang
* @DATE: 2022/7/21 16:26
*/
public class MRGraphEntity {
private List<XbarRPoint> list;
private Double xbarbar;
private Double Rbar;
public MRGraphEntity(List<XbarRPoint> list, Double xbarbar,Double Rbar){
this.list = list;
this.Rbar = Rbar;
this.xbarbar = xbarbar;
}
public List<XbarRPoint> getList() {
return list;
}
public void setList(List<XbarRPoint> list) {
this.list = list;
}
public Double getXbarbar() {
return xbarbar;
}
public void setXbarbar(Double xbarbar) {
this.xbarbar = xbarbar;
}
public Double getRbar() {
return Rbar;
}
public void setRbar(Double rbar) {
Rbar = rbar;
}
}

290
ym-process-inspection/src/main/java/com/cnbm/processInspection/graphAnalyzed/mr/MeanRGraph.java Normal file
View File

@@ -0,0 +1,290 @@
package com.cnbm.processInspection.graphAnalyzed.mr;
import com.cnbm.basic.entity.ProductFeatures;
import com.cnbm.common.spc.math.Math;
import com.cnbm.common.spc.math.StandardDiviation;
import com.cnbm.common.spc.util.DataUtils;
import com.cnbm.influx.config.InfluxClient;
import com.cnbm.influx.constant.Constant;
import com.cnbm.influx.param.QueryDataParam;
import com.cnbm.influx.param.Range;
import com.cnbm.influx.param.Tag;
import com.cnbm.processInspection.controlCoefficientConstant.XBarRCoefficients;
import com.cnbm.processInspection.dto.InterpretationListArg;
import com.cnbm.qualityPlanning.common.StatisticalControlledTest;
import com.cnbm.qualityPlanning.entity.*;
import com.influxdb.query.FluxRecord;
import com.influxdb.query.FluxTable;
import lombok.Data;
import java.time.Instant;
import java.util.ArrayList;
import java.util.List;
/**
* @Desc: "均值标准差 控制图 , 计算类"
* @Author: caixiang
* @DATE: 2022/7/20 14:26
* 使用方式：① 先new MeanStandardDeviationGraph 对象 ；② 再initialData 初始化数据；③ 再get 控制限
*
* 步骤：
* ① 先读mysql表，查询 product_features 表，先读到 sample_size（样本量）
* ② 再依据 influx.argName == mysql.product_feature.name && 时间段 查询所有的 参数数据
* ③ 拿到参数数据后，分组 整合成List<Point>,
* 计算控制限
* 计算 母体 的 \sigma 、 bar{x} 。。。
* 计算CPK 、CPU 、CPL这些
* ④ 如果配置了判读方案，还要 调用 StatisticalControlledTest Function 检验。
* ⑤
*/
@Data
public class MeanRGraph {
//母体的 μ = xbarbar
private Double miu;
//母体的 xigma ，全局计算的;
private Double xigma;
//母体的 xigma ，全局计算的;
private Double totalXigma;
//需要去检验的 判读方案 eg:[1,2,3,4,6]
private List<InterpretationListArg> interpretationScheme;
//这里 是特有参数 -- begin
private Double rbar;
private Double xbarbar;
private Double a2;
private Double d3;
private Double d4;
//这里 是特有参数 -- end
private Integer sampleSize;
private String argName;
private SpecificationLimit specificationLimit;
private MRGraphEntity mrGraphEntity;
public MeanRGraph(ProductFeatures productFeatures) throws Exception {
if(productFeatures.getSampleSize()==null || productFeatures.getName()==null){
throw new Exception("ProductFeatures 参数异常");
}
this.argName = productFeatures.getName();
this.sampleSize = productFeatures.getSampleSize().intValue();
this.a2 = XBarRCoefficients.getA2(sampleSize);
this.d3 = XBarRCoefficients.getD3(sampleSize);
this.d4 = XBarRCoefficients.getD4(sampleSize);
this.specificationLimit = new SpecificationLimit(
productFeatures.getUsl()==null?null:productFeatures.getUsl(),
productFeatures.getSl()==null?null:productFeatures.getSl(),
productFeatures.getUsl()==null?null:productFeatures.getUsl()
);
}
private Double[] toDoubleArray(Object[] o){
Double[] res= new Double[o.length];
for(int i=0;i<o.length;i++){
res[i] = (Double) o[i];
}
return res;
}
public void isNeedInterpretation(List<InterpretationListArg> list){
this.interpretationScheme = list;
}
/**
* name : 初始化数据函数
* desc : 从influxdb 里面读取数据，然后 加工处理成 我需要的
* 步骤：
* ①
* */
public void initialDate(QueryDataParam queryDataParam){
queryDataParam.setBucket(Constant.bucket);
List<String> dropNames = new ArrayList<>();
dropNames.add("transationId");
dropNames.add("inspectionSheet");
queryDataParam.setDropedTagNames(dropNames);
queryDataParam.setTag(new Tag("argName",argName));
queryDataParam.setRange(new Range(DataUtils.getBeforeDate(10).toInstant(), Instant.now()));
List<FluxTable> query = InfluxClient.Client.query(queryDataParam);
//1. 先从fluxdb 里面提取原始数据
List<Double> originData = new ArrayList<>();
for (FluxTable fluxTable : query) {
List<FluxRecord> records = fluxTable.getRecords();
for (FluxRecord fluxRecord : records) {
//因为 传进去的就是Double 类型，所以取出来，自然而然就是Double
originData.add(Double.parseDouble(fluxRecord.getValueByKey("_value").toString()));
}
}
//2. convert to XbarSPoint
List<List<Double>> doubleListList = DataUtils.fixedGroup(originData, sampleSize);
List<XbarRPoint> list = new ArrayList<>();
List<Double> xbarArray = new ArrayList<>();
List<Double> rArray = new ArrayList<>();
for(int i=0;i<doubleListList.size();i++){
Double[] doubleList = toDoubleArray(doubleListList.get(i).toArray());
Double xbar = Math.getMean(doubleList);
Double r = Math.range(doubleList);
xbarArray.add(xbar);
rArray.add(r);
list.add(new XbarRPoint(
i,
new Double(0),
xbar,
r,
doubleListList.get(i)
)
);
}
//3. 以及一些 控制图的全局参数。
Double xbarbar = Math.getMean(toDoubleArray(xbarArray.toArray()));
this.miu = xbarbar;
StandardDiviation standardDiviation =Math.StandardDiviation(toDoubleArray(originData.toArray()));
this.xigma = standardDiviation.getNormal();
this.totalXigma = standardDiviation.getTotality();
Double rbar = Math.getMean(toDoubleArray(rArray.toArray()));
this.mrGraphEntity = new MRGraphEntity(list,xbarbar,rbar);
//4.判读方案 校验
// 开始
if(this.interpretationScheme != null){
List<Point> forXBar = new ArrayList<>();
List<Point> forR = new ArrayList<>();
for(XbarRPoint x:list){
forXBar.add(new Point(x.getPosition(),x.getXbar()));
forR.add(new Point(x.getPosition(),x.getR()));
}
ControlLimit xbarcl = getXbarCL();
ControlLimit rcl = getRCL();
for(InterpretationListArg arg : this.interpretationScheme){
switch (arg.getNumber())//值必须是整型或者字符型
{
case 1:
StatisticalControlledTest.rule1(forXBar,xbarcl);
StatisticalControlledTest.rule1(forR,rcl);
break;
case 2:
StatisticalControlledTest.rule2(forXBar,xbarcl,arg.getArg1());
StatisticalControlledTest.rule2(forR,rcl,arg.getArg1());
break;
case 3:
StatisticalControlledTest.rule3(forXBar,arg.getArg1());
StatisticalControlledTest.rule3(forR,arg.getArg1());
break;
case 4:
StatisticalControlledTest.rule4(forXBar,arg.getArg1());
StatisticalControlledTest.rule4(forR,arg.getArg1());
break;
case 5:
StatisticalControlledTest.rule5(forXBar,xbarcl,arg.getArg1(), arg.getArg2());
StatisticalControlledTest.rule5(forR,rcl,arg.getArg1(), arg.getArg2());
break;
case 6:
StatisticalControlledTest.rule6(forXBar,xbarcl,arg.getArg1(), arg.getArg2());
StatisticalControlledTest.rule6(forR,rcl,arg.getArg1(), arg.getArg2());
break;
case 7:
StatisticalControlledTest.rule7(forXBar,xbarcl,arg.getArg1());
StatisticalControlledTest.rule7(forR,rcl, arg.getArg1());
break;
case 8:
StatisticalControlledTest.rule8(forXBar,xbarcl, arg.getArg1());
StatisticalControlledTest.rule8(forR,rcl, arg.getArg1());
break;
default:
break;
}
}
for(int i=0;i<list.size();i++){
list.get(i).setrUnsatisfiedRules(forR.get(i).getUnsatisfiedRules());
list.get(i).setXbarUnsatisfiedRules(forXBar.get(i).getUnsatisfiedRules());
}
}
// 结束
}
/**
* desc: get Xbar控制图 的控制限
* 注意：此函数 要在 initialDate()函数执行之后
* */
public ControlLimit getXbarCL(){
return new ControlLimit(
(this.mrGraphEntity.getXbarbar() + a2*this.mrGraphEntity.getRbar()),
this.mrGraphEntity.getXbarbar(),
(this.mrGraphEntity.getXbarbar() - a2*this.mrGraphEntity.getRbar())
);
}
/**
* desc: get R控制图 的控制限
* 注意：此函数 要在 initialDate()函数执行之后
* */
public ControlLimit getRCL(){
return new ControlLimit(
( d4 * this.mrGraphEntity.getRbar() ) ,
this.mrGraphEntity.getRbar() ,
( d3 * this.mrGraphEntity.getRbar())
);
}
/**
* desc: get 工序能力
* 注意：此函数 要在 initialDate()函数执行之后
* */
public ProcessCapability getProcessCapacity(){
SpecificationLimit sp = this.specificationLimit;
List<String> warming = new ArrayList<>();
Float usl = sp.getUSL();
Float lsl = sp.getLSL();
if(usl==null || lsl==null){
if(usl==null && lsl!=null){
if(this.miu<lsl){
warming.add("平均值已不能满足规范要求, 有相当多数据均为不合格, 则表示完全没有工序能力");
return new ProcessCapability(null,null,null,(double)0,warming);
}
Double cpl = ( (this.miu-lsl) / (3 * this.xigma) );
return new ProcessCapability(null,null,null,cpl,warming);
} else if (lsl == null && usl!=null) {
if(this.miu>usl){
warming.add("平均值已不能满足规范要求, 有相当多数据均为不合格, 则表示完全没有工序能力");
return new ProcessCapability(null,null,(double)0,null,warming);
}
Double cpu = ( (usl-this.miu) / (3 * this.xigma) );
return new ProcessCapability(null,null,cpu,null,warming);
}else {
return new ProcessCapability(null,null,null,null,warming);
}
}
//下面就是 usl 和 lsl 都部位null的情况
Float k1 = usl-lsl;
Float k2 = usl+lsl;
Double cp = ( k1 / (6*this.xigma) );
Double k = ( java.lang.Math.abs( (this.miu-k2/2) ) / (k1/2) );
if(k>1){
warming.add("此工艺参数K 大于1, 说明均值μ已经偏离到规范范围意外，即相当多的工艺参数均不满足规范要求，这说明工艺加工结果很差，该设备根本不适用于批量生产");
}
Double cpk = cp * ( 1- k);
if(cpk<0){
warming.add("Cpk=0，该工序完全没有能力 批量生产");
cpk = (double) 0;
}
Double cpu = ( (usl-this.miu) / (3 * this.xigma) );
Double cpl = ( (this.miu-lsl) / (3 * this.xigma) );
return new ProcessCapability(cp,cpk,cpu,cpl,warming);
}
}

4
ym-process-inspection/src/main/java/com/cnbm/processInspection/graphAnalyzed/MSDGraphEntity.java → ym-process-inspection/src/main/java/com/cnbm/processInspection/graphAnalyzed/ms/MSDGraphEntity.java
View File

@@ -1,4 +1,4 @@
package com.cnbm.processInspection.graphAnalyzed;
package com.cnbm.processInspection.graphAnalyzed.ms;
@@ -8,7 +8,7 @@ import java.util.Date;
import java.util.List;
/**
* @Desc: ""
* @Desc: "xbar-s 结果类"
* @Author: caixiang
* @DATE: 2022/7/21 16:26
*/

12
ym-process-inspection/src/main/java/com/cnbm/processInspection/graphAnalyzed/MeanStandardDeviationGraph.java → ym-process-inspection/src/main/java/com/cnbm/processInspection/graphAnalyzed/ms/MeanStandardDeviationGraph.java
View File

@@ -1,4 +1,4 @@
package com.cnbm.processInspection.graphAnalyzed;
package com.cnbm.processInspection.graphAnalyzed.ms;
import com.cnbm.basic.entity.ProductFeatures;
import com.cnbm.common.spc.math.Math;
@@ -20,7 +20,6 @@ import lombok.Data;
import java.time.Instant;
import java.util.ArrayList;
import java.util.Arrays;
import java.util.List;
/**
@@ -99,7 +98,10 @@ public class MeanStandardDeviationGraph {
* */
public void initialDate(QueryDataParam queryDataParam){
queryDataParam.setBucket(Constant.bucket);
queryDataParam.setDropedTagName("transationId");
List<String> dropNames = new ArrayList<>();
dropNames.add("transationId");
dropNames.add("inspectionSheet");
queryDataParam.setDropedTagNames(dropNames);
queryDataParam.setTag(new Tag("argName",argName));
queryDataParam.setRange(new Range(DataUtils.getBeforeDate(10).toInstant(), Instant.now()));
@@ -141,12 +143,12 @@ public class MeanStandardDeviationGraph {
}
//3. 以及一些 控制图的全局参数。
Double xbarbar = Math.getMean(toDoubleArray(xbarArray.toArray()));
this.xbarbar = Math.getMean(toDoubleArray(xbarArray.toArray()));
this.miu = xbarbar;
StandardDiviation standardDiviation =Math.StandardDiviation(toDoubleArray(originData.toArray()));
this.xigma = standardDiviation.getNormal();
this.totalXigma = standardDiviation.getTotality();
Double sbar = Math.getMean(toDoubleArray(sArray.toArray()));
this.sbar = Math.getMean(toDoubleArray(sArray.toArray()));
Double rbar = Math.getMean(toDoubleArray(rArray.toArray()));
this.msdGraphEntity = new MSDGraphEntity(list,xbarbar,sbar,rbar);