Print

Please login or register.
1 Hour 1 Day 1 Week 1 Month Forever
Login with username, password and session length

[Helge A Bentsen]

A long time ago(possible in a galaxy far far away) someone on this board referenced a simulation software where you could visualize interference pattern from source spacing.

I remember it as pretty simple but it got the point across to people not that fluent in simulation software reading.
Does anybody remember what this software was/is?

[Scott Helmke]

You mean like Meyer's "MAPP" software? 

[Helge A Bentsen]

You mean like Meyer's "MAPP" software?

Sort of, just a lot simpler visually.

[David Sturzenbecher]

Sort of, just a lot simpler visually.

There is an iOS App called "ripple tank" that might be close to what you are thinking about.   

[Frank Koenig]

A long time ago(possible in a galaxy far far away) someone on this board referenced a simulation software where you could visualize interference pattern from source spacing.

I remember it as pretty simple but it got the point across to people not that fluent in simulation software reading.
Does anybody remember what this software was/is?

Helge,

I don't know if this solves your problem but I couldn't resist putting the entire source code of a program in a forum reply. This is in R (under Windows) and you would need to download and install the R interpreter -- easily done.

https://www.r-project.org/

The parameters are edited right into the code and are fairly self-explanatory. If this doesn't look like fun you can download free SPL modeling programs (2D and 3D) from Danley. But maybe I don't understand what you're looking for.

--Frank

#---------------------------------------------------------------------------------------------
rm(list = ls())
X <- 100         #length of prediction (ft)
Y <- 100         #width of prediction (ft)
res <- 4         #spatial resolution (points/ft)
M <- X * res
N <- Y * res
x <- (1:M - M/2) / res
y <- (1:N - N/2) / res
f <- 62.5         #frequency (Hz)
omega <- 2 * pi * f   #frequency (rad/s)
c <- 1130         #velocity (ft/s)
k <- omega / c      #wave number (1/ft)
I <- 4         #number of sources used in this run
S <- matrix(nrow = 10, ncol = 4)   #source parameters -- max of 10 sources allocated
S[1,1] <- 4.52      #x coord (ft)
S[1,2] <- 3         #y coord (ft)
S[1,3] <- 1         #sound pressure (relative)
S[1,4] <- .004      #delay (s)
S[2,1] <- 0
S[2,2] <- 1
S[2,3] <- 1
S[2,4] <- 0
S[3,1] <- 0
S[3,2] <- -1
S[3,3] <- 1
S[3,4] <- 0
S[4,1] <- 4.52
S[4,2] <- -3
S[4,3] <- 1
S[4,4] <- .004
d <- matrix(nrow = M, ncol = N)               #distance to source at each point(ft)
p <- matrix(data = complex(1, 0, 0), nrow = M, ncol = N)   #complex sound pressure at each point
for (i in 1:I){      #for each source
  for (m in 1:M){      #for each x value
    d[m, 1:N] <- ((x[m] - S[i,1])^2 + (y[1:N] - S[i, 2])^2)^.5
    p[m, 1:N] <- p[m, 1:N] + complex(1, modulus = S[i,3] / d[m, 1:N], argument = k * d[m, 1:N]  + S[i,4] * omega)
  }
}
windows(7, 7)      #plot magnitude
image(x, y, -20*log(Mod(p), 10), zlim = c(0, 42), col = rainbow(256, start = 0, end = .)
windows(7, 7)      #plot phase
image(x, y, Arg(p), zlim = c(-pi, pi), col = rainbow(256, start = 0, end = .)
windows(7,5)      #plot magnitude cross section
plot.default(x, -20*log(Mod(p[1:M, N/2]), 10), "l")
grid()
#cal <- matrix(nrow = 1, ncol =    #plot color test strip
#cal[1, 1:8] <- 6 * 0:7
#windows(2,7)
#image(1, 1:8, cal, zlim = c(0, 42), col = rainbow(256, start = 0, end = .)
#---------------------------------------------------------------------------------

I see the forum software helpfully replaces all occurrences of a the numeral 8 followed by a closing parenthesis with a smiley with sunglasses. Aren't computers wonderful?

[Tim McCulloch]

I see the forum software helpfully replaces all occurrences of a the numeral 8 followed by a closing parenthesis with a smiley with sunglasses. Aren't computers wonderful?

There is an option to "turn off BB Code" which will disable the smilies.

[Dave Garoutte]

Who allowed these geeks on here?

[Scott Holtzman]

Helge,

I don't know if this solves your problem but I couldn't resist putting the entire source code of a program in a forum reply. This is in R (under Windows) and you would need to download and install the R interpreter -- easily done.

https://www.r-project.org/

The parameters are edited right into the code and are fairly self-explanatory. If this doesn't look like fun you can download free SPL modeling programs (2D and 3D) from Danley. But maybe I don't understand what you're looking for.

--Frank

#---------------------------------------------------------------------------------------------
rm(list = ls())
X <- 100#length of prediction (ft)
Y <- 100#width of prediction (ft)
res <- 4#spatial resolution (points/ft)
M <- X * res
N <- Y * res
x <- (1:M - M/2) / res
y <- (1:N - N/2) / res
f <- 62.5#frequency (Hz)
omega <- 2 * pi * f#frequency (rad/s)
c <- 1130#velocity (ft/s)
k <- omega / c#wave number (1/ft)
I <- 4#number of sources used in this run
S <- matrix(nrow = 10, ncol = 4)#source parameters -- max of 10 sources allocated
S[1,1] <- 4.52#x coord (ft)
S[1,2] <- 3#y coord (ft)
S[1,3] <- 1#sound pressure (relative)
S[1,4] <- .004#delay (s)
S[2,1] <- 0
S[2,2] <- 1
S[2,3] <- 1
S[2,4] <- 0
S[3,1] <- 0
S[3,2] <- -1
S[3,3] <- 1
S[3,4] <- 0
S[4,1] <- 4.52
S[4,2] <- -3
S[4,3] <- 1
S[4,4] <- .004
d <- matrix(nrow = M, ncol = N)#distance to source at each point(ft)
p <- matrix(data = complex(1, 0, 0), nrow = M, ncol = N)#complex sound pressure at each point
for (i in 1:I){#for each source
  for (m in 1:M){#for each x value
    d[m, 1:N] <- ((x[m] - S[i,1])^2 + (y[1:N] - S[i, 2])^2)^.5
    p[m, 1:N] <- p[m, 1:N] + complex(1, modulus = S[i,3] / d[m, 1:N], argument = k * d[m, 1:N]  + S[i,4] * omega)
  }
}
windows(7, 7)#plot magnitude
image(x, y, -20*log(Mod(p), 10), zlim = c(0, 42), col = rainbow(256, start = 0, end = .)
windows(7, 7)#plot phase
image(x, y, Arg(p), zlim = c(-pi, pi), col = rainbow(256, start = 0, end = .)
windows(7,5)#plot magnitude cross section
plot.default(x, -20*log(Mod(p[1:M, N/2]), 10), "l")
grid()
#cal <- matrix(nrow = 1, ncol = #plot color test strip
#cal[1, 1:8] <- 6 * 0:7
#windows(2,7)
#image(1, 1:8, cal, zlim = c(0, 42), col = rainbow(256, start = 0, end = .)
#---------------------------------------------------------------------------------

I see the forum software helpfully replaces all occurrences of a the numeral 8 followed by a closing parenthesis with a smiley with sunglasses. Aren't computers wonderful?

Frank there is a code tag in BBCODE that won't parse anything and uses a fixed pitch font to preserve formatting.

Sent from my VS996 using Tapatalk

[Luke Geis]

I think he is thinking of Ray-end? which as of this moment I believe is no longer available?

[Scott Holtzman]

Frank there is a code tag in BBCODE that won't parse anything and uses a fixed pitch font to preserve formatting.

Sent from my VS996 using Tapatalk

I should have posted code

Code: [Select]

[code]test[/code]

Sent from my VS996 using Tapatalk

[Andrew Broughton]

https://forums.prosoundweb.com/index.php/topic,119435.msg1119306.html#msg1119306

[Frank Koenig]

https://forums.prosoundweb.com/index.php/topic,119435.msg1119306.html#msg1119306

Right. The Danley software is GPA. I believe Danley retained Sabastian to write and maintain it. He was present at the Danley demo room at INFOCOM a few years ago. -F

[Frank Koenig]

I should have posted code

Code: [Select]

[code]test[/code]

Sent from my VS996 using Tapatalk

Thanks for the tips on killing the smileys. As it turns out if you cut and paste the code from the browser, as you would need to to actually run it, everything comes out normal. The Smileys only render as such when viewed in the browser. I'm always amused by stuff like that.

As for the program itself, I threw it together about 6 years ago mostly as a fun R exercise. I played with it a little yesterday and, while fundamentally working, there are some issues that I may try to fix. There are, of course, endless embellishments that could be made. I mainly put it up here to demonstrate how a very few lines of code sometimes are sufficient to model a complex physical situation and give you the ability to interact, experiment, and learn. Maybe I'll get a few R converts, although, frankly, Python might be a better bet these days.

--Frank

[Mal Brown]

R is a pretty wonderful tool if you are a data geek.  In the last iteration of my day gig I was analyzing 2500 samples per second of aeronautics telemetry.  R was a godsend...  I’ve moved onto a much tidier data set so my R skills have slipped...  now it”s SQL and Ms-PowerBI. Which is a pretty fun toolset as well.

[Frank Koenig]

The sample grid is now symmetrical about the origin and a sample occurs exactly at (0,0). The colors are prettier and I did some testing. The example output below is for a two point-source end-fire array set for max cancellation at 100 Hz, evaluated at 100 Hz.





I'll put up the new source if anyone is interested. -F

[Scott Holtzman]

R is a pretty wonderful tool if you are a data geek.  In the last iteration of my day gig I was analyzing 2500 samples per second of aeronautics telemetry.  R was a godsend...  I’ve moved onto a much tidier data set so my R skills have slipped...  now it”s SQL and Ms-PowerBI. Which is a pretty fun toolset as well.

With 7 cloud data centers we collect a ton of data too

SQLhas a ton of limitations.  Mal Brown have you tried a time series database (TSDB)?  I am having great success with Influx.  The command line too will be familiar to anyone with SQL experience.

We are starting to us more modern tools to access including some cool elasticsearch and AI analysis of data sets is the next frontier.

Grafana is a great visualizer.

Since moving to a data driven network operations center we resolve 97% resolution prior to service impairment.  Even better more than half of troubles are resolved by automated incident response software (Pagerduty and Swimlane).


Have fun.

Sent from my VS996 using Tapatalk

[ProSoundWeb Community]