c5_1_1.m

来自「包含了大量的DSP应用子程序」· M 代码 · 共 13 行

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));