Re: Linear Interpolation Inversion for N even?
- From: "Roger Stafford" <ellieandrogerxyzzy@xxxxxxxxxxxxxxxxxxxxxx>
- Date: Wed, 7 May 2008 21:04:03 +0000 (UTC)
Greg Heath <heath@xxxxxxxxxxxxxxxx> wrote in message <8c13e682-
e924-4aea-8d41-d5d78bfd4d50@xxxxxxxxxxxxxxxxxxxxxxxxxxx>...
On May 6, 11:54 pm, "Roger Stafford"work for any N.
<ellieandrogerxy...@xxxxxxxxxxxxxxxxxxxxxx> wrote:
Greg Heath <he...@xxxxxxxxxxxxxxxx> wrote in message
<6ca9feef-9014-4491-
b1d2-60acf4d08...@xxxxxxxxxxxxxxxxxxxxxxxxxxxx>...> This doesn't
(2),...,xmGreg
-------
Greg, my understanding of your problem is that from xm(1),xm
suffice to(n-1), and some s, you want to recover the original x's, where
xm(1) = (x(1)+x(2))/2
xm(2) = (x(2)+x(3))/2
...
xm(n-1) = (x(n-1)+x(n))/2
I proposed s = x(1)-x(2)+x(3)-x(4)+x(5)....
but you don't know x you only know xm.
I admit you'll have to do a little work, but such an s would actually
Let nrecover the original x's. This can be shown by the following example.
way= 4. Define matrix A by
A = [.5 .5 0 0;
0 .5 .5 0;
0 0 .5 .5;
1 -1 1 -1];
b = [xm;s];
but s(x) is unknown.
where xm is a column vector. The equation
A*x = b;
corresponds to the original equalities above, and solving it for the x's by
formulaof x = A\b would recover their original values uniquely. The general
neverfor the determinant of A when its size is n x n is (-.5)^(n-1)*n, so A can
alsobe singular and these linear equations always have a unique solution.
Of course there is an infinitude of other ways of defining s that would
evenwork, but I thought this one fit in with the spirit of your original sum,
thethough that one didn't actually work for even values of n. As you see,
-----------------above s works for all n, even or odd.
I see now that
2*sum(j=1,N-1){ xm(j)} = 2*sum(j=1,N){ x(j) } - ( x(1) + x(N) )
= - x(1) - x(N)% for mean(x) = 0
and
2*sum(j=1,N-1){(-1)^(j-1) * xm(j) } = x(1) + (-1)^(N-2) * x(N)
= x(1) - x(N)% for N odd
= x(1) + x(N)% for N even
Which leads to the simultaneous solution of x(1) and x(N) when
N is odd and an undetermined solution when N is even.
Hope this helps.
Greg
Greg, I confess I'm feeling a bit exasperated at this point. Originally, you
asked the question, and I quote, "Is there a way to recover x if no original
points are known and N is even?" (in reference to only knowing the midpoint
values xm.) The answer is a flat NO for both N even or N odd. There is no
way the x's can be recovered without some additional information.
You also stated, "Under certain conditions the sampled function x(1:N) can
be reconstructed from the linearly intepolated midpoints." That circumstance
was that N be odd and the sum(x) be known, and this is true. It is not true if
N is even.
I proposed that instead of knowing sum(x), suppose you knew
s = sum((-1).^(0:N-1).*x)
which is just x(1)-x(2)+x(3)-x(4)+... In that case the x's can be recovered if s
is known, no matter whether N is even or odd.
Your response was to complain that you didn't know s because you didn't
know x. Why in the odd N case doesn't it also bother you that you don't know
sum(x) because you don't know x? What is there so special about sum(x), as
opposed to the sum of the x's with alternating sign reversals?
You will never be able to reconstruct the x values knowing only the xm
values. That is a fact. You must know some one other additional bit of
information about them. Answering a hypothesis about such additional
knowledge with the statement that "but you don't know x you only know xm"
is not going to bring you any closer to a resolution to your problem! You
must make up your mind whether or not you are going to permit anything
else to be known about these x's other than their midpoints, xm. If you
aren't, there is no solution to the problem. If you are, it isn't reasonable to
object to such additional knowledge simply because "but you don't know x."
In your last paragraph you appear to be tentatively trying out the hypothesis
of knowing the alternate-signs sum, but your attempt to do the
reconstruction is entirely incorrect. That is not the way it is done. It must be
done as I described in the earlier article, solving the equation A*x = b which I
defined there.
Roger Stafford
.
- Follow-Ups:
- Re: Linear Interpolation Inversion for N even?
- From: Greg Heath
- Re: Linear Interpolation Inversion for N even?
- References:
- Re: Linear Interpolation Inversion for N even?
- From: Greg Heath
- Re: Linear Interpolation Inversion for N even?
- Prev by Date: Re: gammainc function
- Next by Date: GMSK Modulator (.m code)
- Previous by thread: Re: Linear Interpolation Inversion for N even?
- Next by thread: Re: Linear Interpolation Inversion for N even?
- Index(es):
Relevant Pages
|
Loading