3
$\begingroup$

Is it possible to write the following function as $H(x)=$'some expresssion` ? $D(x) = H(x) + H(x-1)$

Edit: Hey everyone, thanks for all the great responses, and just to clarify H(x) and D(x) are always going to be polynomials, I wasn't sure if that made too big of a difference so I didn't mention it before.

Edit 2: I'm sincerely sorry if I was too general, I've only really used polynomial equations in my education thus far, so I didn't realize they might not seam just regular to everyone else. Once again I'm very sorry If I didn't give everyone an opportunity to answer the question because of my vagueness. I sincerely appreciate all of your answers and time.

  • 2
    @Ratz: Restricting to polynomials is a very different situation than "just two regular ol' functions". Now you've had everybody chasing wild geese for weeks.2011-11-13

7 Answers 7

1

I'm not sure if I'm interpreting the question correctly but I am reading it as follows.

Suppose that $D:\mathbb{R}\rightarrow\mathbb{R}$ is a polynomial with real coefficients. Does there exist a polynomial $H:\mathbb{R}\rightarrow\mathbb{R}$ such that $D(x)=H(x)+H(x-1)$.

Let $D(x)=\sum_{i=0}^na_ix^i$ ($a_n\neq0$). It's pretty clear that we need $H(x)$ to be of degree $n$ also. Let $H(x)=\sum_{i=0}^nb_ix^i$. Now

$H(x)+H(x-1)=\sum_{i=0}^n(b_ix^i+b_i(x-1)^i)$ $=\sum_{i=0}^n\left(b_ix^i+b_i\left(\sum_{j=0}^i{i\choose j}x^{i-j}(-1)^j\right)\right)$

Now we can ask the question what is the coefficient of $x^m$? Now the indices $i$, $j$ satisfy $0\leq j\leq i\leq n$. Hence for the $x^{i-j}$, we need $i-j=m$. This can only happen for $(i,j)\in{(m,0),(m+1,1),\dots,(m+(n-m),n-m)}$. Introduce the index $l=0,1,\dots,n-m$. Using this we can write the coefficient of $x^m$ and hence we solve the system of equations

$a_m=b_m+\sum_{l=0}^{m-n}b_{m+l}{m+l\choose l}(-1)^{m+1}$ for $m=0,1,\dots,n$.

These are equations of the form

$a_m=\sum_{k=m}^n\alpha^m_{k}b_k$. $(*)_m$

Hence the equation $a_n=\alpha^m_nb_n$ determines $b_n$ and we can solve sequentially for $b_{n-1},b_{n-2},\dots,b_0$ so that when we solve the $(*)_m$ the $b_k$ for $k>m$ have already been found. This implies that $H$ is in fact unique.

I hope this is what you were looking for and sorry it's a bit messy --- it's 2 a.m. in Ireland!

4

Yes, although maybe not the prettiest one: $ D(x) - D(x - 1) + \cdots = H(x) + H(x - 1) - (H(x - 1) + H(x - 2)) + \cdots $ Using summation notation, we have $ \sum_{n=0}^\infty (-1)^nD(x-n) = \sum_{n=0}^\infty (-1)^n(H(x-n) + H(x - n - 1)) = H(x) $ And there you go.

  • 5
    You're making some assumption about convergence, no?2011-10-17
3

$H(x)$ could be defined as anything on $[0,1)$, then the relation $ D(x) = H(x) + H(x-1)\tag{1} $ would define $H(x)$ on the rest of $\mathbb{R}$. For example, for $x\ge0$, $ H(x)=(-1)^{\lfloor x\rfloor}H(x-\lfloor x\rfloor)+\sum_{k=0}^{\lfloor x\rfloor-1}(-1)^kD(x-k)\tag{2} $ and for $x<0$, $ H(x)=(-1)^{\lfloor x\rfloor}H(x-\lfloor x\rfloor)+\sum_{k=\lfloor x\rfloor}^{-1}(-1)^kD(x-k)\tag{3} $ Polynomial solutions:

Some interest was shown for polynomial solutions when $D(x)$ is a polynomial. As was done in another recent answer, we can use Taylor's Theorem to show that $ e^{-\mathcal{D}}P(x)=P(x-1)\tag{4} $ where $\mathcal{D}=\dfrac{\mathrm{d}}{\mathrm{d}x}$, at least for polynomial $P$. Using $(4)$, $(1)$ becomes $ D(x)=\left(I+e^{-\mathcal{D}}\right)H(x)\tag{5} $ Noting that $ (1+e^{-x})^{-1}=\tfrac{1}{2}+\tfrac{1}{4}x-\tfrac{1}{48}x^3+\tfrac{1}{480}x^5+\dots\tag{6} $ we can get a polynomial solution of $(1)$ with $ H(x)=\left(\tfrac{1}{2}I+\tfrac{1}{4}\mathcal{D}-\tfrac{1}{48}\mathcal{D}^3+\tfrac{1}{480}\mathcal{D}^5+\dots\right)D(x)\tag{7} $ Polynomial example:

For example, if $D(x)=x^2$, then using $(7)$, $H(x)=\tfrac{1}{2}x^2+\tfrac{1}{2}x$. Check: $\left(\tfrac{1}{2}x^2+\tfrac{1}{2}x\right)+\left(\tfrac{1}{2}(x-1)^2+\tfrac{1}{2}(x-1)\right)=x^2$

  • 0
    @Ratz: The polynomial solutions are unique and are given explicitly above. For higher order polynomials, you have to compute more terms of the Taylor expansion for $(1+e^{-x})^{-1}$ and plug them into $(7)$. I had added the polynomial solutions before you had narrowed the scope of your question to polynomials, but since the answer is so nice for polynomials, I added them, too.2011-11-13
2

I'll use lower case, so $d(x) = h(x) + h(x-1)$.

Let $D(x) = \sum \limits_{n=0}^{\infty} d(n) x^n/n!$

and $H(x) = \sum \limits_{n=0}^{\infty} h(n) x^n/n!$.

H'(x) = \sum \limits_{n=1}^{\infty} h(n) x^{n-1}/(n-1)! = \sum \limits_{n=0}^{\infty} h(n+1) x^n/n! so H(x) + H'(x) = \sum \limits_{n=0}^{\infty} (h(n)+h(n+1)) x^n/n! = \sum \limits_{n=0}^{\infty} d(n+1) x^n/n!.

Since D'(x) = \sum \limits_{n=0}^{\infty} d(n+1) x^n/n!, H(x) + H'(x) = D'(x).

The usual integrating factor is $e^x$: e^x D'(x) = e^x(H(x) + H'(x)) =(e^x H(x))', so e^x H(x) = \int e^x D'(x) dx or H(x) = e^{-x} \int e^x D'(x) dx . If you can do the integral (and I haven't made too many mistakes) , you can get $H$.

You can integrate by parts, but you still have to deal with $\int e^x D(x) dx$.

Now I'll see what was entered while I was entering this ...

  • 0
    Hey this looks promising but it didn't work for the original case I intended. If you look at the simple situation H(x) = x then D(x) = 2x-1 and D'(x)= 2. But then plug into the equation H(x) = e^-x*(|(e^x*2)dx) and the answer comes out to be 2.2011-11-02
2

For $x$ an integer, up to a constant, $H(x)$ is the alternating sum of $D(x)$ from 0 up to $x$ (no reason to start at $-\infty$ like Arthur, although this might sometimes be preferable) since

$((-1)^x D(0)- (-1)^x D(1)+\dots + D(x))+ ((-1)^{x-1}D(0) - (-1)^{x-1}D(1) + \dots + D(x-1)) = D(x) $

and for $H(x)=0$ the equation gives a recurrence that fixes $D$ given $D(0)$.

Therefore, your problem includes the calculations of all finite sums and you cannot expect a reasonable answer without additional assumptions on your functions.

For example, if your $D$ is a polynomial, then you know that the answer is also a polynomial that can easily be found by an interpolation formula. Clearly, the polynomial equation will then hold for all $x$. For a hypergeometric expression, the Gosper algorithm can be adapted.

For an integrable function, the Euler-Maclaurin formula can sometimes be of use.

But there is no general answer. Sums do not always have simple answers and even if they do, there is no simple way to find them.

2

I will assume $x$ is an integer, and the functions are $\mathbb N \to \mathbb R$.

If you have a function $H$ such that $D(x)=H(x)+H(x-1)$ then you can add to $H$ the function $c (-1)^x$ for any constant $c$, and it will still be a solution. For example, while $0=0+0$, also $0=(-1)^x+(-1)^{x+1}$.

So the solution is not uniquely determined. How do you expect me to guess which $H$ did you mean when $D=0$? I will instead find all $H$ satisfying the equation.

Suppose

$D(x)=H(x)+H(x-1)$

this means

$H(x) = D(x) - H(x-1)$

This is a recurrence relation. If you know $H(0)$ you can compute all values of $H$ by going back, for example:

$H(2) = D(2) - H(1) = D(2) - (D(1) - H(0)) = D(2) - D(1) + H(0)$

and generally, the signs alternate:

$H(x) = \sum \limits_{i=1}^{x} (-1)^{x-i} D(i) + (-1)^x H(0)$

(easy to write formula for negative numbers too)

General formula:

$H(x) = \sum \limits_{i=1}^{x} (-1)^{x-i} D(i) + (-1)^x c$

where $c$ is any constant. The equation will be satisfied for any choice of $c$. (For more about that stuff, read about affine spaces or nonhomogeneous equations.)

  • 0
    @Ratz: So when I wrote half an hour earlier that $H$ is the alternating sum, it was not what you were looking for?2011-11-13
1

You mentioned an example in one of your comments: if $H(x)=x$, then $D(x)=2x-1$. So let's suppose $D(x)=2x-1$ (for, say, $x=0,1,2,\dots$), and try to recover $H(x)$.

We have $1=D(1)=H(0)+H(1)$. What is $H(0)$? Well, it could be anything. You could take $H(0)=17$; then for $H(x)$ you would get the sequence $17,-16,19,-14,21,-12,23,-10,\dots$ given by the formula $H(x)=x+17(-1)^x$. The general solution is not just $H(x)=x$ but $H(x)=x+H(0)(-1)^x$ In short, you want a single function $H(x)$, but what we are all telling you in various ways is that for any given $D(x)$ there are infinitely many different $H(x)$.