ʵÑé3 FSK (ASK)µ÷ÖÆ½âµ÷ʵÑé

ͨÐŹ¤³ÌʵÑé½ÌѧÖÐÐÄ Í¨ÐÅϵͳԭÀíʵÑ鱨¸æ

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

ÁªÏµ¿Í·þ£º779662525#qq.com(#Ìæ»»Îª@) ËÕICP±¸20003344ºÅ-4