zlb_fskcode.m

用数字调制技术中的fsk调制方法

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信号的时域图形');