samp5_19.m

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

%Samp5_19
N=5;  %滤波器阶数
Wn=1000;  %边界频率
[b,a]=besself(N,Wn); %设计Bessel滤波器
figure(1)
[H,w]=freqs(b,a,512);   %绘出滤波器的复数频率特性
magH=abs(H);phaH=unwrap(angle(H));  %求幅频响应和相频响应
subplot(2,1,1),plot(w/(2*pi),20*log10(magH));  %以频率为横轴绘幅频响应
grid on;xlabel('频率/Hz');ylabel('振幅/dB');
subplot(2,1,2),plot(w/2/pi,angle(H)*180/pi);%绘制相频响应
grid on;xlabel('频率/Hz');ylabel('相位/^o')
figure(2)    %新的绘图框
dt=1/10000;   %信号采样间隔
f1=100;f2=1000;  %输入两个频率成分
t=0:dt:0.1;
x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t); %输入信号
H=[tf(b,a)];   %滤波器在MATLAB系统中的表示
[y,t1]=lsim(H,x,t);  %模拟滤波器的输出
subplot(2,1,1),plot(t,x),title('输入信号')  %绘出输入信号
subplot(2,1,2),plot(t1,y)  %绘出输出信号
title('输出信号'),xlabel('时间/s')