Fs=2000;
fp=700;
fs=500;
wp=2*pi*fp/Fs;%通帶截止頻率%
ws=2*pi*fs/Fs;%阻帶截止頻率%
Bt=wp-ws;%設計過渡帶帶寬%
NO=ceil（6.2*pi/Bt）;%計算所需h（n）長度NO，ceil取大于等于x的最小整數(shù)%
wc=（wp+ws）/2/pi;%計算理想帶通濾波器截止頻率%
N=NO+mod（NO+12）-18;%確保h（n）長度N是奇數(shù)%
n=0:N-1;
hn=fir1（N-1wchanning（N））;%調(diào)用firl計算帶通FIR數(shù)字濾波器的h（n）%
[Hw]=freqz（hn11024）;%計算頻率響應函數(shù)h（n）%
figure（1）;
magH=20*log10（abs（H）/max（abs（H）））;%計算幅度%
subplot（311）;
stem（nhn‘.‘）;
xlabel（‘n‘）;
ylabel（‘h（n）‘）;
title（‘漢寧窗FIR數(shù)字低通濾波器的單位脈沖響應‘）;
subplot（312）;
plot（w/pi*Fs/2magH）;%繪制幅度特性%
title（‘漢寧窗FIR數(shù)字低通濾波器的幅度特性‘）;
xlabel（‘頻率/Hz‘）;
ylabel（‘幅度magH‘）;
grid?on;
subplot（313）;
plot（w/pi*Fs/2unwrap（angle（H）））;
title（‘漢寧窗FIR數(shù)字低通濾波器的相位特性‘）;
n=0:1023;
dt=1/Fs;
t=n*dt;
Tp=1024*dt;
f1=250;
f2=800;
x=sin（2*pi*f1*t）+sin（2*pi*f2*t）;
y=filter（hn1x）;
%=====================原始信號
f1=sin（2*pi*f1*t）;
f2=sin（2*pi*f2*t）;
N=500;
F1=abs（fft（f1N））;
F2=abs（fft（f2N））;
X=abs（fft（xN））;
Y=abs（fft（yN））;
figure（2）;
subplot（411）;
plot（tf1）;
title（‘原始信號f1（t）‘）;