做一个最简单的2FSK调制。
废话不多说,直接上代码。
clear;
clf;
close all
timeLimit = 0.1;
sampleNum = 800;
fs = 6000;
t = linspace(0,timeLimit-1/fs,timeLimit*fs);
fcarrier1 = 200;
fcarrier2 = 600;
sourceRate = 100;
sourceBit = randi([0 1],1,timeLimit*sourceRate);
sourceSample = reshape(repmat(sourceBit,fs/sourceRate,1),1,[]);
carrierSapme1 = sin(2*pi*fcarrier1*t);
carrierSapme2 = sin(2*pi*fcarrier2*t);
fskSample = (sourceSample==1).*carrierSapme2+(sourceSample==0).*carrierSapme1;
h=figure;
% set(gcf,'Position',get(0,'ScreenSize'));
set(gcf,'Position',[30 30 1200 300])
for i = 1:4:length(t)-4
plot(t(i:i+4),fskSample(i:i+4),'r','LineWidth',2);
axis([0,timeLimit,-1.2,1.2])
hold on
plot(t(i:i+4),sourceSample(i:i+4),'b','LineWidth',1);
legend("\fontsize{8}2FSK信号","\fontsize{8}原始信号")
title("\fontsize{8}2FSK调制信号,fc1=200Hz, fc2=600Hz")
hold on
drawnow
f=getframe(gcf);
imind=frame2im(f);
[imind,cm] = rgb2ind(imind,256);
if i == 1
imwrite(imind,cm,'FSKModulation.gif','GIF', 'Loopcount',inf,'DelayTime',0.2);
else
imwrite(imind,cm,'FSKModulation.gif','GIF','WriteMode','append','DelayTime',0.2);
end
end
运行之后,效果是这样的: