clear?all;close?all;clc
fs=10000;??%采樣頻率
t=0:0.005:1*pi;??%采樣時(shí)間
x=4*sin（2*200*pi*t）.*（sin（2*4500*pi*t））+25*（sin（2*4500*pi*t））;??%原信號(hào)x
N=length（x）;?%計(jì)算原信號(hào)的長(zhǎng)度
f=fs*（0:N-1）/N;??%頻率分布

y=fft（x）;??%對(duì)原時(shí)域信號(hào)x進(jìn)行fft，得到頻域信號(hào)y
subplot（311）
plot（fabs（y））??%繪制原信號(hào)的頻譜y
title（‘原信號(hào)頻譜‘）

%設(shè)計(jì)濾波器FIR濾波器
f1=128/fs;??%上截止頻率（要除以采樣頻率，歸一化）
f2=4096/fs;?%下截止頻率（要除以采樣頻率，歸一化）
b=fir1（48[f1?f2]）;??%設(shè)計(jì)帶通濾波器
c=freqz（b1N）;???%頻率特性
subplot（312）???
plot（fabs（c））;??%繪制帶通濾波器
title（‘濾波器‘）
xlabel（‘f（Hz）‘）

%頻域?yàn)V波
y1=y.*c‘;???%對(duì)原信號(hào)頻域?yàn)V波，得到濾波之后的頻域信號(hào)y1（原頻域信號(hào)y*濾波器c）
subplot（313）?
pl

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

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

-----------?---------??----------?-----??----

?????????????????1436????????????????????1