highpasslvbo.m

来自「实现数据的高通、带阻、低通的数字滤波。 非常适合初学者使用」· M 代码 · 共 35 行

function y=highpasslvbo(fs,fc,x);
%实现数据的高通滤波
%FIR的窗口法：汉明窗
%48阶，通常用
%fs：数据x的采样频率；fc：3dB截止频率。
% load sda52;
% x=sda52;
% fs=4*10^4;

X=x;
L=length(X);
f=fs;
l=L./2;

a=1;
wn=(fc*2)./fs;
% wn=0.3;
b=fir1(48,wn,'high');
% freqz(b,1,2048);
% hold on;

y=filter(b,a,X);
Y=fft(y);
Pyy = Y.* conj(Y)/L;
F = (f./1000)*(0:l)/L;
plot(F,Pyy(1:l+1));
hold on;
% axis([0 0.5 0 15000]);
title('Frequency content of Y')
xlabel('Frequency /KHz)')
ylabel('Amplitude:W/Hz')

y=y;%filtered data
%