q612b.m

matlab在电子信息中的应用和实例中的源代码

%《MATLAB及在电子信息课程中的应用》第六章例612程序q612b
% 方波经过三阶巴特沃斯低通滤波器后波形
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
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	% 画出滤波后的频谱
set(gcf,'color','w')                % 设置图形背景色为白色
y = Y*exp(j*w'*t)/pi*dw; 			% 对中段频谱求付利叶逆变换
subplot(1,2,2)
plot(t,f,t,y,'linewidth',1.5),grid	% 画出原波形及滤波后的波形