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

q611g							% 調(diào)用上題的程序,但數(shù)據(jù)已確定
% 截止頻率為wc=10的三階巴特沃斯低通濾波器轉(zhuǎn)移函數(shù)為
% H(s)=500/(s^3+20s^2+200s+1000
H = freqs(1000,[1,20,200,1000],w);
% 如用信號處理工具箱，語句為 
% [b,a]==butter(wn,階數(shù);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	% 畫出原波形及濾波后的波形