%《數字信號處理教程——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');                % 置圖形背景色為白