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