·¢²¼Ê±¼ä : ÐÇÆÚÈý ÎÄÕ»ùÓÚlmsË㷨ʱÓò¾ùºâÆ÷µÄÉè¼Æ£¨ÖÕ£© ±¾¿Æ±ÏÒµÂÛÎĸüÐÂÍê±Ï¿ªÊ¼ÔĶÁ
н®´óѧ¿Æѧ¼¼ÊõѧԺ±ÏÒµÂÛÎÄ£¨Éè¼Æ£©
plot(s); %ÐźÅsʱÓò²¨ÐÎ title('ÐźÅsʱÓò²¨ÐÎ'); xlabel('n');
axis([0,N,-a-1,a+1]);
xn=awgn(s,5); %¼ÓÈë¾ùֵΪÁãµÄ¸ß˹°×ÔëÉù %ÉèÖóõÖµ
y=zeros(1,N); %Êä³öÐźÅy
y(1:k)=xn(1:k); %½«ÊäÈëÐźÅxnµÄÇ°k¸öÖµ×÷ΪÊä³öyµÄÇ°k¸öÖµ w=zeros(1,k); %ÉèÖóéÍ·¼ÓȨ³õÖµ e=zeros(1,N); %Îó²îÐźŠ%ÓÃLMSËã·¨µü´úÂ˲¨ for i=(k+1):N
XN=xn((i-k+1):(i)); y(i)=w*XN'; e(i)=s(i)-y(i); w=w+u*e(i)*XN;
end
pp(q,:)=(e((k+1):N)).^2; end subplot(312)
plot(xn); %ÐźÅsʱÓò²¨ÐÎ title('ÐźÅs¼ÓÔëÉùºóµÄʱÓò²¨ÐÎ'); subplot(313)
plot(y); %ÐźÅsʱÓò²¨ÐÎ title('×ÔÊÊÓ¦Â˲¨ºóµÄÊä³öʱÓò²¨ÐÎ'); for b=1:N-k
bi(b)=sum(pp(:,b))/g; %ÇóÎó²îͳ¼Æƽ¾ù end
figure(2); %Ëã·¨ÊÕÁ²ÇúÏß
13
н®´óѧ¿Æѧ¼¼ÊõѧԺ±ÏÒµÂÛÎÄ£¨Éè¼Æ£©
t=1:N-k; plot(t,bi,'r'); title('ÊÕÁ²ÇúÏß');
hold on %½«Ã¿´ÎÑ»·µÄͼÐÎÏÔʾ½á¹û±£´æÏÂÀ´
ÐźÅsʱÓò²¨ÐÎ20-20100200300500600700nÐźÅs¼ÓÔëÉùºóµÄʱÓò²¨ÐÎ400800900100050-501002003004005006007008009001000×ÔÊÊÓ¦Â˲¨ºóµÄÊä³öʱÓò²¨ÐÎ20-201002003004005006007008009001000
ͼ4.1 ¦Ì=0.00026ϵͳÊäÈë¡¢ÔëÉù¡¢Êä³öͼ
ÊÕÁ²ÇúÏß10.5001002003004005006007008009001000ͼ4.1 ¦Ì=0.00026ÊÕÁ²ÇúÏß
14
н®´óѧ¿Æѧ¼¼ÊõѧԺ±ÏÒµÂÛÎÄ£¨Éè¼Æ£©
ÐźÅsʱÓò²¨ÐÎ20-202004006008001000nÐźÅs¼ÓÔëÉùºóµÄʱÓò²¨ÐÎ50-502004006008001000×ÔÊÊÓ¦Â˲¨ºóµÄÊä³öʱÓò²¨ÐÎ50-502004006008001000
ͼ4.2 ¦Ì=0.00058ϵͳÊäÈë¡¢ÔëÉù¡¢Êä³öͼ
ÊÕÁ²ÇúÏß10.5002004006008001000
ͼ4.2 ¦Ì=0.00058ÊÕÁ²ÇúÏß
15
н®´óѧ¿Æѧ¼¼ÊõѧԺ±ÏÒµÂÛÎÄ£¨Éè¼Æ£©
ÐźÅsʱÓò²¨ÐÎ20-202004006008001000nÐźÅs¼ÓÔëÉùºóµÄʱÓò²¨ÐÎ50-502004006008001000×ÔÊÊÓ¦Â˲¨ºóµÄÊä³öʱÓò²¨ÐÎ50-502004006008001000
ͼ4.3¦Ì=0.00116ϵͳÊäÈë¡¢ÔëÉù¡¢Êä³öͼ
ÊÕÁ²ÇúÏß10.5002004006008001000ͼ4.3 ¦Ì=0.00116ÊÕÁ²ÇúÏß
16