exm046_3.m

北京航空航天大学出版社张志涌、徐彦琴主编的matlab教程附带的mfile例程压缩包

%exm046_3.m
clear,randn('state',1)
ws=1000;					 %采样频率
t=0:1/ws:0.4;	
x=sin(2*pi*10*t)+cos(2*pi*100*t)+0.2*randn(size(t));	    %生成带噪声的多频率信号
wn=ws/2;					 %Nyquest频率
[B,A]=butter(10,30/wn);	%截止频率为30/wn的10阶ButterWorth低通滤波器
y=filter(B,A,x);			  %进行（初值为0的）滤波处理
clf,shg
plot(t,x,'b-',t,y,'r.','MarkerSize',10)
legend('Input','Output',0)