Re: Optimal sweep duration for BRIR measurements (Brian FG Katz )


Subject: Re: Optimal sweep duration for BRIR measurements
From:    Brian FG Katz  <brian.katz@xxxxxxxx>
Date:    Mon, 27 Jul 2015 18:46:05 +0200
List-Archive:<http://lists.mcgill.ca/scripts/wa.exe?LIST=AUDITORY>

This is a multipart message in MIME format. ------=_NextPart_000_0066_01D0C89C.7E7A6BF0 Content-Type: multipart/related; boundary="----=_NextPart_001_0067_01D0C89C.7E7A6BF0" ------=_NextPart_001_0067_01D0C89C.7E7A6BF0 Content-Type: multipart/alternative; boundary="----=_NextPart_002_0068_01D0C89C.7E7A6BF0" ------=_NextPart_002_0068_01D0C89C.7E7A6BF0 Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable Dear John, =20 As others have pointed out, increasing the length of the sweep increases your signal-to-noise ratio. For large room acoustics, we typically use sweeps of 20 to 40 seconds. This of course depends on the size of your = room, its reverberation time, and the power of your source. The weaker the = source, the longer the sweep. You should basically do a test and see what SNR = you get. If you need more, and your measurement chain is at its limit, the = only option is longer sweep. Longer sweeps will not help much with impulsive interruptions, while averaging will.With a single sweep, basically that freq-range during the noise is lost. If you are just measuring from the = RIR, this may not be an issue, as measurement parameters are often in octave-bands, but this is not true for auralization usage; corrupt data = is corrupt data, though I haven=92t gone through a thorough study of this = actual case. =20 =20 The sweep should definitely be longer than the RT, for room = measurements. Then, don=92t forget that you need to record the sweep length PLUS the = RT or more if your SNR is better than 60dB! =20 We have also compared averaging repeated sweeps vs. longer sweeps. = Avoiding the recent developments in overlapping sweep processing, this basic repetition approach =91requires=92 the decay of the first sweep to = finish before you launch the second sweep. As such, 3x20 second sweeps take longer = than 1 x 60 second sweep, due to the additional pauses. Of course, without repetitions, you have no backup in case something goes wrong, like a = door slam or something. So, I tend to use repeated sweeps and take the best = 1.=20 =20 We do all our processing in MatLab, and have never had an issue (within = the last 20 years) of processing long sweeps in real halls.=20 =20 As you are considering BRIR, and not anechoic HRTFs, you are subject to = the same conditions. If you want to convolve the BRIR directly, you will = need to ensure that the SNR is sufficiently high that the noise-floor is not = audible as a late reverb part of the BRIR. Some noise extension methods exist = from older studies on basic auralization and scale model RIR auralizations. I cannot imagine a 2s sweep for BRIR unless you are measuring an office of other room with <1 sec RT. I think there is a fault in your correlation-based analysis for reliability, and there are too many = factors to consider in comparing BRIR of different lengths. First, examine your = SNR. =20 Regarding distortion of the source, this should be an issue for the processing element (different from burning out you speaker). This is = because one of the strengths of the sweep method (when done correctly) is that = any harmonic distortion components of a higher frequency that the excitation signal at the time or folded back to BEFORE the direct sound after deconvolution by the excitation signal (as that frequency has yet to be generated). We presented a through work on this feature, extending = Farina=92s earlier works, to general conditions: =20 M. R=E9billat, R. Hennequin, E. Corteel, and B. Katz, =93Identification = of cascade of Hammerstein models for the description of nonlinearities in vibrating devices,=94 J. Sound and Vibration, vol. 330, pp. 1018=961038, = 2011, (doi:10.1016/j.jsv.2010.09.012). =20 This method lets you actually extract and analyse the different harmonic distortions (THD etc.) as well as allowing for the modelling of = non-linear responses. =20 =20 Cheers, Brian=20 -- Brian FG Katz, Ph.D, HDR Resp. Groupe Audio & Acoustique LIMSI - CNRS Rue John von Neumann Campus Universitaire d'Orsay, B=E2t 508 91405 Orsay cedex=20 France Phone. + 33 (0)1 69 85 80 67 - Fax. + 33 (0)1 69 85 80 88 <http://www.limsi.fr/> http://www.limsi.fr web_group: <http://www.limsi.fr/Scientifique/aa/> = http://www.limsi.fr/Scientifique/aa/ web_theme: <http://www.limsi.fr/Scientifique/aa/thmsonesp/> http://www.limsi.fr/Scientifique/aa/thmsonesp/ =20 De : AUDITORY - Research in Auditory Perception [mailto:AUDITORY@xxxxxxxx De la part de Anders Tornvig = Christensen Envoy=E9 : lundi 27 juillet 2015 09:07 =C0 : AUDITORY@xxxxxxxx Objet : Re: Optimal sweep duration for BRIR measurements =20 Hello John, The sweep method at a single frequency is an approximation to a = steady-state measurement with a pure tone. Longer sweeps give higher signal to noise ratio per sweep because it spends more time per frequency.=20 Short repeated sweeps (but not shorter than the length of the impulse response) are good, if time-varying or sudden noise that doesn't average = out is likely to contaminate the measurement. Sweep duration (and "rate" in general) also matters if the system (room = in your case) is nonlinear, time-variant, or both, but that's another discussion. Something is wrong with your implementation if the temporal offset of = the impulse responses you measure depends on the sweep duration. You should = be able to check this by connecting the output of your sound card directly = to its input. Also note that wrongly measured or wrongly computed impulse responses may be very reproducible in terms of correlation. Best, Anders PhD student in acoustics Aalborg University, Denmark _____ =20 From: AUDITORY - Research in Auditory Perception = [AUDITORY@xxxxxxxx on behalf of John Culling [CullingJ@xxxxxxxx Sent: Friday, July 24, 2015 5:25 PM To: AUDITORY@xxxxxxxx Subject: Optimal sweep duration for BRIR measurements Dear all, =20 Basic Q=85 =20 Does anyone have insight into the optimum sweep duration using Farina's method for measuring room impulses responses? =20 More detailed background=85 =20 We are planning to make an extensive series of measurements, and in preparation have been testing the method using different sweep durations. One way to = check the method is to correlate the impulses respones from repeated measurements or = those generated=20 with different durations. To our surprise short sweeps (1-2 seconds) = appear to give more=20 reliable results (repeated sweeps correlate, r>0.98) than longer ones. Comparing sweeps of different durations is a little trickier, because we find a temporal offset that reduces the correlation and can only be partially overcome by using cross-correlation. Nonetheless, it is apparent that durations from 1 second upwards correlate well, = while going below one second leads to reliable IRs, but ones that are inaccurate when compared with those from longer sweep durations.=20 =20 Our surprising conclusion is that ~2s should be fine, but Farina refers = to an ISO standard that=20 recommends very long sweeps (Farina has an example of 50s) to help = overcome noise. This seems an unintuitive rationale to us, since longer sweeps should increase both the=20 signal energy captured and the noise energy, and the method does not = involve averaging as far as I understand. Longer durations should help address brief interupting sounds, but=20 I am unsure if that it what was the idea. In the presence of continuous noise, we did not=20 notice any improvement in the IRs produced by longer sweeps. =20 The nascent plan is to take >1 short sweep for each measurement and = reject IRs that that don't correlate well with another. =20 Any insights/advice appreciated, =20 John. =20 Prof. John Culling School of Psychology, Cardiff University Tel: +44 (0)29 2087 4556 Yr Athro John Culling Yr Ysgol Seicoleg, Prifysgol Caerdydd=20 Ff=F4n : +44 (0)29 2087 4556 =20 =20 ------=_NextPart_002_0068_01D0C89C.7E7A6BF0 Content-Type: text/html; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable <META HTTP-EQUIV=3D"Content-Type" CONTENT=3D"text/html; = charset=3Diso-8859-1"> <html xmlns:v=3D"urn:schemas-microsoft-com:vml" = xmlns:o=3D"urn:schemas-microsoft-com:office:office" = xmlns:w=3D"urn:schemas-microsoft-com:office:word" = xmlns:m=3D"http://schemas.microsoft.com/office/2004/12/omml" = xmlns=3D"http://www.w3.org/TR/REC-html40"><head><meta name=3DGenerator = content=3D"Microsoft Word 14 (filtered medium)"><!--[if = !mso]><style>v\:* {behavior:url(#default#VML);} o\:* {behavior:url(#default#VML);} w\:* {behavior:url(#default#VML);} .shape {behavior:url(#default#VML);} </style><![endif]--><style><!-- /* Font Definitions */ @xxxxxxxx {font-family:Calibri; panose-1:2 15 5 2 2 2 4 3 2 4;} @xxxxxxxx {font-family:Tahoma; panose-1:2 11 6 4 3 5 4 4 2 4;} /* Style Definitions */ p.MsoNormal, li.MsoNormal, div.MsoNormal {margin:0cm; margin-bottom:.0001pt; font-size:11.0pt; font-family:"Calibri","sans-serif";} a:link, span.MsoHyperlink {mso-style-priority:99; color:#0563C1; text-decoration:underline;} a:visited, span.MsoHyperlinkFollowed {mso-style-priority:99; color:#954F72; text-decoration:underline;} p {mso-style-priority:99; margin:0cm; margin-bottom:.0001pt; font-size:12.0pt; font-family:"Times New Roman","serif";} p.MsoAcetate, li.MsoAcetate, div.MsoAcetate {mso-style-priority:99; mso-style-link:"Texte de bulles Car"; margin:0cm; margin-bottom:.0001pt; font-size:8.0pt; font-family:"Tahoma","sans-serif";} p.msochpdefault, li.msochpdefault, div.msochpdefault {mso-style-name:msochpdefault; margin:0cm; margin-bottom:.0001pt; font-size:12.0pt; font-family:"Calibri","sans-serif";} span.emailstyle17 {mso-style-name:emailstyle17; font-family:"Calibri","sans-serif"; color:windowtext;} span.EmailStyle20 {mso-style-type:personal-reply; font-family:"Calibri","sans-serif"; color:#1F497D;} span.TextedebullesCar {mso-style-name:"Texte de bulles Car"; mso-style-priority:99; mso-style-link:"Texte de bulles"; font-family:"Tahoma","sans-serif";} .MsoChpDefault {mso-style-type:export-only; font-size:10.0pt;} @xxxxxxxx WordSection1 {size:612.0pt 792.0pt; margin:70.85pt 70.85pt 70.85pt 70.85pt;} div.WordSection1 {page:WordSection1;} --></style><!--[if gte mso 9]><xml> <o:shapedefaults v:ext=3D"edit" spidmax=3D"1027" /> </xml><![endif]--><!--[if gte mso 9]><xml> <o:shapelayout v:ext=3D"edit"> <o:idmap v:ext=3D"edit" data=3D"1" /> </o:shapelayout></xml><![endif]--></head><body lang=3DFR = link=3D"#0563C1" vlink=3D"#954F72"><div class=3DWordSection1><p = class=3DMsoNormal><span style=3D'color:#1F497D'>Dear = John,<o:p></o:p></span></p><p class=3DMsoNormal><span = style=3D'color:#1F497D'><o:p>&nbsp;</o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:#1F497D'>As others = have pointed out, increasing the length of the sweep increases your = signal-to-noise ratio. For large room acoustics, =A0we typically use = sweeps of 20 to 40 seconds. This of course depends on the size of your = room, its reverberation time, and the power of your source. The weaker = the source, the longer the sweep. You should basically do a test and see = what SNR you get. If you need more, and your measurement chain is at its = limit, the only option is longer sweep. Longer sweeps will not help much = with impulsive interruptions, while averaging will.With a single sweep, = basically that freq-range during the noise is lost. If you are just = measuring from the RIR, this may not be an issue, as measurement = parameters are often in octave-bands, but this is not true for = auralization usage; corrupt data is corrupt data, though I haven&#8217;t = gone through a thorough study of this actual case. = =A0<o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:#1F497D'><o:p>&nbsp;</o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:#1F497D'>The sweep = should definitely be longer than the RT, for room measurements. Then, = don&#8217;t forget that you need to record the sweep length PLUS the RT = or more if your SNR is better than 60dB!<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:#1F497D'><o:p>&nbsp;</o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:#1F497D'>We have = also compared averaging repeated sweeps vs. longer sweeps. Avoiding the = recent developments in overlapping sweep processing, this basic = repetition approach &#8216;requires&#8217; the decay of the first sweep = to finish before you launch the second sweep. As such, 3x20 second = sweeps take longer than 1 x 60 second sweep, due to the additional = pauses. Of course, without repetitions, you have no backup in case = something goes wrong, like a door slam or something. So, I tend to use = repeated sweeps and take the best 1. <o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:#1F497D'><o:p>&nbsp;</o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:#1F497D'>We do all = our processing in MatLab, and have never had an issue (within the last = 20 years) of processing long sweeps in real halls. = <o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:#1F497D'><o:p>&nbsp;</o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:#1F497D'>As you are = considering BRIR, and not anechoic HRTFs, you are subject to the same = conditions. If you want to convolve the BRIR directly, you will need to = ensure that the SNR is sufficiently high that the noise-floor is not = audible as a late reverb part of the BRIR. Some noise extension methods = exist from older studies on basic auralization and scale model RIR = auralizations. I cannot imagine a 2s sweep for BRIR unless you are = measuring an office of other room with &lt;1 sec RT. I think there is a = fault in your correlation-based analysis for reliability, and there are = too many factors to consider in comparing BRIR of different lengths. = First, examine your SNR. <o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:#1F497D'><o:p>&nbsp;</o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:#1F497D'>Regarding = distortion of the source, this should be an issue for the processing = element (different from burning out you speaker). This is because one of = the strengths of the sweep method (when done correctly) is that any = harmonic distortion components of a higher frequency that the excitation = signal at the time or folded back to BEFORE the direct sound after = deconvolution by the excitation signal (as that frequency has yet to be = generated). We presented a through work on this feature, extending = Farina&#8217;s earlier works, to general = conditions:<o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:#1F497D'><o:p>&nbsp;</o:p></span></p><p class=3DMsoNormal = style=3D'margin-left:35.4pt'><span lang=3DEN-GB = style=3D'color:#1F497D'>M. R=E9billat, R. Hennequin, E. Corteel, and B. = Katz, &#8220;Identification of cascade of Hammerstein models for the = description of nonlinearities in vibrating devices,&#8221; J. Sound and = Vibration, vol. 330, pp. 1018&#8211;1038, 2011, = (doi:10.1016/j.jsv.2010.09.012).<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:#1F497D'><o:p>&nbsp;</o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:#1F497D'>This method = lets you actually extract and analyse the different harmonic distortions = (THD etc.) as well as allowing for the modelling of non-linear = responses. =A0<o:p></o:p></span></p><p class=3DMsoNormal><span = lang=3DEN-GB style=3D'color:#1F497D'><o:p>&nbsp;</o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:#1F497D'>Cheers,<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:#1F497D'>Brian = <o:p></o:p></span></p><div><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:#1F497D'>--<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB = style=3D'font-size:10.0pt;color:#1F497D'>Brian FG Katz, Ph.D, = HDR<o:p></o:p></span></p><p class=3DMsoNormal><span = style=3D'font-size:10.0pt;color:#1F497D'>Resp. Groupe Audio &amp; = Acoustique<o:p></o:p></span></p><p class=3DMsoNormal><!--[if gte vml = 1]><v:shapetype id=3D"_x0000_t75" coordsize=3D"21600,21600" o:spt=3D"75" = o:preferrelative=3D"t" path=3D"m@xxxxxxxx@xxxxxxxx@xxxxxxxx@xxxxxxxx@xxxxxxxx@xxxxxxxx@xxxxxxxx@xxxxxxxx" filled=3D"f" = stroked=3D"f"> <v:stroke joinstyle=3D"miter" /> <v:formulas> <v:f eqn=3D"if lineDrawn pixelLineWidth 0" /> <v:f eqn=3D"sum @xxxxxxxx 1 0" /> <v:f eqn=3D"sum 0 0 @xxxxxxxx" /> <v:f eqn=3D"prod @xxxxxxxx 1 2" /> <v:f eqn=3D"prod @xxxxxxxx 21600 pixelWidth" /> <v:f eqn=3D"prod @xxxxxxxx 21600 pixelHeight" /> <v:f eqn=3D"sum @xxxxxxxx 0 1" /> <v:f eqn=3D"prod @xxxxxxxx 1 2" /> <v:f eqn=3D"prod @xxxxxxxx 21600 pixelWidth" /> <v:f eqn=3D"sum @xxxxxxxx 21600 0" /> <v:f eqn=3D"prod @xxxxxxxx 21600 pixelHeight" /> <v:f eqn=3D"sum @xxxxxxxx 21600 0" /> </v:formulas> <v:path o:extrusionok=3D"f" gradientshapeok=3D"t" o:connecttype=3D"rect" = /> <o:lock v:ext=3D"edit" aspectratio=3D"t" /> </v:shapetype><v:shape id=3D"Image_x0020_2" o:spid=3D"_x0000_s1026" = type=3D"#_x0000_t75" = style=3D'position:absolute;margin-left:-.25pt;margin-top:.45pt;width:66.5= 5pt;height:34.6pt;z-index:251659264;visibility:visible;mso-wrap-style:squ= are;mso-width-percent:0;mso-height-percent:0;mso-wrap-distance-left:9pt;m= so-wrap-distance-top:0;mso-wrap-distance-right:9pt;mso-wrap-distance-bott= om:0;mso-position-horizontal:absolute;mso-position-horizontal-relative:te= xt;mso-position-vertical:absolute;mso-position-vertical-relative:text;mso= -width-percent:0;mso-height-percent:0;mso-width-relative:page;mso-height-= relative:page'> <v:imagedata src=3D"cid:image001.png@xxxxxxxx" o:title=3D"" = croptop=3D"22737f" cropbottom=3D"8693f" /> <w:wrap type=3D"square"/> </v:shape><![endif]--><![if !vml]><img width=3D89 height=3D46 = src=3D"cid:image002.jpg@xxxxxxxx" align=3Dleft hspace=3D12 = v:shapes=3D"Image_x0020_2"><![endif]><span = style=3D'font-size:10.0pt;color:#1F497D'>LIMSI - CNRS<br></span><span = style=3D'font-size:10.0pt;color:#1F497D'>Rue John von Neumann<br>Campus = Universitaire d'Orsay, B=E2t 508<br>91405 Orsay cedex</span><span = style=3D'font-size:10.0pt;color:#1F497D'> <o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-US = style=3D'font-size:10.0pt;color:#1F497D'>France<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-US = style=3D'font-size:10.0pt;color:#1F497D'>Phone.&nbsp;</span><span = lang=3DEN-GB style=3D'font-size:10.0pt;color:#1F497D'>+&nbsp; 33 (0)1 69 = 85 80 67 - Fax.&nbsp; +&nbsp; 33 (0)1 69 85 80 88</span><span = lang=3DEN-GB style=3D'color:#1F497D'><br></span><span lang=3DEN-US = style=3D'color:#1F497D'><a href=3D"http://www.limsi.fr/" = target=3D"_blank"><span lang=3DEN-GB = style=3D'font-size:10.0pt;color:navy'>http://www.limsi.fr</span></a>=A0=A0= =A0 </span><span lang=3DEN-GB = style=3D'font-size:10.0pt;color:#1F497D'>web_group: </span><span = style=3D'font-size:10.0pt;color:#1F497D'><a = href=3D"http://www.limsi.fr/Scientifique/aa/"><span lang=3DEN-GB = style=3D'color:blue'>http://www.limsi.fr/Scientifique/aa/</span></a></spa= n><span lang=3DEN-GB = style=3D'font-size:10.0pt;color:#1F497D'>=A0=A0=A0=A0 web_theme: = </span><span style=3D'font-size:10.0pt;color:#1F497D'><a = href=3D"http://www.limsi.fr/Scientifique/aa/thmsonesp/"><span = lang=3DEN-GB = style=3D'color:blue'>http://www.limsi.fr/Scientifique/aa/thmsonesp/</span= ></a></span><span lang=3DEN-GB = style=3D'font-size:10.0pt;color:#1F497D'><o:p></o:p></span></p></div><p = class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:#1F497D'><o:p>&nbsp;</o:p></span></p><div><div = style=3D'border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm = 0cm 0cm'><p class=3DMsoNormal><b><span = style=3D'font-size:10.0pt;font-family:"Tahoma","sans-serif"'>De&nbsp;:</s= pan></b><span = style=3D'font-size:10.0pt;font-family:"Tahoma","sans-serif"'> AUDITORY - = Research in Auditory Perception [mailto:AUDITORY@xxxxxxxx <b>De = la part de</b> Anders Tornvig Christensen<br><b>Envoy=E9&nbsp;:</b> = lundi 27 juillet 2015 09:07<br><b>=C0&nbsp;:</b> = AUDITORY@xxxxxxxx<br><b>Objet&nbsp;:</b> Re: Optimal sweep = duration for BRIR measurements<o:p></o:p></span></p></div></div><p = class=3DMsoNormal><o:p>&nbsp;</o:p></p><div><div><p = class=3DMsoNormal><span lang=3DEN-GB = style=3D'font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black'>= Hello John,<br><br>The sweep method at a single frequency is an = approximation to a steady-state measurement with a pure tone. Longer = sweeps give higher signal to noise ratio per sweep because it spends = more time per frequency. <br><br>Short repeated sweeps (but not shorter = than the length of the impulse response) are good, if time-varying or = sudden noise that doesn't average out is likely to contaminate the = measurement.<br><br>Sweep duration (and &quot;rate&quot; in general) = also matters if the system (room in your case) is nonlinear, = time-variant, or both, but that's another discussion.<br><br>Something = is wrong with your implementation if the temporal offset of the impulse = responses you measure depends on the sweep duration. You should be able = to check this by connecting the output of your sound card directly to = its input. Also note that wrongly measured or wrongly computed impulse = responses may be very reproducible in terms of = correlation.<br><br>Best,<br>Anders<br><br>PhD student in = acoustics<br>Aalborg University, Denmark<o:p></o:p></span></p><div><div = class=3DMsoNormal align=3Dcenter style=3D'text-align:center'><span = lang=3DEN-GB style=3D'font-size:12.0pt;font-family:"Times New = Roman","serif";color:black'><hr size=3D2 width=3D"100%" = align=3Dcenter></span></div><div id=3DdivRpF437927><p class=3DMsoNormal = style=3D'margin-bottom:12.0pt'><b><span lang=3DEN-GB = style=3D'font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black'>= From:</span></b><span lang=3DEN-GB = style=3D'font-size:10.0pt;font-family:"Tahoma","sans-serif";color:black'>= AUDITORY - Research in Auditory Perception [AUDITORY@xxxxxxxx = on behalf of John Culling [CullingJ@xxxxxxxx<br><b>Sent:</b> = Friday, July 24, 2015 5:25 PM<br><b>To:</b> <a = href=3D"mailto:AUDITORY@xxxxxxxx">AUDITORY@xxxxxxxx</a><br>= <b>Subject:</b> Optimal sweep duration for BRIR measurements</span><span = lang=3DEN-GB style=3D'font-size:12.0pt;font-family:"Times New = Roman","serif";color:black'><o:p></o:p></span></p></div><div><div><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:black'>Dear = all,<o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>&nbsp;<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:black'>Basic = Q&#8230;<o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>&nbsp;<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:black'>Does anyone = have insight into the optimum sweep duration using Farina's = method<o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>for measuring room impulses = responses?<o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>&nbsp;<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:black'>More detailed = background&#8230;<o:p></o:p></span></p><p class=3DMsoNormal><span = lang=3DEN-GB style=3D'color:black'>&nbsp;<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:black'>We are = planning to make an extensive series of measurements, and in preparation = have<o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>been testing the method using different sweep = durations. One way to check the method<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:black'>is to = correlate the impulses respones from repeated measurements or those = generated <o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>with different durations. &nbsp;To our surprise = short sweeps (1-2 seconds) appear to give more <o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:black'>reliable = results (repeated sweeps correlate, r&gt;0.98) than longer ones. = Comparing sweeps<o:p></o:p></span></p><p class=3DMsoNormal><span = lang=3DEN-GB style=3D'color:black'>of different durations is a little = trickier, because we find a temporal offset that = reduces<o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>the correlation and can only be partially overcome = by using cross-correlation. Nonetheless,<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:black'>it is = apparent that durations from 1 second upwards correlate well, while = going below one<o:p></o:p></span></p><p class=3DMsoNormal><span = lang=3DEN-GB style=3D'color:black'>second leads to reliable IRs, but = ones that are inaccurate when compared with those = from<o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>longer sweep durations. <o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>&nbsp;<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:black'>Our = surprising conclusion is that ~2s should be fine, but Farina refers to = an ISO standard that <o:p></o:p></span></p><p class=3DMsoNormal><span = lang=3DEN-GB style=3D'color:black'>recommends very long sweeps (Farina = has an example of 50s) to help overcome noise.<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:black'>This seems an = unintuitive rationale to us, since longer sweeps should increase both = the <o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>signal energy captured and the noise energy, and = the method does not involve averaging<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:black'>as far as I = understand. Longer durations should help address brief interupting = sounds, but <o:p></o:p></span></p><p class=3DMsoNormal><span = lang=3DEN-GB style=3D'color:black'>I am unsure if that it what was the = idea. In the presence of continuous noise, we did not = <o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>notice any improvement in the IRs produced by = longer sweeps.<o:p></o:p></span></p><p class=3DMsoNormal><span = lang=3DEN-GB style=3D'color:black'>&nbsp;<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:black'>The nascent = plan is to take &gt;1 short sweep for each measurement and reject IRs = that<o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>that don't correlate well with = another.<o:p></o:p></span></p><p class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>&nbsp;<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB style=3D'color:black'>Any = insights/advice appreciated,<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>&nbsp;<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>John.<o:p></o:p></span></p><p = class=3DMsoNormal><span lang=3DEN-GB = style=3D'color:black'>&nbsp;<o:p></o:p></span></p><table = class=3DMsoNormalTable border=3D0 cellspacing=3D0 cellpadding=3D0 = style=3D'border-collapse:collapse'><tr><td width=3D301 valign=3Dtop = style=3D'width:225.4pt;border:solid #BFBFBF 1.0pt;border-right:solid = windowtext 1.0pt;padding:0cm 5.4pt 0cm 5.4pt'><p class=3DMsoNormal><span = style=3D'font-size:12.0pt;font-family:"Times New Roman","serif"'>Prof. = John Culling</span><o:p></o:p></p><p class=3DMsoNormal><span = style=3D'font-size:12.0pt;font-family:"Times New Roman","serif"'>School = of Psychology, Cardiff University</span><o:p></o:p></p><p = class=3DMsoNormal><span style=3D'font-size:12.0pt;font-family:"Times New = Roman","serif"'>Tel: +44 (0)29 2087 4556</span><o:p></o:p></p></td><td = width=3D301 valign=3Dtop style=3D'width:225.4pt;border:solid #BFBFBF = 1.0pt;border-left:none;padding:0cm 5.4pt 0cm 5.4pt'><p = class=3DMsoNormal><span style=3D'font-size:12.0pt;font-family:"Times New = Roman","serif"'>Yr Athro John Culling</span><o:p></o:p></p><p = class=3DMsoNormal><span style=3D'font-size:12.0pt;font-family:"Times New = Roman","serif"'>Yr Ysgol Seicoleg, Prifysgol Caerdydd = </span><o:p></o:p></p><p class=3DMsoNormal><span = style=3D'font-size:12.0pt;font-family:"Times New = Roman","serif"'>Ff=F4n&nbsp;: +44 (0)29 2087 = 4556</span><o:p></o:p></p></td></tr></table><p class=3DMsoNormal><span = lang=3DEN-GB style=3D'font-size:12.0pt;font-family:"Times New = Roman","serif";color:black'>&nbsp;</span><span lang=3DEN-GB = style=3D'color:black'><o:p></o:p></span></p><p class=3DMsoNormal><span = lang=3DEN-GB = style=3D'color:black'>&nbsp;<o:p></o:p></span></p></div></div></div></div= ></div></div></body></html> ------=_NextPart_002_0068_01D0C89C.7E7A6BF0-- ------=_NextPart_001_0067_01D0C89C.7E7A6BF0 Content-Type: image/png; name="image001.png" Content-Transfer-Encoding: base64 Content-ID: <image001.png@xxxxxxxx> iVBORw0KGgoAAAANSUhEUgAAAMsAAADLCAIAAABPmfBkAAAAAXNSR0IArs4c6QAAAARnQU1BAACx jwv8YQUAAAAJcEhZcwAAIdUAACHVAQSctJ0AABerSURBVHhe7Z2JW1NX2sDnf4F2phWoLOKCrdpa 7XSZts902tk6nZnO0na+URaxaq0du8x0l0XFDdQiaN0XZF/DvmYFwpaEkJA9uUs2EPt87829xAAq CeYohPf3HHkw5567nPzyvufce3P52dTU1G0EIcbP4N9PCEIMNAwhCxqGkAUNQ8iChiFkQcMQsqBh CFnQMIQsaBhCFjQMIQsahpAFDUPIgoYhZEHDELKgYQhZ0DCELGgYQhY0DCELGoaQBQ1DyIKGIWRB wxCyoGEIWdAwhCxoGEIWNAwhCxqGkAUNQ8iChiFkQcMQsqBhCFnQMIQsaBhCFjQMIQsahpAFDUPI goYhZEHDELKgYQhZ0DCELGgYQhY0DCELGoaQBQ1DyIKGIWRBwxCyoGEIWdAwhCxoGEIWNAwhCxqG kAUNQ8iChiFkQcMQsqBhCFnQMIQsaBhCFjQMIQsahpAFDUPIgoYhZEHDELKgYQhZ0DCELGgYQhY0 DCELGoaQBQ1DyIKGIWRBwxCyoGEIWdAwhCxoGEIWNAwhCxqGkAUNQ8iChiFkQcMQsqBhCFnQMIQs aBhClkgzDA5nysekj4mJCfh569YteEVYYlnCd4u/Q3geTrdElGFwLNBr0IlOl0sqHs799nLO15f6 FGqX2+2dmIDejLCPU5Dw3eL1eqFbOlr7s766eCTr2siwDrqF90xYjgyRY5i/HxmGvn65edOaDxNi UqE8t25XRWk7wzBQtQwjGd8tHq5b2KLT1esTd/DdsmX97rYWBcuy8IEk2i2RYxifBVjWOT5uevn5 TxLi0hOe8pW4tLde/cJstrBOJ/RmxBxv8HCfOpZV9qvXxUOHpE13S/q7f/zOYrU5nS7oN3Ld8ggM Yxk3w7ghZQn/DxNgmMfjsdsdyoGRDat2Cv3oK5vW7urrG3RQFCywDMMYZEOb3d4o6k56cltgt2x5 Zs/wsJqiaVAwYgy7/eVnF7as3wPlV5v/U1Uu/il8Gwd1XC6X2WKRy8UvPc19RoWujEt/4+VP+5WD EMagr5ebYXC8TqfTZLa0tHSvXQkd4o9haW+/+d/BoWGr1er2eCLEsMmJW1ue/ggOL953qO//JVeo CAe8YUajSSEXHzxw+pnknYlx6VA2rM48mneht18JvQxD3WVoGIwcDAajTN775WcFaxN28N3y3Nqd xUUl/QODFosVPngRYhhDuTat2S18hp5Kf/f3WUJFOOANA42UA0N1DY0nTpxN++Cb1Pe/Lig419DY BC9C1bKNYUaTua9/oKa2Ie9w4bZ/fp2x7buiMxcbm9sGh0YskRTDGNq1aW2AYX84IFSEA+hK6CkY h6k12m6xtKqm7tqNm9eu34RfeiQyzajW7nAs23GY1WZXqTWd3eKKqhrolus3SmvqGiRS+ah2DIan kTMOI2oYHAhMFSEjwHhreEQtlSk6OrvbO7vglxGV2mQ2w3xqec4l+VMVMH4YGh4Bq6BbOrq6ZYpe tWbUbLHCFDty5pJEDQP8J34sZqtOp9eMjkHo0un1MCcHvaAqxLOL0DOBnTNvR/mXD6kVcfhuoWkG Pntj092iHx+3+rqFPxctLEqAiDIMeorvTRh5NInkBcfKC09W9vdqYIDPn269V1d6vZM2K9Ov0Irq 5RfONRYcqThxuDw/ryI/D36Ww3/PFjZUlfX0dA6ajQ5YmdBsGpfTYzTYuzuHq8t7LpwVFRZU8w3z D5cfP8Q1v3KhuUXUNzSgoxzOW7dCeDsnvJPQxKi3KXu1PZ1DbS3KxgZFk0gB+zkyPC4sNB9w1BCl uLNiDFt+s73gWNm5olrtqHHebgkLEWUYAIcD/aWQqdfGpyfEcievU5Iyx/XWe3XixMTklfPNf/nt d+sTM1fHpq1akcq1ivOd7AgsMalJK1KTY1LXrcz4zUuf5R+pgGOB5jYbk/vt1Vdf+BReh9qkFWnc GfO4NH/DeDhSmLvFbE9ekbomLu2Z1bt2/OuYuGuI3/p9gB27er7pH388AE3WrUyHtrB+2D1/ef7p PRCXhaXng++Wa5daV8VshwOE/dmy6WOW4WbW9+qZcBFphvm4feZkbSK8zb6tgBkdrUqhJgA48GZR 75uv/DeJd2J6r+YvcWnQ5KVn9+3bXfh8yp7EGPAp2ObxT3ELpyRkZPzfMYiawq7Mobtj8I2XvkgS 1nyPlcel3bjcJjQIjr07TiXETndLTKp+zCpUkCQiDfvpzMkav2GJMantzWDY7MOsLOt+ZvWHibFp vpNzoZX4lRnwBidCwJtTFUzhthiT+s5b31nMlLA3Acglqs3rPwIV7r9jcICF+dVCm+D4OMCwVWjY gxBoGMQwn2EzcLu8m1P2cCksoHABJjaNy4axqZAxITH5C5dAIdTd/10H5yCRxXDJdHVAw+SY7fC6 LxrNXD4m9fN9xcIOBfC3d3ICF46HkLmCy9GzyvqkzEGlTmgTHGhY2JjPsNv7dp6e8S6uzEhcsf2F jXs/3vlD0cm6mkpxW1O/pHsYxnNQxJ1DME6/eK5pT1oBJJdZXnIlLn3Viu0wZsr55ur5IhGM99tb lNCws1XZ2tRfcqU9L7vkvT9lc57NdPSZ5A9NRvtPU3dmuJOTtzas/tC/AITJD/6cU1hQC1ONwFJ0 uq6rAwZzoQ2h0LCwcX/DKMrJ7QaMxH0LwFuekphxvqjBbrvnwIgHpoFdbYPrEnb423IlLu1Xz3/a K9NMTtzvVAiooxkxvv/nXG4mMd0WcvTJoxWBogz0jSXHZwi1K1JPHCn3eiaFugcGDQsb9zdMqzFD 8vLvBmSr/XuKfEP/+bti6tbUe3/JATP8zZNi0y6ebeTq5m9+W6e1bFyzy982ISYNpgtCpQ8w9Y5h T25vEvUKFeEADQsb9zdMPWKEpObfDXgjCwtqhLog+PzjYjDD33x94o5e2ahQFwQvbPzY3xbmoen/ Ph74FiikmtUBhtVUSYSKcICGhY15DUuKmWHYD6FMyjjDAmLYhuRM9bBBqAuCrRv2+tuCYRnbTgS+ BYZx29rETP8Cr27dD1FNqHtg0LCw8ZANUw0Fe3oduL9hwEubPvEvAIO8lKTMt371v23/Ov7ZJ2e/ /PzC919fuWvJ+f5a8enajhYldLKwojmgYWFjSRv2xf5zgbMBocDhwEahwIT0rmUFlO2Jsamvb91/ 8lil03mX0/1oWNhY0oaBHL9//avATYRUYGoMR/TX333HMm5hjdMsdsNu3eK++QRAk+BbBYKGAfMa BlgtdPoHx7idhOAUeGYk+BKb+vl/fpx1JWPxGgbLTExMuHVS5/V/ukfb4PfJyYWcpEHDgGAMA7hz b+0D335x4bevfbk5ZQ/MWNc8lb76HiUp1ncJKyY18HTuxrW7DeN2YXU+FqlhvF5Ol4sVfeXMimYq 97IQx70LuVkUDQOCNCwQm4WWiVWNdQpRrXxuaaiRVZX1XDzbuCutICVgHnrnRN00i9QwSIsut5um GUdZJpsbbf/xHYeN+0oFaCcsETRzDAvnffqBRJhhwXP0UCl3bYpff0zqwe+vCxU+Fqlhvq+5shaL 1XL5fefBKMexlHHdKEXTHk+wNyf5cbu9z6Vw3zXiyx9+/bVQEW6Wq2G3q8vFd+74iEnN/e6aUONj kRrm9d2Aqx832C78zXUois5+TDsksdns7tANg229vPkT/5Xj9at2QvDna3w/w8ayNazsRmei/3LF UjEMsqGDorRjOuv5d8EwNita3VNpXuh3D/8ReGtKXNovN+3r6hiiHE6hOkwsZcMWbhvMQFPfP+q/ ZgrpMv9IuVDnYzEaBrUQw+wOh2ZUa/oRDItmc6JV9QUGo4llue+oCMsFTfGZ+hk37sWlJz6x/dmU 3f/++6ED/7tUmF/z45kGGJ9W3OyqKOkMrnTNvSdi6Ro2dfs2RbHVpcEffueNy62FBTV7M09yY9yA C6bJT6U3N/QJ6/WxeA2z2R0ajVZ3/j32YLTrYJT+StqYTg+pcyJ0wyiH69Wt+/29IBTulvYQ7kUO LIlx6W++8oWw9mmWdJZ8+41v/DsfQuE6cMYrv3nlv7cmZySZRZoluceNQAzTaMeuZTpzo1wHo60F W9RqtYNy8N9UEZYLlttyifq5dbugR+53v2goZW3CDv2YTVi9j5ANO1El1AVBeA2bdW/FxMTkM8l3 7kBcWIlfmbE6Lq2ytFtY6TQfZ57y59DFNdLns+RIyRdsNhfDqNwn1JI6q83ucrtvLej5ZhqV8d3f Z62J37Hg0BVYVsdn9PeOCavmmCo+XRdoWGfbgFDjA7Y+6+6dM6dqhbog+HL/j4HJaGPyTtWwUagL ghef3edvC4btSM0PfAtYxrM+acaDg0IrcWlw4BtW7Tx/pkFYYwD7PvzB/9lIik0b1834WBIi2JH+ qHZsqOqI8wDEsChndpS2Oot/0MjCDONR9mp3bjux+emP1q7MSF6RmvQklO0LKC9s2CuscRp4y59O yuRrN6zZZbXM+MIFy7hfenafv3lyzHblDEHnoat9cH1CxnTz1Hfe/FaoCI6dqfn+TUOkuVDcOOtm 6Ne27vcvEHxZ9ST3FbcNyTuhS2mKFdY1k4Ya2Zq4NH751178zOsJ+YzmAggqhlEUPaYb72++yn7P xTBXbpTxzG+1Y3qaYRZ2+cgHv90pJ+uGcD0yNN6n0DbWy2urJCGVmkqJTmvxrWoGEKjqqqQ1FWKd FnLBjLdwauq2zUq3iHprKyWNDYpRtSnUXA8r57ZeKenuHKYc7Lx96AeWBL/lMg3sdk2luE8xOrep QW+DPQ88xnkLLN/TNawaNtznG3I8ffLR2iqpqFbuWzLY3X4QgophYNK4wdDX3Ujl/IIzDBLl4VXq AandbneTfDAQEgHMbxjkQZZlDUZTr1yiP/6i02cYmwWJ8gD/6MoHSZRIxBOcYdxD9MzKgaGRi+mQ IvlEaTu+QT+mgfC2gAuUyPJhfsNgAZgzWqy2oWGVvOoUm/U4nyjZnF+MSqptdgfR55shS52gDPOd EuOmk+K2Bkd2LG8YpEvj2bdNRsPCTu4jy4T5DQNgOgnZUKcfl0ml+hOvCImSG41F62TVDorCMIbc i6AMgxDFPc/YZFYODstvfMvm+s5Z+Ir54t+NJpPTxT2TXVgaQQIIyjBYxs09C9SmUo+2t7faDq/y G8bkPqmX3OSfyU76SVTIUiQow4Dp8656iUyhKnzHnyidOVHW4t+ZzCb+caDC0ggyTbCGgT3gkNFo Ug4OdVeeYbJhRunLlYdgUvn4eHuxw2GHCUGQa0OWD8EaBou5XC6b3a7WaDq7urQFb4BbQhg7FMUc jDdoR/iHPWOuRAIJwTDvxAR/O7W8t6/nZp7zwJ3xvjMn2nxlm81mW/DdFkikEqxhADejdLksFuvw iKqjrWns6FYhUfoKnf2Ysb2IxiE/MpMQDANvYKTloCgY78sUvR2lBXT2z/2GQaHy1phUUv5PFga/ WiSyCcEwYMI33jdxf3RD1d7WPnb8RVdutNMfyXKj7D+8zv0lJu5yOM4rEY7QDBPCmEMIY81VV+y5 K12H7uRKNucx6/UMu93q9uCADOEIzTAAMiCEMbOFC2Mdnd19xdtd2YJefHHmPG5uPk5TVOh/5QWJ QEI2DKRxezwURcOkUtHbL6qv1eVt5m8a8xfm+2iL+Cp/Yw+O+pc5IRsGgGQul8tqs6s1oz1iafON fDZrhmFQ6IMxlv4aGPWT/sNMyCJnIYaBMfzdFuMGo3JgqLW1vfvcR76z/HcMg+E/dXSNdaiFv5qE ki1bFmIYwIUxt9vucGjHdIrePpGoXn38VZhXBo76nbnR9OEE62CT0+mESShKtjxZoGEAdwLW6bRY bZArxRJZdWXZ+OGN/ivi/kLlxtoGRBDJ8CTZ8mThhkFDyJX8H18dGBzu7OppKCm0ZsfNMswF6fJw vK2vEsdky5OFGwZAGIN5pYOidPrxvv6BlraO5ks59IHH5kgWxeSttHafZRgaZ5fLjQcyDFwByVwu t9Vq4wdkjc2trRe+oXKemGUYZE9nVrSt8TBD2TweD54nWz48kGE8EJacThdIphnVyuSKhsbmtrP/ mXXJki8sSHZzp8Nu9njcU1PcM62FVSCRSxgMgzVwJ/pZ1mS2qNQaiVRe19DYUryfOQBTy9mSOXMf cxS+YjeqYSqKGXM5EAbDABAFch/NcJKNqNRiqbRe1Nh6/ksq6y5jMsiYdF6KXXqZn2BixoxswmMY AKJ4vF6KZowm8/CISiyR1Yua6i8dMuckzT2FAbGNyf6Fveoz2qIDNScmuLNlmDQjkrAZBnCSQSTz SQaRTCKTi5paKq4Wjx9KmSuZ82A0mx1N5W+0a6QuFxfMMGNGJOE0DAiUDMZkMnlvU3Nrecll5fE3 YZg/SzKu5EIwi7GX7+WDGTTHSBZhhNkwgE+X/JhMrRmVK3pb2jqqKso7Cj+0Z6+4S8aEeJYdRRc8 7+gsYlgG2k5OTqJnEUP4DQN4yRiWtdpso9qxvv6Bjq7umtr62ktHdIc3sTnRs+724STLBc8eo4vf oJTVLG3z+jzDvBkBEDEMADnAEpgtWq22MZ1+YHCoWywVNTbfLCvrPvUv+oDwpLvZ5VAUeyDaceY1 eqTRydITE/ilkiUPKcMATjLfn9xyOCiD0QTDMqlcwWXM6pqyy0WDR3/NjczmnDDjSm4Um/uE49QW qu04bTPA+Azj2dKFoGEArJy7rOR2UzRttljVGi1kzK4eSUNjc3nZTdHZrzR5zztzou86OIMXoYo9 sZ4qzWT6ylyUxeOB1Mn94UGwjehuI2GErGE84ITb42FZp81m1+sNwyNqGP63dXTWNTSWlJY2F+0f O7yJzv6572H9swuM2Jy50S6YCpxIocv30IoS1qB00Wavxw0Bkp978ggbQxYZD8Mw3gDwDDImRdEW ixWG/8qBIbFUBkmzpq7h2o3S6nNZI3m/ZL5/nItnd02dvhdhlAaFhsB2M43qOkWrWliHEUZsTpZx u51utwvMm/B6fcUziefYFgEPwzA/EHK4YOZ02h0OGJlpRrXKgUHwrLm1vbZeVFp6s/zSqe7CdN3B dXc/eRZYcrjizHmSPRTPHl3N5G9gCl9mzr3FnH+bufIuFPbCn9iK3bA5vDD1aHmohgEQVOAt9yfN cYMRPOtXgmfy1vbO+oam8sqaspJrtecPSk6/N3pkqy03lsmOhgHZPYdrQib1/QI/efOgSRZ39ZOi uMfMQvgUNo88dB62YQBsEYIKeAYzAJphrDbbdDwbgvFZZ3dPc0tbXX0jqFZSWl52pbjp7OeKU++q 8l40HVprz4l35MQw2T9ns7lHl3HO+ayC/7LZj9PZK6jcBPvRdYZTr+uv7x4X37RabZCaYVvCtpGH ziMwjCfQM4Zh7XaHyWTWjulGVJp+5YBEJu/o7G5pbYdZJwzUKqo4266VlJWUlFReKWy8erTl0gHJ 1a9llz5VXN6vuP6V+EZ2V1l+e83F9pYmaCiRyiH/grVmi4Xx3b0tbBV56Dwyw/zADoABHt/4jJsH WCGkGWEqMKJSgyVyRZ9YIuvo6oEc2tTcKmpqgRkoX2DoBj/rRU3wIlTBpKG9swtCICyv6O0fGh4Z 0+mtNjt/j5CwMeSh8+gNA2AfYHwGIc3r9bohpLGsw0GBakaTWacfB9uGR9SDQ8MQ2xS9fVK5QiKV 9Uik3WIoEig9YilYJZMrQKy+/gFYUqXWgF5ms4WiaRyHPVoWhWGBgGcghMfrhfETzAZAEbuDguGU yWw2GIz6cYN2TA/OqTWjKrWvqDTwE/4LL0IViAWzB1gYphEwyOPvpMW55CNk0Rnmh49q4AcENhDF JxwLIzbIpA6KsjscNvuMAi9CFQ1aMSwsDKELL58vBhavYX5gDwE+tvHCwaBtLvA6AAsAsDCItfgP bTmwBAybi0+5uyMsgSwalqRhyBICDUPIgoYhZEHDELKgYQhZ0DCELGgYQhY0DCELGoaQBQ1DyIKG IWRBwxCyoGEIWdAwhCxoGEIWNAwhCxqGkAUNQ8iChiFkQcMQsqBhCFnQMIQsaBhCFjQMIQsahpAF DUPIgoYhZEHDELKgYQhZ0DCELGgYQhY0DCELGoaQBQ1DyIKGIWRBwxCS/PTT/wPb9nLFn1yCPgAA AABJRU5ErkJggg== ------=_NextPart_001_0067_01D0C89C.7E7A6BF0 Content-Type: image/jpeg; name="image002.jpg" Content-Transfer-Encoding: base64 Content-ID: <image002.jpg@xxxxxxxx> /9j/4AAQSkZJRgABAQEAYABgAAD/2wBDAAoHBwkHBgoJCAkLCwoMDxkQDw4ODx4WFxIZJCAmJSMg IyIoLTkwKCo2KyIjMkQyNjs9QEBAJjBGS0U+Sjk/QD3/2wBDAQsLCw8NDx0QEB09KSMpPT09PT09 PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT3/wAARCAAuAFkDASIA AhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQA AAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3 ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWm p6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEA AwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSEx BhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElK U1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3 uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwD2Wiii gCOGeK5j8yCVJUyRuRgwyOvIpLm4jtLWW4nbbFEhd2xnAAya5j4aDHhBMf8APzP/AOjDXRao5i0m 8kXG5IHYZGRkKTVSjaVgDTNSttY06G+snL28y7kYqQcfQ1arzxPEWrNofhWy06S3gvNWjJe4aIFU 2jJwo45p3iG78WeGIrSQ63bXgu51tgsloE2M3RuD2q/ZO9rgeg0Vweuy+LPC+kvq02uW17HAVMtu 1oEDAsBgEc96TxJ4qvbfWoLVdUt9FtJLRZ0uJrfzfNYnlQegxSVJvZgd4zBFLMcKoyTVLR9Zs9ds BeafIZICxXcVKnIODwapaJ4l0vV4VtbbVra+u0iHm7BtLkDlgvpmq/gO8l1Dwpb3Nxs815JASiBR w5A4HsKTi0ncDo6KKKgChrZddJmaPO5cNx7HNW7eZbi3jlQ5V1BFLLGJYmjYZVhg1g6ZeHSLttOv DiItmFz0Ge1YTn7OopPZ6ej6febQjzwaW61/zM2y8O+J9AWW10XUNNewMrSRrdRNvTcckZXr1rSt LTxLcRXkGsT6Y0M1u8cf2ZHDByMAnPbrXQ0V1Obe5icda+DbyB/CzNcQH+x0ZZsZ+ckY+Xj+daXi 3w/P4gt9Pjt5oozbXkdw3mZ5VeoGO9b9FHtJXuBjeLtFm8Q+GrvTbaSOOWbbtaTO0YYHnH0qpqen eIX+zwac+jtaxQqpF3E7NvHBIxxjpXSUUlNpWA4/S/C2qnxJbarq02mp9ljdIo7GEpuLDB3E1r+E 9Fm8P+H4bC4kjkkjd2LR5x8zEjr9a2aKcpuWjAKKKKgAqlqWmQ6lDtkGHH3WHUVdopSipKz2Gm4u 6Ocin1PRf3csRurcdCOoH1q9F4itJANyTo390xk/yrVxnrTPJjznYufpWMaLhpCWnZ6/8E1lVU9Z R18tBysGUMOhGaWiitzEKKKKACiiigAooooA/9k= ------=_NextPart_001_0067_01D0C89C.7E7A6BF0-- ------=_NextPart_000_0066_01D0C89C.7E7A6BF0 Content-Type: application/pkcs7-signature; name="smime.p7s" Content-Transfer-Encoding: base64 Content-Disposition: attachment; filename="smime.p7s" MIAGCSqGSIb3DQEHAqCAMIACAQExCzAJBgUrDgMCGgUAMIAGCSqGSIb3DQEHAQAAoIIL2DCCA20w ggJVoAMCAQICAQAwDQYJKoZIhvcNAQEFBQAwLDELMAkGA1UEBhMCRlIxDTALBgNVBAoTBENOUlMx DjAMBgNVBAMTBUNOUlMyMB4XDTA5MDEyMTA4NTExM1oXDTI5MDEyMTA4NTExM1owLDELMAkGA1UE BhMCRlIxDTALBgNVBAoTBENOUlMxDjAMBgNVBAMTBUNOUlMyMIIBIjANBgkqhkiG9w0BAQEFAAOC AQ8AMIIBCgKCAQEAuLp3q1edEiUZFUSKI1XaYw0GbVnrwQSgF3kS8yeCSOqNvjzLNjV4GwTVKQO1 aUgp+6ena8Tc0ikQ8+JOlaUArUl1whzrIjjNmQh/8vEkqvEqTQxss6Iidiv7d7ZgPD2IV6nrrxal kGFRgP/DfLmug70lXaA9lLlsNpzITT1j3yfOxA1ktUIlkTYVnJNLIRk09ikTQc4Pu5LvlqxvVDrz wCsEGdH+vNzu8eFOQ47mE9wn2rsAcl2zfy4wmGjGPQUsAgdKEjbN2r7PDGsKTBoK69sB8uGcrELE klxXJrbW9dSKvWscDDMaLu2mhIQ2Jsyk7/4k+Rtfri1ZA1AuDwhx9QIDAQABo4GZMIGWMA8GA1Ud EwEB/wQFMAMBAf8wHQYDVR0OBBYEFFCXtg33rDMXr/EdRjxrO/8AoOXlMFQGA1UdIwRNMEuAFFCX tg33rDMXr/EdRjxrO/8AoOXloTCkLjAsMQswCQYDVQQGEwJGUjENMAsGA1UEChMEQ05SUzEOMAwG A1UEAxMFQ05SUzKCAQAwDgYDVR0PAQH/BAQDAgEGMA0GCSqGSIb3DQEBBQUAA4IBAQA0MWtAIb+H 9V58oaX1/xjSJOZQvGeYrjsMGMeICK2T0QYJkGr+qBTxkA1rSHDZ0Gd6hZRqn/BpndZ4/9e0/OoI 7c5aCFpTLq/18t47spk5/X2x68xUdAUXEGZfGs429Uvhn5Ox+CuKiBEx/PNqo40d1gpd6vlM1uJC 5jwLLuxiTGkOLCaepo3nN69faiEMRekv80m3G/VfIVIgJAowymK9AYJmTEHyy0L9Yq9EeylpGRWl gUTow958lDNmnhOPUYPPU/gN5668TvqByB8sKWTda2MjnB/H93c64jLqaSH1nTnCd/oGtNIIoZBH DKxRaBwsD32n1qDRVfU+k/SqjRnNMIIDtjCCAp6gAwIBAgIBAzANBgkqhkiG9w0BAQUFADAsMQsw CQYDVQQGEwJGUjENMAsGA1UEChMEQ05SUzEOMAwGA1UEAxMFQ05SUzIwHhcNMDkwMTIxMDkwMzUy WhcNMjkwMTIwMDkwMzUyWjA1MQswCQYDVQQGEwJGUjENMAsGA1UEChMEQ05SUzEXMBUGA1UEAxMO Q05SUzItU3RhbmRhcmQwggEiMA0GCSqGSIb3DQEBAQUAA4IBDwAwggEKAoIBAQCcqWRqtAcjGcO+ CAjE4ipd7dQo3v4/oO2SRGsWRzUTN410qfVNg+QFG48IYUzed1vXqf+7ZZ3JrDwP8Yl7kyKgRw9f Br6OeKALeddCZWFyY5JAiltnEPrWJz7Yd00eedtvSxIAUlji2EIFsjM3d+CNbB28+PbylFgl94LZ htlc89guzf1wiQSH7znRKtiSjGkvHFJl+mCPEH/Mj/lcLhTtTIvgnTe40jyUylD9QXre9jzsruy9 XBbdsPBzmXqB/zsz0gCeTYzGWzguM8N2hspyVNst4hD9klBok4Ttd6mJaVGF6+CUo/mXCBGMJPAP uji0LNbXbE5j/rQxn++PplqlAgMBAAGjgdkwgdYwDwYDVR0TAQH/BAUwAwEB/zAdBgNVHQ4EFgQU EePZ0VJHG1mxPBt4Zmv0oYjtClswVAYDVR0jBE0wS4AUUJe2DfesMxev8R1GPGs7/wCg5eWhMKQu MCwxCzAJBgNVBAYTAkZSMQ0wCwYDVQQKEwRDTlJTMQ4wDAYDVQQDEwVDTlJTMoIBADAOBgNVHQ8B Af8EBAMCAQYwPgYDVR0fBDcwNTAzoDGgL4YtaHR0cDovL2NybHMuc2VydmljZXMuY25ycy5mci9D TlJTMi9nZXRkZXIuY3JsMA0GCSqGSIb3DQEBBQUAA4IBAQBP6eMX5kwn9RdqUFXq7s9Mqr3uLN12 PlmgLwtoidpD0B4dLGsS+WcTnmUrpITSMc+CoMrvIugQbOnPF3ysrYeFTIrh1XtM5tkOsqHkKZ4g rw8+e7l/JEwy4Ij+etqsFmChda6JnIpmnZS6fqsOICPt98eEXYJIFN7E1LDFTiyH1I0DES821Slv rplnGQyNbVog39GfxsuvVe4LyrFbD+9cdPTzRoljR34szr/IMA0ARvYY5mw/MH6QbhCNHzlJz5Np 68iewmdzr0qL9nxCKN2l3e2vyVUbZD9UvGEmEOuR1en37tD2Cwv6ppDdYPRAnQunnPjOLyE4/qsG v9oNd/pnMIIEqTCCA5GgAwIBAgICdy0wDQYJKoZIhvcNAQEFBQAwNTELMAkGA1UEBhMCRlIxDTAL BgNVBAoTBENOUlMxFzAVBgNVBAMTDkNOUlMyLVN0YW5kYXJkMB4XDTEzMTIxNjE0NTgyNVoXDTE1 MTIxNjE0NTgyNVowZzELMAkGA1UEBhMCRlIxDTALBgNVBAoTBENOUlMxEDAOBgNVBAsTB1VQUjMy NTExEzARBgNVBAMTCkJyaWFuIEthdHoxIjAgBgkqhkiG9w0BCQEWE2JyaWFuLmthdHpAbGltc2ku ZnIwggEiMA0GCSqGSIb3DQEBAQUAA4IBDwAwggEKAoIBAQDogOd9JNpls4A+x5qzEterjFSY2bGu 0BX2r3dsVhY0DB2cg9ovZmTFR/Ro+QjGPWZHDn6eY3O5HTSMLR/Mh7zVR+o9dD+IGKr5W12HQShD XduRDdZgFOwAX1skYhKZ1C7TdYiShnxl/hAiHY+IXvpHaGJlpIYHjD0+pMJ8y1uTDjuskNHOBC3l XhvMI8hpE+NleJxJLZsmGZwW41IrDnsuhiagv3nOGum4dEY2w2MfJzS8NHuQApY74NseIkrx/wfC 2VnnoWkQoyC0C3qBem90IqVl0Nnu2kvqZKFG1/rJFQHSr76z9el02WGxaC1pP982kUenPd1CT8/5 aYPRH1j5AgMBAAGjggGPMIIBizAMBgNVHRMBAf8EAjAAMBEGCWCGSAGG+EIBAQQEAwIEsDAOBgNV HQ8BAf8EBAMCBeAwegYJYIZIAYb4QgENBG0Wa0NlcnRpZmljYXQgQ05SUzItU3RhbmRhcmQuIFBv dXIgdG91dGUgaW5mb3JtYXRpb24gc2UgcmVwb3J0ZXIg4CBodHRwOi8vaWdjLnNlcnZpY2VzLmNu cnMuZnIvQ05SUzItU3RhbmRhcmQvMB0GA1UdDgQWBBS0jqVYl1rSaoQ9UONWYrVIZ4tZrTBUBgNV HSMETTBLgBQR49nRUkcbWbE8G3hma/ShiO0KW6EwpC4wLDELMAkGA1UEBhMCRlIxDTALBgNVBAoT BENOUlMxDjAMBgNVBAMTBUNOUlMyggEDMB4GA1UdEQQXMBWBE2JyaWFuLmthdHpAbGltc2kuZnIw RwYDVR0fBEAwPjA8oDqgOIY2aHR0cDovL2NybHMuc2VydmljZXMuY25ycy5mci9DTlJTMi1TdGFu ZGFyZC9nZXRkZXIuY3JsMA0GCSqGSIb3DQEBBQUAA4IBAQCHKbX7H/CDt3W23zz+LXBrb0JX6LGF C6AeCRblVwh2vXj+HE08GCuh6BW+QDk6POnHTg7oWgINtmorqIw0yKO8UGOwsKvBmsbq/Z3qfLti ewnz3AKYoS6PmAHmViXooij9YItH8QKyorpmzw0NLLMdQ5KI4aKj1XNQAQgYvuo9ML4O4gK3vv/K LnIwrT3YjPLloYHMEaCOkfl9HnWrXYMPg9HsUG1gty+JBTHQmYOS9u9QbKZAjyq9j2jSIL+cpU7i U6r3k3sETHOaMwpffTfm2+gpdQ8KO/UCtyY9PAlnJHz71w/6iULRvnLS14SjwOaGA99iT6neEd3j 9gkLWXW6MYIDCzCCAwcCAQEwOzA1MQswCQYDVQQGEwJGUjENMAsGA1UEChMEQ05SUzEXMBUGA1UE AxMOQ05SUzItU3RhbmRhcmQCAnctMAkGBSsOAwIaBQCgggGlMBgGCSqGSIb3DQEJAzELBgkqhkiG 9w0BBwEwHAYJKoZIhvcNAQkFMQ8XDTE1MDcyNzE2NDYwNVowIwYJKoZIhvcNAQkEMRYEFKVlgAQO 53RPE6CHZtBAeqZAAYObMEoGCSsGAQQBgjcQBDE9MDswNTELMAkGA1UEBhMCRlIxDTALBgNVBAoT BENOUlMxFzAVBgNVBAMTDkNOUlMyLVN0YW5kYXJkAgJ3LTBMBgsqhkiG9w0BCRACCzE9oDswNTEL MAkGA1UEBhMCRlIxDTALBgNVBAoTBENOUlMxFzAVBgNVBAMTDkNOUlMyLVN0YW5kYXJkAgJ3LTCB qwYJKoZIhvcNAQkPMYGdMIGaMAsGCWCGSAFlAwQBKjALBglghkgBZQMEARYwCgYIKoZIhvcNAwcw CwYJYIZIAWUDBAECMA4GCCqGSIb3DQMCAgIAgDAHBgUrDgMCBzANBggqhkiG9w0DAgIBQDANBggq hkiG9w0DAgIBKDAHBgUrDgMCGjALBglghkgBZQMEAgMwCwYJYIZIAWUDBAICMAsGCWCGSAFlAwQC ATANBgkqhkiG9w0BAQEFAASCAQBnABxw4hatQj8cEsU3NsSysq2xpM7Hx7Iwp7raHGRT5VMJIqjU 6zluqiTOqwQHtVlRtc2gLyd5yZbqRD6V83CicLKEe5f7QBJTbTJiRXtCBKr91IytlfrMJLOCtp0i aErY5SItIIxFaf9wpMklcv+P68LTueNlC/88dBMn0Egkv1942Fq6xvZgCapPh7vJZfHUnWNnNmT2 Q2X37aFJkpxPf+EaafPi+ir/T9eR7oogCw+2e0M6eLzK5DqY7CfxXja18ZJ6D1jiCVPGlaPOKmMa qARuDoYcqMEN8LIlFoX5NZS217ifZAT0I/6MdwKu0cUdYOjdIYs2EclKJbRLqs9gAAAAAAAA ------=_NextPart_000_0066_01D0C89C.7E7A6BF0--


This message came from the mail archive
/var/www/postings/2015/
maintained by:
DAn Ellis <dpwe@ee.columbia.edu>
Electrical Engineering Dept., Columbia University