带通滤波器的matlab程序设计

fs=500;                                                 % 采样率
n=length(data);                                         % 滤波数据长度
t = 0:1/fs:(n-1)/fs;                                    % 横坐标刻度值
lpass=100;
hpass=200;                                              % 带通滤波的带宽
filterorder =2;                                         % 滤波器阶数
filtercutoff = [2*lpass/fs 2*hpass/fs];                 % 频率归一化 对应频率/nyquist频率，nyquist频率是采样率的一半 
[f_b, f_a] = butter(filterorder,filtercutoff);          % 调用butter()函数返回巴特沃斯滤波器传递函数的零极点系数和增益
sig= filtfilt(f_b,f_a,data);                            % 零相位数字滤波 通过x在正向和反向两个方向上处理输入数据来执行零相位数字滤波 。在向前过滤数据后，反转过滤后的序列，然后将其运回过滤器。
noise=data-sig;                                         % 噪声为原始信号减去滤波后信号的剩余信号
SNR=10*log((sig.^2)/(noise.^2));                        % 10lg(Ps/Pn)，其中Ps和Pn分别代表信号和噪声的有效功率，也可以换算成电压幅值的比率关系：20Lg(Vs/Vn)，Vs和Vn分别代表信号和噪声电压的“有效值”。用dB表示。