Tài liệu Đề tài " Radon inversion on Grassmannians via G˚ardingGindikin fractional integrals " pptx

LINK DOWNLOAD MIỄN PHÍ TÀI LIỆU "Tài liệu Đề tài " Radon inversion on Grassmannians via G˚ardingGindikin fractional integrals " pptx": http://123doc.vn/document/1036260-tai-lieu-de-tai-radon-inversion-on-grassmannians-via-g-ardinggindikin-fractional-integrals-pptx.htm

786 ERIC L. GRINBERG AND BORIS RUBIN
k × k matrix. Each function f on the Grassmannian G
n,k
can be identified
with the relevant function f(x)onV
n,k
which is O(k) right-invariant, i.e.,
f(xγ)=f(x) ∀γ ∈ O(k) (the group of orthogonal k × k matrices). The
right O(k) invariance of a function on the Stiefel manifold simply means that
the function is invariant under change of basis within the span of a given
frame, and hence “drops” to a well-defined function on the Grassmannian.
The aforementioned identification enables us to reach numerous important
statements and to achieve better understanding of the matter by working with
functions of a matrix argument.
Definition 1.1. Given η ∈ G
n,k
and y ∈ V
n,
,≤ k, we define
(1.7) Cos
2
(η, y)=y

Pr
η
y, Sin
2
(η, y)=y

Pr
η
⊥
y,
where η
⊥
denotes the (n − k)-subspace orthogonal to η.
Both quantities represent positive semidefinite  ×  matrices. This can
be readily seen if we replace the linear operator Pr
η
by its matrix xx

where
x =[x
1
, ,x
k
] ∈ V
n,k
is an orthonormal basis of η. Clearly,
Cos
2
(η, y) + Sin
2
(η, y)=I

.
We introduce the following mean value operators
(1.8)
(M
r
f)(ξ)=

Cos
2
(ξ,x)=r
f(x)dm
ξ
(x), (M
∗
r
ϕ)(x)=

Cos
2
(ξ,x)=r
ϕ(ξ)dm
x
(ξ),
x ∈ V
n,k
,ξ∈ G
n,k

,r∈P
k
; dm
ξ
(x) and dm
x
(ξ) are the relevant induced
measures. A precise definition of these integrals is given in Section 3. According
to this definition, (M
∗
r
ϕ)(x) is well defined as a function of η ∈ G
n,k
, and (up
to abuse of notation) one can write (M
∗
r
ϕ)(x) ≡ (M
∗
r
ϕ)(η). Operators (1.8)
are matrix generalizations of the relevant Helgason transforms for k = 1 (cf.
formula (35) in [H2, p. 96]). The mean value M
∗
r
ϕ with the matrix-valued
averaging parameter r ∈P
k
serves as a substitute for M
∗
v
ϕ in (1.4). For
r = I
k
, operators (1.8) coincide with the Radon transform (1.1) and its dual,
respectively (see §4).
Theorem 1.2. Let f ∈ L
p
(G
n,k
), 1 ≤ p<∞. Suppose that ϕ(ξ)=
(Rf)(ξ),ξ∈ G
n,k

, 1 ≤ k<k

≤ n − 1,k+ k

≤ n, and denote
(1.9) α =(k

− k)/2, ˆϕ
η
(r) = (det(r))
α−1/2
(M
∗
r
ϕ)(η),c=
Γ
k
(k/2)
Γ
k
(k

/2)
.
Then for any integer m>(k

− 1)/2,
(1.10) f(η)=c
(L
p
)
lim
r→I
k
(D
m
+
I
m−α
+
ˆϕ
η
)(r),
RADON INVERSION ON GRASSMANNIANS
787
the differentiation being understood in the sense of distributions. In particular,
for k

− k =2,  ∈ N,
(1.11) f = c
(L
p
)
lim
r→I
k
(D

+
ˆϕ
η
)(r).
If f is a continuous function on G
n,k
, then the limit in (1.10) and (1.11) can
be treated in the sup-norm.
This theorem gives a family of inversion formulae parametrized by the
integer m. They generalize (1.4) to the higher rank case and f ∈ L
p
. The
equality (1.10) coincides with (1.4), if k =1,m = k

, and has the same
structure. Moreover, (1.10) covers the full range (1.3), including even and odd
cases for k

− k. A simple structure of the formula (1.10) is based on the fact
that the analytic family {I
α
+
} includes the identity operator, namely, I
0
+
= I.
Here one should take into account that I
α
+
w for Re α ≤ (k − 1)/2 is defined by
analytic continuation (for sufficiently good w) or in the sense of distributions;
see Section 2.2 and [Gi].
As in the classical Funk-Radon theory, Theorem 1.2 is preceded by a
similar one for zonal functions. The results for this important special case are
as follows.
Definition 1.3 (-zonal functions). Let O(n) be the group of orthogonal
n × n matrices. Fix  so that 1 ≤  ≤ n − 1. Given ρ ∈ O(n − ), let
g
ρ
=

ρ 0
0 I


∈ O(n).
A function f(η)onG
n,k
is called -zonal if f(g
ρ
η)=f(η) for all ρ ∈ O(n − ).
If  = k = 1 then an -zonal function depends only on one variable,
sometimes called height.
In the following theorems we employ the notion of rank of a symmetric
space. This can be defined in various equivalent ways, e.g., using Lie algebras,
maximal totally geodesic flat subspaces or invariant differential operators [H3].
The rank of G
n,k
can be computed: rank G
n,k
= min (k, n−k). Rank comes up
in the harmonic analysis of functions on Grassmannians, and the injectivity
dimension criterion (1.3) can be motivated by means of rank considerations
[Gr3]. Here we do not use the intrinsic definition of rank explicitly, but it
surfaces autonomously in the analysis.
Theorem 1.4. Choose  so that 1 ≤  ≤ min (k, n − k)(= rank G
n,k
),
and let f(η) be an integrable -zonal function on G
n,k
.
(i) There is a function f
0
(s) on P

so that
f(η)
a.e.
= f
0
(s),s= Cos
2
(η, σ

),σ

=

0
I


∈ V
n,
,
788 ERIC L. GRINBERG AND BORIS RUBIN
and
(1.12)

G
n,k
f(η)dη =
Γ

(n/2)
Γ

(k/2) Γ

((n − k)/2)
I


0
f
0
(s)dµ(s),
(1.13) dµ(s) = (det(s))
(k−−1)/2
(det(I

− s))
(n−k−−1)/2
ds.
(ii) If  ≤ k

− k, 1 ≤ k<k

≤ n − 1, then the Radon transform
(Rf)(ξ),ξ∈ G
n,k

, is represented by the G˚arding-Gindikin fractional integral
as follows:
(1.14) (Rf)(ξ)=c (det(S))
−(k

−−1)/2
(I
α
+
˜
f
0
)(S),
where
˜
f
0
(s) = (det(s))
(k−−1)/2
f
0
(s),
α =(k

− k)/2,S= Cos
2
(ξ,σ

) ∈P

,c=Γ

(k

/2)/Γ

(k/2).
Let us comment on this theorem. The identity (1.12) gives precise infor-
mation about the weighted L
1
space to which f
0
(s) belongs. This information
is needed to keep convergence of numerous integrals which arise in the analysis
below under control. The condition 1 ≤  ≤ rank G
n,k
is natural. It reflects
the geometric fact that G
n,k
is isomorphic to G
n,n−k
and is necessary to en-
sure absolute convergence of the integral in the right-hand side of (1.12). The
additional condition  ≤ k

− k in (ii) is necessary for absolute convergence of
the fractional integral in (1.14), but it is not needed for (Rf)(ξ) because the
latter exists pointwise almost everywhere for any integrable f. This obvious
gap can be reduced if we restrict ourselves to the case when (Rf )(ξ), as well
as f, is a function on the cone P

. To this end we impose the extra condition
1 ≤  ≤ rank G
n,k

and get
(1.15) 1 ≤  ≤ min(rank G
n,k
, rank G
n,k

) = min (k,n − k

).
This condition does not imply  ≤ k

− k. Hence we need a substitute for
(1.14) which holds for  satisfying (1.15) and enables us to invert Rf.
Theorem 1.5. Let  satisfy 1 ≤  ≤ min(k,n − k

), and suppose that
ϕ(ξ)=(Rf)(ξ),ξ∈ G
n,k

, where f(η) is an integrable -zonal function on
G
n,k
.
(i) There exist functions f
0
(s) and F
0
(S) so that
f(η)
a.e.
= f
0
(s),s= Cos
2
(η, σ

),ϕ(ξ)
a.e.
= F
0
(S),S= Cos
2
(ξ,σ

).
If
ˆ
f
0
(s) = (det(s))
(k−−1)/2
f
0
(s) and
ˆ
F
0
(S) = (det(S))
(k

−−1)/2
F
0
(S) then
(1.16) I
(n−k

)/2
+
ˆ
F
0
= cI
(n−k)/2
+
ˆ
f
0
,c=Γ

(k

/2)/Γ

(k/2).
RADON INVERSION ON GRASSMANNIANS
789
(ii) The function f
0
(s) can be recovered by the formula
(1.17) f
0
(s)=c
−1
(det(s))
−(k−−1)/2
(D
m
+
I
m−α
+
ˆ
F
0
)(s),
α =(k

− k)/2,m∈ N,m>(k

− 1)/2,
where D
m
+
is understood in the sense of distributions.
Natural analogs of Theorems 1.4 and 1.5 hold for the dual Radon trans-
form. For k = 1, these results were obtained in [Ru2]. Unlike the case
k = 1 (where pointwise differentiation is possible), we cannot do the same
for k>1. The treatment of D
m
+
in the sense of distributions is unavoidable
in the framework of the method (even for smooth f), because of convergence
restrictions. The latter are intimately connected with the complicated struc-
ture of the boundary of P
k
(or P

). It is important to note that in the -zonal
case inversion formulae for the Radon transform and its dual hold without the
assumption k + k

≤ n.
A few words about technical tools are in order. We were inspired by
the papers of Herz [Herz] and Petrov [P2] (unfortunately the latter was not
translated into English). The key role in our argument belongs to Lemma 2.2
which extends the notion of bispherical coordinates [VK, pp. 12, 22] to Stiefel
manifolds and generalizes Lemma 3.7 from [Herz, p. 495].
The paper is organized as follows. Section 2 contains preliminaries and
derivation of basic integral formulae. In the rank-one case these formulae are
known to every analyst working on the sphere. We need their extension to
Stiefel and Grassmann manifolds. In Section 2 we also prove part (i) of The-
orem 1.4 (see Corollary 2.9). In Section 3 we introduce mean value operators,
which can be regarded as matrix analogs of geodesic spherical means on S
n−1
,
and which play a key role in our consideration. In Section 4 we complete the
proof of the main theorems. Theorem 4.6 covers part (ii) of Theorem 1.4, and
a similar statement holds for the dual Radon transform R
∗
. Theorem 4.10 im-
plies (1.16) and the corresponding equality for R
∗
. Inversion formulae (1.10),
(1.11), (1.17), and an analog of (1.17) for R
∗
are proved at the end of the
section.
Acknowledgements. The work was started in Summer 2000 when B.
Rubin was visiting Temple University in Philadelphia. He expresses gratitude
to his co-author, Professor Eric Grinberg, for the hospitality. Both authors
are grateful to the referee for his comments and valuable suggestions owing to
which the original text of the paper was essentially improved.
2. Preliminaries
2.1. Notation, matrix spaces and Siegel gamma functions. The main
references for the following are [M, Ch. 2 and Appendix], [T, Ch. 4], [Herz]. We
790 ERIC L. GRINBERG AND BORIS RUBIN
recall some basic facts and definitions. Let
M
n,k
be the space of real matrices
having n rows and k columns. One can identify
M
n,k
with the real Euclidean
space R
nk
so that for x =(x
i,j
) the volume element is dx =

n
i=1

k
j=1
dx
i,j
.
In the following x

denotes the transpose of x, 0 (sometimes with subscripts)
denotes zero entries; I
k
is the identity k ×k matrix; e
1
, ,e
n
are the canonical
coordinate unit vectors in R
n
.
Let S
k
be the space of k × k real symmetric matrices r =(r
i,j
),r
i,j
= r
j,i
.
A matrix r ∈ S
k
is called positive definite (positive semidefinite) if a

ra > 0
(a

ra ≥ 0) for all vectors a =0inR
k
; this is commonly expressed as r>0
( r ≥ 0). Given r
1
,r
2
∈ S
k
, the inequality r
1
>r
2
means r
1
− r
2
∈P
k
. The
following facts are well known; see, e.g., [M], [T]:
(i) If r>0 then r
−1
> 0.
(ii) For any matrix x ∈
M
n,k
,x

x ≥ 0.
(iii) If r ≥ 0 then r is nonsingular if and only if r>0.
(iv) If r>0,s>0,r− s>0 then s
−1
− r
−1
> 0 and det(r) > det(s).
(v) A symmetric matrix is positive definite (positive semidefinite) if and
only if all its eigenvalues are positive (nonnegative).
(vi) If r ∈ S
k
then there exists an orthogonal matrix γ ∈ O(k) such that
γ

rγ = diag(λ
1
, ,λ
k
) where each λ
j
is real and equal to the j
th
eigenvalue
of r.
(vii) If r is a positive semidefinite k × k matrix then there exists a positive
semidefinite k × k matrix, written as r
1/2
, such that r = r
1/2
r
1/2
.
We hope that, with these properties in mind, the reader will find more
transparent the numerous calculations with functions of a matrix variables that
occur throughout the paper.
The set S
k
of symmetric k × k matrices is a vector space of dimension
k(k +1)/2 and is a measure space isomorphic to R
k(k+1)/2
with the volume
element dr =

i≤j
dr
i,j
.Forr ≥ 0 we shall use the notation |r| = det(r).
Given positive semidefinite matrices r and R in S
k
, the symbol

R
r
f(s)ds
denotes integration over the set
{s : s ∈P
k
,r<s<R}.
For Ω ⊂P
k
, the function space L
p
(Ω) is defined in the usual way with respect
to the measure dr. The set P
k
is a convex cone in S
k
. It is a symmetric
space of the group GL(k, R) of non-singular k × k real matrices. The action
of g ∈ GL(k, R)onr ∈P
k
is given by r → g

rg. This action is transitive (but
not simply transitive). The relevant invariant measure on P
k
has the form
(2.1) dµ(r)=|r|
−d

1≤i≤j≤k
dr
i,j
,d=(k +1)/2,
RADON INVERSION ON GRASSMANNIANS
791
[T, p. 18]. Let T
k
be the group of upper triangular matrices t of the form
(2.2) t =






t
1
.t
i,j
.
.
0 t
k






,t
i
> 0,t
i,j
∈ R.
Each r ∈P
k
has a unique representation r = t

t, t ∈ T
k
, so that
(2.3)

P
k
f(r)dr =
∞

0
t
k
1
dt
1
∞

0
t
k−1
2
dt
2

∞

0
t
k
˜
f(t
1
, ,t
k
) dt
k
,
˜
f(t
1
, ,t
k
)=2
k
∞

−∞

∞

−∞
f(t

t)

i<j
dt
i,j
[T, p. 22], [M, p. 592]. In this last integration the diagonal entries of the matrix
t are given by the arguments of
˜
f, and the strictly upper triangular entries of
t are variables of integration.
To the cone P
k
one can associate the Siegel gamma function
(2.4) Γ
k
(α)=

P
k
e
−tr(r)
|r|
α−d
dr, tr(r) = trace of r.
By (2.3), it is easy to check [M, p. 62] that this integral converges absolutely
for Re α>d− 1, and represents the product of the usual Γ-functions:
(2.5) Γ
k
(α)=π
k(k−1)/4
Γ(α)Γ(α −
1
2
) Γ(α −
k − 1
2
).
For the corresponding Beta function we have [Herz, p. 480]
(2.6)
R

0
|r|
α−d
|R − r|
β−d
dr = B
k
(α, β)|R|
α+β−d
,
B
k
(α, β)=
Γ
k
(α)Γ
k
(β)
Γ
k
(α + β)
;Reα, Re β>d− 1; R ∈P
k
.
2.2. G˚arding-Gindikin fractional integrals. Let
Q = {r ∈P
k
:0<r<I
k
}
be the “unit interval” in P
k
. Let f be a function in L
1
(Q). The G˚arding-
Gindikin fractional integrals of f of order α are defined by
(2.7) (I
α
+
f)(r)=
1
Γ
k
(α)
r

0
f(s)|r − s|
α−d
ds,
792 ERIC L. GRINBERG AND BORIS RUBIN
(2.8) (I
α
−
f)(r)=
1
Γ
k
(α)
I
k

r
f(s)|s − r|
α−d
ds,
where r ∈ Q, d =(k +1)/2, Re α>d− 1. Both integrals are finite for
almost all r ∈ Q. To see this it suffices to show that the integrals

I
k
0
(I
α
±
f)(r)dr
are finite for any nonnegative f ∈ L
1
(Q). By changing the order of integration,
and evaluating inner integrals according to (2.6), we get
I
k

0
(I
α
+
f)(r)dr = c
I
k

0
f(s)|I
k
− s|
α
ds,
I
k

0
(I
α
−
f)(r)dr = c
I
k

0
f(s)|s|
α
ds,
c =Γ
k
(d)/Γ
k
(α + d). Since the right-hand sides of these equalities are ma-
jorized by const ||f||
L
1
(Q)
, the statement follows.
The equality (2.6) also implies the semigroup property
(2.9) I
α
±
I
β
±
f = I
α+β
±
f, f ∈ L
1
(Q), Re α, Re β>d− 1.
For s =(s
i,j
) ∈P
k
, we define the following differential operators in the
s-variable:
D
+
= det

η
i,j
∂
∂s
i,j

,η
i,j
=

1ifi = j
1/2ifi = j,
D
−
=(−1)
k
2
D
+
.(2.10)
If f is sufficiently good, then
(2.11) D
m
±
I
α
±
f = I
α−m
±
f, m ∈ N, Re α>m+ d − 1,
(see, e.g., [G˚a]). Let D(Q) be the space of infinitely differentiable functions
supported in Q.Forw ∈D(Q), the integrals I
α
±
w can be extended to all α ∈ C
as entire functions of α, so that I
0
±
w = w, I
α
±
I
β
±
w = I
α+β
±
w and D
m
±
I
α
±
w =
I
α
±
D
m
±
w = I
α−m
±
w for all α, β ∈ C and all m ∈ N [Gi]. This enables us to
define I
α
±
f for f ∈ L
1
(Q) and Re α ≤ d − 1 in the sense of distributions by
setting
(I
α
±
f,w)=

Q
(I
α
±
f)(r)w(r)dr =(f,I
α
∓
w),w∈D(Q).
Note that explicit construction of the analytic continuation of I
α
±
w is rather
complicated if w does not vanish identically on the boundary of Q (cf. [G˚a],
[R1], [R2]). In order to invert ϕ = I
α
+
f for f ∈ L
1
(Q) and Re α>d− 1 in the
sense of distributions, let m ∈ N,m− Re α>d− 1. By (2.9), I
m
+
f = I
m−α
+
ϕ,
and therefore
(f,w) ≡

Q
f(r)w(r)dr =(D
m
+
I
m−α
+
ϕ, w) ≡ (ϕ, I
m−α
−
D
m
−
w).
2.3. Stiefel manifolds. Let V
n,k
= {x ∈ M
n,k
: x

x = I
k
} be the Stiefel
manifold of orthonormal k-frames in R
n
,n ≥ k.Forn = k, V
n,n
=O(n)
RADON INVERSION ON GRASSMANNIANS
793
represents the orthogonal group in R
n
. The Stiefel manifold is a homogeneous
space with respect to the action V
n,k
 x → γx ∈ V
n,k
,γ∈ O(n), so that
V
n,k
=O(n)/O(n − k). The group O(n) acts on V
n,k
transitively. The same
is true for the group SO(n)={γ ∈ O(n) : det (γ)=1} provided n>k.Itis
known that dim V
n,k
= k(2n − k − 1)/2. We fix invariant measures dx on V
n,k
and dγ on SO(n) normalized by
(2.12) σ
n,k
≡

V
n,k
dx =
2
k
π
nk/2
Γ
k
(n/2)
and

SO(n)
dγ = 1 [M, p. 70], [J, p. 57].
Lemma 2.1 (polar decomposition). Almost all x ∈
M
n,k
,n≥ k (specif-
ically, all matrices x ∈
M
n,k
of rank k), can be decomposed as
x = vr
1/2
,v∈ V
n,k
,r= x

x ∈P
k
so that dx =2
−k
|r|
(n−k−1)/2
drdv.
This statement can be found in [Herz, p. 482], [GK, p. 93], [M, pp. 66,
591].
Lemma 2.2 (bi-Stiefel decomposition). Let k and  be arbitrary integers
satisfying 1 ≤ k ≤  ≤ n − 1,k+  ≤ n. Almost all x ∈ V
n,k
can be represented
in the form
(2.13) x =

ur
1/2
v(I
k
− r)
1/2

,u∈ V
,k
,v∈ V
n−,k
,r∈P
k
,
so that
(2.14)

V
n,k
f(x)dx =
I
k

0
dν(r)

V
,k
du

V
n−,k
f

ur
1/2
v(I
k
− r)
1/2

dv,
(2.15) dν(r)=2
−k
|r|
γ
|I
k
− r|
δ
dr, γ =
 − k − 1
2
,δ=
n −  − k − 1
2
.
Proof.Fork = 1, this statement is well known and represents bispherical
decomposition on the unit sphere; cf. [VK, pp. 12, 22]. For the general case
related to Stiefel manifolds the proof is essentially the same as that of the
slightly less general Lemma 3.7 from [Herz, p. 495]. For convenience of the
reader we sketch this proof.
Let us check (2.13). If x =

a
b

∈ V
n,k
,a∈ M
,k
,b∈ M
n−,k
, then
I
k
= x

x = a

a + b

b. By Lemma 2.1 for almost all a we have a = ur
1/2
. Hence
b

b = I
k
− r, and therefore b = v(I
k
− r)
1/2
. This gives (2.13). The explicit
meaning of “almost all” in Lemma 2.2 becomes clear from Lemma 2.1 having
been applied to the matrices a and b.
794 ERIC L. GRINBERG AND BORIS RUBIN
In order to prove (2.14) we write it in the form
(2.16)

V
n,k
f(x)dx =

0<a

a<I
k
|I
k
− a

a|
δ
da

V
n−,k
f

a
v(I
k
− a

a)
1/2

dv
and show the coincidence of the two measures, dx and
˜
dx = |I
k
− a

a|
δ
dadv.
Following [Herz], we consider the Fourier transforms
F
1
(s)=

V
n,k
etr(is

x)dx and F
2
(s)=

V
n,k
etr(is

x)
˜
dx,
where s ∈
M
n,k
, etr(Λ) = e
tr
(Λ)
, and show that F
1
= F
2
. To this end we
employ the Bessel functions A
λ
(r) of Herz for which
(2.17)

V
n,k
etr(is

x)dx =2
k
π
nk/2
A
(n−k−1)/2

1
4
s

s

.
Let
s =

s
1
s
2

,x=

a
v(I
k
− a

a)
1/2

; s
1
∈
M
,k
,s
2
∈
M
n−,k
; a ∈
M
,k
.
Then s

x = s

1
a + s

2
v(I
k
− a

a)
1/2
, and we have
F
2
(s)=

a

a<I
k
etr(is

1
a)|I
k
− a

a|
δ
da

V
n−,k
etr(is

2
v(I
k
− a

a)
1/2
)dv.
By (2.17) (use the equality tr(is

2
vR) = tr(iR
−1
Rs

2
vR) = tr(iRs

2
v) with
R =(I
k
− a

a)
1/2
) the inner integral is evaluated as
2
k
π
(n−)k/2
A
δ

1
4
Rs

2
s
2
R

=2
k
π
(n−)k/2
A
δ

1
4
s

2
s
2
R
2

(the last equality holds because of the invariance property A
δ
(R
−1
rR)=
A
δ
(r)). Thus
F
2
(s)=

a

a<I
k
etr(is

1
a)ϕ(a

a)da,
ϕ(r)=2
k
π
(n−)k/2
|I
k
− a

a|
δ
A
δ

1
4
s

2
s
2
(I
k
− r)

.
The function F
2
(s) can be transformed by the generalized Bochner formula

M
,k
etr(iy

a)ϕ(a

a)da = π
k/2
g

1
4
y

y

,
g(Λ) =

P
k
A
γ
(Λr)|r|
γ
ϕ(r)dr, γ =
 − k − 1
2
,y∈
M
,k
RADON INVERSION ON GRASSMANNIANS
795
[Herz, p. 493], that yields
F
2
(s)=2
k
π
nk/2
I
k

0
A
γ

1
4
s

1
s
1
r

|r|
γ
A
δ

1
4
s

2
s
2
(I
k
− r)

|I
k
− r|
δ
dr,
γ, δ being defined by (2.15). This integral can be evaluated using the formula
(2.6) from [Herz, p. 487]. The result is
F
2
(s)=2
k
π
nk/2
A
(n−k−1)/2

1
4
(s

1
s
1
+ s

2
s
2
)

=2
k
π
nk/2
A
(n−k−1)/2

1
4
s

s

.
By (2.17), the latter coincides with F
1
(s).
Remark 2.3. The assumptions k +  ≤ n and k ≤  in Lemma 2.2 are
necessary for absolute convergence of the integral

I
k
0
in the right-hand side of
(2.14). It would be interesting to prove this lemma directly, without using the
Fourier transform. Such a proof would be helpful in transferring Lemma 2.2
and many other results of the paper to the hyperbolic space (cf. [VK, pp. 12,
23], [BR], [Ru5] for the rank-one case).
Lemma 2.4. Let x ∈ V
n,k
,y∈ V
n,
;1≤ k, ≤ n.Iff is a function of
 × k matrices then
(2.18)
1
σ
n,k

V
n,k
f(y

x)dx =
1
σ
n,

V
n,
f(y

x)dy.
Proof. We should observe that formally the left-hand side is a function of
y, while the right-hand side is a function of x. In fact, both are constant. To
prove (2.18) let G = SO(n),g∈ G, g
1
= g
−1
. The left-hand side is

G
f(y

gx)dg =

G
f((g
1
y)

x)dg
1
which equals the right-hand side.
We shall need a “lower-dimensional” representation of integrals of the form
(2.19) I
f
=

V
n,k
f(A

x)dx, A ∈ M
n,
;0<k<n, 0 <<n.
For k =  = 1 such a representation is well known. In the following lemma we
do not specify assumptions for the function f. For our purposes it suffices to
assume only that the integral (2.19) is absolutely convergent and therefore well
defined for all or almost all A. This enables us to give a proof which consists, in
fact, of a number of applications of the Fubini theorem. Furthermore, for our
purposes it suffices to consider matrices A for which A

A is positive definite.