Show Posts

This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.


Messages - Brian Bi

Pages: [1] 2
1
Ch 9 / Re: 9.3 problem 18
« on: March 25, 2013, 01:28:57 PM »
So it is a spiral point but I didn't zoom in closely enough?

2
Ch 9 / 9.3 problem 18
« on: March 25, 2013, 12:32:43 AM »
I'm having some trouble getting this problem to work out. There are four critical points: (0,0), (2, 1), (-2, 1), and (-2, -4). At the critical point (-2, -4), the Jacobian is \begin{pmatrix} 10 & -5 \\ 6 & 0 \end{pmatrix} with eigenvalues $5 \pm i\sqrt{5}$. Therefore it looks like it should be an unstable spiral point. However, when I plotted it, it looked like a node. Has anyone else done this problem?

http://www.math.psu.edu/melvin/phase/newphase.html

3
MidTerm / Re: MT Problem 3
« on: March 07, 2013, 12:19:03 AM »
I wrote that undetermined coefficients does not apply because the ODE does not have constant coefficients.

4
MidTerm / Re: MT Problem 2b
« on: March 06, 2013, 11:10:43 PM »
By inspection, $y = 1/2$ is a solution.


Darn, why did I not see this...
Just goes to show that slowing down during a test and looking at the question with a calm mind can do wonders
I also did not see this during the test. Is there anyone who did?

5
MidTerm / Re: MT Problem 2b
« on: March 06, 2013, 10:01:33 PM »
By inspection, $y = 1/2$ is a solution.

6
Reading Week Challenge / Re: Reading Week Bonus problem 2
« on: February 16, 2013, 05:10:49 PM »
Assume that $z(x) > 0$ on $(x_0, x_1)$. We will show a contradiction.

Since $y(x) > 0$ and $z(x) > 0$ on $(x_0, x_1)$, recall the result of Bonus Problem 1:
\begin{equation}
W'(x) = (q(x) - Q(x))y(x)z(x)
\end{equation}
We know that $W'(x) \leq 0$ on $(x_0, x_1)$. By the mean value theorem, $W(x_0) \geq W(x_1)$. Consider what happens if we have $W(x_0) = W(x_1)$. Then if for some $x \in (x_0, x_1)$ we had $W(x) \neq W(x_0)$, then by the mean value theorem $W'$ would have to be positive on either $(x_0, x)$ or $(x, x_1)$, a contradiction. Conclude that if $W(x_0) = W(x_1)$, then $W$ is constant on $(x_0, x_1)$, implying that $q - Q$ vanishes identically on $(x_0, x_1)$. We proceed under the assumption that this is not the case, so that $W(x_0) > W(x_1)$.

Writing out the Wronskians explicitly gives
\begin{align}
\begin{vmatrix} y(x_0) & z(x_0) \\ y'(x_0) & z'(x_0) \end{vmatrix} &> \begin{vmatrix} y(x_1) & z(x_1) \\ y'(x_1) & z'(x_1)  \end{vmatrix} \\
\begin{vmatrix} 0 & z(x_0) \\ y'(x_0) & z'(x_0) \end{vmatrix} &> \begin{vmatrix} 0 & z(x_1) \\ y'(x_1) & z'(x_1) \end{vmatrix} \\
-z(x_0) y'(x_0) &> -z(x_1) y'(x_1) \label{eq:foo}
\end{align}
Now, because $y(x_0) = 0$ and $y(x_0) > 0$ on $(x_0, x_1)$, we must have $y'(x_0) \geq 0$. Likewise $y'(x_1) \leq 0$. We thus consider the following:
Case 1: If $z(x_0) > 0$ and $z(x_1) > 0$, then $-z(x_0) y'(x_0) \leq 0$ and $-z(x_1) y'(x_1) \geq 0$; together with $(\ref{eq:foo})$ we obtain the desired contradiction.
Case 2: If either $z(x_0) < 0$ or $z(x_1) < 0$, this contradicts the assumption that $z > 0$ on $(x_0, x_1)$ since $z$ is continuous.
Case 3: If $z(x_0) = 0$, then the RHS of $(\ref{eq:foo})$ must be negative, implying $z(x_1) < 0$. As in Case 2, find a contradiction.
Case 4: If $z(x_1) = 0$, then the LHS of $(\ref{eq:foo})$ must be positive, implying $z(x_0) < 0$. As in Case 2, find a contradiction.

So it is impossible for $z$ to be strictly positive on $(x_0, x_1)$; this implies that $z$ either has a zero, changes sign, or is strictly negative on $(x_0, x_1)$. The latter is impossible by an analysis analogous to the foregoing; and a sign change implies the presence of a zero anyway by the intermediate value theorem, so we are done.

7
Reading Week Challenge / Re: Reading Week Bonus problem 1
« on: February 16, 2013, 02:41:29 PM »
Is the second equation supposed to read: $z'' + Q(x)z = 0$ ?

Write out the Wronskian explicitly:
\begin{equation}
W[y,z](x) = \begin{vmatrix} y & z \\ y' & z' \end{vmatrix} = yz' - y'z
\end{equation}

Differentiate w.r.t. $x$:
\begin{equation}
W' = \frac{d}{dx} (yz' - y'z) = y'z' + yz'' - y''z - y'z' = yz'' - y''z
\end{equation}

From the equations (1) and (2) we obtain $y'' = -qy$ and $z'' = -Qz$ so
\begin{equation}
W' = y(-Qz) - (-qy)z = (q-Q)yz \leq 0
\end{equation}
where the inequality follows from $Q > q$ and $y, z > 0$.

8
Ch 4 / Re: Bonus problem for week 5b
« on: February 08, 2013, 01:31:32 PM »
The general solution to a homogeneous linear ODE with constant coefficients is given by
\begin{equation}
y = \sum_{i=1}^n \sum_{j=0}^{p_i - 1} A_{ij} t^j e^{r_i t} \label{eqn:general}
\end{equation}
where $r_1, r_2, ..., r_n$ are the distinct roots of the characteristic polynomial and $p_i$ is the multiplicity of the $i$th distinct root.

We see that a term of the form $Ae^t$ is given by $j = 0$ and $r_i = 1$ in ($\ref{eqn:general}$), and that a term of the form $Bte^{-t}$ is likewise given by $j = 1$ and $r_i = -1$. Since $j = 1$, the multiplicity of the root $r_i = -1$ must be at least 2. So the minimal characteristic polynomial ought to be $$(x - 1)(x + 1)^2 = x^3 + x^2 - x - 1$$Conclude that the desired ODE is $$y''' + y'' - y' - y = 0$$

9
Ch 4 / Re: Bonus problem for week 5a
« on: February 08, 2013, 01:23:31 PM »
We can tackle this by "brute force":
\begin{align}
0 &= y_1'' + py_1' + qy_1 \\
&= (e^t)'' + p(e^t)' + q(e^t) \\
&= e^t + pe^t + qe^t \\
&= e^t(1 + p + q)
\end{align}
and dividing through by $e^t$ we obtain $1 + p + q = 0$.
Also
\begin{align}
0 &= y_2'' + py_2' + qy_2 \\
&= (te^{-t})'' + p(te^{-t})' + q(te^{-t}) \\
&= (e^{-t} - te^{-t})' + p(e^{-t} - te^{-t}) + qte^{-t} \\
&= (-e^{-t} - e^{-t} + te^{-t}) + pe^{-t} - pte^{-t} + qte^{-t} \\
&= e^{-t}((1-t)p + tq + (t-2)
\end{align}
and dividing through by $e^{-t}$ we obtain $(1-t)p + tq + (t-2) = 0$.

Solving this linear system for $p$ and $q$ we arrive at
\begin{align}
p &= -\frac{2}{2t-1} \\
q &= \frac{3-2t}{2t-1}
\end{align}

so we conclude that the desired ODE is $$(2t-1)y'' - 2y' - (2t-3)y = 0$$

10
Technical Questions / Re: Preview
« on: February 03, 2013, 07:53:50 PM »
However I don't think you are correct about Chrome: Preview works for me in Chrome and it is Webkit. I just checked with Opera--Preview works, and with Maxthon--Preview works too!

I didn't say that it doesn't work---I said that it reloads the page when you click "Preview". If everything were working properly then it wouldn't have to reload the page, but would insert the preview directly into the current page. This is evidence that none of the browsers really do what they're supposed to do. Only Firefox has a weird bug that causes the preview not to work at all. Even still, fixing the XML would improve the user experience on all browsers.

11
Technical Questions / Re: Preview
« on: February 03, 2013, 01:38:36 AM »
After two hours of debugging, I've figured out the problem. The async endpoint is not sending back valid XML. Instead it sends something like this:

Code: [Select]
<script type="text/javascript" src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS_HTML">
    MathJax.Hub.Config({
        tex2jax: {
            inlineMath: [ ['$','$'] ],
            displayMath: [ ['$$','$$'] ],},
        TeX: {equationNumbers: {
            autoNumber: "AMS"
            }, Augment: {  Definitions: {
           macros: {
             overbracket:  ['UnderOver','23B4',1],
             underbracket: ['UnderOver','23B5',1],
           }
         }}},
    });
</script>
<script type="text/javascript" src="http://code.jquery.com/jquery-latest.min.js"></script> <link href="http://aimath.org/knowlstyle.css" rel="stylesheet" type="text/css" /> <script type="text/javascript" src="http://aimath.org/knowl.js"></script><?xml version="1.0" encoding="UTF-8"?>
<smf>
<preview>
<subject><![CDATA[Re: Previous years resources]]></subject>
<body><![CDATA[testing]]></body>
</preview>
<errors serious="0" topic_locked="0">
<caption name="guestname" class="" />
<caption name="email" class="" />
<caption name="evtitle" class="" />
<caption name="subject" class="" />
<caption name="question" class="" />
</errors>
<last_msg>1111</last_msg>
</smf>

This is not valid XML because it contains multiple root nodes. Only the <smf> node is actually desired in the response, I believe.

Actually, both Chrome and Firefox are "broken" because of this. Chrome sees that the XML is invalid and falls back to reloading the page (which it shouldn't have to do if the XML was correct, I think). On the other hand, Firefox thinks that XML was returned, and chokes on trying to process it.

12
Quiz 2 / Re: Day Section, Question 1
« on: February 01, 2013, 05:43:49 PM »
I wonder, isn't that when 0<α<1, e(α-1)t tend to 0 as t tend to infinity, and eαt tend to infinity? As a result, isn't that all nonzero solutions become unbounded as t tend to zero when α>0 instead of α>1?
The coefficient on $e^{\alpha t}$ might be zero, so you can have solutions that are just $A e^{(\alpha-1)t}$. These will decay to zero as $t \to \infty$.

13
Technical Questions / Preview
« on: February 01, 2013, 10:44:54 AM »
The preview post feature doesn't work at all. It just gets stuck at "Fetching preview..."

14
Ch 3 / Re: Problem of the week 4b
« on: February 01, 2013, 12:13:33 AM »
Add and subtract the first and second equations to obtain:
\begin{align}
(y+z)'' + K(y+z) &= 0 \label{added} \\
(y-z)'' + (K+2L)(y-z) &= 0 \label{subtracted}
\end{align}
Since $K, L > 0$, both equations are of the form $u'' + \omega^2 u = 0$, with general solution $u = A \cos (\omega t) + B \sin (\omega t)$. So the general solution to $(\ref{added})$ is
\begin{equation}
y+z = A \cos (\sqrt{K} t) + B \sin (\sqrt{K} t)
\end{equation}
and the general solution to $(\ref{subtracted})$ is
\begin{equation}
y - z = C \cos (\sqrt{K+2L} t) + D \sin(\sqrt{K+2L}t)
\end{equation}
Using the identities $y = \frac{1}{2}((y+z)+(y-z))$ and $z = \frac{1}{2}((y+z)-(y-z))$ we obtain the general solution to $(\ref{eq-1})$:
\begin{align}
y &= A' \cos(\omega_1 t) + B' \sin(\omega_1 t) + C' \cos(\omega_2 t) + D' \sin(\omega_2 t) \\
z &= A' \cos(\omega_1 t) + B' \sin(\omega_1 t) - C' \cos(\omega_2 t) - D' \sin(\omega_2 t)
\end{align}
where $A' = A/2$ and so on, and the frequencies are $\omega_1 = \sqrt{K}$, and $\omega_2 = \sqrt{K+2L}$.

15
Ch 3 / Re: Problem of the week 4a
« on: January 31, 2013, 07:40:19 PM »
The amplitude $\sqrt{A^2 + B^2}$ is then:
\begin{align}
\sqrt{A^2+ B^2} &= \left[\left(\frac{\omega^2 - \beta^2}{(\omega^2-\beta^2)^2 + (\alpha \beta)^2}\right)^2 + \left(\frac{\alpha\beta}{(\omega^2 -\beta^2)^2 + (\alpha \beta)^2}\right)^2\right]^{1/2} \\
&= \left[\frac{(\omega^2 - \beta^2)^2 + (\alpha\beta)^2}{((\omega^2-\beta^2)^2+(\alpha\beta)^2)^2}\right]^{1/2} \\
&= ((\omega^2 - \beta^2)^2 + (\alpha \beta)^2)^{-1/2}
\end{align}
This will be maximized if we minimize the expression $(\omega^2 - \beta^2)^2 + (\alpha \beta)^2$, which we expand to give
$\omega^4 + (\alpha^2 - 2\omega^2) \beta^2 + \beta^4$. This is a quadratic in $\beta^2$ with positive leading coefficient, so it has
a global minimum at $\beta^2 = -\frac{1}{2}(\alpha^2 - 2\omega^2)$, or $\beta = \sqrt{\omega^2 - \alpha^2/2}$. This value of $\beta$, the resonance frequency, is real and nonzero only if $\alpha < \sqrt{2} \omega$, so resonance occurs only in underdamped systems (and only those that are sufficiently underdamped).

Pages: [1] 2