Soubor:MATLABSChartForPairedXBarAndSChart.png
Obsah stránky není podporován v jiných jazycích.
Z Wikiknih
MATLABSChartForPairedXBarAndSChart.png (560 × 420 pixelů, velikost souboru: 4 KB, MIME typ: image/png)
Tento soubor pochází z Wikimedia Commons a mohou ho používat ostatní projekty. Níže jsou zobrazeny informace, které obsahuje jeho tamější stránka s popisem souboru.
Popis
PopisMATLABSChartForPairedXBarAndSChart.png |
English: A en:MATLAB-generated s en:Control chart for a process that experienced a 1.5σ drift starting at midnight. This chart is paired with an xbar chart to form an en:Xbar and s chart. |
Datum | |
Zdroj | Vlastní dílo |
Autor | DanielPenfield |
Licence
Já, držitel autorských práv k tomuto dílu, ho tímto zveřejňuji za podmínek následující licence:
Tento soubor podléhá licenci Creative Commons Uveďte autora-Zachovejte licenci 3.0 Unported
- Dílo smíte:
- šířit – kopírovat, distribuovat a sdělovat veřejnosti
- upravovat – pozměňovat, doplňovat, využívat celé nebo částečně v jiných dílech
- Za těchto podmínek:
- uveďte autora – Máte povinnost uvést autorství, poskytnout odkaz na licenci a uvést, pokud jste provedli změny. Toho můžete docílit jakýmkoli rozumným způsobem, avšak ne způsobem naznačujícím, že by poskytovatel licence schvaloval nebo podporoval vás nebo vaše užití díla.
- zachovejte licenci – Pokud tento materiál jakkoliv upravíte, přepracujete nebo použijete ve svém díle, musíte své příspěvky šířit pod stejnou nebo slučitelnou licencí jako originál.
Source code
en:Perl
#!/usr/bin/perl -w
use strict;
use Math::Random;
my %shiftSchedule = (
"first" => { "start" => 6.00, "end" => 14.00 },
"second" => { "start" => 14.00, "end" => 22.00 },
"third" => { "start" => 22.00, "end" => 6.00 }
);
my $shift = "third"; # shift to monitor
my $inspectionRate = 1 / 2; # every 1/2 hour
my $drift = 1.5; # drift to simulate
my $m = 25; # samples in control chart setup
my $n = 12; # observations per sample
my $target = 100.0; # quality characteristic target
my $hour;
my $i;
my $j;
my $minute;
my $observation;
my $setupM = $m;
#
# simulate control chart setup
#
open(SETUPCSV, ">xs_setup.csv") || die "! can't open \"xs_setup.csv\" ($!)\n";
for ($i = 1; $i <= $m; $i++) {
for ($j = 0; $j < $n; $j++) {
$observation = $target + random_normal();
if ($j > 0) {
print SETUPCSV ",";
}
printf SETUPCSV "%7.3f", $observation;
}
printf SETUPCSV "\r\n";
}
close(SETUPCSV);
#
# simulate control chart monitoring
#
open(MONITORINGCSV, ">xs_monitoring.csv") || die "! can't open \"xs_monitoring.csv\" ($!)\n";
$m = $shiftSchedule{$shift}{"end"} - $shiftSchedule{$shift}{"start"};
if ($m < 0) {
$m += 24;
}
$m /= $inspectionRate;
for ($i = 1; $i <= $m; $i++) {
$hour = int($i * $inspectionRate + $shiftSchedule{$shift}{"start"});
if ($hour >= 24) {
$hour -= 24;
}
$minute = ($i & 0x1) ? (60 * $inspectionRate) : 0;
printf MONITORINGCSV "'%d:%02d'", $hour, $minute;
for ($j = 0; $j < $n; $j++) {
$observation = $target + random_normal();
if ($i >= (0.25 * $m)) {
if ($i < (0.75 * $m)) {
$observation += ($drift / (0.5 * $m)) * ($i - (0.25 * $m));
} else {
$observation += $drift;
}
}
printf MONITORINGCSV ",%7.3f", $observation;
}
printf MONITORINGCSV "\r\n";
}
close(MONITORINGCSV);
en:MATLAB
%
% display an s control chart in MATLAB
%
clear
%
% Phase I
%
% compute the control chart center line and control limits based on a
% process that is simulated to be in a state of statistical control
%
setupobservations = csvread('xs_setup.csv');
setupstats = controlchart(setupobservations, 'charttype', 's');
%
% Phase II
%
% read in the process observations representing the monitoring phase
%
observations = importdata('xs_monitoring.csv');
%
% first column is the time of the observation (24 hour clock)
%
halfhourlylabel = observations.rowheaders;
%
% second column consists of the observations (counts of
% nonconformances per rational subgroup)
%
monitoringobservations = observations.data;
%
% just display labels on the "on the hour" ticks
%
emptylabel = cell(size(monitoringobservations,1) - size(halfhourlylabel,1), 1);
emptylabel(:) = {''};
hourlylabel = vertcat(halfhourlylabel(2:2:end), emptylabel);
%
% plot the control chart for the monitoring phase observations
%
monitoringstats = controlchart(monitoringobservations, ...
'charttype', 's', ...
'label', halfhourlylabel.', ...
'mean', setupstats.mu, ...
'sigma', setupstats.sigma);
title('s chart for quality characteristic XXX')
xlabel('Sample')
ylabel('Standard deviation (units)')
%
% the labels supplied to controlchart() only appear when the user
% selects a plotted point with her mouse--we have to explicitly
% set labels in the X axis if we want them
%
set(gca,'XTickLabel', hourlylabel)
Položky vyobrazené v tomto souboru
zobrazuje
Nějaká hodnota bez položky na Wikidatech
22. 6. 2013
Historie souboru
Kliknutím na datum a čas se zobrazí tehdejší verze souboru.
Datum a čas | Náhled | Rozměry | Uživatel | Komentář | |
---|---|---|---|---|---|
současná | 22. 6. 2013, 14:05 | 560 × 420 (4 KB) | DanielPenfield | User created page with UploadWizard |
Využití souboru
Tento soubor nepoužívá žádná stránka.
Globální využití souboru
Tento soubor využívají následující wiki:
- Využití na en.wikibooks.org
Metadata
Tento soubor obsahuje dodatečné informace, poskytnuté zřejmě digitálním fotoaparátem nebo scannerem, kterým byl pořízen. Pokud byl soubor od té doby změněn, některé údaje mohou být neplatné.
Autor | |
---|---|
Držitel autorských práv |
|
Název obrázku |
|
Rozlišení obrázku na šířku | 37,8 dpc |
Rozlišení obrázku na výšku | 37,8 dpc |
Datum a čas digitalizace |
|
Použitý software |
|