Questions tagged [inverse-problem]
The inverse-problem tag has no usage guidance.
81
questions
1
vote
1
answer
75
views
Solving inverse problem using black box implementation of the kernel
My question is related to Solving regularized least squares problem using black-box computation of $\mathbf{A}\mathbf{x}$ and $\mathbf{A}^T\mathbf{x}$.
In case, the problem is formulated as:
\begin{...
1
vote
1
answer
53
views
Matrix-vector multiplication representation of Total Variation function
I'm reading a paper - Total Variation Superiorized Conjugate Gradient Method for Image Reconstruction on total variation regularization and conjugate gradients. In page $3$, the authors define the ...
3
votes
1
answer
102
views
Solving regularized least squares problem using black-box computation of $\mathbf{A}\mathbf{x}$ and $\mathbf{A}^T\mathbf{x}$
Let $\mathbf{A} \in \mathbb{R}^{n \times n}$. I'm working in a problem where I have a black-box algorithmic solution to compute the products $\mathbf{A}\mathbf{x}$ and $\mathbf{A}^T \mathbf{x}$ given ...
8
votes
1
answer
124
views
The Different Solutions for Filter Coefficients Estimation for Periodic Convolution and Full Convolution
As a continuation of the question Least Squares Solution Using the DFT vs Wiener Hopf Equations raised by Dan Boschen.
The question is, given the model:
$$ \boldsymbol{y} = \boldsymbol{h} * \...
5
votes
1
answer
71
views
Academic reference for a specific type of regularized inverse filtering
Let $y(t) = (h * x)(t) + n(t)$ be some observed signal where $h(t)$ is some filter / impulse response, $x(t)$ is some input signal we are interested in, and $n(t)$ is noise.
In order to recover $x$ ...
0
votes
2
answers
75
views
Get the inverse transfer function from the measured response
What is a numerically stable way to obtain the inverse transfer function from a measured response?
I have a system that shows a low-pass behavior. I would like to increase the bandwidth by some form ...
0
votes
1
answer
106
views
Reconstructing a signal from a Nyquist plot
I have a system which is like a blackbox which has just one input which could be a sinusoidal wave which is a sum of a range of frequencies, now the problem is that I dont have the time-domain output ...
-1
votes
1
answer
140
views
How to compute the inverse Z-transform
How to compute the inverse Z-transform of the form $$ G(z)=\frac{z^{2n}}{a(z^{2n})+b(z^n)+c} $$
I started by taking
$$ F(z^n)=G(z) $$
so
$$ \frac{F(z)}{z}=\frac{z}{az^2+bz+c}$$
This can be solved ...
2
votes
1
answer
90
views
Is the system $y\left(t\right)=\int _{t^3-1}^{t^3}x\left(s\right)ds\:$ invertible? [duplicate]
I have the following system: $$y\left(t\right)=\int _{t^3-1}^{t^3}x\left(s\right)ds\:$$
I was told to determine if it's invertible system, casual system, memoryless system and linear system. I was ...
1
vote
2
answers
710
views
How to use inverse 2D Fourier transform to reconstruct the original image?
I have managed to get the forward Fourier transform of an image to the frequency space like so:
But I cannot for the life of me reconstruct the original image from the inverse Fourier transform of ...
10
votes
1
answer
152
views
Image / Video Upscaling (Super Resolution) Algorithm Explanation (Image and Video Upscaling from Local Self Examples)
So, I'm trying to implement the classical algorithm described in this paper Image and Video Upscaling from Local Self-Examples and this presentation to serve as a baseline for comparison with AI/NN-...
3
votes
1
answer
59
views
Blind source separation for asynchronously observed mixture channels
Given your practical and theoretical expertise: Does ICA work reliably when applied to a multidimensional mixture (observation) $X = (X^1, \cdots, X^d)$ if the different channels $X^i$ of the ...
2
votes
2
answers
356
views
Is a neural network an adaptive filter?
I am confused as to the difference between neural networks and adaptive filters: As far as I understand it, "neural networks" are largely used for solving inverse problems, where an unknown ...
0
votes
1
answer
63
views
Identifying a signal by its power spectrum
Background: There is a method in optics for determining the electric field of light (both intensity and phase) via a three step process:
Add a known phase shift to the light (called the "...
5
votes
1
answer
604
views
Reconstruction of a Signal from Sub Sampled Spectrum by Compressed Sensing
Context
Attempting to reproduce an illustrative example of compressive sampling from Candes-Wakin 2008. Specifically, the L1 recovery of a sparse signal shown on pg 5 in Fig. 2.
Using my code (below), ...
7
votes
1
answer
78
views
Estimate the Image Using Multi Realizations of Its Convolution with a Known Filters Using Wiener Filter
Suppose we have a corrupted image $Y = H*X + \epsilon$ formed by taking an image $X$, convolving it with a point-spread function $H$, and adding gaussian noise $\epsilon$. Then we know that the Wiener ...
2
votes
2
answers
77
views
Cocktail Party Problem with a Single Signal of Data (Single Mic)
I have been doing some multimodal signal analysis, and sometimes ICA is used for detecting statistically independent components.
From my understanding, say if you have 2 sources and 2 receivers/...
7
votes
1
answer
302
views
Solve Efficiently the 1D Total Variation Regularized Least Squares Problem (Denoising / Deblurring)
How to solve a 1D Least Squares with Total Variation Regularization?
I know gradient based methods, I wonder how much faster / efficient I can get.
5
votes
1
answer
95
views
How Could One Accelerate the Convergence of the Least Mean Squares (LMS) Filter?
How can the convergence of an LMS filter be accelerated?
Can we do better than the vanilla algorithm?
0
votes
1
answer
49
views
Sparse FFT of ramp using zero padding
I would like to sparsely represent a linear function/ramp in the Fourier domain. In an attempt to improve the sparsity, I have zero padded it. With this padding, it is possible in the example I tried ...
6
votes
1
answer
282
views
Tikhonov Regularization for Complex Matrices
Tikhonov regularization is used to regularize ill-posed inverse problems if the matrix $A \in \mathbb{R}^{n,m}$ to be inversed has a high condition number. For example
$$
A=\begin{bmatrix}1&1\\
1&...
1
vote
3
answers
170
views
Inverse Fourier transform: where am I going wrong?
I am studying a course in signal processing, currently we are examining Fourier transforms. I got stuck on an exercise with an inverse Fourier transform.
I am supposed to find the inverse Fourier ...
4
votes
1
answer
153
views
Regularization for inverse filter design
Given a $2 \times 2$ matrix, $C$, suppose I want to compute a filter matrix $H = C^{-1}$ and that I need to add regularization for practical purposes (e.g., for an audio filter, regularization is ...
1
vote
0
answers
175
views
Explain the relationship between Tikhonov regularization, SVD , least squres, and the Wiener filter
I found in the
Wikipedia site for Tikhonov regularization
that the SVD for a Tikhonov regularized problem take us to the least squares regularized solution:
\begin{eqnarray}
\hat{x}= V D U^T b,
\end{...
7
votes
1
answer
73
views
Stochastic Methods for Image Deconvolution Problem
If we convolve an image with a point spread function and from the resulting image to find the input image, can we use any stochastic approaches? I feel like we will not be able to. A single image ...
2
votes
2
answers
440
views
Inverse system of sinc?
I'm doing some self-study for an important exam I'll have in late March and came across the following question:
So, using the convolution properties, if I want to find an identity system so that the ...
11
votes
2
answers
774
views
Can Principal Component Analysis (PCA) Solve the Cocktail Party Problem?
I'm looking into the cocktail party problem and trying to figure out whether something like Principal Component Analysis is enough to separate out all the various voices at the cocktail party into its ...
8
votes
3
answers
1k
views
Deconvolution of a 1D Time Domain Wave Signal Convolved with Series of Rect Signals
I have a synthesized signal (the bottom of the following figure), which is the convolution of the input signal (at the top) and the objective function (in the middle). The intention is to retrieve the ...
2
votes
0
answers
84
views
Adapting Richardson Lucy (RL) Deconvolution for Shot Noise Limited Coherent Imaging
I am an experimental physicist who is collecting a series of coherent imaging of trapped gas. If you are familiar with phase contrast imaging, you may understand what I mean by coherent imaging. The ...
2
votes
1
answer
329
views
An invertible system with memory
Suppose $\mathcal{L}$ be invertible system with memory. Does $\mathcal{L}^{-1}$ have memory necessarily?
Intuitively I think the answer is "yes". There are many examples showing that. For ...
2
votes
1
answer
92
views
In computed tomography (CT), why is 'Inverse Problem of Radon transform' studied?
As I know, the Radon transform is a very important tool in CT by Beer's law. Thus, finding $f(x)$ of Radon transform $Rf(L):=\int_{L} f(x)dl(x)$ is helpful in CT. Nowaday, the Filtered back-projection ...
5
votes
1
answer
83
views
Regularization for Inverse Problems using the Singular Value Decomposition (SVD)
I am reading these lecture notes on Optimisation and Inverse Problems in Imaging, and I have difficulties understanding how figures on page 20 (Figure 3.2) or page 21 (Figure 3.3).
Precisely, I don't ...
6
votes
3
answers
247
views
Deconvolving a 1d Signal Using a Lookup Table
assuming I measure a signal that has different PSFs per position in time.
for example:
...
0
votes
1
answer
756
views
Is Differentiation as a system, is an invertible system?
is the following system invertible?
as I understand it, invertible means finding an inverse function which should return back the original input from an output of the given system.
if so I ...
7
votes
2
answers
172
views
Estimating Convolution Input Under the Assumption of Sparsity and Constant Non Zero Values Using Compressive Sensing Approach
I was wondering about if there is compressive sensing algorithm to estimate the sparse vector where the number of non-zeros values and amplitude of every non-zeros value are known. For example, assume ...
0
votes
0
answers
17
views
Whats the Correct Approach to Estimate the PSF of a Moving Detector
In lab, I did a bunch of scans using a radiation source and a detector. My source emits a gaussian beam (whose dimensions I know), and my detector is a uniform 7mmx7mm square. These are stationary and ...
8
votes
3
answers
259
views
Estimating the Signal by Deconvolution with a Prior on the Filter Coefficients and the Signal Samples
Assume I have signal $y[n]$ which is a result of convolution between channel $h[n]$ and signal $x[n]$. which means:
$$y[n] = h[n] \ast x[n]$$ where $\ast$ is the convolution operation
The signal $...
7
votes
2
answers
267
views
Deconvolution of an Image Acquired by a Square Uniform Detector
So, I acquired some images by scanning a radiation source with a square detector like in the following gif.
Where the dashed grid represents reality, the 3x3 square my detector, and the 4x4 my ...
7
votes
4
answers
4k
views
Are all LTI systems invertible? If not, what is a good counterexample?
I have been trying to figure this out for a while now. Everywhere I have looked I could easily find examples of invertible LTI systems, but I could not find any counterexamples. Can anybody shed some ...
0
votes
0
answers
59
views
Solving equation with convolution
I have the measured signal $y(t)$ that can be modeled in the frequency domain as $Y(f)$:
$$Y(f) = X(f)\cdot A(f) - [X(f)\cdot B(f)] \ast C(f)$$
where $\ast$ is the convolution.
I know $A(f)$, $B(f)$,...
2
votes
0
answers
39
views
Extracting 'structure' post permutation
I have particle activity as shown in the left pane of animation below. The activity is clustered and it moves slowly. Sometimes these clusters merges together. On the right side of it, its shuffled ...
5
votes
2
answers
428
views
What Is the Relation Between Deblurring and Deconvolution in Computer Vision and Image Processing?
The deblurring problem can be modelled as follows
$$
f = \phi u + \epsilon, \; \epsilon \sim N(0, \sigma)
$$
where $\phi$ is a filter (e.g. a low-pass filter) and $\epsilon$ is a Gaussian noise.
In ...
0
votes
2
answers
138
views
Partial Fractions
Attached is image with solution and my attempt.
I am trying to calculate the coefficients for partial fractions expansion of the following:
$$ H(e^{j\omega}) = \frac{ \frac{1}{3} e^{-j2\omega} }{(1-\...
9
votes
3
answers
1k
views
Deconvolution of Synthetic 1D Signals - How To?
I convolved a square wave with a Gaussian wave using linear convolution. Can I get the original square wave back by deconvolving my output with the Gaussian function?
I took the FFT of both signals, ...
7
votes
1
answer
1k
views
Deconvolution of a 1D Signal with Known Kernel (Square Wave)
I have a signal measured from a radiation detector in a narrow beam of radiation. The peaks I get are quasi-gaussian in shape, see attached picture.
The signal is not a function of time, rather a ...
6
votes
1
answer
2k
views
What Are the Differences between Super Resolution, Denoising and Deblurring?
In the fields of computer vision and image processing, what are the differences between Super Resolution, Denoising and Deblurring?
6
votes
3
answers
2k
views
Deconvolution Using Complex Division in The Frequency Domain
Consider these two signals:
a = [1 1 0 0 0 0 0 0]
b = [1 0 1 0 0 0 0 0]
their convolution is
c = a * b = [1 1 1 1 0 0 0 0]
...
4
votes
2
answers
1k
views
Auto-correlation function, an inverse problem
$x[n]$ is a complex function $n=0,1,2,\cdots,L-1 $
we assume $x[n]$ is periodic in its index: $x[n+L]=x[n]$
Its auto-correlation function $C[n]$ is uniquely defined as:
$$
C[n]=\sum_{i=0}^{L-1} x[i+...
6
votes
1
answer
636
views
Intuitive Meaning of Regularization in Imaging Inverse Problems
Hello Every one I have been trying to understand the intuitive meaning of using a regularizer in images. Specifically what does the Total Variation regularizer do in images and how is it able to ...
0
votes
1
answer
157
views
Inverse FFT - Synch the Phase
Is there any way to synchronize phase of output of inverse DFT in each buffer? When I send the output of inverse DFT to the ...