hc561.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第四章例5.6.1程序hc561
% 混合信号通过一阶数字滤波器的输出
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
% 程序 5.6.1
%
t=0:0.005:3; x1=sin(10*t);  % 原始低频信号x1
x=x1+sin(300*t);            % 低频信号x1加高频干扰
b=[1 1];a=[20 -18];         % 滤波器系数
y=filter(b,a,x);            % 求滤波器输出y
subplot(1,2,1),plot(t,x)      % 画原始合成信号
title('低频信号加高频干扰')
subplot(1,2,2),plot(t,y,t,x1,':') % 画原低频及输出信号
title('原低频信号及滤波输出')
[H,w]=freqz(b,a,[0.05,0.15,1.5]);    % 求出三个频点的频率特性
abs(H),angle(H)*180/pi          % 显示幅度和以度为单位的相角
set(gcf,'color','w');                % 置图形背景色为白