MATLAB希尔伯特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） 
plot（fabs（y1））  %绘制滤波之后的频域信号y1
title（‘原信号频域滤波之后‘）

%反变换y1，得到滤波之后的时域信号x1
x1=ifft（y1）; %反变换y1，得到滤波之后的时域信号x1
figure
subplot（211）
plot（tx）  %绘制时域原信号x
title（‘时域原信号‘）  
subplot（212）
plot（treal（x1））   %绘制滤波之后的时域信号x1
title（‘滤波之

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

     文件       1436  2014-03-15 22:55  MATLAB希尔伯特Hilbert变换求包络谱 源程序代码\chengxu.m

     文件        331  2014-01-07 20:20  MATLAB希尔伯特Hilbert变换求包络谱 源程序代码\Read Me.url

     文件        294  2014-02-06 11:28  MATLAB希尔伯特Hilbert变换求包络谱 源程序代码\联系我们.txt

     目录          0  2014-03-15 22:55  MATLAB希尔伯特Hilbert变换求包络谱 源程序代码

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

                 2061                    4