# Inversion of Fourier Transforms by Means of Scale-Frequency Series

Authors: Nassar H. S. Haidar
Journal: International Journal of Engineering Mathematics
 (2014)

---

## Abstract

We report on inversion of the Fourier transform when the frequency variable can be scaled in a variety of different ways that improve the resolution of certain parts of the frequency domain. The corresponding inverse Fourier transform is shown to exist in the form of two dual scale-frequency series. Upon discretization of the continuous scale factor, this Fourier transform series inverse becomes a certain nonharmonic double series, a discretized scale-frequency (DSF) series. The DSF series is also demonstrated, theoretically and practically, to be rate-optimizable with respect to its two free parameters, when it satisfies, as an entropy maximizer, a pertaining recursive nonlinear programming problem incorporating the entropy-based uncertainty principle.


---

## Body

## 1. Introduction

We revisit the classical problem of recovering a Fourier transformable signal(1)
f
(
t
)
=
1
2
π
∫
R
‍
f
^
(
μ
)
e
i
μ
t
d
μ
=
F
-
1
{
f
^
(
μ
)
}
in F
(
R
), R
=
(
-
∞
,
∞
), from its Fourier transform; see, for example, [1, 2], image f
^
(
μ
). The space F is the set of all real functions f for which the integral f
^
(
μ
)
=
∫
R
‍
f
(
t
)
e
-
i
μ
t
d
t exists.Obviously, the function to be integrated in (1) is the product
(2)
f
^
(
μ
)
e
i
μ
t
=
f
^
(
μ
)
[
cos
⁡
μ
t
+
i
sin
μ
t
]
,
t
∈
R
,
and when the variable |
t
|
→
∞, the function e
i
μ
t will oscillate extremely rapidly. Then in order to follow the variations of the product f
^
(
μ
)
e
i
μ
t meaningfully in a μ-quadrature formula, see, for example, [3], for (1), there is a need for a large number of μ-points, even for slowly varying f
^
(
μ
). Consequently, the faster the f
^
(
μ
) decreases as |
μ
|
→
∞, the more the tractable computationally the integral (1) will be, and the more accurate are its numerical evaluations.Bandlimiting of signals has been a practical way for easing the previous problem. A fact that has so far been motivating the wide engineering interest in bandlimited signals, for which(3)
F
{
f
(
t
)
}
=
{
g
^
(
μ
)
,
|
μ
|
≤
2
π
B
0
,
otherwise
=
f
^
(
μ
)
,
where B is the band width of f
(
t
) with an image f
^
(
μ
) having the [
-
2
π
B
,
2
π
B
] compact support of g
^
(
μ
). Inversion of bandlimited signals is currently performed by projections onto convex sets (POCS) algorithms [4], which encompass the Gerchberg-Papoulis algorithm [5].In this work we report on a novel robust semianalytical method for integration of (1) to recover a signal f
(
t
)
∈
F
(
R
) that is not necessarily bandlimited. In this method, the frequency variable can be scaled in a variety of different ways that improve the resolution of certain parts of the frequency domain. The corresponding inverse Fourier transform is shown to exist in the form of two dual scale-frequency series. Upon discretization of the continuous scale factor, this Fourier transform series inverse becomes a certain nonharmonic [6] double series, a discretized scale-frequency (DSF) series. This DSF series is demonstrated, for the first time in this work, to be the proper framework for an entropy-maximizing inverse Fourier transformation. A fact has emerged after theoretically illustrating that the DSF series is rate-optimizable with respect to the scale factor step and the frequency gap size. The robustness of this optimization (which is illustrated by a simple example) reduces to the satisfaction, via entropy-maximization, of an emerging recursive nonlinear programming problem, incorporating the entropy-based uncertainty principle.


## 2. Theory of the Continuous Scale-Frequency Series Inverse

LetZ be the set of integers, Q the set of rational numbers, and Q
c the set of irrational numbers. λ
∈
R
∖
{
0
} will be a scale factor, while the frequency γ can either be m
∈
Z or ν
=
ν
(
m
)
=
(
m
+
ç
)
∈
Q or Q
c (depending on whether ç
∈
Q or Q
c
). The symbol □ means either multiplication or division of γ by λ, and γ
c is the alternative (complementary) symbol to γ; that is, if γ
=
m, then γ
c
=
ν and vise versa. The signal f
(
t
)
∈
F
(
R
) is assumed to be expandable in the variably scaled trigonometric basis:
(4)
q
□
γ
λ
(
t
)
=
{
λ
-
1
[
cos
⁡
(
γ
c
□
λ
)
t
+
sin
(
γ
□
λ
)
t
]
}
.The ensemble ofγ with □ in (4) can be in any of their possible C
2
4
=
6 combinations of 4 quantities, taken two at a time. Accordingly throughout this work we will assume, as in [7], any of
(5)
q
≻
m
λ
(
t
)
=
λ
-
1
[
cos
⁡
ν
λ
t
+
sin
m
λ
t
]
;
□
=
≻
:
γ
□
λ
=
m
λ
,
γ
c
□
λ
=
ν
λ
,
q
≺
m
λ
(
t
)
=
λ
-
1
[
cos
⁡
m
λ
t
+
sin
ν
λ
t
]
;
□
=
≺
:
γ
□
λ
=
ν
λ
,
γ
c
□
λ
=
m
λ
for high frequency resolving scaling (as λ
→
0
), any of
(6)
q
⊢
m
λ
(
t
)
=
λ
-
1
[
cos
⁡
ν
λ
t
+
sin
m
λ
t
]
;
□
=
⊢
:
γ
□
λ
=
m
λ
,
γ
c
□
λ
=
ν
λ
,
q
⊣
m
λ
(
t
)
=
λ
-
1
[
cos
⁡
m
λ
t
+
sin
ν
λ
t
]
;
□
=
⊣
:
γ
□
λ
=
ν
λ
,
γ
c
□
λ
=
m
λ
for low frequency resolving scaling, while any of
(7)
q
⊳
m
λ
(
t
)
=
λ
-
1
[
cos
⁡
m
λ
t
+
sin
m
λ
t
]
;
□
=
⊳
:
γ
□
λ
=
m
λ
,
γ
c
□
λ
=
m
λ
,
q
⊲
m
λ
(
t
)
=
λ
-
1
[
cos
⁡
m
λ
t
+
sin
m
λ
t
]
;
□
=
⊲
:
γ
□
λ
=
m
λ
,
γ
c
□
λ
=
m
λ
for mixed (high and low) frequency resolving scaling.Next consider the sequence (q
□
γ
λ) that supports the bases (4) when R
o
=
(
R
∖
{
0
}
), and let □ be any one of the six symbols ≻, ≺, ⊢, ⊣, ⊳, ⊲, to invoke the following double parameterized Riemann-Lebesgue lemma (RLL).Lemma 1 (see [7]).
ForΓ
□
γ
λ
[
·
]
=
∫
R
‍
[
·
]
q
□
γ
λ
(
t
)
d
t, when m
∈
Z
,
(8)
lim
⁡
|
m
|
→
∞
Γ
□
γ
λ
[
f
]
=
0
,
for all f
∈
F
(
R
), for any scaling factor λ
∈
R
o, and for all six □ symbols.This lemma justifies the existence and guarantees the convergence of the scale-frequency transform(9)
F
□
(
λ
,
m
)
:
=
∫
R
‍
f
(
t
)
q
□
γ
λ
(
t
)
d
t
=
〈
f
,
q
□
γ
λ
(
t
)
〉
,
for all λ
∈
R
o, of f
∈
F
(
R
) with respect to q
□
γ
λ. Make then use of the notation
(10)
{
F
c
(
μ
)
F
s
(
μ
)
}
=
∫
R
‍
f
(
t
)
{
cos
⁡
μ
t
sin
μ
t
}
d
t
in
(11)
f
^
(
μ
)
=
∫
R
‍
f
(
t
)
e
-
i
μ
t
d
t
=
F
{
f
(
t
)
}
to rewrite it as
(12)
f
^
(
μ
)
=
F
c
(
μ
)
-
i
F
s
(
μ
)
.
That is,
(13)
F
c
(
μ
)
=
Re
[
f
^
(
μ
)
]
,
F
s
(
μ
)
=
-
Im
⁡
[
f
^
(
μ
)
]
.Unlike the Fourier transformationf
^
(
μ
), which is a complex function of μ
,
F
□
(
λ
,
m
) is a real function of the scaled frequencies
(14)
ρ
□
=
γ
□
λ
,
μ
□
=
γ
c
□
λ
.Theorem 2.
F
□
-
1
(
λ
,
m
)
=
f
(
t
)
∈
F
(
R
) satisfies
(15)
f
(
t
)
=
1
2
π
∫
R
‍
F
c
(
μ
□
)
cos
⁡
μ
□
t
d
μ
□
λ
+
1
2
π
∫
R
‍
F
s
(
ρ
□
)
sin
ρ
□
t
d
ρ
□
λ
.Proof.
The above result follows directly from the structural linearity of(16)
F
□
(
λ
,
m
)
=
1
λ
F
c
(
μ
□
)
+
1
λ
F
s
(
ρ
□
)
and the linearity of its inverse transformation.The direct inversion formula (15) can be logarithmic in the scaling factor and is computationally of much larger complexity than (1). A simpler indirect inversion is, however, possible, and in a dual process, that is tied to the possible utilization of anyone of the components of f
^
(
μ
) in (13); namely,
(17)
F
□
(
λ
,
m
)
⟺
{
F
□
c
-
1
(
λ
,
m
)
F
□
s
-
1
(
λ
,
m
)
}
=
F
□
-
1
(
λ
,
m
)
=
f
(
t
)
,
when
{
F
c
(
μ
)
F
s
(
μ
)
}
is
utilized
.
Interestingly, this results in 12 distinct scale-frequency forms for the inverse Fourier transform f
(
t
)
.Obviously, the symbols⊢, ⊣ correspond to scaling sines and cosines in the same way, while ≻, ≺ or ⊳, ⊲ correspond to scaling cosines and sines in different ways. It should be noted here that although a Fourier transform is virtually unaffected by a shift of f in time (which mixes the roles of the cosines and sines), the corresponding scale-frequency transform could change drastically. This raises the question on the rather unclear merits of different-way scalings, to the favour of same-way (
⊢
,
⊣
) scalings. It is believed, however, that different-way scalings can particularly reveal their merits when applied to discontinuous f
∈
F
(
R
)
.Theorem 3 (high frequency resolving scaling).
When□
=
≻ in F
□
(
λ
,
m
)
⇔
f
(
t
), the inverse is either
(18)
f
(
t
)
=
F
≻
c
-
1
(
λ
,
m
)
=
1
2
π
M
∑
m
=
0
M
-
1
‍
∫
R
‍
m
F
c
(
m
λ
)
cos
⁡
(
λ
+
1
)
m
t
d
λ
-
1
2
π
M
∑
m
=
0
M
-
1
‍
∫
R
‍
m
{
F
c
(
m
λ
)
+
F
s
(
[
m
+
ç
]
λ
)
-
F
c
(
[
m

+

ç
]
λ
)
}

×

sin
(
λ
+
1
)
m
t
d
λ
,
or the dual representation
(19)
f
(
t
)
=
F
≻
s
-
1
(
λ
,
m
)
=
1
2
π
M
×
∑
m
=
0
M
-
1
‍
∫
R
‍
(
m
+
ç
)
{
F
c
(
m
λ
)
+
F
s
(
[
m
+
ç
]
λ
)
-
F
s
(
m
λ
)
}
×
cos
⁡
(
[
m
+
ç
]
λ
+
m
)
t
d
λ
-
1
2
π
M
×
∑
m
=
0
M
-
1
‍
∫
R
‍
(
m
+
ç
)
F
s
(
[
m
+
ç
]
λ
)
sin
(
[
m
+
ç
]
λ
+
m
)
t
d
λ
.Proof.
By incorporation of(20)
F
≻
(
λ
,
m
)
=
1
λ
F
c
(
μ
≻
)
+
1
λ
F
s
(
ρ
≻
)
in the details of the proofs of Theorem 1 and Corollary 1 of [7].A similar theorem holds also for the case when□
=
≺.Theorem 4 (low frequency resolving scaling).
For□
=
≻ in F
□
(
λ
,
m
)
⇔
f
(
t
), the inverse is either
(21)
f
(
t
)
=
F
⊢
c
-
1
(
λ
,
m
)
=
-
1
2
π
M

×

∑
m
=
0
M
-
1
‍
∫
R
‍
m
1
λ
F
c
(
m
λ
)
cos
⁡
(
1
λ

+

1
)
m
t
d
(
ln
⁡
λ
)
+
1
2
π
M
×
∑
m
=
0
M
-
1
‍
∫
R
‍
m
{
1
λ
F
c
(
m
λ
)
+
1
λ
F
s
(
[
m
+
ç
]
λ
)
-
1
λ
F
c
(
[
m
+
ç
]
λ
)
}
×
sin
(
1
λ
+
1
)
m
t
d
(
ln
⁡
λ
)
,
or the dual representation
(22)
f
(
t
)
=
F
⊢
s
-
1
(
λ
,
m
)
=
1
2
π
M
×
∑
m
=
0
M
-
1
‍
∫
R
‍
(
m
+
ç
)
{
1
λ
F
c
(
m
λ
)
+
1
λ
F
s
(
[
m
+
ç
]
λ
)
-
1
λ
F
s
(
m
λ
)
(
[
m
+
ç
]
λ
)
}
×
cos
⁡
(
[
m
+
ç
]
λ
+
m
)
t
d
(
ln
⁡
λ
)
+
1
2
π
M
∑
m
=
0
M
-
1
‍
∫
R
‍
(
m
+
ç
)
1
λ
F
s
(
[
m
+
ç
]
λ
)
×
sin
(
[
m
+
ç
]
λ
+
m
)
t
d
(
ln
⁡
λ
)
when the logarithmic integrals are conceived in the Principal Value sense.Proof.
By incorporation ofF
⊢
(
λ
,
m
)
=
(
1
/
λ
)
F
c
(
μ
⊢
)
+
(
1
/
λ
)
F
s
(
ρ
⊢
) in the details of the proofs of Theorem 2 and Corollary 2 of [7].A similar theorem holds also for the case when□
=
⊣.


## 3. Discretized Scale-Frequency Series Inverse

Here we discretize theλ axis uniformly over R
o according to λ
k
=
k
δ
∈
[
-
Λ
,
Λ
]
⊂
R, where Λ is the scale factor band width, which is not apriori known. k
∈
[
-
N
,
N
]
∖
{
0
}
=
W
o
⊂
Z, with δ
=
Λ
/
N as the scale sample size, which can be arbitrarily small or large. In the special case when Λ
=
N, δ
=
1
.


### 3.1. Discretized High Frequency Resolving Scaling

The above mentioned dicretization of the scale variable, withk as a scaling factor number, when □
=
≻, leads to
(23)
F
≻
(
k
δ
,
m
)
=
1
k
δ
F
c
(
m
k
δ
)
+
1
k
δ
F
s
(
[
m
+
ç
]
k
δ
)
,
the pertaining phase
(24)
φ
k
,
m
;
δ
=
m
k
δ
+
m
,
and to the directly related discretized notation,
(25)
A
k
,
m
=
δ
2
π
M
m
Re
[
f
^
(
m
k
δ
)
]
,
(26)
B
k
,
m
=
δ
2
π
M
m
[
Im
⁡
[
f
^
(
[
m
+
ç
]
k
δ
)
]
+
Re
[
f
^
(
[
m
+
ç
]
k
δ
)
]
-
Re
[
f
^
(
m
k
δ
)
]
]
.Alternatively, for the dual representation (19), we employ
(27)
φ
k
,
m
;
δ
+
=
(
m
+
ç
)
k
δ
+
m
,
(28)
A
k
,
m
+
=
δ
2
π
M
(
m
+
ç
)
[
Re
[
f
^
(
m
k
δ
)
]
-
Im
⁡
[
f
^
(
[
m
+
ç
]
k
δ
)
]
+
Im
⁡
[
f
^
(
m
k
δ
)
]
]
,
(29)
B
k
,
m
+
=
δ
2
π
M
(
m
+
ç
)
Im
⁡
[
f
^
(
[
m
+
ç
]
k
δ
)
]
to state one of the main results of this work.Theorem 5 (high frequency resolving discrete scaling).
When□
=
≻ in F
□
(
λ
,
m
)
⇔
f
(
t
), the DSF series inverse is either
(30)
f
(
t
)
≈
D
F
≻
c
-
1
(
λ
,
m
)
=
∑
k
∈
W
o
‍
∑
m
=
0
M
-
1
‍
A
k
,
m
cos
⁡
φ
k
,
m
;
δ
t
+
B
k
,
m
sin
φ
k
,
m
;
δ
t
,
or the dual representation
(31)
f
(
t
)
≈
D
F
≻
s
-
1
(
λ
,
m
)
=
∑
k
∈
W
o
‍
∑
m
=
0
M
-
1
‍
A
k
,
m
+
cos
⁡
φ
k
,
m
;
δ
+
t
+
B
k
,
m
+
sin
φ
k
,
m
;
δ
+
t
.Proof.
By considering (13) when following the same arguments of the proofs of Theorem 3 and Corollary 3 of [7].A similar result holds also for the case when□
=
≺.


### 3.2. Discretized Low Frequency Resolving Scaling

Here□
=
⊢ is associated, by discretizing the scale factor, with
(32)
ψ
k
,
m
;
δ
=
m
k
δ
+
m
,
(33)
A
k
,
m
=
-
δ
2
π
M
m
(
1
k
δ
)
2
Re
[
f
^
(
m
k
δ
)
]
,
(34)
B
k
,
m
=
δ
2
π
M
m
(
1
k
δ
)
2
×
{
Re
[
f
^
(
m
k
δ
)
]
-
Im
⁡
[
f
^
(
[
m
+
ç
]
k
δ
)
]
-
Re
[
f
^
(
[
m
+
ç
]
k
δ
)
]
}
.Alternatively, for the dual representation (22), we employ
(35)
ψ
k
,
m
;
δ
+
=
(
m
+
ç
)
k
δ
+
m
,
(36)
A
k
,
m
+
=
δ
2
π
M
(
m
+
ç
)
(
1
k
δ
)
2
×
{
Re
[
f
^
(
m
k
δ
)
]
-
Im
⁡
[
f
^
(
[
m
+
ç
]
k
δ
)
]
+
Im
⁡
[
f
^
(
m
k
δ
)
]
[
f
^
(
[
m
+
ç
]
k
δ
)
]
}
,
(37)
B
k
,
m
+
=
-
δ
2
π
M
(
m
+
ç
)
(
1
k
δ
)
2
Im
⁡
[
f
^
(
[
m
+
ç
]
k
δ
)
]
to state the following main result of this work.Theorem 6 (low frequency resolving discrete scaling).
When□
=
⊢ in F
□
(
λ
,
m
)
⇔
f
(
t
), the DSF series inverse is either
(38)
f
(
t
)
≈
D
F
⊢
c
-
1
(
λ
,
m
)
=
∑
k
∈
W
o
‍
∑
m
=
0
M
-
1
‍
A
k
,
m
cos
⁡
ψ
k
,
m
;
δ
t
+
B
k
,
m
sin
ψ
k
,
m
;
δ
t
or the dual representation
(39)
f
(
t
)
≈
D
F
⊢
s
-
1
(
λ
,
m
)
=
∑
k
∈
W
o
‍
∑
m
=
0
M
-
1
A
k
,
m
+
cos
⁡
ψ
k
,
m
;
δ
+
t
+
B
k
,
m
+
sin
ψ
k
,
m
;
δ
+
t
‍
.Proof.
By considering (13) when following the same arguments of the proofs of Theorem 4 and Corollary 4 of [7].A similar result holds also for the case when□
=
⊣.


## 3.1. Discretized High Frequency Resolving Scaling

The above mentioned dicretization of the scale variable, withk as a scaling factor number, when □
=
≻, leads to
(23)
F
≻
(
k
δ
,
m
)
=
1
k
δ
F
c
(
m
k
δ
)
+
1
k
δ
F
s
(
[
m
+
ç
]
k
δ
)
,
the pertaining phase
(24)
φ
k
,
m
;
δ
=
m
k
δ
+
m
,
and to the directly related discretized notation,
(25)
A
k
,
m
=
δ
2
π
M
m
Re
[
f
^
(
m
k
δ
)
]
,
(26)
B
k
,
m
=
δ
2
π
M
m
[
Im
⁡
[
f
^
(
[
m
+
ç
]
k
δ
)
]
+
Re
[
f
^
(
[
m
+
ç
]
k
δ
)
]
-
Re
[
f
^
(
m
k
δ
)
]
]
.Alternatively, for the dual representation (19), we employ
(27)
φ
k
,
m
;
δ
+
=
(
m
+
ç
)
k
δ
+
m
,
(28)
A
k
,
m
+
=
δ
2
π
M
(
m
+
ç
)
[
Re
[
f
^
(
m
k
δ
)
]
-
Im
⁡
[
f
^
(
[
m
+
ç
]
k
δ
)
]
+
Im
⁡
[
f
^
(
m
k
δ
)
]
]
,
(29)
B
k
,
m
+
=
δ
2
π
M
(
m
+
ç
)
Im
⁡
[
f
^
(
[
m
+
ç
]
k
δ
)
]
to state one of the main results of this work.Theorem 5 (high frequency resolving discrete scaling).
When□
=
≻ in F
□
(
λ
,
m
)
⇔
f
(
t
), the DSF series inverse is either
(30)
f
(
t
)
≈
D
F
≻
c
-
1
(
λ
,
m
)
=
∑
k
∈
W
o
‍
∑
m
=
0
M
-
1
‍
A
k
,
m
cos
⁡
φ
k
,
m
;
δ
t
+
B
k
,
m
sin
φ
k
,
m
;
δ
t
,
or the dual representation
(31)
f
(
t
)
≈
D
F
≻
s
-
1
(
λ
,
m
)
=
∑
k
∈
W
o
‍
∑
m
=
0
M
-
1
‍
A
k
,
m
+
cos
⁡
φ
k
,
m
;
δ
+
t
+
B
k
,
m
+
sin
φ
k
,
m
;
δ
+
t
.Proof.
By considering (13) when following the same arguments of the proofs of Theorem 3 and Corollary 3 of [7].A similar result holds also for the case when□
=
≺.


## 3.2. Discretized Low Frequency Resolving Scaling

Here□
=
⊢ is associated, by discretizing the scale factor, with
(32)
ψ
k
,
m
;
δ
=
m
k
δ
+
m
,
(33)
A
k
,
m
=
-
δ
2
π
M
m
(
1
k
δ
)
2
Re
[
f
^
(
m
k
δ
)
]
,
(34)
B
k
,
m
=
δ
2
π
M
m
(
1
k
δ
)
2
×
{
Re
[
f
^
(
m
k
δ
)
]
-
Im
⁡
[
f
^
(
[
m
+
ç
]
k
δ
)
]
-
Re
[
f
^
(
[
m
+
ç
]
k
δ
)
]
}
.Alternatively, for the dual representation (22), we employ
(35)
ψ
k
,
m
;
δ
+
=
(
m
+
ç
)
k
δ
+
m
,
(36)
A
k
,
m
+
=
δ
2
π
M
(
m
+
ç
)
(
1
k
δ
)
2
×
{
Re
[
f
^
(
m
k
δ
)
]
-
Im
⁡
[
f
^
(
[
m
+
ç
]
k
δ
)
]
+
Im
⁡
[
f
^
(
m
k
δ
)
]
[
f
^
(
[
m
+
ç
]
k
δ
)
]
}
,
(37)
B
k
,
m
+
=
-
δ
2
π
M
(
m
+
ç
)
(
1
k
δ
)
2
Im
⁡
[
f
^
(
[
m
+
ç
]
k
δ
)
]
to state the following main result of this work.Theorem 6 (low frequency resolving discrete scaling).
When□
=
⊢ in F
□
(
λ
,
m
)
⇔
f
(
t
), the DSF series inverse is either
(38)
f
(
t
)
≈
D
F
⊢
c
-
1
(
λ
,
m
)
=
∑
k
∈
W
o
‍
∑
m
=
0
M
-
1
‍
A
k
,
m
cos
⁡
ψ
k
,
m
;
δ
t
+
B
k
,
m
sin
ψ
k
,
m
;
δ
t
or the dual representation
(39)
f
(
t
)
≈
D
F
⊢
s
-
1
(
λ
,
m
)
=
∑
k
∈
W
o
‍
∑
m
=
0
M
-
1
A
k
,
m
+
cos
⁡
ψ
k
,
m
;
δ
+
t
+
B
k
,
m
+
sin
ψ
k
,
m
;
δ
+
t
‍
.Proof.
By considering (13) when following the same arguments of the proofs of Theorem 4 and Corollary 4 of [7].A similar result holds also for the case when□
=
⊣.


## 4. Demonstrative Application

To illustrate the applicability of the previously reported results, we give here three illustrative examples of the high frequency resolving scaling in recovering three well-known simple Fourier preimages.Example 7.
The even real signalf
(
t
)
=
e
-
|
t
| has a real image, f
^
(
μ
)
=
2
/
(
1
+
μ
2
); that is, Im
⁡
[
f
^
(
μ
)
]
=
0. According to (30) it admits the following family of nonharmonic series expansions:
(40)
e
-
|
t
|
≈
D
F
≻
c
-
1
(
λ
,
m
)
=
δ
π
M
∑
k
=
-
N
k
≠
0
N
‍
∑
m
=
0
M
-
1
‍
m
1
+
m
2
k
2
δ
2
cos
⁡
(
k
δ
+
1
)
m
t
+
{
m
1
+
[
m
+
ç
]
2
k
2
δ
2
-
m
1
+
m
2
k
2
δ
2
}
sin
(
k
δ
+
1
)
m
t
,
for various values of ç, δ
,
M, and N. It should be noted here that δ
N defines the unknown scale factor band width Λ
.
Obviously whenç
=
0, the expansions become purely cosinusoidal, and the dual DSF expansions D
F
≻
s
-
1
(
λ
,
m
) can similarly be constructed according to (31).Example 8.
For the odd signalg
(
t
)
=
t
e
-
|
t
| we have a purely imaginary image, g
^
(
μ
)
=
(
-
4
μ
/
(
1
+
μ
2
)
2
)
i; that is, Re
[
g
^
(
μ
)
]
=
0. Consequently, according to (30) it admits the following family of nonharmonic sinusoidal series expansions:
(41)
t
e
-
|
t
|
≈
D
G
≻
c
-
1
(
λ
,
m
)
=
-
2
δ
π
M
×
∑
k
=
-
N
k
≠
0
N
‍
∑
m
=
0
M
-
1
‍
m
[
m
+
ç
]
k
δ
{
1
+
[
m
+
ç
]
2
k
2
δ
2
}
2
sin
(
k
δ
+
1
)
m
t
,
in addition to the dual expansions D
G
≻
s
-
1
(
λ
,
m
)
.Example 9.
The nonsymmetric real signalh
(
t
)
=
(
t
+
1
)
e
-
|
t
| has the complex image h
^
(
μ
)
=
2
/
(
1
+
μ
2
)
-
(
4
μ
/
(
1
+
μ
2
)
2
)
i. Then it admits, by (30), the family of nonharmonic series expansions:
(42)
(
t
+
1
)
e
-
|
t
|
≈
D
H
≻
c
-
1
(
λ
,
m
)
=
δ
π
M
∑
k
=
-
N
k
≠
0
N
‍
∑
m
=
0
M
-
1
‍
m
1
+
m
2
k
2
δ
2
cos
⁡
(
k
δ
+
1
)
m
t
-
{
m
1
+
m
2
k
2
δ
2
+
2
m
[
m
+
ç
]
k
δ
{
1
+
[
m
+
ç
]
2
k
2
δ
2
}
2
-
m
1
+
[
m
+
ç
]
2
k
2
δ
2
m
1
+
m
2
k
2
δ
2
+
2
m
[
m
+
ç
]
k
δ
{
1
+
[
m
+
ç
]
2
k
2
δ
2
}
2
}
sin
(
k
δ
+
1
)
m
t
,
in addition to the dual DSF expansions D
H
≻
s
-
1
(
λ
,
m
).


## 5. Rate-Optimization of the DSF Series Inverse

In DSF series inversion off
^
(
μ
), the number M of frequency points is decided by the need for satisfactory followup of the variation of f
^
(
μ
) over its effective compact support. In contrast, the number of scaling points 2
N
=
2
Λ
/
δ
∈
Z remains rather arbitrary.Besides, these series can be made rate-optimal in the sense that the derivative of a certain norm of them will vanish at a high rate asM and N tend to infinity. Clearly when N
→
∞, δ
→
0 for a given Λ. The two free parameters that can be optimized in this sense are (i) the scale factor step, δ, for a given N, and (ii) the frequency gap size, ç, for a given M.Here we will study this optimization problem for the case of the DSF series inverseD
F
≻
c
-
1
(
λ
,
m
), (30), without loss of generality.Definition 10.
The DSF norm∥
f
∥
□ of f
(
t
)
∈
F
(
R
) is the sum of the amplitudes of all scaled-harmonics appearing in the D
F
≻
c
-
1
(
λ
,
m
) expansion; that is,
(43)
∥
f
∥
□
=
∑
k
=
-
N
k
≠
0
N
‍
∑
m
=
0
M
-
1
‍
A
k
,
m
;
□
2
+
B
k
,
m
;
□
2
,
where A
k
,
m
;
□ and B
k
,
m
;
□ are, respectively, defined by (25)-(26) or (28)-(29) and (33)-(34) or (36)-(37).
The concept of the similar (but different) trigonometric norm∥
f
∥
* was first advanced by Lyaponov in 1954; see, for example, [8], for the standard orthonormal (over [
0,2
π
]) trigonometric system. Our DFS series D
F
≻
c
-
1
(
λ
,
m
) or D
F
≻
s
-
1
(
λ
,
m
) are, however, with respect to the nonharmonic trigonometric systems,
(44)
{
e
i
φ
k
,
m
;
δ
t
}
k
,
m
or
{
e
i
φ
k
,
m
;
δ
+
t
}
k
,
m
;
ç
,
which, for certain δ and ç, may even not be closed in the space L
2
(
R
), and despite the fact that the coefficients A
k
,
m
;
□ and B
k
,
m
;
□ do satisfy the scale-frequency Riemann-Lebesgue lemma [7]. Nevertheless, ∥
f
∥
□ can easily be proved to satisfy the following defining properties of a norm: (i)
∥
c
f
∥
□
=
|
c
|
∥
f
∥
□, c
∈
R (absolute homogeneity),
(ii)
∥
f
+
g
∥
□
≤
∥
f
∥
□
+
∥
g
∥
□ (subadditivity),
(iii)
∥
f
∥
□
=
0, if and only if f
(
t
)
=
0 (separating points).Moreover, as is well known, for piecewise continuous functions, the usual Fourier series (not necessarily (30)) may, in general, not converge absolutely and then the norm ∥
f
∥
* becomes unbounded. Besides, for sufficiently smooth functions, the pertaining Fourier coefficients satisfy the Riemann-Lebesgue lemma and this causes ∥
f
∥
* to approach the uniform norm ∥
f
∥
=
sup
⁡
t
⁡
|
f
(
t
)
| as an upper bound, namely,
(45)
∥
f
∥
≤
∥
f
∥
*
.
Relation (45) is called Lyaponov’s norm inequality.Conjecture 11.
Like∥
f
∥
*, the DSF norm ∥
f
∥
□ also satisfies
(46)
∥
f
∥
≤
∥
f
∥
□
.Proposition 12.
Stationarization of∥
f
∥
□ with respect to δ and ç leads to stabilizing the DSF series inverse.Proof.
Consider, without loss of generality,∥
f
∥
□ for D
F
≻
c
-
1
(
λ
,
m
):
(47)
∥
f
∥
≻
c
=
∑
k
=
-
N
k
≠
0
N
‍
∑
m
=
0
M
-
1
‍
A
k
,
m
2
+
B
k
,
m
2
.
Substitution of (25) and (26) for A
k
,
m
2
+
B
k
,
m
2 in (47), when adopting the simplifying notation
(48)
f
^
δ
=
f
^
(
m
k
δ
)
,
f
^
δ
,
ç
=
f
^
(
[
m
+
ç
]
k
δ
)
,
leads to
(49)
A
k
,
m
2
+
B
k
,
m
2
=
m
2
δ
2
4
π
2
M
2
×
{
〈
Re
[
f
^
δ
,
ç
]
+
Im
⁡
[
f
^
δ
,
ç
]
〉
×
〈
Re
[
f
^
δ
,
ç
]
+
Im
⁡
[
f
^
δ
,
ç
]
-
2
Re
[
f
^
δ
]
〉
+
2
Re
2
[
f
^
δ
]
}
,
then to
(50)
∥
f
∥
≻
c
=
δ
2
π
M
×
∑
k
=
-
N
k
≠
0
N
‍
∑
m
=
0
M
-
1
‍
m
{
〈
Re
[
f
^
δ
,
ç
]
+
Im
⁡
[
f
^
δ
,
ç
]
〉
×
〈
Re
[
f
^
δ
,
ç
]
+
Im
⁡
[
f
^
δ
,
ç
]
-
2
Re
[
f
^
δ
]
〉
+
2
Re
2
[
f
^
δ
]
}
1
/
2
.
Stationarization of∥
f
∥
□ with respect to δ and ç leads, respectively, to
(51)
∂
∂
δ
∥
f
∥
≻
c
=
δ
2
4
π
2
M
2
×
∑
k
=
-
N
k
≠
0
N
‍
∑
m
=
0
M
-
1
‍
m
2
A
k
,
m
2
+
B
k
,
m
2
×
{
∂
δ
〈
Re
[
f
^
δ
,
ç
]
+
Im
⁡
[
f
^
δ
,
ç
]
-
Re
[
f
^
δ
]
〉
×
∂
∂
δ
〈
Re
[
f
^
δ
,
ç
]
+
Im
⁡
[
f
^
δ
,
ç
]
-
Re
[
f
^
δ
]
〉
+
Re
[
f
^
δ
]
∂
∂
δ
Re
[
f
^
δ
]
}
=
0
,
(52)
∂
∂
ç
∥
f
∥
≻
c
=
δ
2
4
π
2
M
2
∑
k
=
-
N
k
≠
0
N
‍
∑
m
=
0
M
-
1
‍
m
2
A
k
,
m
2
+
B
k
,
m
2
×
〈
Re
[
f
^
δ
,
ç
]
+
Im
⁡
[
f
^
δ
,
ç
]
-
Re
[
f
^
δ
]
〉
×
∂
∂
ç
〈
Re
[
f
^
δ
,
ç
]
+
Im
⁡
[
f
^
δ
,
ç
]
〉
=
0
.
Obviously, increasingM and N (or decreasing δ
), which increases the smoothness of D
F
≻
c
-
1
(
λ
,
m
), accelerates the convergence of the series in (51) and (52) to zero. This is accompanied by accelerating the stationarization of (
∂
/
∂
δ
)
∥
f
∥
≻
c and (
∂
/
∂
ç
)
∥
f
∥
≻
c, with a rate ~
M
-
2
N
-
2, and to stabilizing D
F
≻
c
-
1
(
λ
,
m
) with respect to variations in δ and ç. This is a fact that illustrates the consistency of this proposition and its underlying Conjecture 11. The same arguments are also applicable to the rate optimization of every other D
F
□
-
1
(
λ
,
m
)
.Lemma 13.
A rate-optimal(
δ
,
ç
) pair for the DSF series inverse D
F
≻
c
-
1
(
λ
,
m
) exists when the system (51)-(52) is satisfied.Example 14.
To illustrate the previous lemma, let us revisit Example7 with f
^
(
μ
)
=
2
/
(
1
+
μ
2
), to determine the rate-optimal values of δ and ç for it. In the notation of Proposition 12, we have
(53)
Re
[
f
^
δ
]
=
2
1
+
m
2
k
2
δ
2
,
Re
[
f
^
δ
,
ç
]
=
2
1
+
[
m
+
ç
]
2
k
2
δ
2
,
Im
⁡
[
f
^
δ
,
ç
]
=
0
.
First, since(54)
∂
∂
δ
Re
[
f
^
δ
]
=
-
4
m
2
k
2
δ
(
1
+
m
2
k
2
δ
2
)
2
,
∂
∂
δ
Re
[
f
^
δ
,
ç
]
=
-
4
[
m
+
ç
]
2
k
2
δ
(
1
+
[
m
+
ç
]
2
k
2
δ
2
)
2
,
then
(55)
〈
Re
[
f
^
δ
,
ç
]
-
Re
[
f
^
δ
]
〉
=
-
2
ç
[
2
m
+
ç
]
(
1
+
m
2
k
2
δ
2
)
(
1
+
[
m
+
ç
]
2
k
2
δ
2
)
,
∂
∂
δ
〈
Re
[
f
^
δ
,
ç
]
-
Re
[
f
^
δ
]
〉
=
-
4
〈
[
m
+
ç
]
2
(
1
+
[
m
+
ç
]
2
k
2
δ
2
)
2
-
m
2
(
1
+
m
2
k
2
δ
2
)
2
〉
k
2
δ
.
Consider further these relations in (51) to rewrite it as
(56)
∂
∂
δ
∥
f
∥
≻
c
=
2
δ
3
π
2
M
2
∑
k
=
-
N
k
≠
0
N
‍
∑
m
=
0
M
-
1
‍
m
2
k
2
A
k
,
m
2
+
B
k
,
m
2
(
1
+
m
2
k
2
δ
2
)
×
{
〈
[
m
+
ç
]
2
(
1
+
[
m
+
ç
]
2
k
2
δ
2
)
2
-
m
2
(
1
+
m
2
k
2
δ
2
)
2
〉
ç
[
2
m
+
ç
]
(
1
+
[
m
+
ç
]
2
k
2
δ
2
)
×
〈
[
m
+
ç
]
2
(
1
+
[
m
+
ç
]
2
k
2
δ
2
)
2
-
m
2
(
1
+
m
2
k
2
δ
2
)
2
〉
-
m
2
(
1
+
m
2
k
2
δ
2
)
2
〈
[
m
+
ç
]
2
(
1
+
[
m
+
ç
]
2
k
2
δ
2
)
2
-
m
2
(
1
+
m
2
k
2
δ
2
)
2
〉
}
=
0
.
This nonlinear relation inδ can be solved only numerically, for each ç, to yield the rate-optimal δ (i.e., the optimal Λ
) for any given N.
Second, we consider(57)
∂
∂
ç
Re
[
f
^
δ
,
ç
]
=
-
4
k
2
δ
2
[
m
+
ç
]
(
1
+
[
m
+
ç
]
2
k
2
δ
2
)
2
,
∂
∂
ç
Re
[
f
^
δ
]
=
0
,
in (52) to rewrite it as
(58)
∂
∂
ç
∥
f
∥
≻
c
=
2
δ
4
ç
π
2
M
2
×
∑
k
=
-
N
k
≠
0
N
‍
∑
m
=
0
M
-
1
‍
k
2
m
2
A
k
,
m
2
+
B
k
,
m
2
(
1
+
m
2
k
2
δ
2
)
×
[
2
m
+
ç
]
[
m
+
ç
]
(
1
+
[
m
+
ç
]
2
k
2
δ
2
)
3
=
0
.
This can be easily resolved inç and rearranged into the ratio
(59)
ç
=
-
∑
k
=
-
N
k
≠
0
N
‍
∑
m
=
0
M
-
1
‍
[
2
m

+

ç
]
k
2
m
3
A
k
,
m
2

+

B
k
,
m
2
(
1

+

m
2
k
2
δ
2
)
(
1
+
[
m
+
ç
]
2
k
2
δ
2
)
3

×

(
∑
k
=
-
N
k
≠
0
N
‍
∑
m
=
0
M
-
1
‍
[
2
m
+
ç
]
k
2
m
2
A
k
,
m
2
+
B
k
,
m
2
(
1

+

m
2
k
2
δ
2
)
(
1
+
[
m

+

ç
]
2
k
2
δ
2
)
3
)
-
1
for the rate-optimal ç.
Despite its complexity, the previous relation appears as a nonlinearly weighted negative mean ofm, which may always be approximated by some fraction of M, like -
M
/
2, say. Finally, the rate-optimal (
δ
,
ç
) pair for the D
F
≻
c
-
1
(
λ
,
m
) inverse emerges as the solution of the simultaneous system of nonlinear equations (56)-(59), a solution that may not necessarily be unique.


## 6. Recursive Maximum-Entropy Reconstruction of the DSF Series Inverse

The possible nonuniqueness of the rate-optimal(
δ
,
ç
) pair of Lemma 13 calls for imposing an additional weak (inequality) constraint to select a robust value for this pair, which may account for all approximations that are involved in the development of D
F
□
-
1
(
λ
,
m
). The constraint must be the satisfaction of the Heisenberg uncertainty principle (HUP); see, for example, [9], when rate-optimizing δ and ç. The HUP essentially states that f
(
t
)
=
D
F
□
-
1
(
λ
,
m
) and f
^
(
μ
) cannot be both arbitrarily concentrated. In information theory, the definitive measure of concentration of a probability density function, ϕ
(
t
), is the Boltzman entropy
(60)
H
(
ϕ
)
=
-
∫
R
‍
ϕ
(
t
)
ln
⁡
ϕ
(
t
)
d
t
.Hirschman Jr. [10, 11] observed in 1957 that Weyl’s formulation of the HUP is a consequence of a stronger version of the inequality of this principle, namely, the entropy version
(61)
H
(
|
ϕ
|
2
)
H
(
|
ϕ
^
|
2
)
≥
ln
⁡
e
2
,
when ∥
ϕ
∥
2
=
1 and ϕ
^
(
μ
)
≥
0. Hirschman also conjectured that the minimizers of the above sharp inequality were Gaussians, as in the case of the HUP.To state the main result of this section, we adopt the notation(62)
Q
=
∫
R
‍
D
F
≻
c
-
1
(
λ
,
m
)
d
t
,
∥
f
^
∥
2
2
=
∫
R
‍
|
f
^
|
2
d
t
,
the square of the p
=
2 norm, which satisfies Parseval’s identity: ∥
f
^
∥
2
=
∥
f
∥
2.Definition 15.
The DSF Taxi-cab norm1
∥
f
∥
□ of f
(
t
)
∈
F
(
R
) is the sum of the absolute values of all scaled-harmonics appearing in the D
F
≻
c
-
1
(
λ
,
m
) expansion; that is,
(63)
∥
f
∥
□
1
=
∑
k
=
-
N
k
≠
0
N
‍
∑
m
=
0
M
-
1
‍
|
A
k
,
m
;
□
|
+
|
B
k
,
m
;
□
|
,
where A
k
,
m
;
□ and B
k
,
m
;
□ are, respectively, defined by (25)-(26) or (27)-(28) and (33)-(34) or (36)-(37).Lemma 16.
A rate-optimal pair(
δ
,
ç
) for D
F
≻
c
-
1
(
λ
,
m
) with a known Q should satisfy
(64)
1
∥
D
F
≻
c
-
1
(
λ
,
m
)
∥
≻
c
2
[
ln
⁡
∥
D
F
≻
c
-
1
(
λ
,
m
)
∥
≻
c
2
1
-
ln
⁡
Q
2
]
>
Q
4
ln
⁡
⁡
(
e
/
2
)
[
H
(
Q
-
2
|
f
^
|
2
)
-
∥
f
^
∥
2
2
ln
⁡
Q
-
2
]
.Proof.
Conceivef
(
t
)
/
Q as a probability density function ϕ
(
t
). Moreover, if f
(
t
) F
T
↔ f
^
(
μ
), then f
(
t
)
/
Q F
T
↔ f
^
(
μ
)
/
Q. Hence in the notation of our DSF series D
F
≻
c
-
1
(
λ
,
m
), even when it is not necessarily >
0, the previous entropy version for the uncertainty principle writes as
(65)
H
(
Q
-
2
|
D
F
≻
c
-
1
(
λ
,
m
)
|
2
)
≥
ln
⁡
⁡
(
e
/
2
)
H
(
Q
-
2
|
f
^
|
2
)
,
where
(66)
H
(
Q
-
2
|
D
F
≻
c
-
1
(
λ
,
m
)
|
2
)
=
2
∫
R
‍
|
Q
-
2
∑
k
∈
W
o
‍
∑
m
=
0
M
-
1
‍
A
k
,
m
cos
⁡
φ
k
,
m
;
δ
t
+
B
k
,
m
sin
φ
k
,
m
;
δ
t
|
2
×
ln
⁡
|
Q
-
2
∑
k
∈
W
o
‍
∑
m
=
0
M
-
1
‍
A
k
,
m
cos
⁡
φ
k
,
m
;
δ
t
+
B
k
,
m
sin
φ
k
,
m
;
δ
t
|
d
t
,
(67)
H
(
Q
-
2
|
f
^
|
2
)
=
Q
-
2
[
H
(
|
f
^
|
2
)
-
∥
f
^
∥
2
2
ln
⁡
Q
2
]
.
Clearly
(68)
H
(
Q
-
2
|
D
F
≻
c
-
1
(
λ
,
m
)
|
2
)
<
|
Q
-
2
∑
k
∈
W
o
‍
∑
m
=
0
M
-
1
‍
(
A
k
,
m
+
B
k
,
m
)
|
2
×
ln
⁡
|
Q
-
2
∑
k
∈
W
o
‍
(
A
k
,
m
+
B
k
,
m
)
|
2
T
(
k
,
m
;
δ
)
,
where
(69)
T
(
k
,
m
;
δ
)
=
{
2
π
m
(
k
δ
+
1
)
,
m
≠
0
0
,
m
=
0
.
Then replacement of T
(
k
,
m
;
δ
) by 1 in the previous inequality reduces it to
(70)
H
(
Q
-
2
|
D
F
≻
c
-
1
(
λ
,
m
)
|
2
)
<
|
Q
-
2
∑
k
∈
W
o
‍
∑
m
=
0
M
-
1
‍
(
A
k
,
m
+
B
k
,
m
)
|
2
×
ln
⁡
|
Q
-
2
∑
k
∈
W
o
‍
∑
m
=
0
M
-
1
‍
(
A
k
,
m
+
B
k
,
m
)
|
2
.
Now since
(71)
|
Q
-
2
∑
k
∈
W
o
‍
∑
m
=
0
M
-
1
‍
A
k
,
m
+
B
k
,
m
|
2
≤
1
∥
Q
-
1
D
F
≻
c
-
1
(
λ
,
m
)
∥
≻
c
2
=
{
∑
k
=
-
N
k
≠
0
N
‍
∑
m
=
0
M
-
1
‍
|
A
k
,
m
Q
|
+
|
B
k
,
m
Q
|
}
2
,
then
(72)
H
(
Q
-
2
|
D
F
≻
c
-
1
(
λ
,
m
)
|
2
)
<
∥
Q
-
1
D
F
≻
c
-
1
(
λ
,
m
)
∥
≻
c
2
1
ln
⁡
∥
Q
-
1
D
F
≻
c
-
1
(
λ
,
m
)
∥
≻
c
2
1
,
and relation (65) becomes
(73)
∥
Q
-
1
D
F
≻
c
-
1
(
λ
,
m
)
∥
≻
c
2
1
ln
⁡
1
∥
Q
-
1
D
F
≻
c
-
1
(
λ
,
m
)
∥
≻
c
2
>
ln
⁡
⁡
(
e
/
2
)
H
(
Q
-
2
|
f
^
|
2
)
.
Furthermore, since
(74)
∥
Q
-
1
D
F
≻
c
-
1
(
λ
,
m
)
∥
≻
c
2
1
=
Q
-
2
∥
D
F
≻
c
-
1
(
λ
,
m
)
∥
≻
c
2
1
,
then the left-hand side of (72) becomes
(75)
Q
-
2
∥
D
F
≻
c
-
1
(
λ
,
m
)
∥
≻
c
2
1
[
ln
⁡
∥
D
F
≻
c
-
1
(
λ
,
m
)
∥
≻
c
2
1
-
ln
⁡
Q
2
]
.
Finally, consideration of this relation together with (67) in (72) leads to the required result.Since the entropyH
(
Q
-
2
|
D
F
≻
c
-
1
(
λ
,
m
)
|
2
) is a function of the (
δ
,
ç
) free parameters, then variation of these parameters can always be forced towards maximization of this entropy. As is well known, see, for example, [12], a maximum-entropy D
F
≻
c
-
1
(
λ
,
m
) should be a most certain inverse of f
^
(
μ
) and with best achievable smoothness.Theorem 17.
A robust rate-optimal pair(
δ
,
ç
), for D
F
≻
c
-
1
(
λ
,
m
) with a known Q, is a solution to the nonlinear programming problem
(76)
M
i
n
i
m
i
z
e
∥
Q
-
1
D
F
≻
c
-
1
(
λ
,
m
)
∥
≻
c
2
1
ln
⁡
∥
Q
-
1
D
F
≻
c
-
1
(
λ
,
m
)
∥
≻
c
2
1
S
u
b
j
e
c
t
t
o
:
(
64
)
,
(
51
)
-
(
52
)
.Proof.
By consideration of relation (72) of Lemma 16 when maximizing the entropy H
(
Q
-
2
|
D
F
≻
c
-
1
(
λ
,
m
)
|
2
), subject to the inequality constraint (64) and the equality constraints (51) and (52).Q, however, is not apriori known. Hence, we may assume that δ
r
,
ç
r correspond to an rth iteration, Q
r, for Q and D
F
≻
c
-
1
(
λ
,
m
;
δ
r
,
ç
r
) is D
F
≻
c
-
1
(
λ
,
m
) that corresponds to Q
r. These concepts invoke a recursive determination process, indexed by r
=
0,1
,
2,3
,
…
,
L, with L as the termination number, defined, namely, as
(77)
Q
L
+
1
-
Q
L
≤
ϵ
,
where 0
≤
ϵ
≪
1 is an acceptable tolerance.Proposition 18.
A robust rate-optimal pair(
δ
,
ç
)
=
lim
⁡
r
→
L
⁡
(
δ
r
,
ç
r
) for D
F
≻
c
-
1
(
λ
,
m
) exists when (
δ
r
,
ç
r
) is a solution to the recursive nonlinear programming problem
(78)
M
i
n
i
m
i
z
e
∥
Q
r
-
1
D
F
≻
c
-
1
(
λ
,
m
;
δ
r
,
ç
r
)
∥
≻
c
2
1
×
ln
⁡
⁡
∥
Q
r
-
1
D
F
≻
c
-
1
(
λ
,
m
;
δ
r
,
ç
r
)
∥
≻
c
2
1
S
u
b
j
e
c
t
t
o
:
(
64
)
,
(
51
)

-

(
52
)
,
w
i
t
h
D
F
≻
c
-
1

(
λ
,

m
;
δ
r
,
ç
r

)
r
e
p
l
a
c
i
n
g
D
F
≻
c
-
1

(
λ
,

m
)
,
′
Q
r
+
1
=
∫
R
‍
D
F
≻
c
-
1
(
λ
,
m
;
δ
r
,
ç
r
)
d
t
,
and arbitrary Q
0
>
0.Due to convexity of the inequality constraint (64), the nonlinear programming problem of Theorem 17 should be solvable by the method of Kuhn-Tucker multipliers [13]. The rather large computational complexity of the above proposition is comparable to that of the alternative wavelet analysis. It should, however, not represent a difficulty in our present times, when powerful computing resources are available or accessible almost everywhere.


## 7. Conclusion

In this paper, we have demonstrated how evaluation of the inverse Fourier transform can be variably scaled. The reported DSF representation, for, say, high frequency resolving scaling, bypasses the need for quadratures in computing (1). Instead, a repeated evaluation of the real and imaginary parts of f
^, when evaluating the A
k
,
m and B
k
,
m (or the A
k
,
m
+ andB
k
,
m
+) coefficients, of (24)–(26), is needed at a large enough number of scale-frequency points (
m
,
k
), as illustrated in Sections 4 and 5. Low frequency resolving scaling also requires a similar computational effort with the real and imaginary parts of f
^ when evaluating the A
k
,
m and B
k
,
m (or the A
k
,
m
+andB
k
,
m
+) coefficients of (32)–(34).Determination of the rate-optimal scaling factor stepδ and frequency gap number ç, in DSF series has been demonstrated to invoke solving a recursive nonlinear programming problem that maximizes entropy in the framework of the uncertainty principle. It should be emphasized that the scale-frequency approach to inversion of the Fourier transform is a new concept that is not necessariliy fast. In that sense it is not expected to pose a challenge to the FFT algorithm. Its error analysis is, however, of crucial importance and will be a subject for a future publication.



---

*Source: 686785-2014-05-27.xml*