function adaptive1(f) 
length=2000; 
p=100;  %FIR滤波器的阶数p
Ts=0.0001; 
noise=randn(1,length); %生成1*N的随机矩阵，当做宽带信号噪声
k=1:length; 
F=cos(2*pi*f*k*Ts); 
X=noise+F;%输入滤波器的信号
u=0.0001; %μ的取值，u越小收敛越慢效果越好，u较大容易出现震荡现象
Y=zeros(1,length);     %滤波结果矩阵
error=zeros(1,length); %误差函数
Weights=zeros(1,p); %权重函数for n=p:length x=X(n:-1:n-p+1); %滤波器输入值Y(n)=x*Weights';%转置error(n)=F(n)-Y(n); %误差Weights=Weights+2*u.*error(n).*x; %迭代更新w
endsubplot(4,1,1);
plot(F); 
title( '初始信号  F'); 
subplot(4,1,2); 
plot(X); 
title( '混合信号波形  X（X=noise+F）'); 
subplot(4,1,3); 
plot(Y); 
title( '滤波结果 y')
subplot(4,1,4); 
plot(error); 
title( '提取的单频信号波形  F 的误差 ');