%《MATLAB在電子信息課程中的應用》第六章例612程序q612b
% 方波經過三階巴特沃斯低通濾波器后波形
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

q611g							% 調用上題的程序,但數據已確定
% 截止頻率為wc=10的三階巴特沃斯低通濾波器轉移函數為
% H(s)=500/(s^3+20s^2+200s+1000
H = freqs(1000,[1,20,200,1000],w);
% 如用信號處理工具箱，語句為 
% [b,a]==butter(wn,階數;H=freqs(b,a,w),
Y = H.*F;
figure(3),subplot(1,2,1),
plot(w,abs(Y),'linewidth',1.5),grid	% 畫出濾波后的頻譜
y = Y*exp(j*w'*t)/pi*dw; 			% 對中段頻譜求付利葉逆變換
subplot(1,2,2)
plot(t,f,t,y,'linewidth',1.5),grid	% 畫出原波形及濾波后的波形