Hi all!

I am currently playing around and implementing a way to design and test my own skyline diffusers based on this research paper: („The design and application of modular acoustic diffusing elements“ – R. Walker).

I have the following code already (https://github.com/Ohmagar/Skyline_Diffuser/blob/main/Untitled%2015.py) which creates an array of a Primitive Root Diffuser, turns this array into a matrix with given factors (as outlined in the research paper), calculates the individual well depths and the resulting phase changes (convolution / impulse response) of the Diffuser.

Now my questions:

How do I calculate (and plot?) the spatial response?

And how do I calculate how well the diffuser in itself functions? Preferably with a resulting plot in which I can „read“ the distribution, possible phase cancellations happening because of the diffuser, etc?

Any kind of approaches I tried (some with the help of ChatGPT) ended up as me being more confused than before so now I turn to the immediate experts! I know MATLAB well enough that I could translate code from MATLAB into python on my own if that’s your choice of language!

Help is greatly appreciated and have a wonderful Christmas!