高通滤波器设计,高通滤波器的程序设计

高通滤波器的程序设计  

clear

close all

clc

% 系统采样率

fs=8000;

% 原始信号

tt=0:1/fs:0.1;  

y=sin(2*pi*100*tt)+3*sin(2*pi*500*tt)+5*sin(2*pi*1000*tt);

%显示原始信号

subplot(211); plot(tt,y);  

%对原始信号进行谱分析 fy=abs(fftshift(fft(y)));

ww=linspace(-fs/2,fs/2,length(fy));

%显示原始信号的频谱

subplot(212);

plot(ww,fy);

%分离高频，设计高通滤波器  

[N,wc]=buttord(700/4000,900/4000,1,20);

[B,A]=butter(N,wc,'high');

%figure  

freqz(B,A,512);

title('HPF');

% 滤波

figure  

y3=filter(B,A,y);

 % 滤波器启动时的波形会有畸变

y3= y3(50:end);   

tt=0:1/fs:(length(y3)-1)/fs;

%显示滤波后的波形及频谱

subplot(211);plot(tt,y3);

fy3=abs(fftshift(fft(y3)));  

ww=linspace(-fs/2,fs/2,length(fy3));

subplot(212);

plot(ww,fy3);

经过高通滤波器后信号的波形与频谱

原始信号波形及频谱