envelope of rays that propagate into the ellipse for $F=u_x^2 + u_y^2 - 1$












0












$begingroup$


I have the equation $F=u_x^2 + u_y^2 - 1$ and that $u=0$ on the boundary of the ellipse $frac{x^2}{a^2} +frac{y^2}{b^2}=1$



At $u=0$ the ellipse can be parametrised as $x=acos(s), y=bsin(s)$



Charpits eqs give me:



$dx/dt= 2p$



$dy/dt = 2q$



$dp/dt = dq/dt = 0$



$du/dt = 2p^2 + 2q^2$



where $p = frac{partial u}{partial x}$ and $q = frac{partial u}{partial y}$



solving charpits eqs I get:



$p = p_0(s)$



$q = q_0(s)$



$x = 2p_0(s)t + acos(s)$



$y = 2q_0(s)t + bsin(s)$



Since $F(p,q) = p^2 + q^2 - 1$, this implies that $p_0^2 + q_0^2 = 1$



To find the envelope i need to find when the jacobian =0 which is when $2bp_0cos(s)+ 2aq_0sin(s) = 0$



I am confused as to how to draw the envelope of rays that propagate into the ellipse and also how to determine the ridge of discontinuity






ordinary-differential-equations pde







share|cite|improve this question











$endgroup$

















    0












    $begingroup$


    I have the equation $F=u_x^2 + u_y^2 - 1$ and that $u=0$ on the boundary of the ellipse $frac{x^2}{a^2} +frac{y^2}{b^2}=1$



    At $u=0$ the ellipse can be parametrised as $x=acos(s), y=bsin(s)$



    Charpits eqs give me:



    $dx/dt= 2p$



    $dy/dt = 2q$



    $dp/dt = dq/dt = 0$



    $du/dt = 2p^2 + 2q^2$



    where $p = frac{partial u}{partial x}$ and $q = frac{partial u}{partial y}$



    solving charpits eqs I get:



    $p = p_0(s)$



    $q = q_0(s)$



    $x = 2p_0(s)t + acos(s)$



    $y = 2q_0(s)t + bsin(s)$



    Since $F(p,q) = p^2 + q^2 - 1$, this implies that $p_0^2 + q_0^2 = 1$



    To find the envelope i need to find when the jacobian =0 which is when $2bp_0cos(s)+ 2aq_0sin(s) = 0$



    I am confused as to how to draw the envelope of rays that propagate into the ellipse and also how to determine the ridge of discontinuity






    ordinary-differential-equations pde







    share|cite|improve this question











    $endgroup$















      0












      0








      0





      $begingroup$


      I have the equation $F=u_x^2 + u_y^2 - 1$ and that $u=0$ on the boundary of the ellipse $frac{x^2}{a^2} +frac{y^2}{b^2}=1$



      At $u=0$ the ellipse can be parametrised as $x=acos(s), y=bsin(s)$



      Charpits eqs give me:



      $dx/dt= 2p$



      $dy/dt = 2q$



      $dp/dt = dq/dt = 0$



      $du/dt = 2p^2 + 2q^2$



      where $p = frac{partial u}{partial x}$ and $q = frac{partial u}{partial y}$



      solving charpits eqs I get:



      $p = p_0(s)$



      $q = q_0(s)$



      $x = 2p_0(s)t + acos(s)$



      $y = 2q_0(s)t + bsin(s)$



      Since $F(p,q) = p^2 + q^2 - 1$, this implies that $p_0^2 + q_0^2 = 1$



      To find the envelope i need to find when the jacobian =0 which is when $2bp_0cos(s)+ 2aq_0sin(s) = 0$



      I am confused as to how to draw the envelope of rays that propagate into the ellipse and also how to determine the ridge of discontinuity






      ordinary-differential-equations pde







      share|cite|improve this question











      $endgroup$




      I have the equation $F=u_x^2 + u_y^2 - 1$ and that $u=0$ on the boundary of the ellipse $frac{x^2}{a^2} +frac{y^2}{b^2}=1$



      At $u=0$ the ellipse can be parametrised as $x=acos(s), y=bsin(s)$



      Charpits eqs give me:



      $dx/dt= 2p$



      $dy/dt = 2q$



      $dp/dt = dq/dt = 0$



      $du/dt = 2p^2 + 2q^2$



      where $p = frac{partial u}{partial x}$ and $q = frac{partial u}{partial y}$



      solving charpits eqs I get:



      $p = p_0(s)$



      $q = q_0(s)$



      $x = 2p_0(s)t + acos(s)$



      $y = 2q_0(s)t + bsin(s)$



      Since $F(p,q) = p^2 + q^2 - 1$, this implies that $p_0^2 + q_0^2 = 1$



      To find the envelope i need to find when the jacobian =0 which is when $2bp_0cos(s)+ 2aq_0sin(s) = 0$



      I am confused as to how to draw the envelope of rays that propagate into the ellipse and also how to determine the ridge of discontinuity






      ordinary-differential-equations pde




      ordinary-differential-equations pde






      share|cite|improve this question















      share|cite|improve this question













      share|cite|improve this question




      share|cite|improve this question








      edited Dec 10 '18 at 16:58









      Harry49

      6,17331132




      6,17331132










      asked Dec 10 '18 at 12:11









      pablo_mathscobarpablo_mathscobar

      996




      996






















          1 Answer
          1






          active

          oldest

          votes


















          0












          $begingroup$

          I think you might have missed a sign in $dot x_0(s)=-asin(s)$. Then the equation for $p_0,q_0$ reads as
          $$
          0=dot z_0(s)=p_0(s)dot x_0(s)+q_0(s)dot y_0(s)=-ap_0sin(s)+bq_0cos(s).
          $$
          This tells us that the unit vector $(p_0,q_0)$ has to be orthogonal to $(-asin(s), bcos(s))$, thus $$(p_0,q_0)=pmfrac{(bcos(s),asin(s))}{sqrt{b^2cos^2(s)+a^2sin^2(s)}}.$$
          As this should point inwards, chose the negative sign.



          To find the characteristic that some point $(x,y)$ lies on, one has to solve
          $$
          left.begin{aligned}
          x&=x_0+2p_0t=acos(s) - r b cos(s)\
          y&=y_0+2q_0t=bsin(s) - r a sin(s)
          end{aligned}right}
          implies frac{x^2}{(a-rb)^2}+frac{y^2}{(b-ar)^2}=1
          $$
          where $r$ is some positive multiple of $t$.
          The last is equivalent to a 4th degree polynomial equation, by the intermediate value theorem there is at least one solution between $r=0$ and $min(frac ba,frac ab)$, there is another solution for $r>max(frac ba,frac ab)$, but that might be inadmissible.




          Two close-by rays associated to the angles $s$ and $s+ds$ intersect at
          $$
          left.begin{aligned}
          (a - r b) cos(s) = (a - (r+dr) b) cos(s+ds)\
          (b - r a) sin(s) = (b - (r+dr) a) sin(s+ds)
          end{aligned}right}
          implies
          left{begin{aligned}
          0=-b,dr,cos(s)-(a-br)sin(s)ds\
          0=-a,dr,sin(s)+(b-ar)cos(s)ds
          end{aligned}right.
          $$
          so that necessarily $$a(a-br)sin^2(s)+b(b-ar)cos^2(s)=0implies r = frac{a^2sin^2s+b^2cos^2s}{ab}$$



          Plotting that curve against the ellipse gives the picture
          enter image description here






          share|cite|improve this answer









          $endgroup$













          • $begingroup$
            I thought to find the envelope we need to find where the jacobian = 0 so we can find an expression for t? I realise my expression for the jacobian was wrong but the method of solving for t then putting $p_0,q_0,t$ back into the characteristics to find the envelope is still correct , is it not?
            $endgroup$
            – pablo_mathscobar
            Dec 12 '18 at 11:15












          • $begingroup$
            You use $dz=p,dx+q,dy$ along the initial condition curve together with all other previously computed relations to determine $(p_0,q_0)$. The zeros of the determinant of the Jacobian of $$(s,t)mapsto(,x_0(s)+2tp_0(s),, y_0(s)+2tq_0(s),)=(,(a+2tf(s),b)cos(s),,(b+2tf(s),a)sin(s),)$$ tell you where to find the envelope. This is the same as me computing the intersection of neighboring rays with $r=2tf(s)$. The points at $s=kfracpi2$ that are inside the ellipse are the borders of the ridge.
            $endgroup$
            – LutzL
            Dec 12 '18 at 11:47












          • $begingroup$
            I see computing the $frac{partial}{partial s} acos(s) + 2tbcos(s)f(s)$ is quite a messy derivative because of the f(s) term so i guess your way is better? also how did you find the intersection between the 2 rays? what did you do remove the cos(ds) terms
            $endgroup$
            – pablo_mathscobar
            Dec 12 '18 at 12:06












          • $begingroup$
            It does not need to be messy if you leave the derivative at $f'(s)$. But yes, by the chain rule you get $$frac{∂(x,y)}{∂(s,t)}=frac{∂(x,y)}{∂(s,r)}frac{∂(s,r)}{∂(s,t)}text{ with }frac{∂(s,r)}{∂(s,t)}=pmatrix{1&0\2tf'(s)&2f(s)},$$ so that the zeros of the determinants give the same points.
            $endgroup$
            – LutzL
            Dec 12 '18 at 12:10











          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',
          autoActivateHeartbeat: false,
          convertImagesToLinks: true,
          noModals: true,
          showLowRepImageUploadWarning: true,
          reputationToPostImages: 10,
          bindNavPrevention: true,
          postfix: "",
          imageUploader: {
          brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
          contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
          allowUrls: true
          },
          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%2f3033828%2fenvelope-of-rays-that-propagate-into-the-ellipse-for-f-u-x2-u-y2-1%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown

























          1 Answer
          1






          active

          oldest

          votes








          1 Answer
          1






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes









          0












          $begingroup$

          I think you might have missed a sign in $dot x_0(s)=-asin(s)$. Then the equation for $p_0,q_0$ reads as
          $$
          0=dot z_0(s)=p_0(s)dot x_0(s)+q_0(s)dot y_0(s)=-ap_0sin(s)+bq_0cos(s).
          $$
          This tells us that the unit vector $(p_0,q_0)$ has to be orthogonal to $(-asin(s), bcos(s))$, thus $$(p_0,q_0)=pmfrac{(bcos(s),asin(s))}{sqrt{b^2cos^2(s)+a^2sin^2(s)}}.$$
          As this should point inwards, chose the negative sign.



          To find the characteristic that some point $(x,y)$ lies on, one has to solve
          $$
          left.begin{aligned}
          x&=x_0+2p_0t=acos(s) - r b cos(s)\
          y&=y_0+2q_0t=bsin(s) - r a sin(s)
          end{aligned}right}
          implies frac{x^2}{(a-rb)^2}+frac{y^2}{(b-ar)^2}=1
          $$
          where $r$ is some positive multiple of $t$.
          The last is equivalent to a 4th degree polynomial equation, by the intermediate value theorem there is at least one solution between $r=0$ and $min(frac ba,frac ab)$, there is another solution for $r>max(frac ba,frac ab)$, but that might be inadmissible.




          Two close-by rays associated to the angles $s$ and $s+ds$ intersect at
          $$
          left.begin{aligned}
          (a - r b) cos(s) = (a - (r+dr) b) cos(s+ds)\
          (b - r a) sin(s) = (b - (r+dr) a) sin(s+ds)
          end{aligned}right}
          implies
          left{begin{aligned}
          0=-b,dr,cos(s)-(a-br)sin(s)ds\
          0=-a,dr,sin(s)+(b-ar)cos(s)ds
          end{aligned}right.
          $$
          so that necessarily $$a(a-br)sin^2(s)+b(b-ar)cos^2(s)=0implies r = frac{a^2sin^2s+b^2cos^2s}{ab}$$



          Plotting that curve against the ellipse gives the picture
          enter image description here






          share|cite|improve this answer









          $endgroup$













          • $begingroup$
            I thought to find the envelope we need to find where the jacobian = 0 so we can find an expression for t? I realise my expression for the jacobian was wrong but the method of solving for t then putting $p_0,q_0,t$ back into the characteristics to find the envelope is still correct , is it not?
            $endgroup$
            – pablo_mathscobar
            Dec 12 '18 at 11:15












          • $begingroup$
            You use $dz=p,dx+q,dy$ along the initial condition curve together with all other previously computed relations to determine $(p_0,q_0)$. The zeros of the determinant of the Jacobian of $$(s,t)mapsto(,x_0(s)+2tp_0(s),, y_0(s)+2tq_0(s),)=(,(a+2tf(s),b)cos(s),,(b+2tf(s),a)sin(s),)$$ tell you where to find the envelope. This is the same as me computing the intersection of neighboring rays with $r=2tf(s)$. The points at $s=kfracpi2$ that are inside the ellipse are the borders of the ridge.
            $endgroup$
            – LutzL
            Dec 12 '18 at 11:47












          • $begingroup$
            I see computing the $frac{partial}{partial s} acos(s) + 2tbcos(s)f(s)$ is quite a messy derivative because of the f(s) term so i guess your way is better? also how did you find the intersection between the 2 rays? what did you do remove the cos(ds) terms
            $endgroup$
            – pablo_mathscobar
            Dec 12 '18 at 12:06












          • $begingroup$
            It does not need to be messy if you leave the derivative at $f'(s)$. But yes, by the chain rule you get $$frac{∂(x,y)}{∂(s,t)}=frac{∂(x,y)}{∂(s,r)}frac{∂(s,r)}{∂(s,t)}text{ with }frac{∂(s,r)}{∂(s,t)}=pmatrix{1&0\2tf'(s)&2f(s)},$$ so that the zeros of the determinants give the same points.
            $endgroup$
            – LutzL
            Dec 12 '18 at 12:10
















          0












          $begingroup$

          I think you might have missed a sign in $dot x_0(s)=-asin(s)$. Then the equation for $p_0,q_0$ reads as
          $$
          0=dot z_0(s)=p_0(s)dot x_0(s)+q_0(s)dot y_0(s)=-ap_0sin(s)+bq_0cos(s).
          $$
          This tells us that the unit vector $(p_0,q_0)$ has to be orthogonal to $(-asin(s), bcos(s))$, thus $$(p_0,q_0)=pmfrac{(bcos(s),asin(s))}{sqrt{b^2cos^2(s)+a^2sin^2(s)}}.$$
          As this should point inwards, chose the negative sign.



          To find the characteristic that some point $(x,y)$ lies on, one has to solve
          $$
          left.begin{aligned}
          x&=x_0+2p_0t=acos(s) - r b cos(s)\
          y&=y_0+2q_0t=bsin(s) - r a sin(s)
          end{aligned}right}
          implies frac{x^2}{(a-rb)^2}+frac{y^2}{(b-ar)^2}=1
          $$
          where $r$ is some positive multiple of $t$.
          The last is equivalent to a 4th degree polynomial equation, by the intermediate value theorem there is at least one solution between $r=0$ and $min(frac ba,frac ab)$, there is another solution for $r>max(frac ba,frac ab)$, but that might be inadmissible.




          Two close-by rays associated to the angles $s$ and $s+ds$ intersect at
          $$
          left.begin{aligned}
          (a - r b) cos(s) = (a - (r+dr) b) cos(s+ds)\
          (b - r a) sin(s) = (b - (r+dr) a) sin(s+ds)
          end{aligned}right}
          implies
          left{begin{aligned}
          0=-b,dr,cos(s)-(a-br)sin(s)ds\
          0=-a,dr,sin(s)+(b-ar)cos(s)ds
          end{aligned}right.
          $$
          so that necessarily $$a(a-br)sin^2(s)+b(b-ar)cos^2(s)=0implies r = frac{a^2sin^2s+b^2cos^2s}{ab}$$



          Plotting that curve against the ellipse gives the picture
          enter image description here






          share|cite|improve this answer









          $endgroup$













          • $begingroup$
            I thought to find the envelope we need to find where the jacobian = 0 so we can find an expression for t? I realise my expression for the jacobian was wrong but the method of solving for t then putting $p_0,q_0,t$ back into the characteristics to find the envelope is still correct , is it not?
            $endgroup$
            – pablo_mathscobar
            Dec 12 '18 at 11:15












          • $begingroup$
            You use $dz=p,dx+q,dy$ along the initial condition curve together with all other previously computed relations to determine $(p_0,q_0)$. The zeros of the determinant of the Jacobian of $$(s,t)mapsto(,x_0(s)+2tp_0(s),, y_0(s)+2tq_0(s),)=(,(a+2tf(s),b)cos(s),,(b+2tf(s),a)sin(s),)$$ tell you where to find the envelope. This is the same as me computing the intersection of neighboring rays with $r=2tf(s)$. The points at $s=kfracpi2$ that are inside the ellipse are the borders of the ridge.
            $endgroup$
            – LutzL
            Dec 12 '18 at 11:47












          • $begingroup$
            I see computing the $frac{partial}{partial s} acos(s) + 2tbcos(s)f(s)$ is quite a messy derivative because of the f(s) term so i guess your way is better? also how did you find the intersection between the 2 rays? what did you do remove the cos(ds) terms
            $endgroup$
            – pablo_mathscobar
            Dec 12 '18 at 12:06












          • $begingroup$
            It does not need to be messy if you leave the derivative at $f'(s)$. But yes, by the chain rule you get $$frac{∂(x,y)}{∂(s,t)}=frac{∂(x,y)}{∂(s,r)}frac{∂(s,r)}{∂(s,t)}text{ with }frac{∂(s,r)}{∂(s,t)}=pmatrix{1&0\2tf'(s)&2f(s)},$$ so that the zeros of the determinants give the same points.
            $endgroup$
            – LutzL
            Dec 12 '18 at 12:10














          0












          0








          0





          $begingroup$

          I think you might have missed a sign in $dot x_0(s)=-asin(s)$. Then the equation for $p_0,q_0$ reads as
          $$
          0=dot z_0(s)=p_0(s)dot x_0(s)+q_0(s)dot y_0(s)=-ap_0sin(s)+bq_0cos(s).
          $$
          This tells us that the unit vector $(p_0,q_0)$ has to be orthogonal to $(-asin(s), bcos(s))$, thus $$(p_0,q_0)=pmfrac{(bcos(s),asin(s))}{sqrt{b^2cos^2(s)+a^2sin^2(s)}}.$$
          As this should point inwards, chose the negative sign.



          To find the characteristic that some point $(x,y)$ lies on, one has to solve
          $$
          left.begin{aligned}
          x&=x_0+2p_0t=acos(s) - r b cos(s)\
          y&=y_0+2q_0t=bsin(s) - r a sin(s)
          end{aligned}right}
          implies frac{x^2}{(a-rb)^2}+frac{y^2}{(b-ar)^2}=1
          $$
          where $r$ is some positive multiple of $t$.
          The last is equivalent to a 4th degree polynomial equation, by the intermediate value theorem there is at least one solution between $r=0$ and $min(frac ba,frac ab)$, there is another solution for $r>max(frac ba,frac ab)$, but that might be inadmissible.




          Two close-by rays associated to the angles $s$ and $s+ds$ intersect at
          $$
          left.begin{aligned}
          (a - r b) cos(s) = (a - (r+dr) b) cos(s+ds)\
          (b - r a) sin(s) = (b - (r+dr) a) sin(s+ds)
          end{aligned}right}
          implies
          left{begin{aligned}
          0=-b,dr,cos(s)-(a-br)sin(s)ds\
          0=-a,dr,sin(s)+(b-ar)cos(s)ds
          end{aligned}right.
          $$
          so that necessarily $$a(a-br)sin^2(s)+b(b-ar)cos^2(s)=0implies r = frac{a^2sin^2s+b^2cos^2s}{ab}$$



          Plotting that curve against the ellipse gives the picture
          enter image description here






          share|cite|improve this answer









          $endgroup$



          I think you might have missed a sign in $dot x_0(s)=-asin(s)$. Then the equation for $p_0,q_0$ reads as
          $$
          0=dot z_0(s)=p_0(s)dot x_0(s)+q_0(s)dot y_0(s)=-ap_0sin(s)+bq_0cos(s).
          $$
          This tells us that the unit vector $(p_0,q_0)$ has to be orthogonal to $(-asin(s), bcos(s))$, thus $$(p_0,q_0)=pmfrac{(bcos(s),asin(s))}{sqrt{b^2cos^2(s)+a^2sin^2(s)}}.$$
          As this should point inwards, chose the negative sign.



          To find the characteristic that some point $(x,y)$ lies on, one has to solve
          $$
          left.begin{aligned}
          x&=x_0+2p_0t=acos(s) - r b cos(s)\
          y&=y_0+2q_0t=bsin(s) - r a sin(s)
          end{aligned}right}
          implies frac{x^2}{(a-rb)^2}+frac{y^2}{(b-ar)^2}=1
          $$
          where $r$ is some positive multiple of $t$.
          The last is equivalent to a 4th degree polynomial equation, by the intermediate value theorem there is at least one solution between $r=0$ and $min(frac ba,frac ab)$, there is another solution for $r>max(frac ba,frac ab)$, but that might be inadmissible.




          Two close-by rays associated to the angles $s$ and $s+ds$ intersect at
          $$
          left.begin{aligned}
          (a - r b) cos(s) = (a - (r+dr) b) cos(s+ds)\
          (b - r a) sin(s) = (b - (r+dr) a) sin(s+ds)
          end{aligned}right}
          implies
          left{begin{aligned}
          0=-b,dr,cos(s)-(a-br)sin(s)ds\
          0=-a,dr,sin(s)+(b-ar)cos(s)ds
          end{aligned}right.
          $$
          so that necessarily $$a(a-br)sin^2(s)+b(b-ar)cos^2(s)=0implies r = frac{a^2sin^2s+b^2cos^2s}{ab}$$



          Plotting that curve against the ellipse gives the picture
          enter image description here







          share|cite|improve this answer












          share|cite|improve this answer



          share|cite|improve this answer










          answered Dec 10 '18 at 15:55









          LutzLLutzL

          57.6k42054




          57.6k42054












          • $begingroup$
            I thought to find the envelope we need to find where the jacobian = 0 so we can find an expression for t? I realise my expression for the jacobian was wrong but the method of solving for t then putting $p_0,q_0,t$ back into the characteristics to find the envelope is still correct , is it not?
            $endgroup$
            – pablo_mathscobar
            Dec 12 '18 at 11:15












          • $begingroup$
            You use $dz=p,dx+q,dy$ along the initial condition curve together with all other previously computed relations to determine $(p_0,q_0)$. The zeros of the determinant of the Jacobian of $$(s,t)mapsto(,x_0(s)+2tp_0(s),, y_0(s)+2tq_0(s),)=(,(a+2tf(s),b)cos(s),,(b+2tf(s),a)sin(s),)$$ tell you where to find the envelope. This is the same as me computing the intersection of neighboring rays with $r=2tf(s)$. The points at $s=kfracpi2$ that are inside the ellipse are the borders of the ridge.
            $endgroup$
            – LutzL
            Dec 12 '18 at 11:47












          • $begingroup$
            I see computing the $frac{partial}{partial s} acos(s) + 2tbcos(s)f(s)$ is quite a messy derivative because of the f(s) term so i guess your way is better? also how did you find the intersection between the 2 rays? what did you do remove the cos(ds) terms
            $endgroup$
            – pablo_mathscobar
            Dec 12 '18 at 12:06












          • $begingroup$
            It does not need to be messy if you leave the derivative at $f'(s)$. But yes, by the chain rule you get $$frac{∂(x,y)}{∂(s,t)}=frac{∂(x,y)}{∂(s,r)}frac{∂(s,r)}{∂(s,t)}text{ with }frac{∂(s,r)}{∂(s,t)}=pmatrix{1&0\2tf'(s)&2f(s)},$$ so that the zeros of the determinants give the same points.
            $endgroup$
            – LutzL
            Dec 12 '18 at 12:10


















          • $begingroup$
            I thought to find the envelope we need to find where the jacobian = 0 so we can find an expression for t? I realise my expression for the jacobian was wrong but the method of solving for t then putting $p_0,q_0,t$ back into the characteristics to find the envelope is still correct , is it not?
            $endgroup$
            – pablo_mathscobar
            Dec 12 '18 at 11:15












          • $begingroup$
            You use $dz=p,dx+q,dy$ along the initial condition curve together with all other previously computed relations to determine $(p_0,q_0)$. The zeros of the determinant of the Jacobian of $$(s,t)mapsto(,x_0(s)+2tp_0(s),, y_0(s)+2tq_0(s),)=(,(a+2tf(s),b)cos(s),,(b+2tf(s),a)sin(s),)$$ tell you where to find the envelope. This is the same as me computing the intersection of neighboring rays with $r=2tf(s)$. The points at $s=kfracpi2$ that are inside the ellipse are the borders of the ridge.
            $endgroup$
            – LutzL
            Dec 12 '18 at 11:47












          • $begingroup$
            I see computing the $frac{partial}{partial s} acos(s) + 2tbcos(s)f(s)$ is quite a messy derivative because of the f(s) term so i guess your way is better? also how did you find the intersection between the 2 rays? what did you do remove the cos(ds) terms
            $endgroup$
            – pablo_mathscobar
            Dec 12 '18 at 12:06












          • $begingroup$
            It does not need to be messy if you leave the derivative at $f'(s)$. But yes, by the chain rule you get $$frac{∂(x,y)}{∂(s,t)}=frac{∂(x,y)}{∂(s,r)}frac{∂(s,r)}{∂(s,t)}text{ with }frac{∂(s,r)}{∂(s,t)}=pmatrix{1&0\2tf'(s)&2f(s)},$$ so that the zeros of the determinants give the same points.
            $endgroup$
            – LutzL
            Dec 12 '18 at 12:10
















          $begingroup$
          I thought to find the envelope we need to find where the jacobian = 0 so we can find an expression for t? I realise my expression for the jacobian was wrong but the method of solving for t then putting $p_0,q_0,t$ back into the characteristics to find the envelope is still correct , is it not?
          $endgroup$
          – pablo_mathscobar
          Dec 12 '18 at 11:15






          $begingroup$
          I thought to find the envelope we need to find where the jacobian = 0 so we can find an expression for t? I realise my expression for the jacobian was wrong but the method of solving for t then putting $p_0,q_0,t$ back into the characteristics to find the envelope is still correct , is it not?
          $endgroup$
          – pablo_mathscobar
          Dec 12 '18 at 11:15














          $begingroup$
          You use $dz=p,dx+q,dy$ along the initial condition curve together with all other previously computed relations to determine $(p_0,q_0)$. The zeros of the determinant of the Jacobian of $$(s,t)mapsto(,x_0(s)+2tp_0(s),, y_0(s)+2tq_0(s),)=(,(a+2tf(s),b)cos(s),,(b+2tf(s),a)sin(s),)$$ tell you where to find the envelope. This is the same as me computing the intersection of neighboring rays with $r=2tf(s)$. The points at $s=kfracpi2$ that are inside the ellipse are the borders of the ridge.
          $endgroup$
          – LutzL
          Dec 12 '18 at 11:47






          $begingroup$
          You use $dz=p,dx+q,dy$ along the initial condition curve together with all other previously computed relations to determine $(p_0,q_0)$. The zeros of the determinant of the Jacobian of $$(s,t)mapsto(,x_0(s)+2tp_0(s),, y_0(s)+2tq_0(s),)=(,(a+2tf(s),b)cos(s),,(b+2tf(s),a)sin(s),)$$ tell you where to find the envelope. This is the same as me computing the intersection of neighboring rays with $r=2tf(s)$. The points at $s=kfracpi2$ that are inside the ellipse are the borders of the ridge.
          $endgroup$
          – LutzL
          Dec 12 '18 at 11:47














          $begingroup$
          I see computing the $frac{partial}{partial s} acos(s) + 2tbcos(s)f(s)$ is quite a messy derivative because of the f(s) term so i guess your way is better? also how did you find the intersection between the 2 rays? what did you do remove the cos(ds) terms
          $endgroup$
          – pablo_mathscobar
          Dec 12 '18 at 12:06






          $begingroup$
          I see computing the $frac{partial}{partial s} acos(s) + 2tbcos(s)f(s)$ is quite a messy derivative because of the f(s) term so i guess your way is better? also how did you find the intersection between the 2 rays? what did you do remove the cos(ds) terms
          $endgroup$
          – pablo_mathscobar
          Dec 12 '18 at 12:06














          $begingroup$
          It does not need to be messy if you leave the derivative at $f'(s)$. But yes, by the chain rule you get $$frac{∂(x,y)}{∂(s,t)}=frac{∂(x,y)}{∂(s,r)}frac{∂(s,r)}{∂(s,t)}text{ with }frac{∂(s,r)}{∂(s,t)}=pmatrix{1&0\2tf'(s)&2f(s)},$$ so that the zeros of the determinants give the same points.
          $endgroup$
          – LutzL
          Dec 12 '18 at 12:10




          $begingroup$
          It does not need to be messy if you leave the derivative at $f'(s)$. But yes, by the chain rule you get $$frac{∂(x,y)}{∂(s,t)}=frac{∂(x,y)}{∂(s,r)}frac{∂(s,r)}{∂(s,t)}text{ with }frac{∂(s,r)}{∂(s,t)}=pmatrix{1&0\2tf'(s)&2f(s)},$$ so that the zeros of the determinants give the same points.
          $endgroup$
          – LutzL
          Dec 12 '18 at 12:10


















          draft saved

          draft discarded




















































          Thanks for contributing an answer to Mathematics Stack Exchange!


          • Please be sure to answer the question. Provide details and share your research!

          But avoid …



          • Asking for help, clarification, or responding to other answers.

          • Making statements based on opinion; back them up with references or personal experience.


          Use MathJax to format equations. MathJax reference.


          To learn more, see our tips on writing great answers.




          draft saved


          draft discarded














          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3033828%2fenvelope-of-rays-that-propagate-into-the-ellipse-for-f-u-x2-u-y2-1%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown





















































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown

































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown







          -

          Popular posts from this blog

          Wiesbaden

          Image
          Wappen Deutschlandkarte 50.0825 8.24 117 Koordinaten: 50° 5′  N , 8° 14′  O Basisdaten Bundesland: Hessen Regierungsbezirk: Darmstadt Höhe: 117 m ü. NHN Fläche: 203,93 km 2 Einwohner: 278.654 (31. Dez. 2017) [1] Bevölkerungsdichte: 1366 Einwohner je km 2 Postleitzahlen: 65183–65207, 55246, 55252 Vorlage:Infobox Gemeinde in Deutschland/Wartung/PLZ enthält Text Vorwahlen: 0611, 06122, 06127, 06134 Kfz-Kennzeichen: WI Gemeindeschlüssel: 06 4 14 000 LOCODE: DE WIB NUTS: DE714 Stadtgliederung: 26 Ortsbezirke Adresse der Stadtverwaltung: Schlossplatz 6 65183 Wiesbaden Webpräsenz: www.wiesbaden.de Oberbürgermeister: Sven Gerich (SPD) Lage der Landeshauptstadt Wiesbaden in Hessen und im Regierungsbezirk Darmstadt Offizielles Logo der Landeshauptstadt Wiesbaden Flagge der Landeshauptstadt Wiesbaden Wiesbaden is
          Read more

          Marschland

          Image
          Kilometerweite Marsch in den Vier- und Marschlanden am Elbdeich Die Wedeler Marschlandschaft in Schleswig-Holstein Der Hadelner Kanal zur Entwässerung der Marsch Die tief gelegene Marsch und ein Entwässerungskanal Das Schöpfwerk Otterndorf im Land Hadeln zur Entwässerung der Medem verfügte damals bereits über die größte Kreiselpumpe Europas. Das Schöpfwerk in Neuhaus für die Aue und des Neuhaus-Bülkauer Kanals Als Marsch(land) (v. niederdt., altsächs. mersc ) – auch Masch , Mersch oder Schwemmland genannt – bezeichnet man eine nacheiszeitlich entstandene geomorphologische Landform im Gebiet der nordwestdeutschen Küsten und Flüsse sowie vergleichbare Landformen weltweit. Inhaltsverzeichnis 1 Landschaft 2 Entstehung 3 Marschböden 4 Nutzung 5 Marschgebiete in Deutschland 6 Siehe auch 7 Literatur 8 Weblinks 9 Einzelnachweise Landschaft | Marschen sind generell flache Landstriche ohne natürliche Erhebungen. Si
          Read more

          Dieringhausen

          Image
          Dieringhausen Stadt Gummersbach 50.982777777778 7.5316666666667 165 Koordinaten: 50° 58′ 58″  N , 7° 31′ 54″  O Höhe: 165  (160–245)  m Einwohner: 5055  (30. Jun. 2016) Postleitzahl: 51645 Vorwahl: 02261 Lage von Dieringhausen in Gummersbach Im Hohl ca. 1927–1930 Ansicht von Südwesten bei Bünghausen Dieringhausen (im örtlichen Dialekt Dierkesen ) ist ein Ortsteil von Gummersbach im Oberbergischen Kreis, Regierungsbezirk Köln, Nordrhein-Westfalen (Deutschland). Inhaltsverzeichnis 1 Geographie 2 Geschichte 3 Kultur 3.1 Sehenswürdigkeiten 3.2 Regelmäßige Veranstaltungen 4 Infrastruktur und Wirtschaft 4.1 Verkehr 4.1.1 Schienen- und Busverkehr 4.1.2 Straßen 4.2 Öffentliche Einrichtungen 4.2.1 Schulen und Bildungseinrichtungen 4.2.2 Kirchliche Einrichtungen 5 Persönlichkeiten 5.1 In Dieringhausen geboren 5.2 In Dieringhausen gelebt 6 Einzelnachweise 7
          Read more