MATLAB希尔伯特Hilbert变换求包络谱

1、内容简介
略

2、内容说明
略

3、仿真分析

clear all;close all;clc
fs=10000; %采样频率
t=0:0.005:1*pi; %采样时间
x=4*sin(2*200*pi*t).*(sin(2*4500*pi*t))+25*(sin(2*4500*pi*t)); %原信号x
N=length(x); %计算原信号的长度
f=fs*(0:N-1)/N; %频率分布

y=fft(x); %对原时域信号x进行fft，得到频域信号y
subplot(311)
plot(f,abs(y)) %绘制原信号的频谱y
title('原信号频谱')

%设计滤波器FIR滤波器
f1=128/fs; %上截止频率（要除以采样频率，归一化）
f2=4096/fs; %下截止频率（要除以采样频率，归一化）
b=fir1(48,[f1 f2]); %设计带通滤波器
c=freqz(b,1,N); %频率特性
subplot(312)
plot(f,abs(c)); %绘制带通滤波器
title('滤波器')
xlabel('f(Hz)')

%频域滤波
y1=y.*c'; %对原信号频域滤波，得到滤波之后的频域信号y1（原频域信号y*滤波器c）
subplot(313)
plot(f,abs(y1)) %绘制滤波之后的频域信号y1
title('原信号频域滤波之后')

%反变换y1，得到滤波之后的时域信号x1
x1=ifft(y1); %反变换y1，得到滤波之后的时域信号x1
figure
subplot(211)
plot(t,x) %绘制时域原信号x
title('时域原信号')
subplot(212)
plot(t,real(x1)) %绘制滤波之后的时域信号x1
title('滤波之后的时域信号')

%hilbert变换，对x1求包络线
x2=hilbert(real(x1)); %x1的希尔伯特变换x2
x3=abs(x2); %x2取模，得到x3

%对x3进行fft，求包络谱
nfft=1024;
p=abs(fft(x3,nfft)); %将x3补齐到1024个，并fft，得到p，就是包络线的fft---包络谱
figure
plot((0:nfft-1)/nfft*fs,p) %绘制包络谱
title('包络谱')
msgbox('MATLAB编程答疑，请加QQ: 1762016542','MATLAB答疑','help')

4、参考论文
略
答1疑2咨3询4企q鹅q号：1762016542

MATLAB希尔伯特Hilbert变换求包络谱相关推荐

Hilbert变换求信号的包络线及MATLAB代码
一.Hilbert变换表达式 Hilbert变换是信号与的卷积,表达式如下: 由于本质是卷积,因此可以从"线性系统","调幅-调频"等角度思考.进一步可参考如下 ...
VMD分解，matlab代码，包络线，包络谱，中心频率，峭度值，能量熵，近似熵，包络熵，频谱图，希尔伯特变换，包含所有程序MATLAB代码，-西储大学数据集为例
目录 1.选取数据 2.VMD函数-matlab代码 3.采用matlab脚本导入数据并做VMD分解 4. VMD分解图 5.计算中心频率 6.画包络线 7. 画包络谱 8. 计算峭度值 9.计算能量 ...
画幅值matlab,好用的画包络谱和幅值谱matlab函数
使用范例: x=@(t) (1+0.5*cos(9*pi*t)).*cos(200*pi*t+2*cos(10*pi*t))+sin(pi*t).*sin(30*pi*t); t=0:0.01:9.9 ...
MATLAB希尔伯特黄变换HHT
这两天在学习希尔伯特黄变换,也就是HHT,趁着学习的劲赶紧整理整理,用的是MATLAB进行编程,所用到的工具箱便是EMD工具箱,链接如下,请自行下载. 希尔伯特黄变换HHT_HHT-电信代码类资源-C ...
希尔伯特(Hilbert)变换
2.3希尔伯特变换 2.3.1 希氏变换希氏变换是完全在时域中进行的一种特殊的正交变换.也可以看成它是由一种特殊的滤波器完成的. 为了便于理解变换特点,我们首先讨论这种变换在频域中的 ...
matlab画傅立叶变换后相位谱,对一幅图像进行傅里叶变换后，包含频谱（也叫幅度谱）和相位谱两部分，请问那一部分更重要？...
对一幅图像进行傅里叶变换后,包含频谱(也叫幅度谱)和相位谱两部分,请问那一部分更重要? 更多相关问题 [判断题] 抗盐产品切换到中分时,生产线不需要酸洗. [判断题] 机械密封中的载荷系数不是越大越好 ...
经验模式分解（EMD）及希尔伯特-黄变换（HHT）简介及matlab实现
本文介绍过程涉及到两个链接工具包,可以自己网上搜索下载,以下提供了网盘下载的地址,因为作者主要做语音方面工作,所以后面的说明主要以说话人识别为例.(链接:https://pan.baidu.com/s ...
基于EMD-HHT包络谱与堆栈降噪自编码SDAE轴承故障诊断
目录 0.引言 1.方法原理 1.1 EMD-HHT包络谱 1.2 PCA 1.3 SDAE 2.流程 3.具体实现 3.1 数据准备 3.2 EMD-HHT包络谱求取 3.3 PCA降维 3.4 故 ...
频谱、边际谱、包络谱
傅里叶谱(即频谱)表示:某一点频率上的幅值表示在整个信号里和在整个时间范围内,有一个含有此频率的三角函数组分.(横坐标为频率,纵坐标为幅值) 边际谱:作用不同:边际谱可以处理非平稳信号,如果信号中存在 ...
【老生谈算法】matlab实现FFT变换算法源码——FFT变换算法
Matlab实现FFT变换(单边谱及双边谱) 1.文档下载: 本算法已经整理成文档如下,有需要的朋友可以点击进行下载序号文档(点击下载) 本项目文档 [老生谈算法]Matlab实现FFT变换程序源 ...

MATLAB希尔伯特Hilbert变换求包络谱

MATLAB希尔伯特Hilbert变换求包络谱相关推荐

最新文章

热门文章