theta acos(costheta_rnd(n,freq[1/Hz],2,seed)) ""
phi phi_rnd(n,freq[1/Hz],1,seed) ""
phase phase_rnd(n,freq[1/Hz],3,seed) ""
k0 2*pi*freq/c0 "Wave number in air"
kx cos(theta)*k0 "Partial wave vector, x-component"
ky sin(theta)*cos(phi)*k0 "Partial wave vector, y-component"
kz sin(theta)*sin(phi)*k0 "Partial wave vector, z-component"
px_room A/sqrt(2*N)*sum(exp(i*phase)*exp(-i*(kx*x+ky*y+kz*z)),n,1,N) "Room diffuse pressure field, waves in positive x direction"
vx_room A/(acpr.omega*rho0)/sqrt(2*N)*sum(kx*exp(i*phase)*exp(-i*(kx*x+ky*y+kz*z)),n,1,N) "Room diffuse particle velocity field, waves in positive x direction"
p_room_sq 2*abs(px_room)^2 "Room diffuse field, total squared pressure"
p_refl A/sqrt(2*N)*sum(exp(i*phase)*exp(-i*(-kx*x+ky*y+kz*z)),n,1,N) "Pressure field reflected from the concrete surface"
p_wall px_room+p_refl "Total wall pressure (mechanical load)"