<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JSIP</journal-id><journal-title-group><journal-title>Journal of Signal and Information Processing</journal-title></journal-title-group><issn pub-type="epub">2159-4465</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jsip.2014.51003</article-id><article-id pub-id-type="publisher-id">JSIP-42874</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject></subj-group></article-categories><title-group><article-title>
 
 
  The Regularized Low Pass Filter
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>eidong</surname><given-names>Chen</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Math Department, University of Georgia, Athens, GA 30602, USA</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>wchen@math.uga.edu</email></corresp></author-notes><pub-date pub-type="epub"><day>10</day><month>02</month><year>2014</year></pub-date><volume>05</volume><issue>01</issue><fpage>14</fpage><lpage>16</lpage><history><date date-type="received"><day>July</day>	<month>14th,</month>	<year>2013</year></date><date date-type="rev-recd"><day>August</day>	<month>17th,</month>	<year>2013</year>	</date><date date-type="accepted"><day>August</day>	<month>21st,</month>	<year>2013</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   In this paper the low pass filter is discussed in the noisy case. And a regularized low pass filter is presented. The convergence property of the regularized low pass filtering algorithm is proved in theory and tested by numerical results.
      
     
 
</p></abstract><kwd-group><kwd>Noise; Low Pass Filter; Regularization</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Filtering is widely applied in engineering [1-5]. In this paper, the problem of the low pass filtering is analyzed in theory and by examples in detail. A regularized low pass filtering algorithm is presented with the proof of the convergence property and numerical results.</p><p>First, we describe the band-limited signals<img src="3-3400301x\53848be1-41f3-46de-95ce-03bb49bc9c0d.jpg" />. The details can be seen in [<xref ref-type="bibr" rid="scirp.42874-ref6">6</xref>].</p><p>Definition: A function <img src="3-3400301x\1c233cb5-72a3-4d6c-976d-a2572684c264.jpg" /> is said to be <img src="3-3400301x\b5e59711-2e9f-4e3c-af66-6dba839e17fa.jpg" />band-limited if</p><p><img src="3-3400301x\6c2ddfdf-752d-447b-9a13-d68b82e08e0c.jpg" /></p><p>Here <img src="3-3400301x\ad69f34c-68c7-4e6a-915e-03bfd4dfa124.jpg" /> is the Fourier transform of<img src="3-3400301x\3ce8483e-1ed5-429c-8956-f7d085f7f07c.jpg" />:</p><disp-formula id="scirp.42874-formula80706"><label>(1)</label><graphic position="anchor" xlink:href="3-3400301x\fa03e488-ef1e-4ce3-a319-ea16f086b058.jpg"  xlink:type="simple"/></disp-formula><p>We then have the inversion formula:</p><p><img src="3-3400301x\c1208d9e-2719-4447-9e4f-b6293fb80e2e.jpg" /></p><p>In many practical problems, the signal <img src="3-3400301x\b913ab02-272e-42de-8ad3-5e78f2ad721f.jpg" /> is noisy:</p><disp-formula id="scirp.42874-formula80707"><label>(2)</label><graphic position="anchor" xlink:href="3-3400301x\f1609606-e2af-4def-95fc-ba72eddf5a5c.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="3-3400301x\958fd800-c7b3-407f-8f86-8525354f49a5.jpg" /> is the noise</p><disp-formula id="scirp.42874-formula80708"><label>(3)</label><graphic position="anchor" xlink:href="3-3400301x\b7825b60-e831-4f38-980d-eff0ab728a68.jpg"  xlink:type="simple"/></disp-formula><p>and <img src="3-3400301x\5ec45ee5-1e17-4eab-ae39-669a5f6f4436.jpg" /> is the exact band-limited signal.</p><p>In this paper, we will consider the problem low pass filtering:</p><disp-formula id="scirp.42874-formula80709"><label>(4)</label><graphic position="anchor" xlink:href="3-3400301x\cd72bf1c-cea3-4209-9c33-1f7459e5badf.jpg"  xlink:type="simple"/></disp-formula><p>If the signal is noisy however, the filter is not reliable. We will give an example to show that the noise can become very large after the low pass filtering process. So this filter is not reliable in the noisy case. And a regularized low pass filtering algorithm will be presented.</p><p>In section 2 we give the property of the low pass filter. A regularized filtering algorithm and the proof of its convergence are in section 3. The numerical results of some examples are given in section 4. Finally, the conclusion is given in section 5.</p></sec><sec id="s2"><title>2. The Property of the Low Pass Filter</title><p>In this section, we discuss the property of the low pass filter.</p><p>Example. Assume the noise is <img src="3-3400301x\f17a1f72-97a6-4a77-a84b-a1b420919e72.jpg" /> where <img src="3-3400301x\375cbf49-9721-4023-aeca-e4654cdcb955.jpg" /> is a given point in the time domain and <img src="3-3400301x\e3f2e57b-2593-473f-9a05-b640357d4be3.jpg" /> is close to zero. Then the noise signal after the filtering is</p><p><img src="3-3400301x\9d8730ab-e032-48b5-abe2-f0282b5e34a2.jpg" /></p><p>We can see that<img src="3-3400301x\e4d7f4f1-a35d-483a-8441-9a10a4629d41.jpg" />. However, the noise at <img src="3-3400301x\61ea745f-13ff-4082-98a6-3cd642bbd20b.jpg" /> after the filtering is</p><p><img src="3-3400301x\077c9de1-6a2b-4cb9-8425-fa27358f2b8a.jpg" /></p><p>Also at any point<img src="3-3400301x\471fe43c-11a9-4878-afe1-d8ef92cf0f9c.jpg" />, <img src="3-3400301x\76b424d2-4912-4bd3-b673-b9626d138427.jpg" /></p><p><img src="3-3400301x\db76866e-bcf4-48b7-a4bf-8aca1654fe19.jpg" /></p><p>So the error after the filtering becomes<img src="3-3400301x\deeb2f44-5e69-4e17-bca6-da4e2c57308c.jpg" />.</p><p>Remark. This is only an example for analysis. In the section of numerical results we will show that the low pass filter (4) is not very effective for white noise.</p></sec><sec id="s3"><title>3. The Regularized Filtering Algorithm</title><p>First, we consider the regularized Fourier transform [<xref ref-type="bibr" rid="scirp.42874-ref7">7</xref>]:</p><p><img src="3-3400301x\302ef56b-1964-4d90-b972-0990fd33ba43.jpg" /></p><p>where <img src="3-3400301x\33ab266a-96cd-4854-9178-2d6b8da843d3.jpg" /> is the regularization parameter. Here <img src="3-3400301x\822aef0e-1905-4776-bffb-b2f3e1c73b92.jpg" /> is the minimizer of a smoothing functional. We have proved <img src="3-3400301x\099dde9f-f986-4cd0-94a9-ea8555acb9d3.jpg" /> converges to the exact Fourier transform as the error of <img src="3-3400301x\0fd39dcc-1296-4d6a-bfc1-3a368cdcc1aa.jpg" /> approaches to zero. In [<xref ref-type="bibr" rid="scirp.42874-ref7">7</xref>], we have successfully used the regularized Fourier transform in extrapolation. So the weight function</p><p><img src="3-3400301x\c62fc3d0-2acf-4e43-8b8b-c5da4b94671f.jpg" /></p><p>is helpful to solve ill-posed problems.</p><p>Based on the regularized Fourier transform we present the regularized filtering formula:</p><disp-formula id="scirp.42874-formula80710"><label>(5)</label><graphic position="anchor" xlink:href="3-3400301x\8d59a5b6-854d-4c9a-b8fe-cee0309c6f96.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="3-3400301x\894260e2-3ea8-465d-906b-58b8abc066ae.jpg" /> is given in (2).</p><p>The convergence property of this regularized filtering formula is given in the theorem below.</p><p>Theorem 3.1. For<img src="3-3400301x\59193e8a-d098-4efa-8cee-9b4b303b3917.jpg" />, if <img src="3-3400301x\f5b3374f-8309-4cf4-89ab-7d2eb2ea1a39.jpg" /> and <img src="3-3400301x\3e461261-d82a-4ab1-9287-8c0adc05769c.jpg" /> as<img src="3-3400301x\af9a740e-162b-4ca4-b5b5-c761c2c12351.jpg" />, then <img src="3-3400301x\9666ac57-829d-4def-a3b0-e916ad948860.jpg" /> according to the maximum norm as<img src="3-3400301x\37bc7596-11a6-405e-86fb-7e43d64502be.jpg" />.</p><p><img src="3-3400301x\7dda7b92-a36c-4f9b-87d3-8a0ad779f3a7.jpg" />Proof.</p><p><img src="3-3400301x\5b88fcbd-e439-4fc6-b699-5359214ff84b.jpg" /></p><p>where</p><p><img src="3-3400301x\b354952e-f841-4748-9d1a-14b350f1bd63.jpg" /></p><p>For each<img src="3-3400301x\05003842-cae6-4721-924c-542fdd5285e4.jpg" />, there exists <img src="3-3400301x\34f55767-5396-48db-b166-2e7342aadf56.jpg" /> such that</p><p><img src="3-3400301x\4f9b2671-cb30-46d2-be99-7c4ef2dc8483.jpg" /></p><p>Then</p><p><img src="3-3400301x\4fab5b08-3f20-4bb6-9c96-3d278c258910.jpg" /></p><p>where</p><p><img src="3-3400301x\e0229d28-3570-4bf5-a784-bcfae81962bd.jpg" /></p><p>as<img src="3-3400301x\b4ba284a-713d-4cb0-b2fd-3e1bdd9442ba.jpg" />.</p></sec><sec id="s4"><title>4. Numerical Results</title><p>In this section, we give some examples to show that the regularized filtering algorithm (5) is more effective in reducing the noise than the convolution (4).</p><p>Suppose the exact signal in example 1 and 2 is</p><p><img src="3-3400301x\330172b4-0768-440b-988d-c1824162cc6e.jpg" /></p><p>Then construct</p><p><img src="3-3400301x\a971137c-1f37-4b3c-ac5f-d73a8157eb35.jpg" /></p><p>where<img src="3-3400301x\b087ce2f-7675-4bf3-a29d-0ee4dbfa8d5b.jpg" />.</p><p>Example 1. We consider the noise</p><p><img src="3-3400301x\dc7151fe-50de-4383-8f6b-9d0c8dba72b9.jpg" /></p><p>where<img src="3-3400301x\50c56d6e-833a-48dc-b089-3ba09bda5759.jpg" />, <img src="3-3400301x\3687f311-4900-46c1-90a1-546e0c71bf26.jpg" />, and<img src="3-3400301x\c5ff3f29-f48e-422e-ac3c-9b22c964fa93.jpg" />. This noise is used in the analysis of the stability in Section 2.</p><p>The result of (4) and the result of the regularized filtering algorithm with <img src="3-3400301x\0aa9ada0-a0a6-4d6f-8ee7-0d39232cedf8.jpg" /> are in figure 1.</p><p>Example 2. We consider the noise to be white noise that is Gauss distribution whose variance is 0.01. The result of (4) and the result of the regularized sampling algorithm with <img src="3-3400301x\7ec3e29b-3bb9-4076-bc43-eb5b27452188.jpg" /> are in figure 2.</p></sec><sec id="s5"><title>5. Conclusion</title><p>The filter of convolution with sinc function is not stable. For some noises the results of the filtering are not reliable.</p><p>Regularized filtering algorithm is more effective in reducing the noise.</p></sec><sec id="s6"><title>Acknowledgements</title><p>I would like thank University of Georgia for supporting my post doctoral work.</p></sec><sec id="s7"><title>REFERENCES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.42874-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">D. Georgakopoulos and W. Q. Yang, “Circuit Noise Reduction by Analogue Lowpass Filtering and Data Averaging,” Electronics Letters, Vol. 37 No. 19, 2001, pp. 1147-1148. http://dx.doi.org/10.1049/el:20010782</mixed-citation></ref><ref id="scirp.42874-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">P. Daripa and W. Hua, “A Numerical Study of an Ill-Posed Boussinesq Equation Arising in Water Waves and Nonlinear Lattices: Filtering and Regularization Techniques,” Applied Mathematics and Computation, Vol. 101, No. 2-3, 1999, pp. 159-207. http://dx.doi.org/10.1016/S0096-3003(98)10070-X</mixed-citation></ref><ref id="scirp.42874-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">O. Scherzer and J. Weickert, “Relations between Regularization and Diffusion Filtering,” Journal of Mathematical Imaging and Vision, Vol. 12, No. 1, 2000, pp. 43-63. http://dx.doi.org/10.1023/A:1008344608808</mixed-citation></ref><ref id="scirp.42874-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">T. Daboczi and I. Kollarl, “Multiparameter Optimization of Inverse Filtering Algorithms,” IEEE Transactions On Instrumentation and Measurement, Vol. 45, No. 2, 1996, pp. 417-421. http://dx.doi.org/10.1109/19.492758</mixed-citation></ref><ref id="scirp.42874-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">D. Tschumperle and R. Deriche, “Vector-Valued Image Regularization with PDEs: A Common Framework for Different Applications,” IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol. 27, No. 4, 2005, pp. 506-517. http://dx.doi.org/10.1109/TPAMI.2005.87</mixed-citation></ref><ref id="scirp.42874-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">W. Chen, “An Efficient Method for An Ill-Posed Problem —Band-limited Extrapolation by Regularization,” IEEE Transaction on Signal Processing, Vol. 54, No. 12, 2006, pp. 4611-4618. http://dx.doi.org/10.1109/TSP.2006.881255</mixed-citation></ref><ref id="scirp.42874-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">W. Chen, “Band-Limited Extrapolation by Regularization,” LAP Lambert Academic Publishing, 2010.</mixed-citation></ref></ref-list></back></article>