Can we combine convolution and higher powers for locally maximising a function?

The name of the pictureThe name of the pictureThe name of the pictureClash Royale CLAN TAG#URR8PPP











up vote
0
down vote

favorite












Can we somehow find local maximum function value (for strictly positive functions) using a convolution?



My idea is based on the result that $$ lim_pto inftyleft[frac1Nsum_k=1^N (a_k) ^pright]^frac 1 p = max(a_k)$$



Similarly we might be able to use the result that



$$ min(a_k) = max(a_k)-max(max(a_k)-a_k)$$



Could we perhaps use this for building envelope detection? It seems theoretically sound to me, but would it be practically feasible?







share|cite|improve this question























    up vote
    0
    down vote

    favorite












    Can we somehow find local maximum function value (for strictly positive functions) using a convolution?



    My idea is based on the result that $$ lim_pto inftyleft[frac1Nsum_k=1^N (a_k) ^pright]^frac 1 p = max(a_k)$$



    Similarly we might be able to use the result that



    $$ min(a_k) = max(a_k)-max(max(a_k)-a_k)$$



    Could we perhaps use this for building envelope detection? It seems theoretically sound to me, but would it be practically feasible?







    share|cite|improve this question





















      up vote
      0
      down vote

      favorite









      up vote
      0
      down vote

      favorite











      Can we somehow find local maximum function value (for strictly positive functions) using a convolution?



      My idea is based on the result that $$ lim_pto inftyleft[frac1Nsum_k=1^N (a_k) ^pright]^frac 1 p = max(a_k)$$



      Similarly we might be able to use the result that



      $$ min(a_k) = max(a_k)-max(max(a_k)-a_k)$$



      Could we perhaps use this for building envelope detection? It seems theoretically sound to me, but would it be practically feasible?







      share|cite|improve this question











      Can we somehow find local maximum function value (for strictly positive functions) using a convolution?



      My idea is based on the result that $$ lim_pto inftyleft[frac1Nsum_k=1^N (a_k) ^pright]^frac 1 p = max(a_k)$$



      Similarly we might be able to use the result that



      $$ min(a_k) = max(a_k)-max(max(a_k)-a_k)$$



      Could we perhaps use this for building envelope detection? It seems theoretically sound to me, but would it be practically feasible?









      share|cite|improve this question










      share|cite|improve this question




      share|cite|improve this question









      asked 2 days ago









      mathreadler

      13.5k71857




      13.5k71857




















          1 Answer
          1






          active

          oldest

          votes

















          up vote
          0
          down vote













          Yes, we can do this.



          If we consider the function $$f(t)=casest$$



          We can with the convolution



          1. If we are in continuous time:

          $$(g*f)(t)=int_-infty^infty g(tau)cdot f(t-tau) dtau$$



          1. If we are in discrete time:
            $$(g*f)(t)=sum_tau=-infty^infty g(tau)cdot f(t-tau)$$

          Define the non-linear operation $$O(g,p) = (g^p*f)(t)^1/p$$
          Where exponent means multiplicative power.



          For example: if $g(t) = sin(t), O(g,4) = ((sin(x)^4*f))(t)^1/4$



          In practice we usually want much higher $p$ exponents, maybe 32 or 1024.



          We now investigate the test-function



          $$g(t) = sin(60t^2)^2cdot (sin(16t)^2 + cos(8t)^2)$$



          The square is to get same level of minimum values: 0.
          If we analyze $g(t)$, we can see that we have a chirp function modulated by a slower wave. This is to investigate behaviors for different frequencies.



          enter image description here



          We see that for higher frequencies we get reasonable envelope behavior. It is also possible to observe that as we change parameter $N$ above, which frequency range we can calculate maxima for are affected - but also the rectangular shape of local maxima will grow. We may need to smooth these out somehow. But how to do so is outside scope of question.






          share|cite|improve this answer























          • Can you explain what you have plotted here? The red "max-envelope" curve can't be simply $g*f$, since that would give you the average value in a window rather than the maximum.
            – Rahul
            yesterday










          • @Rahul yes you are right, I will correct.
            – mathreadler
            yesterday










          Your Answer




          StackExchange.ifUsing("editor", function ()
          return StackExchange.using("mathjaxEditing", function ()
          StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix)
          StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
          );
          );
          , "mathjax-editing");

          StackExchange.ready(function()
          var channelOptions =
          tags: "".split(" "),
          id: "69"
          ;
          initTagRenderer("".split(" "), "".split(" "), channelOptions);

          StackExchange.using("externalEditor", function()
          // Have to fire editor after snippets, if snippets enabled
          if (StackExchange.settings.snippets.snippetsEnabled)
          StackExchange.using("snippets", function()
          createEditor();
          );

          else
          createEditor();

          );

          function createEditor()
          StackExchange.prepareEditor(
          heartbeatType: 'answer',
          convertImagesToLinks: true,
          noModals: false,
          showLowRepImageUploadWarning: true,
          reputationToPostImages: 10,
          bindNavPrevention: true,
          postfix: "",
          noCode: true, onDemand: true,
          discardSelector: ".discard-answer"
          ,immediatelyShowMarkdownHelp:true
          );



          );








           

          draft saved


          draft discarded


















          StackExchange.ready(
          function ()
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2871877%2fcan-we-combine-convolution-and-higher-powers-for-locally-maximising-a-function%23new-answer', 'question_page');

          );

          Post as a guest






























          1 Answer
          1






          active

          oldest

          votes








          1 Answer
          1






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes








          up vote
          0
          down vote













          Yes, we can do this.



          If we consider the function $$f(t)=casest$$



          We can with the convolution



          1. If we are in continuous time:

          $$(g*f)(t)=int_-infty^infty g(tau)cdot f(t-tau) dtau$$



          1. If we are in discrete time:
            $$(g*f)(t)=sum_tau=-infty^infty g(tau)cdot f(t-tau)$$

          Define the non-linear operation $$O(g,p) = (g^p*f)(t)^1/p$$
          Where exponent means multiplicative power.



          For example: if $g(t) = sin(t), O(g,4) = ((sin(x)^4*f))(t)^1/4$



          In practice we usually want much higher $p$ exponents, maybe 32 or 1024.



          We now investigate the test-function



          $$g(t) = sin(60t^2)^2cdot (sin(16t)^2 + cos(8t)^2)$$



          The square is to get same level of minimum values: 0.
          If we analyze $g(t)$, we can see that we have a chirp function modulated by a slower wave. This is to investigate behaviors for different frequencies.



          enter image description here



          We see that for higher frequencies we get reasonable envelope behavior. It is also possible to observe that as we change parameter $N$ above, which frequency range we can calculate maxima for are affected - but also the rectangular shape of local maxima will grow. We may need to smooth these out somehow. But how to do so is outside scope of question.






          share|cite|improve this answer























          • Can you explain what you have plotted here? The red "max-envelope" curve can't be simply $g*f$, since that would give you the average value in a window rather than the maximum.
            – Rahul
            yesterday










          • @Rahul yes you are right, I will correct.
            – mathreadler
            yesterday














          up vote
          0
          down vote













          Yes, we can do this.



          If we consider the function $$f(t)=casest$$



          We can with the convolution



          1. If we are in continuous time:

          $$(g*f)(t)=int_-infty^infty g(tau)cdot f(t-tau) dtau$$



          1. If we are in discrete time:
            $$(g*f)(t)=sum_tau=-infty^infty g(tau)cdot f(t-tau)$$

          Define the non-linear operation $$O(g,p) = (g^p*f)(t)^1/p$$
          Where exponent means multiplicative power.



          For example: if $g(t) = sin(t), O(g,4) = ((sin(x)^4*f))(t)^1/4$



          In practice we usually want much higher $p$ exponents, maybe 32 or 1024.



          We now investigate the test-function



          $$g(t) = sin(60t^2)^2cdot (sin(16t)^2 + cos(8t)^2)$$



          The square is to get same level of minimum values: 0.
          If we analyze $g(t)$, we can see that we have a chirp function modulated by a slower wave. This is to investigate behaviors for different frequencies.



          enter image description here



          We see that for higher frequencies we get reasonable envelope behavior. It is also possible to observe that as we change parameter $N$ above, which frequency range we can calculate maxima for are affected - but also the rectangular shape of local maxima will grow. We may need to smooth these out somehow. But how to do so is outside scope of question.






          share|cite|improve this answer























          • Can you explain what you have plotted here? The red "max-envelope" curve can't be simply $g*f$, since that would give you the average value in a window rather than the maximum.
            – Rahul
            yesterday










          • @Rahul yes you are right, I will correct.
            – mathreadler
            yesterday












          up vote
          0
          down vote










          up vote
          0
          down vote









          Yes, we can do this.



          If we consider the function $$f(t)=casest$$



          We can with the convolution



          1. If we are in continuous time:

          $$(g*f)(t)=int_-infty^infty g(tau)cdot f(t-tau) dtau$$



          1. If we are in discrete time:
            $$(g*f)(t)=sum_tau=-infty^infty g(tau)cdot f(t-tau)$$

          Define the non-linear operation $$O(g,p) = (g^p*f)(t)^1/p$$
          Where exponent means multiplicative power.



          For example: if $g(t) = sin(t), O(g,4) = ((sin(x)^4*f))(t)^1/4$



          In practice we usually want much higher $p$ exponents, maybe 32 or 1024.



          We now investigate the test-function



          $$g(t) = sin(60t^2)^2cdot (sin(16t)^2 + cos(8t)^2)$$



          The square is to get same level of minimum values: 0.
          If we analyze $g(t)$, we can see that we have a chirp function modulated by a slower wave. This is to investigate behaviors for different frequencies.



          enter image description here



          We see that for higher frequencies we get reasonable envelope behavior. It is also possible to observe that as we change parameter $N$ above, which frequency range we can calculate maxima for are affected - but also the rectangular shape of local maxima will grow. We may need to smooth these out somehow. But how to do so is outside scope of question.






          share|cite|improve this answer















          Yes, we can do this.



          If we consider the function $$f(t)=casest$$



          We can with the convolution



          1. If we are in continuous time:

          $$(g*f)(t)=int_-infty^infty g(tau)cdot f(t-tau) dtau$$



          1. If we are in discrete time:
            $$(g*f)(t)=sum_tau=-infty^infty g(tau)cdot f(t-tau)$$

          Define the non-linear operation $$O(g,p) = (g^p*f)(t)^1/p$$
          Where exponent means multiplicative power.



          For example: if $g(t) = sin(t), O(g,4) = ((sin(x)^4*f))(t)^1/4$



          In practice we usually want much higher $p$ exponents, maybe 32 or 1024.



          We now investigate the test-function



          $$g(t) = sin(60t^2)^2cdot (sin(16t)^2 + cos(8t)^2)$$



          The square is to get same level of minimum values: 0.
          If we analyze $g(t)$, we can see that we have a chirp function modulated by a slower wave. This is to investigate behaviors for different frequencies.



          enter image description here



          We see that for higher frequencies we get reasonable envelope behavior. It is also possible to observe that as we change parameter $N$ above, which frequency range we can calculate maxima for are affected - but also the rectangular shape of local maxima will grow. We may need to smooth these out somehow. But how to do so is outside scope of question.







          share|cite|improve this answer















          share|cite|improve this answer



          share|cite|improve this answer








          edited yesterday


























          answered yesterday









          mathreadler

          13.5k71857




          13.5k71857











          • Can you explain what you have plotted here? The red "max-envelope" curve can't be simply $g*f$, since that would give you the average value in a window rather than the maximum.
            – Rahul
            yesterday










          • @Rahul yes you are right, I will correct.
            – mathreadler
            yesterday
















          • Can you explain what you have plotted here? The red "max-envelope" curve can't be simply $g*f$, since that would give you the average value in a window rather than the maximum.
            – Rahul
            yesterday










          • @Rahul yes you are right, I will correct.
            – mathreadler
            yesterday















          Can you explain what you have plotted here? The red "max-envelope" curve can't be simply $g*f$, since that would give you the average value in a window rather than the maximum.
          – Rahul
          yesterday




          Can you explain what you have plotted here? The red "max-envelope" curve can't be simply $g*f$, since that would give you the average value in a window rather than the maximum.
          – Rahul
          yesterday












          @Rahul yes you are right, I will correct.
          – mathreadler
          yesterday




          @Rahul yes you are right, I will correct.
          – mathreadler
          yesterday












           

          draft saved


          draft discarded


























           


          draft saved


          draft discarded














          StackExchange.ready(
          function ()
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2871877%2fcan-we-combine-convolution-and-higher-powers-for-locally-maximising-a-function%23new-answer', 'question_page');

          );

          Post as a guest













































































          Comments

          Popular posts from this blog

          What is the equation of a 3D cone with generalised tilt?

          Color the edges and diagonals of a regular polygon

          Relationship between determinant of matrix and determinant of adjoint?