ͨÐŹ¤³ÌʵÑé½ÌѧÖÐÐÄ Í¨ÐÅϵͳÔÀíʵÑ鱨¸æ
z1=z+noise; subplot(5,2,5) plot(z1);
title('¼ÓÈëÔëÉùºóÐźŲ¨ÐÎ');
%¶Ô¼ÓÔëÐźŽøÐÐÆµÆ×·ÖÎö ba=fft(z1,512); ba=abs(ba); subplot(5,2,6) plot(ba);
title('¼ÓÔëÐÅºÅÆµÆ×') axis([-200 600 0 150])
%Éè¼ÆÒ»¸öµÍͨÂ˲¨Æ÷ Wp =50/80; Ws = 70/80;
[n,Wn] = buttord(Wp,Ws,1,5) [b,a] = butter(n,Wn);
%¶Ô¼ÓÈëÔëÉùµÄÐźŽøÐÐÂ˲¨ x_fir=filter(b,1,z1);
%¹Û²ìÂ˲¨Ö®ºóµÄÐźŲ¨ÐÎ subplot(5,2,7) plot(x_fir);
title('Â˲¨Ö®ºóµÄÐźÅ') axis([0 3000 -2 2])
%Ïà¸É½âµ÷ x2=x_fir.*y; subplot(5,2,8); plot(x2);
title('ÓëÏà¸ÉÔØ²¨Ïà³Ë²¨ÐÎ') axis([0 3000 -0.5 2])
%¶Ô¼ÓÈëÔëÉùµÄÐźŽøÐÐÂ˲¨ x3=filter(b,1,x2); subplot(5,2,9) plot(x3);
title('ÓëÏà¸ÉÔØ²¨Ïà³ËºóÂ˲¨²¨ÐÎ ') axis([0 3000 -0.5 2.5])
%³éÑùÅоö for i=1:N
if abs(x3(i*100-20))>=0.5;
xn2(i*100-99:i*100)=ones(1,100); else
xn2(i*100-99:i*100)=zeros(1,100); end end
subplot(5,2,10);
plot(xn2);
title('»Ö¸´²¨ÐÎ') axis([0 3000 -1 2])
11