In [75]:

    

import numpy as np
import matplotlib.pyplot as plt


    

In [80]:

    

t =  np.linspace(0.0, 2, num=88001)


    

In [99]:

    

t60m = 6.90776;
tau = 0.9/t60m;
v = 1 - exp(-t/tau);
vt = next(vt for vt in v if vt >= 1-0.001);
index = nonzero(v == vt)[0][0]
(v[index], t[index], index)


    

    
Out[99]:





(0.99990000062946705, 1.3333333333333333, 35200)

In [96]:

    

plot(t,v)
grid('on')


    

In [97]:

    

tau = 0.07;
v = 1 - exp(-t/tau);
(v[4999], t[4999])


    

    
Out[97]:





(0.80270435674990726, 0.11361363636363636)

In [98]:

    

plot(v[0:4999])
grid('on')


    

In [86]:

    

ti = -(4999.0/44100.0)/(np.log(1.0-0.791136)) #500
ti/500


    

    
Out[86]:





0.00014476475183363877

In [87]:

    

ti = -(4999.0/44100.0)/(np.log(1.0-0.926475)) #300
ti/300


    

    
Out[87]:





0.0001447642045839103

In [88]:

    

ti = -(4999.0/44100.0)/(np.log(1.0-0.980063)) #200
ti/200


    

    
Out[88]:





0.00014476482261176465

In [89]:

    

1/0.00014476475183363877


    

    
Out[89]:





6907.758880070359

In [64]:

    

s = 600;
t4 = (4999.0/44100.0)
tau1 = s/6907.0;
vti = 1 - exp(-t4/tau1);
vti


    

    
Out[64]:





0.72880484580304972

In [65]:

    

tau1


    

    
Out[65]:





0.0868683943825105

In [66]:

    

0.6/6.907


    

    
Out[66]:





0.0868683943825105

In [ ]: