Dynamical system 4

Error

There was an error loading the page.

Contributors

Bill Foster

Feedback

From users who are members of Bill's workspace :

Bill Foster

said

Ready to use

8 years ago

History

○

Bill Foster 8 years ago

Gave some feedback: Ready to use

○

Bill Foster 8 years ago

Published this.

○

Bill Foster 11 years, 3 months ago

Created this as a copy of Dynamical system 3.

Name	Status	Author	Last Modified
Dynamical system 4	Ready to use	Bill Foster	16/11/2020 10:02
Dynamical system 6:Centre.	draft	Bill Foster	16/11/2020 10:02
Dynamical system 7: Stable spiral.	Ready to use	Bill Foster	16/11/2020 10:02
Dynamical system 9: Saddle.	Ready to use	Bill Foster	16/11/2020 10:02
Bill's copy of Dynamical system 6:Centre.	Ready to use	Bill Foster	16/11/2020 10:02

There are 7 other versions that do you not have access to.

Question statement

Consider the two-dimensional dynamical system. Complex eigenvalues case , so centres or spirals.

$\begin{align}\dot{x}&=\simplify[std]{{a}*x+{b}*y},\\\dot{y}&=\simplify[std]{{c}*x+{d}*y}.\end{align}$

{ode1()}

The above diagram shows the plot of $(x(t),y(t))$ .

At $t=0$ we have initally $x=-5,\;\;y=5$ . Moving the point gives phase diagrams for the following initial values at $t=0$ :

$x=\;$ $y=\;$

Dotted lines give points at which $\dot{x}=0,\;\dot{y}=0$ .

The plots below are for $x(t),\;x(0)=-1$ and $y(t),\;y(0)=1$ , in black and red respectively where the horizontal axis is $t$ .

{ode(xr,yr)}

Give any introductory information the student needs.

			Name	Type	Generated Value

lam1	number	0.09 - 0.75i
r	integer	1
be	integer	-1
c	number	-1.6908
m	number	-1.2681
al	number	1.12
delta	number	0.5706
natures	list	List of 5 items
vec2	number	1.12 + i
disc	number	-2.25
vec1	number	1.12 - i
d	number	0.93
lam2	number	0.09 + 0.75i
a	number	-0.75
b	number	0.75
tau	number	0.18
revec1	number	1.12
yr	integer	100
rorc	integer	1
s	integer	1
naturetext	string	because $\mathrm{Im}(\lambda_1
xr	integer	10
revec2	number	1.12
nature	string	unstable spiral

Automatically regenerate variables?

Name

Data type

Value

xxxxxxxxxx
 
b*vec1+a

JME function reference

Description

Describe what this variable represents, and list any assumptions made about its value.

Can an exam override the value of this variable?

Generated value: `number`

0.09 - 0.75i

← Depends on:

This variable doesn't seem to be used anywhere.

Parts

Part a) Gap-fill
1. Gap Gap 0. Number entry
  Add
2. Gap Gap 1. Number entry
  Add
3. Gap Gap 2. Number entry
  Add
4. Gap Gap 3. Number entry
  Add
Part b) Gap-fill
1. Gap Gap 0. Number entry
  Add
2. Gap Gap 1. Number entry
  Add
Part c) Gap-fill
1. Gap Gap 0. Mathematical expression
  Add
2. Gap Gap 1. Mathematical expression
  Add
Part d) Gap-fill
1. Gap Gap 0. Mathematical expression
  Add
2. Gap Gap 1. Mathematical expression
  Add
Part e) Choose one from a list

Gap-fill

Ask the student a question, and give any hints about how they should answer this part.

The system can be written in the form $\dot{\boldsymbol{x}}=\mathsf{A}\boldsymbol{x}$, where $\boldsymbol{x}=\pmatrix{x,y}^\mathsf{T}$.

Identify the components of the matrix $\mathsf{A}$.

$\mathsf{A}=\Bigg($	Gap 0.	Gap 1.	$\Bigg)$
	Gap 2.	Gap 3.

Advice

The elements of the matrix $\mathsf{A}$ are the coefficients of $x$ and $y$ in the system of equations, hence

\[\mathsf{A}=\pmatrix{\var{a} & \var{b}\\\var{c} & \var{d}}.\]

The trace $\tau$ and determinant $\delta$ of a matrix

\[\mathsf{A}=\pmatrix{a & b\\c & d}\]

are given by $\tau=a+d$ and $\delta=ad-bc$.

A straight forward substitution for $a=\var{a}$, $b=\var{b}$, $c=\var{c}$, and $d=\var{d}$ reveals that $\tau=\var{tau}$ and $\delta=\var{delta}$.

The eigenvalues $\lambda$ of a matrix $\mathsf{A}$ can be calculated from the expression

\[\mathsf{A}\boldsymbol{x}=\lambda\boldsymbol{x},\]

which is equivalent to

\[\left(\mathsf{A}-\lambda\mathsf{I}\right)\boldsymbol{x}=\boldsymbol{0},\]

where $\mathsf{I}$ is the identity matrix.

This linear system has a solution only when

\[\det\left(\mathsf{A}-\lambda\mathsf{I}\right)=0,\]

i.e. when $(a-\lambda)(d-\lambda)-bc=0$. This leads to $\lambda^2-(a+d)\lambda+(ad-bc)=0$, which is equivalent to

\[\lambda^2-\tau\lambda+\delta=0.\]

The eigenvalues $\lambda_1$ and $\lambda_2$ of the matrix $\mathsf{A}$ are therefore given by

\[\lambda_{1,2}=\frac{-\tau\pm\sqrt{\tau^2-4\delta}}{2}.\]

Following through this calculation for the specific values of $a$, $b$, $c$, and $d$ in this question leads to $\lambda_1=\var{lam1}$ and $\lambda_2=\var{lam2}$.

The eigenvectors corresponding to the eigenvalues $\lambda_1$ and $\lambda_2$ can be found by solving the linear system of equations given by

\[\mathsf{A}\boldsymbol{x}=\lambda\boldsymbol{x},\]

\[\begin{align}ax+by&=\lambda x,\\cx+dy&=\lambda y,\end{align}\]

where $\lambda$ is either $\lambda_1$ or $\lambda_2$.

Because this is an eigenvalue problem, we can arbitrarily choose one of the eigenvector components to be $1$, the $x$-component in this case. The $y$-component can then be found from

\[y=\frac{x(\lambda-a)}{b}\quad\text{or}\quad y=\frac{cx}{\lambda-d},\]

with $x=1$. Both expressions are equivalent, and will lead to the same value for $\lambda$.

Making the necessary substitutions reveals that the $y$-component of the eigenvector corresponding to $\lambda_1$ is $\var{vec1}$, and the $y$-component of the eigenvector corresponding to $\lambda_2$ is $\var{vec2}$.

The nature of the fixed point at the origin can be determined by examining the eigenvalues. If the eigenvalues are real, then the sign of each eigenvalue determines the nature of the fixed point. If the eigenvalues are complex, then the sign of the real part of each eigenvalue determines the nature of the fixed point.

In this case, {naturetext}.

Give a worked solution to the whole question.

Use this tab to check that this question works as expected.

There was an error which means the tests can't run:

Part	Test	Passed?
Gap-fill
Gap-fill		Hasn't run yet
Number entry
Number entry		Hasn't run yet
Number entry
Number entry		Hasn't run yet
Number entry
Number entry		Hasn't run yet
Number entry
Number entry		Hasn't run yet
Gap-fill
Gap-fill		Hasn't run yet
Number entry
Number entry		Hasn't run yet
Number entry
Number entry		Hasn't run yet
Gap-fill
Gap-fill		Hasn't run yet
Mathematical expression
Mathematical expression		Hasn't run yet
Mathematical expression
Mathematical expression		Hasn't run yet
Gap-fill
Gap-fill		Hasn't run yet
Mathematical expression
Mathematical expression		Hasn't run yet
Mathematical expression
Mathematical expression		Hasn't run yet
Choose one from a list
Choose one from a list		Hasn't run yet

This question is used in the following exams:

Loading...

Error

Metadata

Taxonomy: mathcentre

Taxonomy: Kind of activity

Taxonomy: Context

Contributors

Feedback

History

Comment

Checkpoint description

Bill Foster 8 years ago

Bill Foster 8 years ago

Bill Foster 11 years, 3 months ago

Error in variable testing condition

Error

Warning

Generated value: `number`

Parts

Gap-fill

Loading...

Error

Dynamical system 4

Metadata

Taxonomy: mathcentre

Taxonomy: Kind of activity

Taxonomy: Context

Contributors

Feedback

History

Comment

Checkpoint description

Bill Foster 8 years ago

Bill Foster 8 years ago

Bill Foster 11 years, 3 months ago

Error in variable testing condition

Error

Warning

Generated value: number

Parts

Gap-fill

Generated value: `number`