Duffing oscillator












0












$begingroup$


I have to draw the phase space of this modification of Duffing oscillator:
begin{align}
dot x &= y,\
dot y &= x - x^3 - ay - (x^2)y quad text{ when } quad a=1 ; text{ and } ;a=5.
end{align}



I know that critical points are $(0, 0)$, $(-1, 0)$, and $(1, 0),$ and I have demonstrated that when $a =1$ I have a spiral the form $(pm1,0)$ and when $a = 5$ not.



Does someone know how to draw it?



Thank you.










share|cite|improve this question











$endgroup$












  • $begingroup$
    Improve your question, pls.
    $endgroup$
    – Word Shallow
    Jan 15 at 14:33










  • $begingroup$
    Why did you revert back the edit made by Key Flex?
    $endgroup$
    – BigbearZzz
    Jan 15 at 14:50










  • $begingroup$
    Do you mean draw with a computer program or do a sketch by hand if no computer is available?
    $endgroup$
    – LutzL
    Jan 15 at 16:11
















0












$begingroup$


I have to draw the phase space of this modification of Duffing oscillator:
begin{align}
dot x &= y,\
dot y &= x - x^3 - ay - (x^2)y quad text{ when } quad a=1 ; text{ and } ;a=5.
end{align}



I know that critical points are $(0, 0)$, $(-1, 0)$, and $(1, 0),$ and I have demonstrated that when $a =1$ I have a spiral the form $(pm1,0)$ and when $a = 5$ not.



Does someone know how to draw it?



Thank you.










share|cite|improve this question











$endgroup$












  • $begingroup$
    Improve your question, pls.
    $endgroup$
    – Word Shallow
    Jan 15 at 14:33










  • $begingroup$
    Why did you revert back the edit made by Key Flex?
    $endgroup$
    – BigbearZzz
    Jan 15 at 14:50










  • $begingroup$
    Do you mean draw with a computer program or do a sketch by hand if no computer is available?
    $endgroup$
    – LutzL
    Jan 15 at 16:11














0












0








0





$begingroup$


I have to draw the phase space of this modification of Duffing oscillator:
begin{align}
dot x &= y,\
dot y &= x - x^3 - ay - (x^2)y quad text{ when } quad a=1 ; text{ and } ;a=5.
end{align}



I know that critical points are $(0, 0)$, $(-1, 0)$, and $(1, 0),$ and I have demonstrated that when $a =1$ I have a spiral the form $(pm1,0)$ and when $a = 5$ not.



Does someone know how to draw it?



Thank you.










share|cite|improve this question











$endgroup$




I have to draw the phase space of this modification of Duffing oscillator:
begin{align}
dot x &= y,\
dot y &= x - x^3 - ay - (x^2)y quad text{ when } quad a=1 ; text{ and } ;a=5.
end{align}



I know that critical points are $(0, 0)$, $(-1, 0)$, and $(1, 0),$ and I have demonstrated that when $a =1$ I have a spiral the form $(pm1,0)$ and when $a = 5$ not.



Does someone know how to draw it?



Thank you.







ordinary-differential-equations






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited Jan 15 at 16:07









LutzL

60k42057




60k42057










asked Jan 15 at 14:30









CastiCasti

61




61












  • $begingroup$
    Improve your question, pls.
    $endgroup$
    – Word Shallow
    Jan 15 at 14:33










  • $begingroup$
    Why did you revert back the edit made by Key Flex?
    $endgroup$
    – BigbearZzz
    Jan 15 at 14:50










  • $begingroup$
    Do you mean draw with a computer program or do a sketch by hand if no computer is available?
    $endgroup$
    – LutzL
    Jan 15 at 16:11


















  • $begingroup$
    Improve your question, pls.
    $endgroup$
    – Word Shallow
    Jan 15 at 14:33










  • $begingroup$
    Why did you revert back the edit made by Key Flex?
    $endgroup$
    – BigbearZzz
    Jan 15 at 14:50










  • $begingroup$
    Do you mean draw with a computer program or do a sketch by hand if no computer is available?
    $endgroup$
    – LutzL
    Jan 15 at 16:11
















$begingroup$
Improve your question, pls.
$endgroup$
– Word Shallow
Jan 15 at 14:33




$begingroup$
Improve your question, pls.
$endgroup$
– Word Shallow
Jan 15 at 14:33












$begingroup$
Why did you revert back the edit made by Key Flex?
$endgroup$
– BigbearZzz
Jan 15 at 14:50




$begingroup$
Why did you revert back the edit made by Key Flex?
$endgroup$
– BigbearZzz
Jan 15 at 14:50












$begingroup$
Do you mean draw with a computer program or do a sketch by hand if no computer is available?
$endgroup$
– LutzL
Jan 15 at 16:11




$begingroup$
Do you mean draw with a computer program or do a sketch by hand if no computer is available?
$endgroup$
– LutzL
Jan 15 at 16:11










1 Answer
1






active

oldest

votes


















0












$begingroup$

This may help you



enter image description here



generated with python:



import numpy as  np
import matplotlib.pyplot as plt

xmin, xmax = -2.0, 2.0
ymin, ymax = -1.0, 1.0

x, y = np.meshgrid(np.linspace(xmin, xmax, 30, endpoint = True), np.linspace(ymin, ymax, 30, endpoint = True))

a = 1
u = y
v = x - x**3 - a * y - x**2 * y


plt.streamplot(x, y, u, v)
plt.show()





share|cite|improve this answer









$endgroup$













  • $begingroup$
    Thank you very much, your answer has been very helpful.
    $endgroup$
    – Casti
    Jan 15 at 17:16












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%2f3074496%2fduffing-oscillator%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$

This may help you



enter image description here



generated with python:



import numpy as  np
import matplotlib.pyplot as plt

xmin, xmax = -2.0, 2.0
ymin, ymax = -1.0, 1.0

x, y = np.meshgrid(np.linspace(xmin, xmax, 30, endpoint = True), np.linspace(ymin, ymax, 30, endpoint = True))

a = 1
u = y
v = x - x**3 - a * y - x**2 * y


plt.streamplot(x, y, u, v)
plt.show()





share|cite|improve this answer









$endgroup$













  • $begingroup$
    Thank you very much, your answer has been very helpful.
    $endgroup$
    – Casti
    Jan 15 at 17:16
















0












$begingroup$

This may help you



enter image description here



generated with python:



import numpy as  np
import matplotlib.pyplot as plt

xmin, xmax = -2.0, 2.0
ymin, ymax = -1.0, 1.0

x, y = np.meshgrid(np.linspace(xmin, xmax, 30, endpoint = True), np.linspace(ymin, ymax, 30, endpoint = True))

a = 1
u = y
v = x - x**3 - a * y - x**2 * y


plt.streamplot(x, y, u, v)
plt.show()





share|cite|improve this answer









$endgroup$













  • $begingroup$
    Thank you very much, your answer has been very helpful.
    $endgroup$
    – Casti
    Jan 15 at 17:16














0












0








0





$begingroup$

This may help you



enter image description here



generated with python:



import numpy as  np
import matplotlib.pyplot as plt

xmin, xmax = -2.0, 2.0
ymin, ymax = -1.0, 1.0

x, y = np.meshgrid(np.linspace(xmin, xmax, 30, endpoint = True), np.linspace(ymin, ymax, 30, endpoint = True))

a = 1
u = y
v = x - x**3 - a * y - x**2 * y


plt.streamplot(x, y, u, v)
plt.show()





share|cite|improve this answer









$endgroup$



This may help you



enter image description here



generated with python:



import numpy as  np
import matplotlib.pyplot as plt

xmin, xmax = -2.0, 2.0
ymin, ymax = -1.0, 1.0

x, y = np.meshgrid(np.linspace(xmin, xmax, 30, endpoint = True), np.linspace(ymin, ymax, 30, endpoint = True))

a = 1
u = y
v = x - x**3 - a * y - x**2 * y


plt.streamplot(x, y, u, v)
plt.show()






share|cite|improve this answer












share|cite|improve this answer



share|cite|improve this answer










answered Jan 15 at 15:11









caveraccaverac

14.8k31130




14.8k31130












  • $begingroup$
    Thank you very much, your answer has been very helpful.
    $endgroup$
    – Casti
    Jan 15 at 17:16


















  • $begingroup$
    Thank you very much, your answer has been very helpful.
    $endgroup$
    – Casti
    Jan 15 at 17:16
















$begingroup$
Thank you very much, your answer has been very helpful.
$endgroup$
– Casti
Jan 15 at 17:16




$begingroup$
Thank you very much, your answer has been very helpful.
$endgroup$
– Casti
Jan 15 at 17:16


















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%2f3074496%2fduffing-oscillator%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

Human spaceflight

Can not write log (Is /dev/pts mounted?) - openpty in Ubuntu-on-Windows?

張江高科駅