bwyang2.m

来自「完整版的频谱检测算法的MATLAB实现」· M 代码 · 共 27 行

clear
clc
N=500;
M=300;
u=0.2;
n=1:N;
a=[0.98 -0.98]
k=1;
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*xn(i-1)*fn(i);
f2(i)=fn(i)^2;
end;
fm(j,:)=f2;
end;
fw=mean(fm);
semilogy(n,f2,n,fw);