tone_est.py 文件源码

def est_tone_phase(sdata,a,f,sr):
    delta = 120
    min_ang = 0.5
    p = 0
    phase = numpy.zeros(3)
    phase[0] = p+(-delta/180.0)*pi
    phase[1] = p
    phase[2] = p+(delta/180.0)*pi
    while (delta > min_ang):
        (i,p) = find_min_phase(sdata,a,f,sr,phase)
        delta = delta/2.0
        phase[0] = p+(-delta/180.0)*pi
        phase[1] = p
        phase[2] = p+(delta/180.0)*pi
        #print "p = ",(180.0*p/pi),'delta = ',delta

    min_phase = p
    #print "min for phase sweep is ",sum_min,' at offset ',min_index
    return min_phase