Hilbert

MATLAB希尔伯特Hilbert变换求包络谱 源程序代码。

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（fabs（y））  %绘制原信号的频谱y
title（‘原信号频谱‘）

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

%频域滤波
y1=y.*c‘;   %对原信号频域滤波，得到滤波之后的频域信号y1（原频域信号y*滤波器c）
subplot（313） 
pl

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1436  2014-03-15 22:55  Hilbert.m

----------- ---------  ---------- -----  ----

                 1436                    1