<?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">AM</journal-id><journal-title-group><journal-title>Applied Mathematics</journal-title></journal-title-group><issn pub-type="epub">2152-7385</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/am.2014.54063</article-id><article-id pub-id-type="publisher-id">AM-43572</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Continuous Piecewise Linear Approximation of BV Function
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>ua</surname><given-names>Yi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Tao</surname><given-names>Yu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Zhiquan</surname><given-names>Chen</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Jingwen</surname><given-names>Zhu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>School of Mathematics and Physics, Jinggangshan University, Ji’an, China</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>yihua@whu.edu.cn(UY)</email>;<email>yutao@jgsu.edu.cn(TY)</email>;<email>zhujingwen666@163.com(JZ)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>10</day><month>03</month><year>2014</year></pub-date><volume>05</volume><issue>04</issue><fpage>667</fpage><lpage>671</lpage><history><date date-type="received"><day>19</day>	<month>December</month>	<year>2013</year></date><date date-type="rev-recd"><day>19</day>	<month>January</month>	<year>2014</year>	</date><date date-type="accepted"><day>27</day>	<month>January</month>	<year>2014</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>
 
 
   Nonlinear approximation is widely used in signal processing. Real-life signals can be modeled as functions of bounded variation. Thus the variable knot of approximating function could be self- adaptively chosen by balancing the total variation of the target function. In this paper, we adopt continuous piecewise linear approximation instead of the existing piecewise constants approximation. The results of experiments show that this new method is superior to the old one. 
 
</p></abstract><kwd-group><kwd>Nonlinear Approximation; Bounded Variation; Continuous Piecewise Linear Basis Function; Piecewise Constant Basis Function; Compression; Denoising</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The fundamental problem of approximation theory is to resolve a possibly complicated function, called the target function, by simpler, easier to compute functions called the approximants [<xref ref-type="bibr" rid="scirp.43572-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.43572-ref2">2</xref>] . The early methods utilized approximation from finite-dimensional linear spaces. In the beginning, these were typically spaces of polynomials, both algebraic and trigonometric. It was noted shortly thereafter that there were some advantages to be gained by not limiting the approximations to come from linear spaces [<xref ref-type="bibr" rid="scirp.43572-ref2">2</xref>] . In nonlinear approximation, the approximating function is not restricted to come from spaces of piece wise polynomials with a fixed partition; rather, the partition was allowed to depend on the target function. In principle, the idea was simple: we should use a finer mesh where the target function is not very smooth (singular) and a coarser mesh where it is smooth. Thus, an important question in nonlinear approximation is how we should measure this smoothness in order to obtain definitive results.</p><p>Kahane (1961) [<xref ref-type="bibr" rid="scirp.43572-ref3">3</xref>] gave a result which was concerned with how to obtain the variable knot of function in bounded variation space. Zhang (2009) [<xref ref-type="bibr" rid="scirp.43572-ref4">4</xref>] gave the corresponding numerical experiments by using the piecewise constants function. Moreover, the advantages of Zhang’s method over other denoising methods, such as Visushrink [<xref ref-type="bibr" rid="scirp.43572-ref5">5</xref>] and SureShrinkage [<xref ref-type="bibr" rid="scirp.43572-ref6">6</xref>] were also analyzed [<xref ref-type="bibr" rid="scirp.43572-ref4">4</xref>] .</p><p>In this paper, we use piecewise linear basis functions instead of the piecewise constants basis functions adopted by Kahane, which is motivated by continuous piecewise linear approximation studied by Y. Shi (2010) [<xref ref-type="bibr" rid="scirp.43572-ref7">7</xref>] . The experimental results of this new method are superior to those of the old one.</p></sec><sec id="s2"><title>2. The Method for the Selection of Variable Knots and Continuous Piecewise Linear Approximation</title><p>We can consider that a real-life signal f is defined on <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\406bf40a-6f7e-4f8d-8ec2-8f7eb8bfc187.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\6b50a5fe-c98e-4b4a-9f36-1f5e27168b3d.png" xlink:type="simple"/></inline-formula>. In order to self-adaptively select the knots according to the target function, we can balance the total variation (TV) of this function [<xref ref-type="bibr" rid="scirp.43572-ref4">4</xref>] . The TV of f is given by</p><disp-formula id="scirp.43572-formula152071"><label>(1.1)</label><graphic position="anchor" xlink:href="htmlimages\9-7402044x\6f0ae381-1ced-40c8-a2ba-caf2449563d2.png"  xlink:type="simple"/></disp-formula><p>where the supremum is taken over the set of all partitions <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\7d6042ab-9a08-4559-92d2-0a22fe74b1ce.png" xlink:type="simple"/></inline-formula> of<inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\40d8e7f7-0243-44ee-97cb-e0c1e169e4b4.png" xlink:type="simple"/></inline-formula>. The TV function <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\03ddace8-eae9-423d-b7ce-2675a1c892f4.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\bd40a1da-64f9-4cd7-b825-a5ff0ab1c27d.png" xlink:type="simple"/></inline-formula> is non-negative monotone continuous function on <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\06a44e2e-9f86-48f9-818a-eb2aa3593526.png" xlink:type="simple"/></inline-formula> with<inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\7843ca61-c5b6-47c8-9c8b-08e6d9a7793f.png" xlink:type="simple"/></inline-formula>. We divide the range of <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\8c2a041b-53fd-4dd4-9a72-f2362cec5bb5.png" xlink:type="simple"/></inline-formula> (which is the interval<inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\63088c87-755f-4cc7-aac1-7d8caca5f4e7.png" xlink:type="simple"/></inline-formula>) on the y-axis into n pieces corresponding to the y values<inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\e5811332-1a02-4ca9-befa-656519d54b36.png" xlink:type="simple"/></inline-formula>. We denote the preimage of these points by<inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\0615cc8d-b9a2-4f48-9268-1dbb3f3d9311.png" xlink:type="simple"/></inline-formula>, we have</p><disp-formula id="scirp.43572-formula152072"><label>(1.2)</label><graphic position="anchor" xlink:href="htmlimages\9-7402044x\3e0811bc-685a-4e37-b0e1-8f242f7c84a5.png"  xlink:type="simple"/></disp-formula><p>Hence, <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\b6f8200f-c117-4a12-b795-b6cda19e2b95.png" xlink:type="simple"/></inline-formula>forms the partition balancing the TV of <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\220004f7-27c1-476d-ac60-49ae2bc29ece.png" xlink:type="simple"/></inline-formula> over each interval<inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\c121d149-31c5-4bbb-b72e-022742aaf528.png" xlink:type="simple"/></inline-formula>. Then the target function can be approximated by piecewise constants [<xref ref-type="bibr" rid="scirp.43572-ref4">4</xref>] . That is,</p><disp-formula id="scirp.43572-formula152073"><label>(1.3)</label><graphic position="anchor" xlink:href="htmlimages\9-7402044x\9aeb3f3c-1454-4843-8523-c768c534e630.png"  xlink:type="simple"/></disp-formula><p>is the approximating function of<inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\db5c5739-123b-4c4f-8605-e2d0ba426b83.png" xlink:type="simple"/></inline-formula>. In this case, the<inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\118813b7-1044-4c3a-bb0b-d275e00d2230.png" xlink:type="simple"/></inline-formula>—error to <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\4e8e0592-b281-49c5-8723-d6a526181b5b.png" xlink:type="simple"/></inline-formula> is</p><disp-formula id="scirp.43572-formula152074"><label>(1.4)</label><graphic position="anchor" xlink:href="htmlimages\9-7402044x\a0365da5-5928-4ad0-adc2-bc6b304e1325.png"  xlink:type="simple"/></disp-formula><p>However, <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\c5d5a586-cc4c-4e93-b245-f76a09227cdb.png" xlink:type="simple"/></inline-formula>is not continuous. In order to approximate the target function by continuous function, we construct the continuous piecewise linear basis function as follows [<xref ref-type="bibr" rid="scirp.43572-ref7">7</xref>] :</p><disp-formula id="scirp.43572-formula152075"><label>(1.5)</label><graphic position="anchor" xlink:href="htmlimages\9-7402044x\55992c4d-48a2-4c59-805a-e7546255fb84.png"  xlink:type="simple"/></disp-formula><p>We can substitute <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\629c933b-e11f-4dd4-b0f5-0f49f8e27cd7.png" xlink:type="simple"/></inline-formula> for Haar scaling function(the piecewise constant function) to achieve continuous piecewise linear approximation, and the corresponding error analysis has been thoroughly studied in [<xref ref-type="bibr" rid="scirp.43572-ref7">7</xref>] .</p><p>We define</p><disp-formula id="scirp.43572-formula152076"><label>(1.6)</label><graphic position="anchor" xlink:href="htmlimages\9-7402044x\0f947a7a-812e-4af8-80b1-0392e5fb585b.png"  xlink:type="simple"/></disp-formula><p>where<inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\2ad9c671-1707-471e-9f63-051682197923.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\0bb99540-2760-4449-a703-920d63506230.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\d16bb89e-18d2-4b8a-a20c-ec7755d024f2.png" xlink:type="simple"/></inline-formula>,<inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\c7a39080-6ab5-4d28-a350-8fcab9cdba5b.png" xlink:type="simple"/></inline-formula>. Then, <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\42598006-af12-4db2-a406-57ed7b489dd2.png" xlink:type="simple"/></inline-formula>achieves the continuous piecewise linear approximation for the target function.</p></sec><sec id="s3"><title>3. Numerical Experiments</title><p>In this section, we give comparisons of methods in signal compression and denoising by two examples. The signal leleccum is composed of 1000 points, presented in <xref ref-type="fig" rid="fig1">Figure 1</xref>(a). TV function of the original function is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>(b). By dividing the range of the TV function on the y-axis into n pieces to get <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\5a8765c2-9531-4adf-ab6b-3fb3700e7525.png" xlink:type="simple"/></inline-formula> y-values, the preimage of these y-values form the self-adaptively knots of the target function. In this experiment of this paper, we let<inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\3c29ded0-e104-4ca9-acfd-b4a17aefbb10.png" xlink:type="simple"/></inline-formula>. So the original signal can be represented by a compressed version with 500 points. <xref ref-type="fig" rid="fig1">Figure 1</xref>(e) provides the compressed signal by Zhang’s method, i.e., using the piecewise constants approximation. <xref ref-type="fig" rid="fig1">Figure 1</xref>(c) presents the compressed signal by using continuous piecewise linear approximation. We use <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\21a9eaf6-14c3-440b-bd3f-4bf2a54b79e0.png" xlink:type="simple"/></inline-formula> norm to measure the error of compression. That is,</p><disp-formula id="scirp.43572-formula152077"><label>(1.7)</label><graphic position="anchor" xlink:href="htmlimages\9-7402044x\b602b5b4-951d-4a35-a712-c7941ba51fda.png"  xlink:type="simple"/></disp-formula><p>We see that the proposed method is superior to the old one as long as we notice that <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\d321974d-40b9-4446-8f46-310bda23eaa8.png" xlink:type="simple"/></inline-formula> for our method is less than <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\ef32251d-ebad-4ac8-b057-b6e8c43dc314.png" xlink:type="simple"/></inline-formula> for the old one. The compressed signal with <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\cf1b8075-38c1-453c-a81c-ffaaa028aef2.png" xlink:type="simple"/></inline-formula> by VisuShrink is also presented in <xref ref-type="fig" rid="fig1">Figure 1</xref>(d).</p><p>The second example is related to a natural noisy signal which is shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>(a). The restored signals by using our method and by using the old one are respectively shown in Figures 2(d) and (c). The advantages of our method over the old one consist of two aspects. Firstly, <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\80b73431-73f3-41fe-b4dc-8cc78dce8db1.png" xlink:type="simple"/></inline-formula>of the denoised signal in <xref ref-type="fig" rid="fig2">Figure 2</xref>(d) is larger than <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\dcb85f3a-33f9-48b2-b994-53a28bad5029.png" xlink:type="simple"/></inline-formula> in <xref ref-type="fig" rid="fig2">Figure 2</xref>(c). At this time, SNR of the denoised signal <inline-formula><inline-graphic xlink:href="tmlimages\9-7402044x\5dd02ff2-ccb0-4d55-b5f8-f77a2e0d9d0f.png" xlink:type="simple"/></inline-formula> is defined as</p><disp-formula id="scirp.43572-formula152078"><label>(1.8)</label><graphic position="anchor" xlink:href="htmlimages\9-7402044x\bc85fcc8-8e65-46d3-8e4b-777ef77d1ff4.png"  xlink:type="simple"/></disp-formula><p>Secondly, the signal in <xref ref-type="fig" rid="fig2">Figure 2</xref>(d) is more smooth. See the rectangle parts of Figures 2(c) and (d) for details. The sawtooth nature of denoised signal in <xref ref-type="fig" rid="fig2">Figure 2</xref>(c) may be caused by the shape of piecewise constant basis function.</p></sec><sec id="s4"><title>4. Conclusion</title><p>Nonlinear approximation is the theoretical foundation of compression and denoising of signals. Zhang has presented a general partition method to obtain the variable knots according to the target function. However, the drawback of Zhang’s method is that the smoothness of the restored signal is bad when n is small. In this paper, we use continuous piecewise linear basis function as a substitute for the piecewise constant basis function adopted by Zhang. Our method performs better than the old one both in terms of SNR and vision.</p></sec><sec id="s5"><title>Funding</title><p>The work was supported by jskjz [<xref ref-type="bibr" rid="scirp.43572-ref2012">2012</xref>] 32-8, [<xref ref-type="bibr" rid="scirp.43572-ref2012">2012</xref>] 32-7, National Natural Science Foundation of China under Grant No. 11347210, the Doctoral Starting up Foundation of Jinggangshan University (Grant No. JZB1304), the Natural Science Foundation of Jiangxi Province, China (Grant No. 20132BAB211018).</p></sec><sec id="s6"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.43572-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">DeVore, R.A. and Lorentz, G.G. (1993) Constructive Approximation. Volume 303. Springer, Berlin.</mixed-citation></ref><ref id="scirp.43572-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">DeVore, R.A. (1998) Nonlinear Approximation. In Acta Numerica, Cambridge University Press, Cambridge, 51-150.</mixed-citation></ref><ref id="scirp.43572-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Kahane, J.P. (1961) Teora Constructiva de Funciones, Volume 5. Universidad de Buenos Aires, Buenos Aires.</mixed-citation></ref><ref id="scirp.43572-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Zhang, T., Fan, Q. and Shu, H. (2009) Approximation of BV Function by Piecewise Constants and Its Application in Signal Denoising and Compression. 2nd International Congress on Image and Signal Processing, Tianjin, 17-19 October 2009, 1-5.</mixed-citation></ref><ref id="scirp.43572-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Donoho, D.L. and Johnstone, J.M. (1994) Ideal Spatial Adaptation by Wavelet Shrinkage. Biometrika, 81, 425-455. http://dx.doi.org/10.1093/biomet/81.3.425</mixed-citation></ref><ref id="scirp.43572-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Donoho, D.L. and Johnstone, I.M. (1995) Adapting to Unknown Smoothness via Wavelet Shrinkage. Journal of the American Statistical Association, 90, 1200-1224. http://dx.doi.org/10.1080/01621459.1995.10476626</mixed-citation></ref><ref id="scirp.43572-ref7"><label>7</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Shi Y.Y. and Li</surname><given-names> X.Y. </given-names></name>,<etal>et al</etal>. (<year>2010</year>)<article-title>Continuous Piecewise Linear Approximation Based on Haar Scale Transform</article-title><source> Journal of Liaoning Technical University (Natural Science)</source><volume> 3</volume>,<fpage> 521</fpage>-<lpage>524</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref></ref-list></back></article>