## Linear Filtering

Let be an -by-1 output vector from a linear filtering operation

where is an -by- filter matrix and is a real -by-1 input vector. Let be the -by-1 vector of squared values of the elements of . For example, let

 (8.10)

Then, we may write the elements of as

where is an -by- matrix of all zeros except a single one on the main diagonal in location .

We now test the SPA solution against a Gaussian mixture approximation. We used a filter length of 3 as shown in equation (8.10) and a feature size of . The filter coefficients were . A Gaussian mixture approximation of under the WGN assumption was obtained using 5000 samples and 20 mixture components. Next, 100 new samples were generated and the log PDF from the SPA was compared with the mixture approximation.

The results of the experiment are shown in Figure 8.2. On the X-axis is the output of software/pdf_quadspa.m, and on the Y-axis is the Gaussian mixture approximation. Light-blue points are computed using usual saddlepoint. Black points are comouted using saddlepoint plus correction term (next term in Taylor series). See software/module_quadspa_test.m with task='filter'. The correction term can be invoked with level = 1.

Baggenstoss 2017-05-19