ÏÖ´úÊý×ÖÐźŴ¦ÀíʵÑ鱨¸æ ÏÂÔØ±¾ÎÄ

hold off;

figure(6)

plot(w/pi,10*log10(PP2/50)); hold on; title('Э·½²î·¨¹¦ÂÊÆ×¹À¼Æ±È½Ï'); ylabel('P(dB)'); xlabel('w/pi');

plot(w/pi,10*log10(P5),'r'); legend('Э·½²î·¨','ÕæÊµÆµÆ×'); grid on; hold off;

figure(7)

plot(w/pi,10*log10(PP3/50)); hold on; title('ÐÞÕýЭ·½²î·¨¹¦ÂÊÆ×¹À¼Æ±È½Ï'); ylabel('P(dB)'); xlabel('w/pi');

plot(w/pi,10*log10(P5),'r'); legend('ÐÞÕýЭ·½²î·¨','ÕæÊµÆµÆ×'); grid on; hold off;

figure(8)

plot(w/pi,10*log10(PP4/50)); hold on; title('Burg·¨¹¦ÂÊÆ×¹À¼Æ±È½Ï'); ylabel('P(dB)'); xlabel('w/pi');

plot(w/pi,10*log10(P5),'r'); legend('Burg·¨','ÕæÊµÆ×'); grid on; hold off;

3¡¢¼ÆËã»úÁ·Ï°3£ºÎ¬ÄÉÔëÉùÒÖÖÆ(Àý6.6µÄÀ©Õ¹ÊµÑé)¡£¼Ù¶¨Í¼6.8ÖÐËùÐèµÄÐźÅd(n)ÊÇÒ»¸öÕýÏÒÐòÁÐd(n)?sin(n?0??),

?0?0.02?, ÔëÉùÐòÁÐv1(n)ºÍv2(n)¶¼ÊÇAR(1) ¹ý³Ì£¬·Ö

v1(n)?0.8v1(n?1)?g(n)

±ðÓÉÈçϵÄÒ»½×²î·Ö·½³Ì²úÉú£º

v2(n)??0.6v2(n?1)?g(n)

ÆäÖÐg(n)ÊÇÁã¾ùÖµ¡¢µ¥Î»·½²îµÄ°×ÔëÉù£¬Óëd(n)²»Ïà¹Ø¡£

(a) ÊÔÓÃMatlab³ÌÐò²úÉúx(n)ºÍv2(n)µÄ500¸öÑùµã£¬»­³ö²¨ÐÎͼ¡£

(b)»ùÓÚx(n)ºÍv2(n)µÄ500¸öÑùµã£¬Éè¼Æp½×µÄ×îÓÅFIRάÄÉÂ˲¨Æ÷£¬ÓÉv2(n)¹À¼Æv1(n)£¬

½ø¶ø¹À¼Æ³öd(n)£¬ÆäÖн×Êýp·Ö±ðȡΪp=3£¬6£¬9£¬12£¬ÊÔ¼ÆËã¸÷ÖÖÇé¿öϹÀ¼Æv1(n)ʱµÄƽ¾ùƽ·½Îó²î(¾ù·½Îó²îµÄÑù±¾¹À¼Æ£¬ÒªÐðÊö¹À¼Æ·½°¸)£¬²¢»­³ö¶Ôd(n)¹À¼ÆµÄ½á¹û¡£ (c)ÓÐʱ¸¨Öú¹Û²âÊý¾ÝÖÐÒ²»á©ÈëһЩd(n)Ðźţ¬¼´¸¨Öú¹Û²âÐźŲ»½öÊÇv2(n)£¬¶øÊÇ

v0(n)?v2(n)??d(n)

ÊÔÕë¶Ôp=12µÄÇé¿ö£¬·Ö±ðÈ¡¼¸¸ö²»Í¬µÄ?Öµ(Èç0.1, 0.5, 1.0)£¬Ñо¿ÕâʱµÄÔëÉùÒÖÖÆÐÔÄÜ¡£ (d)ÈôÖ»ÓÐһ·¹Û²âx(n)?d(n)?v1(n)µÄ1000¸öÑùµã£¬ÄãÄÜÏë°ì·¨½üËÆÍê³É¶ÔÔëÉùv1(n)µÄ

ÓÐЧÒÖÖÆÂð£¿ÊÔ½âÊÍÄãµÄ·½·¨µÄ»ù±¾Ô­Àí£¬ÐðÊöÄãµÄʵÏÖ·½°¸¡£

ͼ6.8 Óи¨¹Û²âÊý¾ÝµÄάÄÉÔëÉùÒÖÖÆÆ÷µÄÔ­Àíͼ

·ÂÕæ½á¹û£º

£¨a£©

ͼ3.1v2(n)µÄ²¨ÐÎ

ͼ3.2 x(n)µÄ500¸öÑùµãµÄ²¨ÐÎ

£¨b£©

»ùÓÚX(n)ºÍV2(n)µÄ500¸öÑùµã£¬¿ÉÒԵõ½

1N?11N?1?v2(k)??V2(n)V2(n?k)¡¢r?xv2(k)??x(n)V2(n?k) rNn?0Nn?0Çó½âWiener-HopfRv2w?rxv2·½³Ì¿ÉÒԵõ½×îÓÅFIRάÄÉÂ˲¨Æ÷¡£

1N?1??1(n)?v1(n))2¼ÆËãµÃµ±p=3¡¢6¡¢9¡¢12ʱ£¬¹À¾ù·½Îó²îµÄÑù±¾¹À¼Æ¿ÉÒÔÓÃ???(vNn?0¼Æv1(n)ʱµÄƽ¾ùƽ·½Îó²î·Ö±ðΪ0.7849¡¢0.2173¡¢0.0747¡¢0.0453¡£

ͼ3.3Â˲¨Æ÷½×Êýp=3ʱµÄ¹À¼ÆÖµÓëÕæÊµÖµ¶Ô±È

ͼ3 .4Â˲¨Æ÷½×Êýp=6ʱµÄ¹À¼ÆÖµÓëÕæÊµÖµ¶Ô±È