»ùÓÚlmsË㷨ʱÓò¾ùºâÆ÷µÄÉè¼Æ£¨ÖÕ£© ±¾¿Æ±ÏÒµÂÛÎÄ

·¢²¼Ê±¼ä : ÐÇÆÚÈý ÎÄÕ»ùÓÚ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

ÁªÏµºÏͬ·¶ÎÄ¿Í·þ£ºxxxxx#qq.com(#Ì滻Ϊ@)