用MATLAB產生FSK調制信號
1．實驗要求：

（1）隨機產生一個包含十個元素的數組，該數組中的元素非0即1，用作FSK信號的輸入數據。

（2）對該數組中的二進制數進行FSK調制，0調制在5000Hz上，1調制在8000Hz上，碼元速率為200碼元/秒。

（3）畫出所產生的FSK信號的時域圖形和頻域圖形。

2．實驗源程序

%隨機產生一個包含十個元素的數組，該數組中的元素非0即1，用作FSK信號的輸入數據。

clear all;

x=randn(10,1)>0;           %Randn正態分布的隨機數矩陣

fl=5000;fh=8000;fs=100000;

ts=1/200;                  %碼元速率200波特

tt=(0:1/fs:ts);

t=[tt;tt+ts;tt+2*ts;tt+3*ts;tt+4*ts;tt+5*ts;tt+6*ts;tt+7*ts;tt+8*ts;tt+9*ts];

y=zeros(10,length(tt)); %Zeros  零矩陣

i=1;

%對該輸入信號FSK調制

while i<=10

y(i,:)=x(i)*sin(2*pi*fh*t(i,:))+~x(i)*sin(2*pi*fl*t(i,:));

i=i+1;

end

t=reshape(t',length(tt)*10,1);

y=reshape(y',length(tt)*10,1);

subplot(211);plot(t,y);

title('FSK信號的時域圖形');

%該輸入信號的頻域圖形

n=length(y);

r=fft(y)/n;r=fftshift(r);

f=linspace(-fs/2,fs/2,n);

subplot(212);

plot(f,abs(r));

set(gca,'XTick',-fs/2:5000:fs/2);

title('FSK信號的頻譜圖');