LBF.zip_DSP 滤波_DSP，滤波

clear all close all clc %%初始值设置 fs=10^5; %采样频率 f=200; %正弦信号频率 %%信号产生 t=0:1/fs:0.05; %信号采样的时间值 original_signal=sin(2*pi*f*t); %原始信号 mix_signal=original_signal+0.1*randn(1, length(t)); %加噪声后的信号 %%滤波器设计 fp=200; %通带边缘频率 fstop=550; %阻带边缘频率 deltaw=2*pi*(fstop-fp)/fs; %数字频率过渡带 fc=(fp+fstop)/2; %截止频率 wc=2*pi*fc/fs; %数字截止频率 N=ceil(6.6*pi/deltaw); %窗口长度 n=1: N; %窗口长度取样 tao=floor(N/2); %冲击响应延时 hd=sin((n-tao)*wc)./((n-tao)*pi);%理想滤波器的冲击响应 %手动计算0/0式的极限 if(tao<=N &&tao>0) hd(tao)=wc/pi; end w=0.54-0.46*cos(2*pi*n/(N-1)); %汉明窗函数 h=w.*hd; %实际滤波器的冲击响应 %%滤波 y=filter(h, 1, mix_signal); %y为滤波后的信号 %%画图 figure(2) plot(mix_signal, 'r.-'); hold on plot(original_signal, 'g.-'); plot(y, 'b.-');