a guest Jan 12th, 2017 74 in 8 hours
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
- def W_rad(msq_vel=msq_flr_plate, sigma=radiationRatio(), p=p0, c=c0, S=flr_area2):
- """flr_area2 is the section of the floor with lenght 10m"
- return p*c*S*sigma*msq_vel
- def SPL_r(W, r, S=flr_area2, p=p0, c=c0):
- """returns the SPL by a source in dB at a distance r, given a known input
- sound power and the area it's vibrating from, W/S = I"""
- p_sqr_rms = p*c*W/S
- Lp = 10*np.log10(p_sqr_rms/(spl_ref**2))
- return Lp
RAW Paste Data