emd信号分解

1.什么是HHT？

HHT就是先将信号进行经验模态分解（EMD分解），然后将分解后的每个IMF分量进行Hilbert变换，得到信号的时频属性的一种时频分析方法。

2.EMD分解的步骤。

EMD分解的流程图如下：

3.实例演示。

给定频率分别为10Hz和35Hz的两个正弦信号相叠加的复合信号，采样频率fs=2048Hz的信号，表达式如下：y=5sin(2*pi*10t)+5*sin(2*pi*35t)

(1)为了对比，先用fft对求上述信号的幅频和相频曲线。

1.

2.

3.

4.

5.

6.

7.

8. function fftfenxi clear;clc; N=2048; %fft默认计算的信号是从0开始的 t=linspace(1,2,N);deta=t(2)-t(1);1/deta x=5*sin(2*pi*10*t)+5*sin(2*pi*35*t); % N1=256;N2=512;w1=0.2*2*pi;w2=0.3*2*pi;w3=0.4*2*pi; %

x=(t>=-200&t<=-200+N1*deta).*sin(w1*t)+(t>-200+N1*deta&t<=-200+N2*d

eta).*sin(w2*t)

1.什么是HHT？

HHT就是先将信号进行经验模态分解（EMD分解），然后将分解后的每个IMF分量进行Hilbert变换，得到信号的时频属性的一种时频分析方法。

2.EMD分解的步骤。

EMD分解的流程图如下：

3.实例演示。

给定频率分别为10Hz和35Hz的两个正弦信号相叠加的复合信号，采样频率fs=2048Hz的信号，表达式如下：y=5sin(2*pi*10t)+5*sin(2*pi*35t)

(1)为了对比，先用fft对求上述信号的幅频和相频曲线。

1. 2. 3. 4. 5. 6. 7. 8.

function fftfenxi clear;clc; N=2048;

?t默认计算的信号是从0开始的

t=linspace(1,2,N);deta=t(2)-t(1);1/deta x=5*sin(2*pi*10*t)+5*sin(2*pi*35*t);

% N1=256;N2=512;w1=0.2*2*pi;w2=0.3*2*pi;w3=0.4*2*pi; %

x=(t>=-200&t<=-200+N1*deta).*sin(w1*t)+(t>-200+N1*deta&t<=-200+N2*deta).*sin(w2*t)+(t>

1.什么是HHT？

HHT就是先将信号进行经验模态分解（EMD分解），然后将分解后的每个IMF分量进行Hilbert变换，得到信号的时频属性的一种时频分析方法。

2.EMD分解的步骤。

EMD分解的流程图如下：

3.实例演示。

给定频率分别为10Hz和35Hz的两个正弦信号相叠加的复合信号，采样频率fs=2048Hz的信号，表达式如下：y=5sin(2*pi*10t)+5*sin(2*pi*35t)

(1)为了对比，先用fft对求上述信号的幅频和相频曲线。

1. 2. 3. 4. 5. 6. 7. 8.

function fftfenxi clear;clc; N=2048;

?t默认计算的信号是从0开始的

t=linspace(1,2,N);deta=t(2)-t(1);1/deta x=5*sin(2*pi*10*t)+5*sin(2*pi*35*t);

% N1=256;N2=512;w1=0.2*2*pi;w2=0.3*2*pi;w3=0.4*2*pi; %

x=(t>=-200&t<=-200+N1*deta).*sin(w1*t)+(t>-200+N1*deta&t<=-200+N2*deta).*sin(w2*t)+(t>

1.什么是HHT？

HHT就是先将信号进行经验模态分解（EMD分解），然后将分解后的每个IMF分量进行Hilbert变换，得到信号的时频属性的一种时频分析方法。

2.EMD分解的步骤。

EMD分解的流程图如下：

3.实例演示。

给定频率分别为10Hz和35Hz的两个正弦信号相叠加的复合信号，采样频率fs=2048Hz的信号，表达式如下：y=5sin(2*pi*10t)+5*sin(2*pi*35t)

(1)为了对比，先用fft对求上述信号的幅频和相频曲线。

1. 2. 3. 4. 5. 6. 7. 8.

function fftfenxi clear;clc; N=2048;

?t默认计算的信号是从0开始的

t=linspace(1,2,N);deta=t(2)-t(1);1/deta x=5*sin(2*pi*10*t)+5*sin(2*pi*35*t);

% N1=256;N2=512;w1=0.2*2*pi;w2=0.3*2*pi;w3=0.4*2*pi; %

x=(t>=-200&t<=-200+N1*deta).*sin(w1*t)+(t>-200+N1*deta&t<=-200+N2*deta).*sin(w2*t)+(t>

西 南 交 通 大 学 本科毕业设计（论文）

EMD信号分析方法端点效应分析

年 级:2005级 学 号:20051198 姓 名:郭云喜

专 业:测控技术与仪器 指导老师:宁静

2009年6月

西南交通大学本科生毕业设计（论文） 第I页

院 系 机械工程学院 专 业 测控技术与仪器 年 级 2005级 姓 名 郭云喜

题 目 EMD信号分析方法端点效应分析

指导教师

评 语

西 南 交 通 大 学 本科毕业设计（论文）

EMD信号分析方法端点效应分析

年 级:2005级 学 号:20051198 姓 名:郭云喜

专 业:测控技术与仪器 指导老师:宁静

2009年6月

西南交通大学本科生毕业设计（论文） 第I页

院 系 机械工程学院 专 业 测控技术与仪器 年 级 2005级 姓 名 郭云喜

题 目 EMD信号分析方法端点效应分析

指导教师

评 语

西 南 交 通 大 学 本科毕业设计（论文）

EMD信号分析方法端点效应分析

年 级:2005级 学 号:20051198 姓 名:郭云喜

专 业:测控技术与仪器 指导老师:宁静

2009年6月

西南交通大学本科生毕业设计（论文） 第I页

院 系 机械工程学院 专 业 测控技术与仪器 年 级 2005级 姓 名 郭云喜

题 目 EMD信号分析方法端点效应分析

指导教师

评 语

西 南 交 通 大 学 本科毕业设计（论文）

EMD信号分析方法端点效应分析

年 级:2005级 学 号:20051198 姓 名:郭云喜

专 业:测控技术与仪器 指导老师:宁静

2009年6月

西南交通大学本科生毕业设计（论文） 第I页

院 系 机械工程学院 专 业 测控技术与仪器 年 级 2005级 姓 名 郭云喜

题 目 EMD信号分析方法端点效应分析

指导教师

评 语

实验四 方波信号的分解与合成

任何电信号都是由各种不同频率、幅度和初相的正弦波迭加而成的。1822年法国数学家傅里叶在研究热传导理论时提出并证明了将周期函数展开为正弦级数的原理。奠定了傅里叶级数的理论基础、揭示了周期信号的本质，即任何周期信号（正弦信号除外）都可以看作是由无数不同频率、不同幅度的正弦波信号叠加而成的，就像物质都是由分子或者原子构成一样。周期信号的基本单元信号是正弦谐波信号。

一、实验目的

1、通过对周期方波信号进行分解，验证周期信号可以展开成正弦无穷级数的基本原理，了解周期方波信号的组成原理。

2、测量各次谐波的频率与幅度，分析方波信号的频谱。 3、观察基波与不同谐波合成时的变化规律。

4、通过方波信号合成的实验，了解数字通信中利用窄带通信系统传输数字信号（方波信号）的本质原理。

二、实验原理

1、一般周期信号的正弦傅里叶级数

按照傅里叶级数原理，任何周期信号在满足狄利克雷条件时都可以展开成如式2-3-1所示的无穷级数

?a0?A0? f(t)???ancos(n?t)??bnsin(n?t)???Ancos(n?t??n) （2-4-1）

2n?12n?1n?1其中Ancos(n?t??n)称为周期信号的n谐波分量，n次谐波

实验四 方波信号的分解与合成

任何电信号都是由各种不同频率、幅度和初相的正弦波迭加而成的。1822年法国数学家傅里叶在研究热传导理论时提出并证明了将周期函数展开为正弦级数的原理。奠定了傅里叶级数的理论基础、揭示了周期信号的本质，即任何周期信号（正弦信号除外）都可以看作是由无数不同频率、不同幅度的正弦波信号叠加而成的，就像物质都是由分子或者原子构成一样。周期信号的基本单元信号是正弦谐波信号。

一、实验目的

1、通过对周期方波信号进行分解，验证周期信号可以展开成正弦无穷级数的基本原理，了解周期方波信号的组成原理。

2、测量各次谐波的频率与幅度，分析方波信号的频谱。 3、观察基波与不同谐波合成时的变化规律。

4、通过方波信号合成的实验，了解数字通信中利用窄带通信系统传输数字信号（方波信号）的本质原理。

二、实验原理

1、一般周期信号的正弦傅里叶级数

按照傅里叶级数原理，任何周期信号在满足狄利克雷条件时都可以展开成如式2-3-1所示的无穷级数

?a0?A0? f(t)???ancos(n?t)??bnsin(n?t)???Ancos(n?t??n) （2-4-1）

2n?12n?1n?1其中Ancos(n?t??n)称为周期信号的n谐波分量，n次谐波