Introduction - If you have any usage issues, please Google them yourself
clear all
clc
t=0:1/1000:10-1/1000
s=sin(2*pi*t)
snr=20
s_power=var(s) varience of s
linear_snr=10^(snr/10)
factor=sqrt(s_power/linear_snr)
noise=randn(1,length(s))*factor
x=s+noise Ó É SNR¼ Æ Ë ã Ë æ » úÔ ë É ù
x1=noise Ô ë É ùÔ ´ Ê ä È ë
x2=noise
w1=0 È ¨Ï µ Ê ý ³ õ Ö µ
w2=0
e=zeros(1,length(x))
y=0
u=0.05
for i=1:10000 LMSË ã ·¨
y(i)=w1*x1(i)+w2*x2(i)
e(i)=x(i)-y(i)
w1=w1+u*e(i)*x1(i)
w2=w2+u*e(i)*x2(i)
end
figure(1)
subplot(4,1,1)