小波变换系数计算实例_小波变换原理与基本案例分析

科研训练报告

一：科研训练目标及意义

在深入理解经典

Fourier

变换的基础上，通过系统介绍加窗

Fourier

变换

(

Windowed Fourier Transfrom, WFT

)

,

小波变换(

Wavelet trsnaform,WT

)等定

义，

理解时间分辨率，

频率分辨率及，

时频局部化性质等基本概念；

通过编程

实现基本的时频分析方法，理解时变信号中时频结构信息。

二：科研训练任务

1.

理解加窗

Fourier

变换的意义

2.

编程实现加窗

Fourier

分析几种合成信号

3.

利用加窗

Fourier(

小波变换

)

分析几种典型的合成信号

4.

观察分析窗函数，窗长变化，频率，分辨率的影响

5.

可尝试利用加窗

Fourier

(小波变换)正变换分析实测语音等信号

三

:

报告结构

四：加窗

Fourier

分析

v.1.Fourier

正变换

v.1.a

离散

Fourier

正变换数学描述

有限长序列通过离散傅里叶变换，并通过快速傅里叶变换算法

FFT

即得：

𝑋

(

𝑘

)

=

∑

𝑥

(

𝑗

)

𝑊

ℵ

(𝑗−1)(𝑅−1)

ℵ

𝑗=1

其中：

𝑊

ℵ

=

𝑒

−2∗𝜋∗𝑗/𝑁

v.1.b.Fourier

正变换

Matlab

实现

被分析信号为

N

维列向量；

x

=

[

x

(

0

)

,

x

(

1

)

,

x

(

2

)

,

x

(

3

)

⋯

x

(

n

−

1

)]

𝑇

变换因子为：

[

𝑊

0

⋯

𝑊

0

⋮

⋱

⋮

𝑊

0

⋯

𝑊

(𝑁−1)(𝑁−1)

]

DFT

结果为

X

=

[

𝑋

(

0

)

,

X

(

1

)

,

X

(

2

)

,

X

(

3

)

⋯

X

(

n

−

1

)]

𝑇

以矩阵形式表示

DFT

：

𝑋

𝑁

=

𝑊

𝑁

𝑥

𝑁

(具体代码参见代码

code erst

)

v.1.c.Fourier

正变案例分析总结

step1

检测代码的正确性



当输入一个离散的冲激函数的

Delta

函数，

其

DFT

变换的幅度谱为一

等幅，幅度为

1

的白色谱。

(参见图

a.3.1

)

Step2

分析信号特性



当输入一个正弦函数时，

运行结果如下：

幅度谱为一对称的离散冲激

信号，

冲激的横坐标为其频率，纵坐标为其幅度。由观测可看出在主

频的附近有明显的频谱泄露。

(参见图

a.3.2

)



当输入等幅的正弦波叠加则运行结果显示，

两个正弦波的幅值并不相

同，存在明显的栅栏现象。

(参见图

a.3.3

)



当输入一分时频信号时，该信号表达式为：

𝑦

=

{

𝑡

<

0.5,

(

sin

(

2

∗

pi

∗

50

∗

t1

))

𝑡

≥

0.5,

(

sin

(

2

∗

pi

∗

120

∗

t1

))

其快速傅里叶变换幅频特性显示如图

a.3.4

，与图

a.3.3

相比，其信号频率仍然

为

50

，

120hertz

。并无变化。

图

v.a.3.1

输入信号为序列长度为

32

的离散冲激信号

Xot(1)

与

Xot(2)

的区别在于其相位不同，

后面是经过

DFT

变换的幅度谱和相位

谱

横坐标为点数，纵坐标为幅值和相位。

图

v.1.c.1

输入信号为

sin(2*pi*50*t)

，横坐标为模拟频率，纵坐标为其幅度值。

图

v.1.c.2

输入信号为

sin(2*pi*50*t)+sin(2*pi*120*t)

横坐标为模拟频率，纵坐标为其幅度值。