bwyang3.m

完整版的频谱检测算法的MATLAB实现

clear
clc
N=500;
M=300;
u=[0.1 0.05 0.02];
n=1:N;
a=[-0.99 0.99]
k=1;
for l=1:3
for j=1:M
vn=randn(1,N);
a1=[1 a(k)];
b1=1;
xn=filter(b1,a1,vn);
d=var(xn);
xn=xn/sqrt(d);
fn(1)=xn(1);
wn(1)=0;
wn(2)=wn(1);
for i=2:length(xn)
fn(i)=xn(i)-wn(i)*xn(i-1);
wn(i+1)=wn(i)+u(l)*xn(i-1)*fn(i);
f2(i)=fn(i)^2;
end;
fm(j,:)=f2;
end;
fw(l,:)=mean(fm);
end;
semilogy(n,fw(1,:),n,fw(2,:),n,fw(3,:));