c5_1_1.m

dsp入门与实践一书的源代码

clear all;                                   
clf;                                         
signal = sin(2*pi*0.055*[0:1000-1]');		
noise=randn(1,1000);				
nfilt=fir1(11,0.4); 				
fnoise=filter(nfilt,1,noise); 		     
d=signal.'+fnoise;  				
w0 = nfilt.' -0.01; 				
mu = 0.05; 					
s=initse(w0,mu); 				
[y,e,s] = adaptse(noise,d,s);
plot(0:199,signal(1:200),0:199,e(1:200));