samp7_8.m

选取滤波器（低通滤波器和带阻滤波器）采样点之间信号并与对应的输入信号进行比较的源程序。

%Samp7_8
clf;N=300;Fs=1000;         %数据总数和采样频率
Order=100;    %滤波器阶数
n=[0:N-1];t=n/Fs;     %时间序列
wn=[100 200]/(Fs/2); %边界频率转换为归一化频率，据6-19式
b=fir1(Order,wn,'stop');          % 设计 100 阶的带阻滤波器
figure(1)
[H,f]=freqz(b,1,512,Fs);
subplot(2,1,1),plot(f,20*log10(abs(H)))
xlabel('频率/Hz');ylabel('振幅/dB');grid on;
subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))
xlabel('频率/Hz');ylabel('相位/^o');grid on;
f1=50;f2=150;  %输入信号频率
x=sin(2*pi*f1*t)+sin(2*pi*f2*t);       %输入信号
y=fftfilt(b,x,256);% 对数据进行滤波
figure(2)
subplot(2,1,1);plot(t,x); grid on;                %绘制输入信号
title('输入信号');
subplot(2,1,2);plot(t,y);                   %绘制输出信号
hold on;plot(Order/2/Fs*ones(1,2),ylim, 'r')    %绘制延迟到的时刻
title('输出信号');
grid on;  xlabel('时间/s')