WEB开发网      婵犵數濞€濞佳囧磹婵犳艾鐤炬い鎰堕檮閸嬬喐銇勯弽銊с€掗梻鍕閺岋箑螣娓氼垱笑闂佽姘﹂褔婀佸┑鐘诧工妤犲憡绂嶉崜褏纾奸弶鍫涘妼缁楁岸鏌熷畡鐗堝殗闁诡喒鏅犲畷褰掝敃閵堝棙顔忔繝鐢靛仦閸ㄥ爼骞愰幘顔肩;闁规崘绉ぐ鎺撳亹闁绘垶锕╁Λ鍕⒑閹肩偛濡奸悗娑掓櫇缁顓兼径妯绘櫇闂佹寧绻傞弻濠囨晝閸屾稓鍘甸柣搴㈢⊕閿氶柣蹇ョ稻缁绘繃绻濋崘銊т紝闂佽鍨伴崯鏉戠暦閻旂⒈鏁傞柛鈾€鏅欑槐妯衡攽閻愬樊鍤熷┑顔藉劤铻為柛鏇ㄥ墯閸欏繘鏌嶉崫鍕櫣缂佲偓婢跺绠鹃柟瀛樼箘閿涘秵顨ラ悙顏勭伈闁诡喖缍婂畷鎯邦槻婵℃彃顭烽弻娑㈠Ω閵夈儺鍔夌紓浣稿€哥粔褰掑极閹剧粯鏅搁柨鐕傛嫹 ---闂傚倷鐒︾€笛兠洪埡鍛闁跨噦鎷�
开发学院软件开发Python 正态分布概率计算 阅读

正态分布概率计算

 2012-05-13 12:21:13 来源:WEB开发网 闂傚倷绶氬ḿ褍螞閹绢喖绠柨鐕傛嫹闂傚倷绀侀幉锟犲垂閻㈠灚宕查柟鎵閸庡秵銇勯幒鎴濃偓鐢稿磻閹炬枼妲堟繛鍡楃С濞岊亞绱撻崒姘扁枌闁瑰嚖鎷�婵犵數濮幏鍐川椤撴繄鎹曢梻渚€娼уú銈吤洪妸鈺佺劦妞ゆ帊鑳堕埊鏇㈡煏閸モ晛浠х紒杈╁仱閺佹捇鏁撻敓锟�闂傚倷绶氬ḿ褍螞閹绢喖绠柨鐕傛嫹  闂傚倷鑳舵灙缂佺粯顨呴埢宥夊即閵忕姵鐎梺缁樺姇閹碱偆鐥閺屾洘绻濊箛鏇犵獥闂佺厧澹婃禍婊堚€旈崘顏佸亾閿濆骸澧ù鐘欏嫮绠鹃柛顐ゅ枎閻忋儲銇勯弴妯哄姦妞ゃ垺鐟╅幃鈩冩償閵囧府鎷�
核心提示: def st_norm(u): '''标准正态分布''' import math x=abs(u)/math.sqrt(2) T=(0.0705230784,0.0422820123,0.0092705272, 0.000152014

 def st_norm(u):

    '''标准正态分布'''
    import math
    x=abs(u)/math.sqrt(2)
    T=(0.0705230784,0.0422820123,0.0092705272,
       0.0001520143,0.0002765672,0.0000430638)
    E=1-pow((1+sum([a*pow(x,(i+1))
                    for i,a in enumerate(T)])),-16)
    p=0.5-0.5*E if u<0 else 0.5+0.5*E
    return(p)
 
def norm(a,sigma,x):
    '''一般正态分布'''
    u=(x-a)/sigma
    return(st_norm(u))
 
while 1:
    '''输入一个数时默认为标准正态分布
    输入三个数(空格隔开)时分别为期望、方差、x
    输入 stop 停止'''
    S=input('please input the parameters:\n')
    if S=='stop':break
    try:
        L=[float(s) for s in S.split()]
    except:
        print('Input error!')
        continue
    if len(L)==1:
        print('f(x)=%.5f'%st_norm(L[0]))
    elif len(L)==3:
        print('f(x)=%.5f'%norm(L[0],L[1],L[2]))
    else:
        print('Input error!')
 

Tags:正态分布 概率 计算

编辑录入:爽爽 [复制链接] [打 印]
赞助商链接