How to Extract Data from an Integral

  • Context: Undergrad 
  • Thread starter Thread starter ecastro
  • Start date Start date
  • Tags Tags
    Data Integral
Click For Summary

Discussion Overview

The discussion revolves around the extraction of a function ##B(x)## from an integral equation of the form ##A = \int B(x) C(x) dx##, where participants explore whether this can be done numerically or analytically given known values for ##A## and ##C(x)##. The conversation touches on theoretical aspects, potential methods, and specific applications, particularly in image reconstruction.

Discussion Character

  • Exploratory
  • Technical explanation
  • Debate/contested
  • Mathematical reasoning
  • Experimental/applied

Main Points Raised

  • Some participants question whether ##A## is a definite or indefinite integral, which affects the nature of the problem.
  • It is noted that for any function ##B(x)##, there are infinitely many functions ##C(x)## that can yield the same value for ##A##.
  • Participants discuss the analogy to discrete cases, where knowing ##A## and ##C_j## leads to underdetermined systems for ##B_j##.
  • Some suggest that methods exist for solving integral equations, but uniqueness of the solution for ##B(x)## is not guaranteed.
  • One participant raises the possibility of using Singular Value Decomposition (SVD) as a method, questioning its applicability to the problem.
  • There is mention of a specific application in image filtering, where ##C(x)## represents a filter and ##A## is a filtered image, leading to inquiries about reconstructing the unfiltered image ##B(x)##.
  • Concerns are expressed that the mathematical problem cannot be solved uniquely for ##B(x)## without additional information.
  • Another participant introduces Hilbert space concepts, indicating that the problem can be framed in terms of scalar products, suggesting that multiple solutions exist unless further constraints are applied.
  • Some propose that if ##B(x)## is assumed to be a specific type of function (e.g., polynomial or exponential), it may be possible to extract it with sufficient equations.

Areas of Agreement / Disagreement

Participants generally agree that the problem is underdetermined and that multiple solutions exist. There is no consensus on a specific method for extracting ##B(x)##, and various approaches are debated without resolution.

Contextual Notes

Limitations include the dependence on the nature of the integral (definite vs. indefinite), the assumptions about the forms of ##B(x)## and ##C(x)##, and the lack of unique solutions without additional constraints or information.

ecastro
Messages
249
Reaction score
8
Consider this form:

##A = \int B\left(x\right) C\left(x\right) dx##

I have the values for ##A## and ##C\left(x\right)## (a value of ##C## per value of ##x##), is there a way that I can extract ##B\left(x\right)## numerically or analytically?

Thank you in advance.
 
Physics news on Phys.org
A question, is ##A## a number or a function ##A=A(x)##, in other words the integration is definite on an interval or is indefinite?
 
Given any B(x) there are an infinite number of different functions, C(x), that will give the same A.
 
ecastro said:
Consider this form:

##A = \int B\left(x\right) C\left(x\right) dx##

I have the values for ##A## and ##C\left(x\right)## (a value of ##C## per value of ##x##), is there a way that I can extract ##B\left(x\right)## numerically or analytically?

Thank you in advance.

Yes, equations of this for are called integral equations. There are numerous methods that have been developed to solve integral equations. However, it's not always possible to find a unique solution.
 
I will assume that the integral in the original post is a definite integral, or else the questions won't make sense.

This is analogous to the discrete case, where you have number Bj and Cj, for 1 ≤ j ≤ n and a number A such that

A = ∑ Bj Cj

where the sum is over j = 1, 2, 3, ..., n.

Consider the simplest non-trivial case, where n = 2:

A = B1 C1 + B2 C2,​

where you know A and the Cj, and are trying to determine the B1 and the B2.

This gives you one linear equation in two unknowns. As we know from algebra, we need two linear equations (and independent ones at that) to determine two unknowns. It does, of course, present a constraint on the B1 and B2 — but that is not enough to solve for them.

This not only gives an idea of why the original integral equation cannot be solved for B(x), but this can even be translated into a corresponding integral equation, with the C(x) function defined as a step function, constant on each of two consecutive unit-length intervals. This discrete case shows that even if B(x) were likewise required to be a step function constant on each of those unit intervals, you still wouldn't be able to solve for it.
 
Ssnow said:
A question, is ##A## a number or a function ##A=A(x)##, in other words the integration is definite on an interval or is indefinite?

The variable ##A## is a number for this case. I apologise if I haven't mentioned that the integral is definite.

the_wolfman said:
Yes, equations of this for are called integral equations. There are numerous methods that have been developed to solve integral equations. However, it's not always possible to find a unique solution.

I'm trying to search these methods with Google, but I always find sites on how to solve a particular integral equation, not to extract a function from it.

zinq said:
I will assume that the integral in the original post is a definite integral, or else the questions won't make sense.

This is analogous to the discrete case, where you have number Bj and Cj, for 1 ≤ j ≤ n and a number A such that

A = ∑ Bj Cj

where the sum is over j = 1, 2, 3, ..., n.

Consider the simplest non-trivial case, where n = 2:

A = B1 C1 + B2 C2,​

where you know A and the Cj, and are trying to determine the B1 and the B2.

This gives you one linear equation in two unknowns. As we know from algebra, we need two linear equations (and independent ones at that) to determine two unknowns. It does, of course, present a constraint on the B1 and B2 — but that is not enough to solve for them.

This not only gives an idea of why the original integral equation cannot be solved for B(x), but this can even be translated into a corresponding integral equation, with the C(x) function defined as a step function, constant on each of two consecutive unit-length intervals. This discrete case shows that even if B(x) were likewise required to be a step function constant on each of those unit intervals, you still wouldn't be able to solve for it.

Is Singular Value Decomposition method valid for this case? I just happen to read on this technique and I do not know if it is applicable for this problem.

In point of fact, I have encountered this form of equation in image filters, where ##C\left(x\right)## is the filter and ##x## is the wavelength. I have a filtered image, which is ##A##, and I want to extract the unfiltered image which is ##B\left(x\right)##. Do you have any thoughts on how to solve this kind of problem?
 
I'm no expert on image reconstruction, but I can say that the math problem as you have described it cannot be solved uniquely for the function B(x) — not even close.

Possibly if you have more information about B(x) and C(x), something more intelligent can be said.
 
zinq said:
I'm no expert on image reconstruction, but I can say that the math problem as you have described it cannot be solved uniquely for the function B(x) — not even close.
Yes, and "not even close" is an understatement.
If the equation you are trying to solve is $$A=\int_{a}^{b} B(x)C(x) dx \ \ (1)$$
any function ##D## satisfying ##\displaystyle \int_{a}^{b} D(x)C(x) dx \neq 0## is a solution up to a constant, as $$B(x)=\frac{A}{\int_{a}^{b} D(x)C(x) dx} D(x)$$ will be a solution to ##(1)##.
 
In Hilbert space terminology, your problem can be written as [itex]A=<B, C>[/itex], where < , > is the scalar product of that particular space (in this case the scalar product is the integral). This means that, for example, [itex]<B+H, C>=<B, C>+<H, C>[/itex]. From this simple fact follows that whenever [itex]<H, C>=0[/itex], [itex]A = <B, C>+<H, C> = <B+H, C>[/itex]. Therefore, there can be no single solution for [itex]B[/itex].

An interesting question might be: find the "minimal" solution for B. This implies that B is some scalar multiple of C, say λ⋅C. This again means that [itex]A=<\lambda\cdot C, C>=\lambda<C, C>[/itex]. Since you know C(x), you can calculate [itex]<C, C>=\int_{a}^{b}C^{2}(x)dx[/itex] and then finding λ is trivial.
 
  • #10
If you assume it is a type of function, such as a polynomial or exponential, then maybe it is possible with enough equations for the unknowns.
 

Similar threads

  • · Replies 6 ·
Replies
6
Views
3K
  • · Replies 8 ·
Replies
8
Views
4K
  • · Replies 11 ·
Replies
11
Views
2K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 14 ·
Replies
14
Views
3K
  • · Replies 12 ·
Replies
12
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
3K