bwyang1.m

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

clear
clc
N=500;
M=300;
u=0.1;
n=1:(N+1);
a=[-0.98 0.98]
for k=1:2
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);
end;
wm(j,:)=wn;
end;
w(k,:)=wn;
ew(k,:)=mean(wm);
end;
plot(n,w(1,:),n,ew(1,:),n,w(2,:),n,ew(2,:));