bsp_5_2_1f.m

% Beispiel für Kap 5.2.1f
%% Beispiel Schraubendurchmesser
al = 0.05;
n = 10;
sigma = 0.05;
mu = 4.95;

fprintf('\nTest auf my bei bekanntem sigma\n');
mu0 = 5;
rng(42);  % For reproducibility
xi = normrnd(mu, sigma, 1, n)
xq = mean(xi)
Tx = sqrt(n)/sigma*abs(xq - 5)
pW = 2*(1-normcdf(Tx));
fprintf('\npW = %8.6f\n', pW)

fprintf('\nTest auf my bei unbekanntem sigma\n');
S = std(xi)
c = tinv(1-al/2, n-1)
Tx = sqrt(n)/S*abs(xq - 5)

fprintf('\nTest auf sigma\n');
c1 = chi2inv(al/2, n-1)
c2 = chi2inv(1-al/2, n-1)
Tx = (n-1)/sigma^2 * S^2

fprintf('\nTest auf gleiche Schraubenlänge\n');
n = 10;
m = 12;
mux = 4.95;
muy = 5.00;
sigma = 0.05;

al = 0.05;
xi = normrnd(mux, sigma, 1, n)
yi = normrnd(muy, sigma, 1, m)
c = norminv(1-al/2)
xq = mean(xi)
yq = mean(yi)
Tx = (xq - yq)/(sigma*sqrt(1/n + 1/m))

fprintf('\nTest auf gleiche Schraubenstreuung\n');
al = 0.01;
n = 10;
m = 10;
mux = 4.95;
muy = 5.00;
sigmax = 0.05;
sigmay = 0.03;
xi = normrnd(mux, sigmax, 1, n)
yi = normrnd(muy, sigmay, 1, m)
c1 = finv(al/2, n-1, m-1)
c2 = finv(1-al/2, n-1, m-1)
sx = std(xi)
sy = std(yi)
Tx = (sx/sy)^2

% bei größerem alpha
al0 = 0.05; 
c1a = finv(al0/2, n-1, m-1);
c2a = finv(1-al0/2, n-1, m-1);
fprintf('\nc-Werte bei al=0.05: %6.4f %6.4f\n', c1a, c2a)