rands.m

正弦干扰陷波LMS算法

function [s,pv]=rands(A,w0)

%产生正弦信号s(n)
%A为幅值 
%w0属于(0,pi)

for n=1:100
    s(n)=A*sin(w0*n);
end

%ps为s(n)的方差
for n=1:100
    p(n)=s(n)*s(n);
end
ps=sum(p(1:100))/100;

%条件ps=10pv
pv=ps/10;