<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">One possible answer to your final observation is that your model is a model of a physical phenomena and if it is a good model it should exhibit much of the same behavior that we see in the physical world. The Fourier analysis I posed should yield results close to what we observe for the continuous case which should be echoed in a good computational model.<div class=""><br class=""></div><div class="">If you have any of my graphics books, they all have an appendix on sampling and aliasing which shows how the sinc and rect functions arise as fundamental to sampling/reconstruction. I don’t use Fourier analysis there. </div><div class=""><br class=""></div><div class="">Here’s some further connections and a tie to Owen.</div><div class=""><br class=""></div><div class="">I joined the EE Dept U Rochester in 1973. I was doing a lot digital image processing then. UR had one of the two Institutes of Optics (the other is at Arizona) and I started working with Brian Thompson, the head of the Institute, on hybrid image processing which got me to learn some about Fourier optics. We taught a joint Optics/EE course on hybrid processing. The optics students didn’t know anything about digital image processing and the EE students knew almost nothing about Fourier analysis in the complex domain. I gave a final exam to the optics students and Brian gave a final to the EE students. We also taught a number of short courses on hybrid processing. We also worked with the image processing group at Xerox Research in Webster. At the time, Owen was in that facility and we worked with the same people but never met. That’s really coming full circle.</div><div class=""><br class=""></div><div class="">A couple of comments about the sinc function. It comes up in sampling and in the continuous world whenever you have a window. There are many consequences of it having negative side lobes. In signal processing, it means that you can’t make a perfect low-pass filter. In optics it comes up in lens design because of the finite size of a lens. Alvy Ray glossed over it in his pixel talk. When you have discrete value and want to produce a finite pixel on an analog device, the value modulates a beam that covers a finite area. Sampling theory tells us that this beam should have the shape of a sinc function. But this is not physically possible because the negative side lobes of the sinc would require the beam to extract energy from the projection surface. So in physical systems such as CRTs, the beam has a Gaussian profile which has no negative values. In digital systems, the equivalent is a profile such as the one Alvy Ray showed that has a finite extent and can be described by polygons. Similar results hold for designing lens coatings, In all cases, these designs are trying minimize the errors, which can be seen in the frequency domain, caused by our inability to a perfect reconstruction.</div><div class=""><br class=""></div><div class="">Ed<br class=""><div class="">
<div>__________<br class=""><br class="">Ed Angel<br class=""><br class="">Founding Director, Art, Research, Technology and Science Laboratory (ARTS Lab)<br class="">Professor Emeritus of Computer Science, University of New Mexico<br class=""><br class="">1017 Sierra Pinon<br class="">Santa Fe, NM 87501<br class="">505-984-0136 (home)<span class="Apple-tab-span" style="white-space: pre;"> </span> <span class="Apple-tab-span" style="white-space: pre;"> </span><a href="mailto:edward.angel@gmail.com" class="">edward.angel@gmail.com</a><br class="">505-453-4944 (cell) <span class="Apple-tab-span" style="white-space: pre;"> </span><span class="Apple-tab-span" style="white-space: pre;"> </span>http://www.cs.unm.edu/~angel<br class=""></div>
</div>
<div><br class=""><blockquote type="cite" class=""><div class="">On Aug 14, 2021, at 9:56 PM, Stephen Guerin <<a href="mailto:stephen.guerin@simtable.com" class="">stephen.guerin@simtable.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div dir="ltr" class=""> Ed, <br class=""><br class="">Yes, that's how I'm seeing it. <br class=""><br class="">For others, Ed's Step function is what I was calling Rect pulse function which is Fourier dual of the Sinc function. Attached is an .mp4 recording of Ed description of the relationship of gap width, frequency, and the amount of spread of the Sinc function which is the diffusion pattern observed. I recorded the mp4 from the illustration/animation" <a href="https://www.olympus-lifescience.com/en/microscope-resource/primer/java/diffraction/" target="_blank" class="">I linked to earlier</a>.<br class=""><br class="">And I think it's pretty cool to think of the gap as a sampler. I suspect this is a well-known idea in optics/physics and old hat to many of you but it's exciting to come onto ideas like this for oneself :-) I can almost hear John Zingale saying, "of course, it's just the<a href="https://en.wikipedia.org/wiki/Convolution_theorem" class=""> convolution theorem</a> applied to a square wave". In the past I would nod and frankly, my eyes glaze if I can't ground it in a microscopic understanding that guides my intuition.<br class=""><br class="">Or if given this amazingly deep statement I came across as I'm searching for connecting sampling and diffraction - "<i class="">the diffraction pattern of an object is the Fourier Transform of the object</i>"<a href="http://www.sci.sdsu.edu/TFrey/Bio750/FourierTransforms.html" class=""> from here</a>, it finally makes sense to me. <br class=""><br class="">And I can practically hear Steve Smith and our dear and late Fred Untersher calling out, "that's how we've been describing holograms to you for 20 years and you always nodded like you understood". </div><div dir="ltr" class=""><br class=""></div><div dir="ltr" class="">Or Ed who brought Pradeep Sen into our world saying what do you think I was showing you with <a href="https://graphics.stanford.edu/papers/dual_photography/" class="">Dual Photography</a>, you idiot?<br class=""><br class="">And Alvy Ray Smith, again Ed bringing into our office, saying that's<a href="https://youtu.be/dvHDXUV7hmQ" class=""> what a Pixel is</a>! it's not a little square nor gaussian point sample, it's Kotelnikov Sampling (Nyquist-Shannon Sampling Theorem),<br class=""><br class="">or potentially worse is Eric Smith and Roger Critchlow shaking their heads saying "you're just confused and making connections that aren't there". :-)</div><div dir="ltr" class=""><br class=""></div><div dir="ltr" class="">--------------------------</div><div dir="ltr" class=""><br class=""></div><div dir="ltr" class="">Now even after having said this, I *still* want to know how the diffraction is happening using only the interaction rules in the model. Obviously, there are no Sinc or Rect functions in the code, nor Fourier transforms explicitly coded. All these wonderful explanations above are emergent properties from the model I would call a macroscopic explanation and description. If nothing else perhaps I learn a better phrase for the level of explanation I'm asking for when you trace an algorithm and understand where the emergent property comes from. (BTW, I think I have a micro answer and will put it in my response to Alex).<br class=""><br class="">-S<br class=""><br class=""><div class=""><div dir="ltr" data-smartmail="gmail_signature" class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class="">_______________________________________________________________________<br class=""><a href="mailto:stephen.guerin@simtable.com" target="_blank" class="">Stephen.Guerin@Simtable.com</a><div class="">CEO, Simtable <a href="http://www.simtable.com/" target="_blank" class="">http://www.simtable.com</a><br class=""><div class="">1600 Lena St #D1, Santa Fe, NM 87505<div class=""><div class="">office: (505)995-0206 <span style="font-size:12.8px" class="">mobile: (505)577-5828</span></div><div class=""><span style="font-size:12.8px" class="">twitter: @simtable</span></div><div class=""><span style="font-size:12.8px" class=""><a href="http://zoom.com/j/5055775828" target="_blank" class="">z</a><a href="http://oom.simtable.com/" target="_blank" class="">oom.simtable.com</a></span></div><div class=""></div></div></div></div></div></div></div></div></div></div></div><br class=""></div><br class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Aug 14, 2021 at 3:57 PM Angel Edward <<a href="mailto:edward.angel@gmail.com" target="_blank" class="">edward.angel@gmail.com</a>> wrote:<br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class=""><div class="">I hope someone can check out the analysis below. </div><div class=""><br class=""></div>If you look at the gap as a sampler, you can do the following analysis using Fourier methods:<div class=""><br class=""></div><div class="">A gap is a window on a continuous function. A perfect gap is a step function multiplying the continuous function. </div><div class=""><br class=""></div><div class="">In the Fourier domain, the Fourier transform of the continuous function on the input side of the gap is convolved with the Fourier transform of gap (the step function).</div><div class=""><br class=""></div><div class="">The Fourier transform of a step function is a sinc (sin(ax)/(ax)) function.</div><div class=""><br class=""></div><div class="">The width of the main lobe of the sinc is inversely proportional to the width of the gap.</div><div class=""><br class=""></div><div class="">Consequently, the smaller the width of the gap, the more a given frequency is distorted because the sinc is wider. Convolution applies the sinc at each frequency of the input function.</div><div class=""><br class=""></div><div class="">I think it gets more complicated when we add in sampling. If we take a number of samples that is proportional to the width of the gap, then as we make the gap smaller there are fewer samples, hence more reconstruction issues which is the second, often overlooked, part of the sampling theorem.</div><div class=""><br class=""></div><div class="">In the limit as the gap goes to zero width, there is no distortion to the continuous function but in the digital world you could have only a single sample.</div><div class=""><br class=""></div><div class="">Ed<br class=""><div class="">
<div class="">__________<br class=""><br class="">Ed Angel<br class=""><br class="">Founding Director, Art, Research, Technology and Science Laboratory (ARTS Lab)<br class="">Professor Emeritus of Computer Science, University of New Mexico<br class=""><br class="">1017 Sierra Pinon<br class="">Santa Fe, NM 87501<br class="">505-984-0136 (home)<span style="white-space:pre-wrap" class=""> </span> <span style="white-space:pre-wrap" class=""> </span><a href="mailto:edward.angel@gmail.com" target="_blank" class="">edward.angel@gmail.com</a><br class="">505-453-4944 (cell) <span style="white-space:pre-wrap" class=""> </span><span style="white-space:pre-wrap" class=""> </span><a href="http://www.cs.unm.edu/~angel" target="_blank" class="">http://www.cs.unm.edu/~angel</a><br class=""></div>
</div>
<div class=""><br class=""><blockquote type="cite" class=""><div class=""><br class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Aug 14, 2021 at 10:17 AM Stephen Guerin <<a href="mailto:stephen.guerin@simtable.com" target="_blank" class="">stephen.guerin@simtable.com</a>> wrote:<br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr" class=""><div dir="ltr" class="">At yesterday's Virtual Friam I asked a question on diffraction and said I would send more background.<br class=""><br class="">The gist of my question is: <br class=""><br class=""><b class="">Even though I completely understand the micro-level rules that generate diffraction in the wave model described below, I still don't have an intuition **how** the gaps in an obstacle have the emergent effect of diffracting waves when wavelengths >= gap width. Can anyone help?</b><br class=""><br class=""><br class="">Background:<br class="">The question arose from my mentoring NM School for the Arts high school students in the <a href="http://nmsupercomputingchallenge.org/" target="_blank" class="">NM Supercomputing Challenge</a> where the students simulated spatial acoustics by appropriating Saint-Venant equations used for shallow water waves to instead model acoustic pressure waves. We wrote a Netogo agent-based model with Python extension for reading / writing the sound files and simulating spatial acoustics.</div><div dir="ltr" class=""><br class=""><span id="m_-6704228213763452108gmail-m_-1165450971226445812cid:ii_ksbvdqhk0" class=""><image.png></span><br class=""><br class=""></div><div dir="ltr" class=""><br class="">The students explored the effects of different room configurations on acoustics. <br class=""><br class="">One configuration of interest was a wall gap illustrated below in the top right under Madelyn's video below. The wall gap is hard to see on right side.<br class=""><br class=""><span id="m_-6704228213763452108gmail-m_-1165450971226445812cid:ii_ksbz1gwy3" class=""><image.png></span><br class=""><br class="">They simulated microphones in Netlogo by recording amplitudes at a patch (red dot below in top-right visualization of room) and simulated speakers (hard-to-see blue dot below red dot on other side of wall) by driving amplitudes at a patch from the time series of amplitudes in .wav files (recordings of a singer and viola performance). They could hear, and through Fourier analysis, see the gap acting as a low-pass filter on the acoustic signal. ie, only the low frequencies were "bending" around the wall to reach the microphone. </div><div dir="ltr" class=""><br class=""></div><div dir="ltr" class="">You can see and listen to this effect and the spectrogram visualization <a href="https://youtu.be/61p97NWJiQ8?t=2117" target="_blank" class="">at time 33:11 in their presentation</a>.<br class=""><br class=""><span id="m_-6704228213763452108gmail-m_-1165450971226445812cid:ii_ksbx6zhp1" class=""><image.png></span><br class=""><br class="">It took me a few weeks after their presentation in the NM Supercomputing Challenge - they got second place - to connect the low pass filter behavior to the concept of diffraction. Had this been a light model and I saw the rainbow effects I would have clued in much faster. Their presentation was a month after finals and they <a href="https://youtu.be/61p97NWJiQ8?t=2761" target="_blank" class="">added this epilogue in the presentation above to identify the effect as diffraction.</a><br class=""><br class="">Their presentation included this <a href="https://youtu.be/BH0NfVUTWG4" target="_blank" class="">physical wavepool video demonstration</a> which was helpful to me to begin to understand the diffraction relationship with frequency and gap width.<br class=""><br class="">Note: my question is not about "describing" the behavior with macroscopic equations or geometric models but fundamentally how does the gap become a point source ala Huygens Principle at the micro-level of the patches interacting with the emergent waves. To help with the distinction, I consider t<a href="https://www.olympus-lifescience.com/en/microscope-resource/primer/java/diffraction/" target="_blank" class="">his interactive model </a> a great macroscopic description of the phenomenon that nicely illustrates the relationship of frequency and gap width but doesn't help me interpret the micro-level interactions giving rise to the diffraction effect in our simple shallow-water model.<br class=""><br class="">The students describe the details of the shallow water model at <a href="https://youtu.be/61p97NWJiQ8?t=870" target="_blank" class="">this point in their presentation</a>:<br class=""><span id="m_-6704228213763452108gmail-m_-1165450971226445812cid:ii_ksbzby8c5" class=""><image.png></span><br class=""><br class=""><br class="">Here is my <a href="https://anysurface.com/sguerin/models/shallowWaterDoubleSlit.html" target="_blank" class="">simplified Netlogo wave model</a> of the same shallow water equations without the acoustics. It's set up to explore double slit but you can change it to single slit and mess with frequency and gap and watch the wave propagations, diffractions and interference patterns<br class=""><a href="https://anysurface.com/sguerin/models/shallowWaterDoubleSlit.html" target="_blank" class="">https://anysurface.com/sguerin/models/shallowWaterDoubleSlit.html</a><br class=""><span id="m_-6704228213763452108gmail-m_-1165450971226445812cid:ii_ksbz5z8m4" class=""><image.png></span><br class=""><br class="">As a related aside, with some follow-up discussions with Ed Angel and Steve Smith I am also trying to understand how the gap might be considered a sampling function on the signal. My intuition is that the diffraction of the wave creates a spreader Sinc function and the gap is Rect function which are Fourier duals. In some way, i see <a href="https://en.wikipedia.org/wiki/Nyquist%E2%80%93Shannon_sampling_theorem" target="_blank" class="">Nyquist-Shannon Sampling Theorem</a> related to the gap. Note that diffraction creates a spreader function on the back wall in single gap experiments and the gap may be considered a Rect pulse when smaller than the wavelength.<br class=""><br class=""><span id="m_-6704228213763452108gmail-m_-1165450971226445812cid:ii_ksbxffke2" class=""><image.png></span><br class=""><br class=""></div><div dir="ltr" class=""><br class=""><br clear="all" class=""><div class=""><div dir="ltr" class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class="">_______________________________________________________________________<br class=""><a href="mailto:stephen.guerin@simtable.com" target="_blank" class="">Stephen.Guerin@Simtable.com</a><div class="">CEO, Simtable <a href="http://www.simtable.com/" target="_blank" class="">http://www.simtable.com</a><br class=""><div class="">1600 Lena St #D1, Santa Fe, NM 87505<div class=""><div class="">office: (505)995-0206 <span style="font-size:12.8px" class="">mobile: (505)577-5828</span></div><div class=""><span style="font-size:12.8px" class="">twitter: @simtable</span></div><div class=""><span style="font-size:12.8px" class=""><a href="http://zoom.com/j/5055775828" target="_blank" class="">z</a><a href="http://oom.simtable.com/" target="_blank" class="">oom.simtable.com</a></span></div></div></div></div></div></div></div></div></div></div></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
</blockquote></div></div>
</blockquote></div>
</div></blockquote></div></div></div><br class="">
</blockquote></div></div>
<span id="cid:f_kscinuuq0"><gapFrequency.mp4></span>_______________________________________________<br class="">Wedtech mailing list<br class=""><a href="mailto:Wedtech@redfish.com" class="">Wedtech@redfish.com</a><br class="">http://redfish.com/mailman/listinfo/wedtech_redfish.com<br class=""></div></blockquote></div><br class=""></div></body></html>