samp7_6.m

MATLAB7.x数字信号处理 光盘内容

%Samp7_6
wp=[0.35 0.65];N=48;    %通带边界频率（归一化频率）和滤波器阶数
Fs=50;
b=fir1(N,wp);           %设计FIR带通滤波器
figure(1)
[H,f]=freqz(b,1,512,Fs);        %以50Hz为采样频率求出滤波器频率响应
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=1;f2=10;f3=20;         %输入信号的三种频率成分
t=0:1/50:3;               %时间序列
x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t)+0.5*sin(2*pi*f3*t);%输入信号
%y=filter(b,1,x);          %可以采用过滤器进行滤波
y=fftfilt(b,x);            %采用 fftfilt 对输入信号滤波
figure(2)
subplot(2,1,1), plot(t,x),title('输入信号')%绘出输入信号波形
subplot(2,1,2),plot(t,y)                  %绘出输出信号波形
hold on;plot(N/2*0.02*ones(1,2),ylim, 'r')    %绘制延迟到的时刻
title('输出信号'),xlabel('时间/s')