mini2

  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)