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) -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)