浙江大学学报(理学版), 2023, 50(3): 266-272 doi: 10.3785/j.issn.1008-9497.2023.03.002

数学与计算机科学

傅里叶正、余弦变换的加权卷积及其应用

向仪,, 冯强,,

延安大学 数学与计算机科学学院,陕西 延安 716000

Weighted convolution of the Fourier sine-cosine transform and its application

XIANG Yi,, FENG Qiang,,

School of Mathematics and Computer Science,Yan'an University,Yan'an 716000,Shaanxi Province,China

通讯作者: ORCID: https://orcid.org/0000-0001-6658-9549,E-mail:yadxfq@yau.edu.cn.

收稿日期: 2022-01-10  

基金资助: 国家自然科学基金资助项目.  62261055.  61861044.  62001193.  11961072
陕西省自然科学基金资助项目.  2022JM-400.  2023-JC-YB-085

Received: 2022-01-10  

作者简介 About authors

向仪(1995—),ORCID:https://orcid.org/0000-0003-1072-5163,男,硕士,主要从事积分变换理论与方法研究. 。

摘要

傅里叶变换是求解积分方程常用的工具。基于傅里叶正弦变换与傅里叶余弦变换,定义了两类傅里叶混合加权卷积,得到了傅里叶正弦变换与傅里叶余弦变换的卷积定理,并研究了这两类卷积运算的性质及Young类不等式,将这两类混合加权卷积应用于求解卷积类积分方程,得到了卷积类积分方程的显式解。

关键词: 傅里叶正-余弦变换 ; 卷积定理 ; Young不等式 ; 积分方程

Abstract

Fourier transform is a powerful tool often used to solve some integral equations. In this paper, two kinds of Fourier sine, cosine mixed weighted convolution are defined based on Fourier sine transform and Fourier cosine transform, the corresponding convolution theorems are obtained. The properties of these two kinds of mixed convolution operation and Young type inequality are also studied. And the explicit solutions of convolution type integral equations are obtained based on these two types of mixed weighted convolution.

Keywords: Fourier sine-cosine transform ; convolution theorem ; Young inequality ; integral equation

PDF (465KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

向仪, 冯强. 傅里叶正、余弦变换的加权卷积及其应用. 浙江大学学报(理学版)[J], 2023, 50(3): 266-272 doi:10.3785/j.issn.1008-9497.2023.03.002

XIANG Yi, FENG Qiang. Weighted convolution of the Fourier sine-cosine transform and its application. Journal of Zhejiang University(Science Edition)[J], 2023, 50(3): 266-272 doi:10.3785/j.issn.1008-9497.2023.03.002

0 引 言

傅里叶变换(Fourier transform,FT)作为一种重要的信号处理工具,在应用数学、光学、图像加密、信号处理等领域具有重要应用1-4。在FT基础上定义的傅里叶正弦变换(Fourier sine transform,FST)与傅里叶余弦变换(Fourier cosine transform,FCT),是求解积分方程常用的工具5-6。用FST与FCT处理奇函数和偶函数的计算复杂度是FT的1/2,因此,用FST、FCT处理奇和偶函数比用FT更有效7

卷积是一种积分变换,常被用于各种应用问题建模。由于新型卷积对积分方程和算子理论有影响,且在展示新特性和应用方面具有较大潜力8-9,吸引了不少研究者的兴趣。近年来,有许多有关傅里叶卷积及其相关理论研究的报道。如文献[10]研究了FST与Kontorovich-Lebedev变换相关的广义卷积,并将卷积定理应用于求解积分微分方程。文献[11]给出了Wiener空间上的卷积运算,建立了Fourier-Feynman变换及其卷积定理,拓展了傅里叶分析理论体系。文献[12]提出了基于矩阵的FST、FCT以及Hartley变换,并将这些理论应用于图像加密。文献[13]研究了Hartley变换、FST、FCT的多重卷积,并将这些卷积应用于求解积分方程,得到了这类方程的显式解。文献[14-17]进一步研究了傅里叶加权卷积运算及其卷积定理,利用所得卷积理论研究了几类卷积积分方程的解。

近年来,作为FST与FCT的广义形式,分数阶正弦变换与分数阶余弦变换引得关注,产生了一些具有代表性的理论成果18-20。但尚未见关于FST与FCT的混合加权卷积成果的报道。混合加权卷积理论不仅是对原有广义卷积理论的拓展,而且可对复杂场景下的实际问题进行建模。因此,研究傅里叶正、余弦变换的混合加权卷积及其相关应用非常必要。

本文在已有文献基础上,研究傅里叶正、余弦变换混合加权卷积及其应用。首先,利用经典卷积以及傅里叶正、余弦变换的性质,给出了傅里叶正、余弦变换的混合加权卷积运算,并推导了相应的卷积定理。其次,研究混合加权卷积运算与已有卷积运算之间的关系,运用混合加权卷积性质得到了Young类不等式。最后,利用混合加权卷积运算,计算了一类卷积类积分方程的解。

1 预备知识

广义卷积定理21

K(f * g)(y)=γ(y)(K1f)(y)(K2g)(y)

其中,函数fgγL1R+K,K1,K2为不同的积分变换。当K,K1,K2为傅里叶变换时,上述广义卷积定理退化为经典的卷积定理:

F(f * g)(y)=(Ff)(y)(Fg)(y)

其中,(f * g)x为经典傅里叶卷积运算,满足

(f * g)(x)=12πRf(x-y)g(y)dy

(Ff)(y)f的傅里叶变换,满足

(Ff)(y)=12πRf(x)e-iyxdx

傅里叶余弦变换(Fcf)(y)与傅里叶正弦变换(Fsf)(y)22分别为:

(Fcf)(y)=2πR+f(x)cos(yx)dx
(Fsf)(y)=2πR+f(x)sin(yx)dx

傅里叶余弦卷积f*0g(x)与傅里叶正弦卷积(f*1g)(x)23分别为:

f*0gx=12πR+f(y)[g(x-y)+g(x+y)]dy,
f*1gx=12πR+f(y)[g(x-y)-g(x+y)]dy

傅里叶余弦与傅里叶正弦的卷积定理分别为:

Fc(f*0g(x))(y)=(Fcf)(y)(Fsg)(y)
Fs((f*1g)(x))(y)=(Fsf)(y)(Fcg)(y)

傅里叶余弦正弦卷积24

(f*2g)(x)=12πR+f(y)[sign(y-x)×g(|y-x|)+g(y+x)]dy,

且满足卷积定理:

Fc(f*2g)(x)(y)=(Fsf)(y)(Fsg)(y)

2 加权卷积及其卷积定理

给出两类新的傅里叶余、正弦加权卷积的定义,并研究这两类卷积的相关性质,推导相应的傅里叶余、正弦加权卷积定理。

定义1 设f(t)L1R,g(t)L1R+,则傅里叶余弦加权卷积f*1γg与傅里叶正弦加权卷积f*2γg分别为:

f*1γgx=14π2πRR+I1+I2+I3+I4×f(u)g(v)dudv,    x>0
f*2γgx=14π2πRR+I3+I4-I1-I2×f(u)g(v)dudv,    x>0

其中,γy=e-ycosy为权函数,且有

I1=1+iu1+iu2+x+v+12,I2=1+iu1+iu2+x+v-12,I3=1+iu1+iu2+x-v+12,I4=1+iu1+iu2+x-v-12

定理1 设f(t)L1R,g(t)L1R+,则有f*1γgf*2γgL1R+,且满足:

Fcf*1γgy=e-ycosyFfyFcgy
Fsf*2γgy=e-ycosyFfyFsgy

证明 由式(13),可得

R+f*1γg(x)dx14π2πRR+2(I1+I2+I3+I4)f(u)g(v)dudvdx

R+1+iu1+iu2+x+v+12dxπ,可得

R+I2dx=R+I3dx=R+I4dxπ

从而

R+f*1γg(x)dx2πRR+f(u)g(v)dudv2πfL1RgL1R+<

f*2γg的证明方法与f*1γg的类似,故有f*1γgf*2γgL1R+

下证式(15)。由式(4)和式(5),可得

e-ycosy(Ff)(y)(Fcg)(y)=1πRR+cosycos(yv)e-(1+iu)y×f(u)g(v)dudv

式(17),可得

1π2πRR+2cos(xy)cosycos(yv)×e-(1+iu)yf(u)g(v)dudvdy=14π2πRR+2[cosy(x+v+1)+cosy(x+v-1)+cosy(x-v+1)+cosy(x-v-1)]e-(1+iu)yf(u)g(v)dudvdy=14π2πRR+(I1+I2+I3+I4)×f(u)g(v)dudv

式(17)和式(18),可得

e-ycosy(Ff)(y)(Fcg)(y)=12π2RR+2cos(xy)(I1+I2+I3+I4)×f(u)g(v)dudvdx

式(5)和式(13),可知式(15)成立。

式(16)的证明方法与式(15)的类似。

定理1证毕。

下面给出傅里叶正、余弦加权卷积与已有卷积的关系。

定理2 设f(t)L1(R),g(t),h(t)L1(R+)

则有

(i) f*0g*1γh=g*1γf*0h

(ii) f*1γg*0h=f*1γg*0h

(iii) f*2γg*2h=f*h2γ*2g

(iv) f*2γg*1h=g*1f*1γh

其中,*0*1*2分别为式(7)、式(8)、式(11)中定义的卷积运算。

证明 先证(iv)。由式(16),可得

Fsf*2γ(g*1h)=e-ycosy(Ff)(y)(Fsg)(y)(Fch)(y)=Fsg*1f*1γh(y)

则有

f*2γ(g*1h)=g*1f*1γh

故式(iv)成立。

(i)~(iii)的证明过程与(iv)的证明类似。

定理2证毕。

定理3 设f(x)Lp(R+,(1+x2) p-1)g(x)LqR+ω(x)LrR+,其中,p,q,r>1,且满足1p+1q+1r=2,则有傅里叶余弦加权卷积不等式

0f*1γgxωxdx22-p2p1πfLpR+,1+x2p-1gLqR+ωLrR+

证明 设p1q1r1>1,满足1p+1p1=11q+1q1=11r+1r1=1,即1p1+1q1+1r1=1。令

Iu,v,x=I1+I2+I3+I4)
Uu,v,x=g(v)qp1ω(x)rp1I1+u21p1
Vu,v,x=f(u)pq1ω(x)rq11+u2p-1q1I1q1
Wu,v,x=f(u)pr1g(v)qr11+u2p-1r1I1r1

由式(21)~式(23),可得

U,V,Wu,v,x=f(u)g(v)ω(x)I(u,v,x)

式(21),可得

ULp1(R+3)p1=R+3g(v)qωxrI(u,v,x)1+u2dudvdx

由Fubini定理,可得

R+2gvqωxrR+I(u,v,x)1+u2dudvdx

由复数的性质及式(20),可得

R+I1(u,v,x)1+u2du=R+1[(x+v+1)2-u2]2+4u2duR+1u2+1duπ2

同理可得

R+I2(u,v,x)1+u2du=R+I3(u,v,x)1+u2du=R+I3(u,v,x)1+u2duπ2

因此,由式(24)~式(27),可得

ULp1(R+3)p1=R+3g(v)qωxr×I(u,v,x)1+u2dudvdx2πgLqR+qωLrR+r

再由式(20),可得

R+I1(u,v,x)dv=R+1+u2[(x+v+1)2+1-u2]2+4u2dvx+11+u2(t-u)2+u2+1dt=arctant-u1+u2x+1π

同理可得

R+Ii(u,v,x)dvπ,i=2,3,4

因此,由式(22)、式(23)、式(29)、式(30),可得

VLq1(R+3)q1=R+3f(u)p1+u2p-1ω(x)r×I(u,v,x)dudvdx4πfLp(R+,(1+x2)p-1)pωLrR+r
WLr1(R+3)r1=R+3f(u)p1+u2p-1g(v)q×I(u,v,x)dudvdx4πR+2f(u)p1+u2p-1g(v)qdudv=4πfLp(R+,(1+x2)p-1)pgLqR+q

由Hölder's不等式,可得

R+f*1γgxωxdx14π2πR+3|fu||gv|×|Iu,v,x||ωx|dudvdx=14π2πR+3Uu,v,xVu,v,x×Wu,v,xdudvdx14π2πULp1(R+3)VLq1(R+3)WLr1(R+3)22-p2p1πfLp(R+,(1+x2)p-1)gLqR+ωLrR+

定理3证毕。

3 卷积方程在积分方程中的应用

卷积方程在许多领域均有非常重要的应用,如辐射能量传播、轴震动等,特别在工程力学和数字信号处理中25,经常会遇到形如式(33)和式(39)的积分方程,求解这些方程是目前研究的热点之一。

fx +λ1-+0+ω1x,u,vϕu ×    gvdudv = hx,λ20+ftω2x,tdt +gx =kx,

其中,

ω1x,u,v=14π2πI1+I2+I3+I4)
ω2x,t=12πft[φ(x-t)+φx+t]

λ1,λ2为复数,ϕkφhL1R+为已知函数,fg为未知函数,I1,I2I3I4见定义1。

定理4 设1-λ1λ2Fcϕ*1γφ(y)0,则式(33)

L1R+上有解,即

fx=hx+(θ*0h)x-λ1ϕ*1γkx-   λ1ϕ*1γk*0θxgx=kx+(k*0θ)x-λ2(φ*0k)x-   λ2((φ*0k)*0θ)x

其中,x>0θL1R+,且满足

Fcθ(y)=λ1λ2Fcϕ*1γφ(y)1-λ1λ2Fcϕ*1γφ(y)

证明 式(33)可改写为

fx+λ1ϕ*1λgx=hxgx+λ2f*0φx=kx

运用卷积定理,得

(Fcf)(y)+λ1e-ycosy(Fϕ)(y)(Fcg)(y)=    Fch(y)(Fcg)(y)+λ2(Fcf)(y)(Fcφ)(y)=    Fck(y)

由Wiener-Levi's定理26,知存在函数θL1R+,使得

Fcθ(y)=λ1λ2Fcϕ*1γφ(y)1-λ1λ2Fcϕ*1γφ(y)

因此,由式(38),得

(Fcf)(y) = [1+Fcθ(y)]Fch(y) -λ1e-ycosy(Fϕ)(y)Fck(y) =[1+Fcθ(y)]Fch(y) -λ1Fcϕ*1γk(y) =Fch(y) + (Fc(θ*0h)(y)) -λ1Fcϕ*1γk(y) -λ1Fcϕ*1γk*0θ(y)

由FCT,有

fx=hx+(θ*0h)x-λ1ϕ*1γkx-λ1ϕ*1γk*0θx

同理可得

gx=kx+(k*0θ)x-λ2(φ*0h)x-λ2((φ*0h)*0θ)x

定理4证毕。

卷积方程在许多场合有重要应用,不同的卷积方程可解决不同的问题。如积分方程:

fx+λ1-+0+ω1x,u,vδufvdudv+0gtμx,tdt=hxgx+λ2-+0+ω3x,u,vτugvdudv+0ftμx,tdt=kx

其中,

μx,t=12πsign(t-x)[μt-x-μt+x]
ω3x,u,v=12π2πI3+I4-I1-I2

λ1,λ2为复数,δ,μ,τ,ω1,ω2L1(R+)为已知函数,f,g为未知函数。

定理5 假设1+H =\ 0,其中,

H=λ1e-ycosyFδ(y)+λ2e-ycosyFτ(y)+λ1λ2e-2ycos2yFτ(y)Fδ(y)-λ1λ2Fsδ(y)(Fsμ)(y),

式(39)具有显式解:

fx=hx+λ2τ*1γhx-λ1(μ*2k)x+(h*0ρ)x+λ2ρ*0τ*1γhx-λ1(ρ*0(μ*2h))xgx=kx+λ1δ*2γkx-λ2(σ*2h)x+(k*1ρ)x+λ1δ*2γh*1ρx-λ2((σ*1h)*1ρ)x

其中,ρL1R+,且满足

(Fcρ)=H1+H

证明 与定理4的证明类似,此证略。

4 结 论

在现有加权傅里叶正、余弦变换卷积的基础上,提出了傅里叶正、余弦变换的混合加权卷积运算,得到了相应的加权卷积定理;研究了混合加权卷积的性质以及Young类不等式;最后利用提出的混合加权卷积,讨论了两类卷积类积分方程的解,得到了相应的显式解。

http://dx.doi.org/10.3785/j.issn.1008-9497.2023.03.002

参考文献

高乾坤刘文清张玉均.

针对湍流噪声的傅里叶光谱数据处理方法

[J]. 光学学报, 20214117): 189-196. DOI:10.3788/AOS202141.1730001

[本文引用: 1]

GAO Q KLIU W QZHANG Y J.

Fourier spectral data processing method for turbulent noise

[J]. Journal of Optics, 20214117): 189-196. DOI:10. 3788/AOS202141.1730001

[本文引用: 1]

徐昭周昕白星.

基于深度学习的相位截断傅里叶变换非对称加密系统攻击方法

[J]. 物理学报, 20217014): 226-232. DOI:10.7498/aps.70. 20202075

XU ZZHOU XBAI Xet al.

Attack method of phase truncated Fourier transform asymmetric encryption system based on deep learning

[J]. Journal of Physics, 20217014): 189-196. DOI:10. 7498/aps.70.20202075

CYCHY SLECHLER SHUANG Z J.

Optimizing the nickel boride layer thickness in aspectroelectrochemical ATR-FTIR thin-film flow cell applied in glycerol oxidation

[J]. Chinese Journal of Catalysis, 20214212): 2206-2215. DOI:10. 1016/S1872-2067(20)63766-4

王娇娇于佳刘惠萍.

基于傅里叶合成全息的彩色全息制作方法

[J]. 中国激光, 2016432): 207-211. DOI:10.3788/CJL201643.0209001

[本文引用: 1]

WANG J JYU JLIU H Pet al.

Fabrication method of color holography based on Fourier synthetic holography

[J]. Chinese Journal of Lasers, 2016432): 207-211. DOI:10.3788/CJL201643.0209001

[本文引用: 1]

LI Y MWEI D YZHANG L N.

Double-encrypted watermarking algorithm based on cosine transform and fractional Fourier transform in invariant wavelet domain

[J]. Information Sciences, 2021551205-227. DOI:10.1016/j.ins.2020.11.020

[本文引用: 1]

GARCIA S RYIH S.

Supercharacters and the discrete Fourier, cosine, and sine transforms

[J]. Communications in Algebra, 2018469): 3745-3765. DOI:10.48550/arXiv.1702.02689

[本文引用: 1]

PEI S CDING J J.

Fractional cosine, sine, and Hartley transforms

[J]. IEEE Transactions on Signal Processing, 2002507): 1661-1680. DOI:10.1109/TSP.2002.1011207

[本文引用: 1]

CASTRO L PGOEL NSILVA A S.

A new convolution operator for the linear canonical transform with applications

[J]. Computational and Applied Mathematics, 2021403): 1-18. DOI:10. 1007/s40314-021-01484-9

[本文引用: 1]

ANH P KCASTRO L PTHAO P Tet al.

Two new convolutions for the fractional Fourier transform

[J]. Wireless Personal Communications, 2017922): 623-637. DOI:10.1007/s11277-016-3567-3

[本文引用: 1]

TUAN T.

On the Fourier-sine and Kontorovich-Lebedev generalized convolution transforms and their applications

[J]. Ukrainian Mathematical Journal, 2020722): 267-279. DOI:10.1007/s11253-020-01782-1

[本文引用: 1]

CHANG S JCHUNG H SCHOI J G.

Generalized Fourier-Feynman transforms and generalized convolution products on Wiener space

[J]. Indagationes Mathematicae New Series, 2017282): 566-579. DOI:10.1016/j.indag.2017.01.004

[本文引用: 1]

LIMA P H E SLIMA J BCAMPELLO de SOUZA R M.

Fractional Fourier, Hartley, cosine and sine number-theoretic transforms based on matrix functions

[J]. Circuits, Systems, and Signal Processing, 2016367): 2893-2916. DOI:10.1007/s00034-016-0447-8

[本文引用: 1]

THAO N XKHOA N MANH P T.

On the polyconvolution for Hartley, Fourier cosine and Fourier sine transforms

[J]. Integral Transforms and Special Functions, 2013247): 517-531. DOI:10. 1080/10652469.2012.714377

[本文引用: 1]

THAO N XKHOA N M.

On the generalized convolution with a weight function for Fourier sine and cosine transforms

[J]. Integral Transforms and Special Function, 2006179): 673-685. DOI:10. 1080/10652460500432071

[本文引用: 1]

THAO N XTUAN V KKHOA N M.

A generalized convolution with a weight function for the Fourier cosine and sine transforms

[J]. Fractional Calculus and Applied Analysis, 200473): 323-337.

KHOA N M.

On the generalized convolution with a weight-function for Fourier cosine, Fourier and Fourier sine transforms

[J]. Southeast Asian Bulletin of Mathematics, 200933285-298.

THAO N XTUAN V KHONG N T.

A Fourier generalized convolution transform and applications to integral equations

[J]. Fractional Calculus and Applied Analysis, 2012153): 493-508. DOI:10. 2478/s13540-012-0035-y

[本文引用: 1]

FENG QLI B Z.

Convolution theorem for fractional cosine-sine transform and its application

[J]. Mathematicl Methods in the Applied Sciences, 2017403651-3665. DOI:10.1002/mma.4251

[本文引用: 1]

FENG QWANG R B.

Fractional convolution associated with a class of integral equations

[J]. IET Signal Processing, 2020141): 15-23. DOI:10. 1049/iet-spr.2019.0140

FENG QYUAN S.

The explicit solutions for a class of fractional Fourier-Laplace convolution equations

[J]. Integral Transforms and Special Functions, 2022342): 128-144. DOI:10.1080/10652469.2022.2093870

[本文引用: 1]

KAKICHEVV ATHAO N X.

On a constructive method for the generalized integral convolution

[J]. Izvestiya Vysshikh Uchebnykh Zavedenii Mathematics, 1998131-40. DOI:10.1007/bfb0066280

[本文引用: 1]

ACHIEZER N I. Lectures on Approximation Theory[M]. MoscowScience Publishing House1965.

[本文引用: 1]

SNEDDON I N. Fourier Transforms[M]. New YorkMcGray-Hill1951.

[本文引用: 1]

THAO N XKAKICHEV V ATUAN V K.

On the generalized convolution for Fourier cosine and sine transforms

[J]. East-West Journal of Mathematics, 1998185-90.

[本文引用: 1]

LIN F RYANG S W.

A weighted H 1 seminorm regularization method for Fredholm integral equations of the first kind

[J]. International Journal of Computer Mathematics, 2014915): 1012-1029. DOI:10.1080/00207160.2013.818137

[本文引用: 1]

NAIMARK M A. Normed Algebras[M]. GroningenWolters-Noordhoff Publisher1972. DOI:10.1007/978-94-009-9260-3_3

[本文引用: 1]

/