时频分布的变形

✍ dations ◷ 2025-04-03 17:17:24 #信号处理

时频分布的变形(Motions on the time-frequency distribution)是指信号经过时频分析(Time-frequency analysis)的结果,在时频平面(time-frequency plane)上图形的各种变形。

对时频分布做各种变形,皆有其对应的物理意义。常见的时频分布的变形有下列几种:水平平移、铅直平移、扩张、斜推、旋转。时频分布的变形在分离信号、滤波器设计、取样定理、调变及多工…等领域上都有相当的帮助。

简单来说,水平平移就是信号在时间轴上的平移。例如一个在0秒至3秒之间大小为1,其他为0的方波信号,被平移成为10秒至13秒之间大小为1,其他为0的方波。若一个信号x经过水平平移t0时间单位后得到的信号为y,则x与y的时频分布关系为:


铅直平移就是信号在频率轴上的平移,在通讯上也被称为调变。将一个基频的讯号经过调变之后变换到射频来传送,可以更加有效的利用频率资源。若一个信号x经过铅直平移f0频率单位后得到y,则


扩张变形对时频分布图形造成的影响是:在时间轴拉伸a倍,而在频率轴压缩a倍。扩张变形的重要特性是时频分布占有的面积不变。若一个信号x经过a倍的扩张变形,得到的结果为y,则

时频分析的扩张

沿着时间轴方向做斜推变形,造成的影响是:时间轴方向的位移量与频率大小成正比。反之,沿着频率轴方向做斜推变形,则会使频率轴方向的位移量与时间大小成正比。沿时间轴的斜推变形,是信号与chirp函数做折积运算的结果;沿频率轴的斜推变形,是信号与chirp函数相乘的结果。若一个信号x沿着时间轴做a倍的斜推变形后得到的结果为y,则

x ( t ) = y ( t ) e j π t 2 a {\displaystyle x(t)=y(t)e^{j\pi {\frac {t^{2}}{a}}}}

W y ( t , f ) = W x ( t a f , f ) {\displaystyle W_{y}(t,\,f)=W_{x}(t-af,\,f)}

x ( t ) = y ( t ) e j π t 2 a {\displaystyle x(t)=y(t)*e^{j\pi {\frac {t^{2}}{a}}}}

W y ( t , f ) = W x ( t , f a t ) {\displaystyle W_{y}(t,\,f)=W_{x}(t,\,f-at)}

时频分析的斜推

旋转变形顾名思义就是把图形对着原点旋转。对信号做傅里叶变换会将图形顺时针方向旋转90度;同样的,做傅利业反变换会将图形逆时钟旋转90度。而分数傅里叶变换可将图形旋转任意的角度。若一个信号x做傅里叶变换后得到的结果为y,则

y = F { x ( t ) } {\displaystyle y=\displaystyle {\mathcal {F}}\left\{x(t)\right\}}

W y ( t , f ) = W x ( f , t ) {\displaystyle W_{y}(t,\,f)=W_{x}(-f,\,t)}

分數傅立葉轉換 

第一种定义:

第二种定义:

ϕ = 0.5 a π {\displaystyle \phi =0.5a\pi } , a {\displaystyle a} 为实数。

a = 1 {\displaystyle a=1} 时 (亦即 ϕ = 0.5 π {\displaystyle \phi =0.5\pi } ),分数傅里叶变换就成了傅里叶变换。

分数傅里叶转换

当一个讯号的时频分布呈现不规则形状时,会有许多不必要的带宽浪费,此时可以将频率与时间的关系近似于多项式,只要将原讯号乘上原讯号的共轭函数,就可以将时频分布拉直,如同直流讯号或是单频讯号的时频分布,这样除了可以降低传输时的带宽之外,因为在时频分布上的面积缩小所以亦可降低总取样点数。

y ( t ) A e j ϕ ( t ) {\displaystyle y(t)\cong Ae^{-j\phi (t)}}

ϕ ( t ) = k = 0 n a k t k {\displaystyle \phi (t)=\sum _{k=0}^{n}a_{k}t^{k}}

f r e q u e n c y = 1 2 π d ϕ ( t ) d t {\displaystyle frequency={\frac {1}{2\pi }}{\frac {\mathrm {d} \phi (t)}{\mathrm {d} t}}}


x ( t ) = y ( t ) e j ϕ ( t ) {\displaystyle x(t)=y(t)e^{j\phi (t)}}

W x ( t , f ) = W y ( t , f 1 2 π d ϕ ( t ) d t ) {\displaystyle W_{x}(t,\,f)=W_{y}(t,\,f-{\frac {1}{2\pi }}{\frac {\mathrm {d} \phi (t)}{\mathrm {d} t}})}


Generalized shearing

利用线性正则变换可以把时频分布做任意的线性变形

定义

线性正则变换有四个参数(a, b, c, d)和一个限制,也就是有三个自由参数,第四个参数则由另外三个所决定。

矩阵 {\displaystyle {\begin{bmatrix}a&b\\c&d\end{bmatrix}}} 的行列式ad - bc = 1


在这里用 F ( a , b , c , d ) ( u ) {\displaystyle F_{(a,b,c,d)}(u)} 表示 f ( u ) {\displaystyle f(u)} 做线性正则变换的结果, O F ( a , b , c , d ) {\displaystyle O_{F}^{(a,b,c,d)}} 表示做参数为{a,b,c,d}的线性正则变换。

线性正则转换

传统的滤波器设计,是在频域对不同频率给定不同的频率响应,借此压抑或是强化某些频率的能量。因为一般而言被处理的信号都是时变的,在加入时频分析工具以及变形之后,可以对信号做更复杂的处理,得到更好的效果,例如:分数傅里叶变换,可以将讯号时频分布旋转至适当角度让讯号及干扰的cut-off line与水平轴垂直,此时讯号再与干扰会在时间轴上分离只要在时间轴上乘上一个适当的window function,就可以分离讯号及噪声。

水平cutoff line

时频分布所占的面积,相当于一个信号所需要取样的次数。同上原因,对信号的时频分布进行适当地旋转,变形及切割之后,可以使取样点数降低或是让取样方法更简单。

讯号时频分割

讯号时频旋转

在通讯系统中,我们为了有效的利用传输资源,会进行分时多工或是分频多工,把时间或频率切割成小的区块再将讯号载在其中。同样的,若利用时频分布的变形,我们有机会把信号调整成更适合小区块的形状再放入,进而节省通讯资源。

利用时频分析处理讯号多工利用时频分析处理讯号多工

相关

  • 杨森制药杨森制药(英语:Janssen Pharmaceutica)是一家比利时制药公司。1953年由比利时化学家保罗·杨森创建,1961年10月24日被美国强生公司收购。1976年保罗·杨森遇见马海德,在几天商
  • 自然 (消歧义)自然可以指:
  • Cavalier-Smith汤玛斯·卡弗利尔-史密斯(Thomas Cavalier-Smith,1942年10月21日-)是一位英国演化生物学家及牛津大学动物学系教授。发表过许多原生生物的分类研究,主要贡献在于提出六大界来区分
  • 南惠提尔南惠提尔(英语:South Whittier)是位于美国加利福尼亚州洛杉矶县的一个人口普查指定地区。南惠提尔的座标为33°56′05″N 118°01′51″W / 33.93472°N 118.03083°W / 33.934
  • 刘易斯·卡斯刘易斯·卡斯(Lewis Cass,1782年10月9日-1866年6月17日),美国军官、政治家,曾任密歇根领地总督(1813年-1831年)、美国战争部长(1831年-1836年)、美国参议员(1845年-1848年、1849年-1857
  • 月部月部,为汉字索引中的部首之一,康熙字典214个部首中的第七十四个(四划的则为第十四个)。就繁体和简体中文中,月部归于四划部首。月部通常是从上、下、左、右方均可为部字。且无其
  • 石板滩街道 (成都市)石板滩街道,原为石板滩镇,是中华人民共和国四川省成都市新都区下辖的一个乡镇级行政单位。2019年12月,撤销木兰镇、石板滩镇,设立石板滩街道,以原木兰镇和原石板滩镇所属行政区域
  • 李旦福李旦福(1971年4月-),中国浙江省玉环市人,中学高级教师,现任玉环市城关第一初级中学校长、党支部书记,玉环市教育局党委委员,中国教育技术协会中小学专业委员会第五届委员会理事。他
  • 霍达 (作家)霍达(1945年11月26日-),北京人,回族,中国女作家。历任第八届全国政协委员,九届全国人大代表,十届、十一届全国政协常委。霍达的家庭是个珠玉世家,自幼酷爱文学艺术。1966年毕业于北京
  • GLAAD媒体奖GLAAD媒体奖(GLAAD Media Award)是由GLAAD授予的奖项,以表彰和奖励媒体的各个分支机构,因为它们在女同性恋、男同性恋、双性恋和变性者(LGBT)社区中的杰出代表以及影响他们生活的