Exemplo n.º 1
0
def BisectionMethod(x0,x1,fcn):
    #stopping Criterion
    tol = 0.0000000000001
    fx = lambda x: fcn(x)[1] - E 
    while abs(x0-x1)> tol:   
        d = (x0+x1)/2.0
        if (fx(x0)*fx(d)) < 0:
            x1 = d
        else:
            x0 = d 
    if sun_angle.sun_azel(d)[0]<200:
        return BisectionMethod(d+tol,24,fcn)    
    return x0
Exemplo n.º 2
0
import matplotlib.pylab as plt
plt.ion()
import sun_angle


# Bisection Method
def BisectionMethod(x0,x1,fcn):
    #stopping Criterion
    tol = 0.0000000000001
    fx = lambda x: fcn(x)[1] - E 
    while abs(x0-x1)> tol:   
        d = (x0+x1)/2.0
        if (fx(x0)*fx(d)) < 0:
            x1 = d
        else:
            x0 = d 
    if sun_angle.sun_azel(d)[0]<200:
        return BisectionMethod(d+tol,24,fcn)    
    return x0


E = 17.23528526
print 'Elevation angle is', E, 'degrees'


h= BisectionMethod(0,24,sun_angle.sun_azel)
print "Hours: ", int(math.modf(h)[1]),', Minutes: ', int(round(math.modf(h)[0]*60))

azel = sun_angle.sun_azel(h)   # h=18 is 6:00pm
print 'Azimuth A =', azel[0], ', Elevation E =', azel[1]