Solve these two coupled first-order differential equations and sketch the flow

Lambda96
Messages
233
Reaction score
77
Homework Statement
solve the two coupled first-order differential equations ##\textbf{f}(\textbf{x}(t))## and sketch the flow ##\phi_t(\textbf{x})##
Relevant Equations
none
Hi,

unfortunately, I have a problem to solve the following task

Bildschirmfoto 2023-07-04 um 20.08.10.png


The equation looks like this:

$$\left(\begin{array}{c} \frac{d}{dt} x(t) \\ \frac{d}{dt} y(t) \end{array}\right)=\left(\begin{array}{c} -a y(t) \\ x(t) \end{array}\right)$$

Since the following is true ##\frac{d}{dt} y(t)=x(t)## I substituted ##x(t)## into the first equation on the left hand side, obtaining a 2nd order differential equation, i.e. ##y''(t)=-ay(t)## to solve this differential equation I then used the Ansatz ##y(t)=e^{\lambda t}## and obtained the following solution.

$$y(t)=c_1 \ e^{i \sqrt{a}t}-c_2 \ e^{-i \sqrt{a}t}$$

I then obtain x(t) using ##\frac{d}{dt} y(t)=x(t)##

$$y(t)=i \ \sqrt{a} \ c_1 \ e^{i \sqrt{a}t}-i \ \sqrt{a} \ c_2 \ e^{-i \sqrt{a}t}$$

Unfortunately, I am now a bit unsure about task b. For the flow ##\phi_t(\textbf{x})## I would now simply draw the vector ##\textbf{f}(\textbf{x}(t))=\left(\begin{array}{c} -a y(t) \\ x(t) \end{array}\right)## with my solution from task a for ##y(t)## and ##x(t)##.

Unfortunately I have problems to draw the vector, because I don't know what the constants ##c_1## and ##c_2## are and unfortunately my solution from task a is complex, so do I have to plot only the real part or is my solution from task a wrong?
 
Physics news on Phys.org
The flow map doesn't assume one particular set of ICs. That would make only one line in the flow map.

Since ##a \in \mathbb R ## your solutions should be real too. This is a physically realizable system. This is a constraint on ##c_1## and ##c_2##. Your ICs will be various values of ##x(0)## and ##y(0)## which will then determine ##c_1## and ##c_2##.

This may help:
http://www.math.sjsu.edu/~simic/Fall05/Math134/flows.pdf

But there's lots of stuff on the web. Steve Strogatz has a nice set of lectures on YouTube in dynamic systems that does flow maps, but it's long.

Look for features like where the derivative(s) are zero to find stable points and asymptotes.

Sorry, it's been a long time since I really did this sort of work, and I'm too lazy to actually solve this system.
 
Last edited:
Thanks DaveE for your help and the links 👍👍
 
  • Like
Likes berkeman and DaveE
##|\Psi|^2=\frac{1}{\sqrt{\pi b^2}}\exp(\frac{-(x-x_0)^2}{b^2}).## ##\braket{x}=\frac{1}{\sqrt{\pi b^2}}\int_{-\infty}^{\infty}dx\,x\exp(-\frac{(x-x_0)^2}{b^2}).## ##y=x-x_0 \quad x=y+x_0 \quad dy=dx.## The boundaries remain infinite, I believe. ##\frac{1}{\sqrt{\pi b^2}}\int_{-\infty}^{\infty}dy(y+x_0)\exp(\frac{-y^2}{b^2}).## ##\frac{2}{\sqrt{\pi b^2}}\int_0^{\infty}dy\,y\exp(\frac{-y^2}{b^2})+\frac{2x_0}{\sqrt{\pi b^2}}\int_0^{\infty}dy\,\exp(-\frac{y^2}{b^2}).## I then resolved the two...
Hello everyone, I’m considering a point charge q that oscillates harmonically about the origin along the z-axis, e.g. $$z_{q}(t)= A\sin(wt)$$ In a strongly simplified / quasi-instantaneous approximation I ignore retardation and take the electric field at the position ##r=(x,y,z)## simply to be the “Coulomb field at the charge’s instantaneous position”: $$E(r,t)=\frac{q}{4\pi\varepsilon_{0}}\frac{r-r_{q}(t)}{||r-r_{q}(t)||^{3}}$$ with $$r_{q}(t)=(0,0,z_{q}(t))$$ (I’m aware this isn’t...
It's given a gas of particles all identical which has T fixed and spin S. Let's ##g(\epsilon)## the density of orbital states and ##g(\epsilon) = g_0## for ##\forall \epsilon \in [\epsilon_0, \epsilon_1]##, zero otherwise. How to compute the number of accessible quantum states of one particle? This is my attempt, and I suspect that is not good. Let S=0 and then bosons in a system. Simply, if we have the density of orbitals we have to integrate ##g(\epsilon)## and we have...
Back
Top