capaplot

Process capability plot

Syntax

p = capaplot(data,specs) p = capaplot(ax,data,specs) [p,h] = capaplot(___)

Description

p = capaplot(data,specs) estimates the mean and variance for the observations in the input vector data, and plots the pdf of the resulting normal distribution. The observations in data are assumed to be normally distributed. The output p is the probability that a new observation from the estimated distribution falls within the range specified by the two-element vector specs. The portion of the distribution between the lower and upper bounds specified in specs is shaded in the plot.

capaplot treats NaN values in data as missing, and ignores them.

p = capaplot(ax,data,specs) plots into the axes specified by ax instead of the current axes (gca). (since R2024a)

[p,h] = capaplot(___) additionally returns handles to the plot elements in h, using any of the input argument combinations in the previous syntaxes.

Examples

collapse all

Create a Process Capability Plot

Open Live Script

Randomly generate sample data from a normal process with a mean of 3 and a standard deviation of 0.005.

rng default; % For reproducibility
data = normrnd(3,0.005,100,1);

Compute capability indices if the process has an upper specification limit of 3.01 and a lower specification limit of 2.99.

S = capability(data,[2.99 3.01])

S = struct with fields:
       mu: 3.0006
    sigma: 0.0058
        P: 0.9129
       Pl: 0.0339
       Pu: 0.0532
       Cp: 0.5735
      Cpl: 0.6088
      Cpu: 0.5382
      Cpk: 0.5382

Visualize the specification and process widths.

capaplot(data,[2.99 3.01]);
grid on

Figure contains an axes object. The axes object with title Probability Between Limits = 0.91292 contains an object of type line. This object represents Normal Distribution.