function [y]=zlb_fskcode(x)
%x=QPSK_rc;
fl=10;%頻率1
fh=20;%頻率2
fd=500;%輸出抽樣頻率（10的7次方以上無法計算）
ts=1/10;%碼元速率
tt=(0:1/fd:ts);%每個碼元用50個點來描述
len_x=length(x);
tr=zeros(len_x,length(tt));
for  m=0:len_x-1
   % tr(m*length(tt)+1:(m+1)*length(tt))=[tt+m*ts];
   tr(m+1,:)=[tt+m*ts];
end
y=zeros(len_x,length(tt));
i=1;
%對該輸入信號FSK調制
while i<=len_x
yr(i,:)=x(i)*sin(2*pi*fh*tr(i,:))+~x(i)*sin(2*pi*fl*tr(i,:));
i=i+1;
end

figure
subplot(211)
tx=(ts:ts:2);%取2s之內的10個輸入序列畫圖（時間范圍可以任意選取）
stem(tx,x(1:length(tx)))
title('用于fsk調制的二進制信道編碼序列');
t=reshape(tr',length(tt)*len_x,1);%reshape():把t'重組為length(t)*10行1列的數組
y=reshape(yr',length(tt)*len_x,1);
subplot(212);
ty=(1/fd:1/fd:2);
plot(ty,y(1:length(ty)))
title('FSK信號的時域圖形');