Ttyjukl

Is there a way to find out the age of climbing ropes?

Practical reasons to have both a large police force and bounty hunting network?

Giving a talk in my old university, how prominently should I tell students my salary?

Can a Mimic (container form) actually hold loot?

ESPP--any reason not to go all in?

What is the oldest European royal house?

Searching for a string that contains the file name

Error in TransformedField

Calculate total length of edges in select Voronoi diagram

Affine transformation of circular arc in 3D

Python function to ask for a multiple-choice answer

Forcing Mathematica's Integrate to give more general answers

Naming Characters after Friends/Family

“I had a flat in the centre of town, but I didn’t like living there, so …”

Is being socially reclusive okay for a graduate student?

Align equations with text before one of them

Named nets not connected in Eagle board design

How to write a chaotic neutral protagonist and prevent my readers from thinking they are evil?

Why doesn't "adolescent" take any articles in "listen to adolescent agonising"?

How spaceships determine each other's mass in space?

Computing the volume of a simplex-like object with constraints

What does it mean when I add a new variable to my linear model and the R^2 stays the same?

Convert an array of objects to array of the objects' values

How do you make a gun that shoots melee weapons and/or swords?

Error in TransformedField

Using spherical derivatives with NDSolveIs there an easy way to transform unit vectors from spherical to Cartesian coordinates?Why does a transform from polar coordinates using `TransformedField` not agree with by-hand calculation?ListDensityPlot in Polar coordinates with a higher efficiencyChanging from Cartesian to polar coordinates in partial differential equationA little confused about coordinate conversionInverse substitution polar-cartesianIntegration using polar coordinatesPlotting a function $psi(rho,theta,phi)$ in spherical coordinatesTransform a field with derivatives

6

1

$begingroup$

I am using TransformedField to convert a system of ODEs from Cartesian to polar coordinates:

TransformedField[

  "Cartesian" -> "Polar",

  {μ x1 - x2 -  σ x1 (x1^2 + x2^2), x1 + μ x2 - σ x2 (x1^2 + x2^2)}, 

  {x1, x2} -> {r, θ}

] // Simplify

and I get the result

{r μ - r^3 σ, r}

but I am pretty sure that the right answer should be

{r μ - r^3 σ, 1}

Where is the error?

coordinate-transformation

share|improve this question

edited 4 hours ago

MarcoB

36.9k556113

asked 4 hours ago

rparpa

856

$endgroup$

add a comment | 

6

1

$begingroup$

I am using TransformedField to convert a system of ODEs from Cartesian to polar coordinates:

TransformedField[

  "Cartesian" -> "Polar",

  {μ x1 - x2 -  σ x1 (x1^2 + x2^2), x1 + μ x2 - σ x2 (x1^2 + x2^2)}, 

  {x1, x2} -> {r, θ}

] // Simplify

and I get the result

{r μ - r^3 σ, r}

but I am pretty sure that the right answer should be

{r μ - r^3 σ, 1}

Where is the error?

coordinate-transformation

share|improve this question

edited 4 hours ago

MarcoB

36.9k556113

asked 4 hours ago

rparpa

856

$endgroup$

add a comment | 

$begingroup$

I am using TransformedField to convert a system of ODEs from Cartesian to polar coordinates:

TransformedField[

  "Cartesian" -> "Polar",

  {μ x1 - x2 -  σ x1 (x1^2 + x2^2), x1 + μ x2 - σ x2 (x1^2 + x2^2)}, 

  {x1, x2} -> {r, θ}

] // Simplify

and I get the result

{r μ - r^3 σ, r}

but I am pretty sure that the right answer should be

{r μ - r^3 σ, 1}

Where is the error?

coordinate-transformation

share|improve this question

edited 4 hours ago

MarcoB

36.9k556113

asked 4 hours ago

rparpa

856

$endgroup$

TransformedField[

  "Cartesian" -> "Polar",

  {μ x1 - x2 -  σ x1 (x1^2 + x2^2), x1 + μ x2 - σ x2 (x1^2 + x2^2)}, 

  {x1, x2} -> {r, θ}

] // Simplify

{r μ - r^3 σ, r}

{r μ - r^3 σ, 1}

share|improve this question

edited 4 hours ago

MarcoB

36.9k556113

asked 4 hours ago

rparpa

856

share|improve this question

edited 4 hours ago

MarcoB

36.9k556113

asked 4 hours ago

rparpa

856

edited 4 hours ago

MarcoB

36.9k556113

edited 4 hours ago

MarcoB

36.9k556113

asked 4 hours ago

rparpa

856

asked 4 hours ago

rparpa

856

1 Answer
1

active

oldest

votes

5

$begingroup$

We can test the built-in TransformedField by defining our own functions.

From $x',y'$ to $r',theta'$, we derive:
$$
r' = left(sqrt{x^2 +y^2} right)'
= frac{(x^2 +y^2)'}{2
sqrt{x^2 +y^2}}=frac{xx' +yy'}{r}
$$
and
$$
theta' = left(arctan frac{y}{x} right)'
= frac{(y/x)'}{1+(y/x)^2} = frac{y' x -x' y}{r^2}.
$$

First, we define

rdot[x1_, x2_] := (x1 (μ x1 - x2 - σ x1 (x1^2 + x2^2)) + x2 (x1 + μ x2 - σ x2 (x1^2 + x2^2)))/r

We now make the substitution and simplify

rdot[r Cos[t], r Sin[t]] // FullSimplify

This yields (matches Mathematica)

$$r' = mu r-r^3 sigma$$

We now do the same for the other

thetadot[x1_,x2_]:=(x1 (x1+μ x2-σ x2 (x1^2+x2^2)) - x2(μ x1-x2-σ x1 (x1^2+x2^2)))/r^2

We now make the substitution and simplify

thetadot[r Cos[t], r Sin[t]] // FullSimplify

This yields (does not match Mathematica - it is as though they are forgetting to divide by $r$)

$$theta'= 1$$

I have asked this question before on this site in two different ways and have never gotten an answer that resolves the matter, but that could just be my denseness!

share|improve this answer

edited 1 hour ago

answered 2 hours ago

MooMoo

7511515

$endgroup$

• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  Have you reported it to the Wolfram tech support?
  $endgroup$
  – Alexey Popkov
  17 mins ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  @AlexeyPopkov: I have not. I have had many issues with it when transforming between different methods. These days, I don't trust it and just create my own transformation rules to do it.
  $endgroup$
  – Moo
  14 mins ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  It is worth to write them about it in order to get it finally fixed. You can even write a short letter to support@wolfram.com with a link to this post.
  $endgroup$
  – Alexey Popkov
  11 mins ago
  
  
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  @AlexeyPopkov: I sent them an email per your suggestion.
  $endgroup$
  – Moo
  5 mins ago
  

  
  
  
  

add a comment | 

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: "387"
};
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: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
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
},
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});


}
});

draft saved

draft discarded

Sign up or log in

StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});

Sign up using Google

Sign up using Facebook

Sign up using Email and Password

Post as a guest

Name

Email

Required, but never shown

StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f192821%2ferror-in-transformedfield%23new-answer', 'question_page');
}
);

Post as a guest

Name

Email

Required, but never shown

5

$begingroup$

We can test the built-in TransformedField by defining our own functions.

From $x',y'$ to $r',theta'$, we derive:
$$
r' = left(sqrt{x^2 +y^2} right)'
= frac{(x^2 +y^2)'}{2
sqrt{x^2 +y^2}}=frac{xx' +yy'}{r}
$$
and
$$
theta' = left(arctan frac{y}{x} right)'
= frac{(y/x)'}{1+(y/x)^2} = frac{y' x -x' y}{r^2}.
$$

First, we define

rdot[x1_, x2_] := (x1 (μ x1 - x2 - σ x1 (x1^2 + x2^2)) + x2 (x1 + μ x2 - σ x2 (x1^2 + x2^2)))/r

We now make the substitution and simplify

rdot[r Cos[t], r Sin[t]] // FullSimplify

This yields (matches Mathematica)

$$r' = mu r-r^3 sigma$$

We now do the same for the other

thetadot[x1_,x2_]:=(x1 (x1+μ x2-σ x2 (x1^2+x2^2)) - x2(μ x1-x2-σ x1 (x1^2+x2^2)))/r^2

We now make the substitution and simplify

thetadot[r Cos[t], r Sin[t]] // FullSimplify

This yields (does not match Mathematica - it is as though they are forgetting to divide by $r$)

$$theta'= 1$$

I have asked this question before on this site in two different ways and have never gotten an answer that resolves the matter, but that could just be my denseness!

share|improve this answer

edited 1 hour ago

answered 2 hours ago

MooMoo

7511515

$endgroup$

• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  Have you reported it to the Wolfram tech support?
  $endgroup$
  – Alexey Popkov
  17 mins ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  @AlexeyPopkov: I have not. I have had many issues with it when transforming between different methods. These days, I don't trust it and just create my own transformation rules to do it.
  $endgroup$
  – Moo
  14 mins ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  It is worth to write them about it in order to get it finally fixed. You can even write a short letter to support@wolfram.com with a link to this post.
  $endgroup$
  – Alexey Popkov
  11 mins ago
  
  
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  @AlexeyPopkov: I sent them an email per your suggestion.
  $endgroup$
  – Moo
  5 mins ago
  

  
  
  
  

add a comment | 

5

$begingroup$

We can test the built-in TransformedField by defining our own functions.

From $x',y'$ to $r',theta'$, we derive:
$$
r' = left(sqrt{x^2 +y^2} right)'
= frac{(x^2 +y^2)'}{2
sqrt{x^2 +y^2}}=frac{xx' +yy'}{r}
$$
and
$$
theta' = left(arctan frac{y}{x} right)'
= frac{(y/x)'}{1+(y/x)^2} = frac{y' x -x' y}{r^2}.
$$

First, we define

rdot[x1_, x2_] := (x1 (μ x1 - x2 - σ x1 (x1^2 + x2^2)) + x2 (x1 + μ x2 - σ x2 (x1^2 + x2^2)))/r

We now make the substitution and simplify

rdot[r Cos[t], r Sin[t]] // FullSimplify

This yields (matches Mathematica)

$$r' = mu r-r^3 sigma$$

We now do the same for the other

thetadot[x1_,x2_]:=(x1 (x1+μ x2-σ x2 (x1^2+x2^2)) - x2(μ x1-x2-σ x1 (x1^2+x2^2)))/r^2

We now make the substitution and simplify

thetadot[r Cos[t], r Sin[t]] // FullSimplify

This yields (does not match Mathematica - it is as though they are forgetting to divide by $r$)

$$theta'= 1$$

I have asked this question before on this site in two different ways and have never gotten an answer that resolves the matter, but that could just be my denseness!

share|improve this answer

edited 1 hour ago

answered 2 hours ago

MooMoo

7511515

$endgroup$

• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  Have you reported it to the Wolfram tech support?
  $endgroup$
  – Alexey Popkov
  17 mins ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  @AlexeyPopkov: I have not. I have had many issues with it when transforming between different methods. These days, I don't trust it and just create my own transformation rules to do it.
  $endgroup$
  – Moo
  14 mins ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  It is worth to write them about it in order to get it finally fixed. You can even write a short letter to support@wolfram.com with a link to this post.
  $endgroup$
  – Alexey Popkov
  11 mins ago
  
  
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  @AlexeyPopkov: I sent them an email per your suggestion.
  $endgroup$
  – Moo
  5 mins ago
  

  
  
  
  

add a comment | 

$begingroup$

We can test the built-in TransformedField by defining our own functions.

From $x',y'$ to $r',theta'$, we derive:
$$
r' = left(sqrt{x^2 +y^2} right)'
= frac{(x^2 +y^2)'}{2
sqrt{x^2 +y^2}}=frac{xx' +yy'}{r}
$$
and
$$
theta' = left(arctan frac{y}{x} right)'
= frac{(y/x)'}{1+(y/x)^2} = frac{y' x -x' y}{r^2}.
$$

First, we define

rdot[x1_, x2_] := (x1 (μ x1 - x2 - σ x1 (x1^2 + x2^2)) + x2 (x1 + μ x2 - σ x2 (x1^2 + x2^2)))/r

We now make the substitution and simplify

rdot[r Cos[t], r Sin[t]] // FullSimplify

This yields (matches Mathematica)

$$r' = mu r-r^3 sigma$$

We now do the same for the other

thetadot[x1_,x2_]:=(x1 (x1+μ x2-σ x2 (x1^2+x2^2)) - x2(μ x1-x2-σ x1 (x1^2+x2^2)))/r^2

We now make the substitution and simplify

thetadot[r Cos[t], r Sin[t]] // FullSimplify

This yields (does not match Mathematica - it is as though they are forgetting to divide by $r$)

$$theta'= 1$$

I have asked this question before on this site in two different ways and have never gotten an answer that resolves the matter, but that could just be my denseness!

share|improve this answer

edited 1 hour ago

answered 2 hours ago

MooMoo

7511515

$endgroup$

rdot[x1_, x2_] := (x1 (μ x1 - x2 - σ x1 (x1^2 + x2^2)) + x2 (x1 + μ x2 - σ x2 (x1^2 + x2^2)))/r

rdot[r Cos[t], r Sin[t]] // FullSimplify

thetadot[x1_,x2_]:=(x1 (x1+μ x2-σ x2 (x1^2+x2^2)) - x2(μ x1-x2-σ x1 (x1^2+x2^2)))/r^2

thetadot[r Cos[t], r Sin[t]] // FullSimplify

share|improve this answer

edited 1 hour ago

answered 2 hours ago

MooMoo

7511515

answered 2 hours ago

MooMoo

7511515

answered 2 hours ago

MooMoo

7511515