| |
|
|
p.specht
| In el Statistik puede ser en pequeño Stückzahlen y Sample-Größen no con el Gauss-schen Normalverteilung trabajo, porque por el "Diskretisierung" gelten hier otro Gesetzlichkeiten.
Wenn una Packung con 100 Schrauben 5 schlechte enthält, y yo entnehme 20: Como wahrscheinlich es entonces, daß en esta 20 : Keine, una, zwei, drei, vier oder todos fünf fehlerhaft son? Hier hilft el Hypergeometrische Verteilung más, de el bekannten Formel "Gesuchte_Fälle / Überhaupt_gleichmögliche_Fälle" entwickelt wurde.
Anbei - sin jede Gewähr - una kleines Programa dazu. Lo enthält todavía no Prüfung, si el Eingabewerte sinnvoll sind; para rumspielen reicht lo aber.
Título de la ventana "Hypergeometrische Verteilung"
Ventana de Estilo 24:set("decimals",15):font 2
Ventana %maxx/10,0 - %maxx/2,%maxy-40
Proc BinCoeff :parámetros N&,k&:var h!=1
'Pascalsches Dreieck: Binärkoeffizient = ´N encima k ´
whileloop k&,1,-1:h!=h!*(N&-&Loop+1)/&Loop
endwhile:volver h!
ENDPROC
Proc HypGFormel :parámetros N&,P!,s&,x&'P(x) = (p x)*(q s-x)/(N s)
Case p!<0:p!=abs(p!/100):Case (p!>=0) and (p!<1):P!=int(p!*N&)
Volver BinCoeff(P!,x&)/BinCoeff(N&,s&) * BinCoeff(N&-P!,s&-x&)
ENDPROC
Proc Hyp0 :parámetros N&,P!,s&'Startwert x=0 el Rekursionsformel
'P!=Stückzahl | Wenn p=0..1 Wahrscheinlichkeitsfaktor, si <0: Prozentanteil
'Wahrscheinlichkeit, en un Sample s de N Elementen con P besonderen
'keines el besonderen Elemente para encontrar
caso p!<0:p!=abs(p!/100):Case (p!>=0) and (p!<1):P!=int(p!*N&)
var H!=1
var Q&=N&-P!
WhileLoop s&,1,-1
H!=H!*(Q&-&Loop+1)/(N&-&Loop+1)
Endwhile
Volver H!
ENDPROC
Proc HypGNext :parámetros N&,P!,s&,x&,f_lastx!'REKURSIONSFORMEL
Case p!<0:p!=abs(p!/100):Case (p!>=0) and (p!<1):P!=int(p!*N&)
caso (x&<0) or (x&>s&):volver 0
var Q&=N&-P!
volver f_lastx!*(P!-x&)/(x&+1)*(s&-x&)/(Q&-(s&-x&)+1)
ENDPROC
'Main
Declarar f!,Addf!,w$,N&,P!,s&,x&
Lup:
CLS
imprimir "\n Anzahl Grundgesamtheit: ";:input w$:caso w$>"":N&=val(w$)
if N&<1:N&=100:P!=5:s&=20:x&=0:CLS:imprimir "\n TESTWERTE 100,5,20":goto "skip":endif
imprimir " Anzahl Besondere Elem.: ";:input w$:caso w$>"":P!=val(w$)
imprimir " Sample-Größe: ";:input w$:caso w$>"":s&=val(w$)
imprimir " Besondere en el Sample: ";:input w$:caso w$>"":x&=val(w$)
Skip:
imprimir "\n Verteilungsdichte"
Imprimir " HyperGVtlg-Rekursionsformel: Check con BinCoeff-Formel:\n"
Whileloop 0,P!
if &Loop : f!=HypGNext(N&,P!,s&,&Loop-1,f!)
más : f!=Hyp0(N&,P!,s&)
endif
Imprimir " x =";&Loop;": P(x)= ";format$("%g",f!);
imprimir tab(44);HypGFormel(N&,P!,s&,&Loop)
Endwhile
Waitinput
cls
imprimir "\n Verteilungsfunktion dazu:\n"
Addf!=0
Whileloop 0,P!
if &Loop : f!=HypGNext(N&,P!,s&,&Loop-1,f!)
más : f!=Hyp0(N&,P!,s&)
endif
Addf!=Addf!+f!
Imprimir " x_max = ";&Loop;": P(x<=";&Loop;")= ";format$("%g",Addf!)
Endwhile
waitinput
Goto "Lup"
|
|
|
| XProfan 11Computer: Gerät, daß es in Mikrosekunden erlaubt, 50.000 Fehler zu machen, zB 'daß' statt 'das'... | 30.05.2021 ▲ |
|
|
|