Skip to main content

differential equations - VariationalD giving the wrong solution?


EDIT: As pointed out in the comments, VariationalD gives a variational derivative (which I don't want), not a derivative with respect to a function (i.e. $\frac{df[x]}{dlog_e(x)}$ as a simple example - this is what I thought it did the first time I read the description). Is it even possible in Mathematica to take the derivative of a function with respect to another function? I know you can use the chain rule to rewrite, using the example just above, $\frac{df[x]}{dlog_e(x)}=x\frac{df[x]}{dx}$. In my case the equivalent would be much more cumbersome and I want to be able to change the function w.r.t. which I'm differentiating. Is it possible to do this in Mathematica?
_


I am using the VariationalMethods` package, specifically the VariationalD command. I am trying to take the derivative of one (complicated) function fitted to experimental data with respect to another (complicated) function.


First I wanted to make sure that the command does what I wanted. So I inputted the example given in VariationalD's "Examples" section:


 VariationalD[y[x] Sqrt[y'[x]], y[x], x]


This gives a result of $$ \frac{d}{dy[x]}(y[x] \cdot y'[x]^{1/2})=\frac{2 y'[x]^2+y[x] y''[x]}{4 y'[x]^{3/2}}. $$


Trying to reproduce this by hand I started with the product rule:


$$ \frac{d}{dy[x]}(y[x] \cdot y'[x]^{1/2})=y[x]\frac{d}{dy[x]}y'[x]^{1/2}+y'[x]^{1/2}\frac{d}{dy[x]}y[x] $$ The second term is just


$$ y'[x]^{1/2}\frac{d}{dy[x]}y[x]=y'[x]^{1/2}. $$


The first term, you use the chain rule on:


$$ y[x]\frac{d}{dy[x]}y'[x]^{1/2}=y[x]\frac{1}{2}y'[x]^{-1/2}\frac{d}{dy[x]}y'[x] $$ $$ \frac{d}{dy[x]}y'[x]=\frac{dx}{dy[x]}\frac{d}{dx}y'[x]=(\frac{dy[x]}{dx})^{-1}y''[x]=\frac{y''[x]}{y'[x]}. $$


Putting all of that together:


$$ \frac{d}{dy[x]}(y[x] \cdot y'[x]^{1/2})=y'[x]^{1/2}+\frac{1}{2}\frac{y[x]y''[x]}{y'[x]^{3/2}}. $$


Simplifying:



$$ \frac{d}{dy[x]}(y[x] \cdot y'[x]^{1/2})=\frac{2 y'[x]^2+y[x]y''[x]}{2 y'[x]^{3/2}}, $$


which is identical to the solution given by VariationalD except for a factor of 2. So I assumed it was a mistake I made somewhere (and it might be) and tried comparing every step I made when doing it by hand with the corresponding step in VariationalD. Eventually I discovered one difference in the step


 VariationalD[Sqrt[y'[x]], y[x], x]

which gives an answer of $\frac{y''[x]}{4y'[x]^{3/2}}$, a factor of 2 different than when I do it. I'm still not sure where the factor of 2 in the $y'[x]^{1/2}\frac{d}{dy[x]}y[x]$ step comes from.


I kept messing around with it and found that


 VariationalD[y'[x],y[x],x]=0

which can't be right. For example, if $y[x]=x^2$ and $y'[x]=2 x$,


$$ \frac{dy'[x]}{dy[x]}=\frac{d(2x)}{dx^2}=2(\frac{dx^2}{dx})^{-1}=2(2x)^{-1}=\frac{1}{x}. $$



So why is VariationalD giving a zero answer? It's even weirder because $\frac{dy'[x]}{dy[x]}$ comes up when you perform the original differentiation by hand as well. It can't be zero there either or the first term would vanish and the answer would just be $y'[x]^{1/2}$.


So am I making a stupid mistake in my differentiation? Even if so, why is VariationalD giving zero for the derivative with respect to a function of its derivative?


I've spent about an hour searching for anything relevant online, and...nothing.


Any help would be appreciated!


EDIT: I've been looking at the description of VariationalD more closely and now I'm not sure it does what I think it does, i.e.


 VariationalD[f[x],g[x],x]

gives $\frac{df[x]}{dg[x]}$ as output. Is this correct? Looking at the description it now seems to me that it would give $f[x]\frac{df[x]}{dg[x]}$. But if that's the case then the output for my original expression is still wrong.



Answer



I think VariationalD[y'[x],y[x],x] indeed equals to zero.



You should not mistake functional derivative with ordinary derivative, where in the former case $y'$ is usually considered an independent variable to $y$. so functional $\mathrm{d} y'/\mathrm{d} y=0$ is just like an ordinary $\mathrm{d} a/\mathrm{d} b=0$.


The same reason you can't write something like this in calculus of variations


$$\frac{\partial J[x,y,y']}{\partial y}=\frac{\partial x}{\partial y}\frac{\partial J[x,y,y']}{\partial x}\,\text{,}$$


just like you can't write this in ordinary calculus:


$$\frac{\partial f(x,y,z)}{\partial y}=\frac{\partial x}{\partial y}\frac{\partial f(x,y,z)}{\partial x}\,\text{.}$$


Comments

Popular posts from this blog

plotting - How to draw lines between specified dots on ListPlot?

I would like to create a plot where I have unconnected dots and some connected. So far, I have figured out how to draw the dots. My code is the following: ListPlot[{{1, 1}, {2, 2}, {3, 3}, {4, 4}, {1, 4}, {2, 5}, {3, 6}, {4, 7}, {1, 7}, {2, 8}, {3, 9}, {4, 10}, {1, 10}, {2, 11}, {3, 12}, {4,13}, {2.5, 7}}, Ticks -> {{1, 2, 3, 4}, None}, AxesStyle -> Thin, TicksStyle -> Directive[Black, Bold, 12], Mesh -> Full] I have thought using ListLinePlot command, but I don't know how to specify to the command to draw only selected lines between the dots. Do have any suggestions/hints on how to do that? Thank you. Answer One possibility would be to use Epilog with Line : ListPlot[ {{1, 1}, {2, 2}, {3, 3}, {4, 4}, {1, 4}, {2, 5}, {3, 6}, {4, 7}, {1, 7}, {2, 8}, {3, 9}, {4, 10}, {1, 10}, {2, 11}, {3, 12}, {4, 13}, {2.5, 7}}, Ticks -> {{1, 2, 3, 4}, None}, AxesStyle -> Thin, TicksStyle -> Directive[Black, Bold, 12], Mesh -> Full, Epilog -> { Line[ ...

dynamic - How can I make a clickable ArrayPlot that returns input?

I would like to create a dynamic ArrayPlot so that the rectangles, when clicked, provide the input. Can I use ArrayPlot for this? Or is there something else I should have to use? Answer ArrayPlot is much more than just a simple array like Grid : it represents a ranged 2D dataset, and its visualization can be finetuned by options like DataReversed and DataRange . These features make it quite complicated to reproduce the same layout and order with Grid . Here I offer AnnotatedArrayPlot which comes in handy when your dataset is more than just a flat 2D array. The dynamic interface allows highlighting individual cells and possibly interacting with them. AnnotatedArrayPlot works the same way as ArrayPlot and accepts the same options plus Enabled , HighlightCoordinates , HighlightStyle and HighlightElementFunction . data = {{Missing["HasSomeMoreData"], GrayLevel[ 1], {RGBColor[0, 1, 1], RGBColor[0, 0, 1], GrayLevel[1]}, RGBColor[0, 1, 0]}, {GrayLevel[0], GrayLevel...

functions - What best practices or performance considerations are there for choosing between Cases, Position, Pick and Select?

Cases , Select , Pick and Position each have different syntaxes and purposes, but there are times when you can express the same calculation equivalently using either of them. So with this input: test = RandomInteger[{-25, 25}, {20, 2}] {{-15, 13}, {-8, 16}, {-8, -19}, {7, 6}, {-21, 9}, {-3, -25}, {21, -18}, {4, 4}, {2, -2}, {-24, 8}, {-17, -8}, {4, -18}, {22, -24}, {-4, -3}, {21, 0}, {19, 18}, {-23, -8}, {23, -25}, {14, -2}, {-1, -13}} You can get the following equivalent results: Cases[test, {_, _?Positive}] {{-15, 13}, {-8, 16}, {7, 6}, {-21, 9}, {4, 4}, {-24, 8}, {19, 18}} Select[test, #[[2]] > 0 &] {{-15, 13}, {-8, 16}, {7, 6}, {-21, 9}, {4, 4}, {-24, 8}, {19, 18}} Pick[test, Sign[test[[All, 2]] ], 1] {{-15, 13}, {-8, 16}, {7, 6}, {-21, 9}, {4, 4}, {-24, 8}, {19, 18}} test[[Flatten@Position[test[[All, 2]], _?Positive] ]] {{-15, 13}, {-8, 16}, {7, 6}, {-21, 9}, {4, 4}, {-24, 8}, {19, 18}} Are there performance or other considerations that should guide which you shou...