Introduction

Linear system to be measured depicts a linear system characterized by an impulse response h(n), driven by an input signal s(n), and producing the output signal r(n). The system identification problem is to estimate h(n) given known input/output signals s(n) and r(n). A practical method for identifying finite impulse responses is the Golay-code measurement technique, described below. Maximum length sequences may be alternatively used for this type of measurement, as they are also noise signals consisting of 1's and -1's, but perfectly inverse filtering the measurement is more computationally intensive.

Golay Code Theory The length L bilevel sequences a(n) and b(n) are Golay if and only if the following condition holds, where ¤ denotes the autocorrelation operator : a ( n ) ¤ a ( n ) + b ( n ) ¤ b ( n ) = 2 L δ ( n ) and δ(n) is the Kronecker delta function. Recall that () can also be written using *, the convolution operator: a ( - n ) * a ( n ) + b ( - n ) * b ( n ) = 2 L δ ( n ) Given that aL(n) and bL(n) are Golay, it turns out that a2L(n)=[aL(n)bL(n)] and b2L(n)=[aL(n)-bL(n)] are also Golay. This means that Golay sequences can be constructed recursively given Golay seed sequences such as a2(n)=[11] and b2(n)=[1-1]. See the MATLAB/Octave source code generate_golay.m for details. Notice also that the resulting bilevel sequences consist of only 1's and -1's. This means that the signal contains the maximum possible power level given that |s(n)|≤1∀n. This property is helpful for minimizing measurement noise. Let ra(n)=a(n)*h(n) be the response due to the Golay code input a(n), and let rb(n)=b(n)*h(n) be the response due to the Golay code input b(n). Due to (), the impulse response h(n) may be determined as follows: h ( n ) = 1 2 L ( a ( n ) ¤ r a ( n ) + b ( n ) ¤ r b ( n ) ) See golay_response.m for more details.

Golay Code Measurement Procedure and Software Generate the Golay codes golayA.wav and golayB.wav using generate_golay.m. Open the pd patch golay.pd, in pd. Ensure that the patch is not in editing mode, and check the “compute audio” box in the main pd window. Adjust the “Output Volume” so that when you click on “Record Response to Golay A”, the system under test is behaving linearly (i.e. not clipping), but so that the input signal to the sound interface is not too noisy. If there is an input volume on the sound interface, adjust it so that the levels approximately match those shown in when you click on “Record Response to Golay A” and “Record Response to Golay B.” If the sound interface has no input volume, then you will need to adjust the “Output Volume” accordingly.

golay.pd after making a measurement with an appropriate input level Once you are satisfied with the results, click the “Write Responses to Disk” button. pd will write the files RespA.wav and RespB.wav to disk. Rename these files so that the names match the measurement you just made. For instance, you might rename them to hpfRespA.wav and hpfRespB.wav if they corresponded to the measurement of a high-pass filter. Run golay_response('hpf') in MATLAB or Octave to analyze the measured response. Plots will be created, and the file hpfImpResp.wav will be written to disk.