Questions tagged [stability]
The study of the propagation of errors in a numerical algorithm.
13
questions
3
votes
1
answer
118
views
Franco Brezzi's didactic paper on the Inf-Sup condition
I am trying to recover the title of a very didactic paper by Franco Brezzi in which he explains the Inf-Sup condition with a very simple case. As I recall, it was a simple Laplacian equation in 1D ...
4
votes
1
answer
161
views
Numerically stable computation of $x^T A x$
We have a large sparse symmetric positive-definite matrix $A \in \mathbb R^{N \times N}$ and a vector $x \in \mathbb R^N$. How do I practically compute the inner product $x^T A x$ when the matrix $A$ ...
0
votes
0
answers
47
views
2nd-order backward difference approximation for temporal terms of Euler–Bernoulli beam
I am trying to solve the Euler–Bernoulli beam numerically in structural dynamics analysis:
$$\frac{\partial^2w(x,t)}{\partial t^2}= \frac {q(x,t)}{\mu(x)}-\frac {1}{\mu(x)}\frac {\partial^2}{\partial ...
2
votes
2
answers
76
views
How to improve and stabilize this code simulating a damped mass-spring oscillator? Runge-Kutta?
I wrote the following function which is simulating a damped mass-spring oscillator. It is being driven by the audio sample input at 44.1 kHz sampling to create the same effect as a resonant bandpass ...
0
votes
2
answers
2k
views
When can I use finite differences for differentiation?
Finite differences are usually used to integrate ODE's and PDE's. However, sometimes they can be used for differentiation which I illustrated simply by using the Matlab code below to differentiate the ...
2
votes
0
answers
80
views
Preventing an Overflow in Exponential Integrating Factor
The following is the Discrete Spectral Vorticity Evolution PDE for incompressible flow:
$$ \frac{\partial \Omega_{pq}}{\partial t} = \nu \left( \frac{\partial^2 \Omega_{pq}}{\partial x^2} + \frac{\...
0
votes
0
answers
45
views
Which numerical method can I use to solve this system of hyperbolic PDEs?
Backround
The mathematical model I am trying to numerically solve models wave propagation inside a cylinder with specific material properties suited for dynamic loading. The cylinder's upper base is ...
1
vote
1
answer
122
views
Step size constraint in Euler backward
I am dealing with an assignment in MATLAB. It has to do with 'self-driving' cars which are driving in-front/behind eachother. Assuming M cars on a single-lane road, each car adjusts its speed based on ...
3
votes
1
answer
188
views
Stability of Euler forward method
I am trying to solve a linear system of ODEs of the form:
$$ \frac{du}{dt} = A u, \quad u(0)=k$$
where $A$ is a 2x2 matrix and $u(t)$ is a 2x1 column vector. I want to solve this numerically, using ...
2
votes
1
answer
146
views
Can this finite difference dispersion be eliminated somehow?
I am trying to solve the wave equation
$$ {\partial ^2u(t,x) \over \partial x^2} = {\partial ^2u(t,x) \over \partial t^2} \tag1 $$
with the following boundary and initial conditions:
$$ {\partial u \...
5
votes
0
answers
90
views
optimization scaling techniques
Consider a convex QP of the form
$$
\min_x \bigl\{\tfrac12 x^\top Q x + q^\top x : Ax\leq b\bigr\}\tag{P}
$$
with dual
$$
\min_y \bigl\{\tfrac12 (A^\top y + q)^\top Q^{\dagger} (A^\top y+q) + b^\top y ...
6
votes
1
answer
1k
views
Why is matrix inversion unstable when svd is stable?
I've heard that matrix is inversion is unstable whereas the SVD is stable.
Now, if $A$ is an invertible matrix, then its SVD is
$$
A = USV^T
$$
Then wouldn't it's inverse just be
$$
A^{-1} = (USV^T)^{...
1
vote
0
answers
144
views
Stability for the 2d diffusion equation
I'm reading this set of notes on numerical stability of the diffusion equation. I can't seem to derive equation (7.17) from the previous equation. I've tried using $\sin^2(\theta)=\frac{1-\cos(2\theta)...