function [H,f]=bp_f(n,f_start,f_cutoff,df1,fs,p)
%帶通濾波器函數    輸入設計的濾波器參數，產生帶通濾波器頻率特性函數H和頻率向量f
%------------------------輸入參數
%n  帶通濾波器的輸入信號長度
%f_start    通帶起始頻率
%f_cutoff   帶通濾波器的截止頻率
%df1        頻率分辨率
%fs         抽樣頻率
%p          濾波器幅度
%----------------------輸出(返回)參數
%H  帶通濾波器頻率響應
%f  頻率向量
n_cutoff = floor(f_cutoff/df1*1.5);      %設計濾波器
n_start = floor(f_start/df1/1.5); 
 
f = [0:df1:df1*(n-1)] -fs/2;    %頻率向量
H = zeros(size(f));
H(n_start+1:n_cutoff) = p*ones(1,n_cutoff-n_start);
H(length(f) - n_cutoff+1:length(f)-n_start) = p*ones(1,n_cutoff-n_start);