Since Fig.6 indicates the existence of fixed and angle-dependent components in the measured impulse responses, and since such angle-independent component is strongly suppressed by baffling in the cabinet enclosure, it is desirable to eliminate this fixed component from the measurements. For this purpose, an iterative algorithm was developed which models the two components separately.
Let denote the number of impulse-response samples in each measured impulse response,and let denote the number of angles (-180:15:180) at which impulse-response measurements were taken. We denote the impulse-response matrix by . Each column of is an impulse response at some horn angle. (Figure 6 can be interpreted as a plot of the transpose of .)
We model as
where is the arrival-time delay, in samples, for the horn output in the th row (the delays clearly visible in Fig.6 as a function of angle). These arrival times are estimated as the location of the peak in the cross-correlation between the th impulse response and the same impulse response after converting it to minimum phase . The diagonal matrix diag denotes a shift operator which delays the th column of by samples. Thus, contains the horn-output impulse response (without the base leakage) shifted to time zero (i.e., the angle-dependent delay is removed). Finally, the error matrix is to be minimized in the least-squares sense.
Each column of the matrix contains a copy of the estimated horn-base leakage impulse-response:
The estimated angle-dependent impulse-responses in are modeled as linear combinations of fixed impulse responses, viewed (loosely) as principal components:
where is the orthonormal matrix of fixed filters (principal components), and is a matrix of weights, found in the usual way by a truncated singular value decomposition (SVD) .