Modulation CPFSK & MSK

%------ Modulation CPFSK & MSK-----

clear;clc;close all;

%---Parametres de simulation------

fe = 200e3;

dt = 1/fe;

Nb = 200;

Db = 5e3;

T = 1/Db;

N=round(T/dt);

modulation = 'cpfsk';

h = 1/2;

fp = 10e3;

df = 100;

f = -10e3 : df : 10e3 ;

t = 0 : dt : Nb*T -dt ;

%-------1)modulateur à phase -----

%---point A: données binaires-----

data=rand(1,Nb)> 0.5 ;

subplot(311);

stem(data);grid on;

xlim([0 20]);

title('donnees binaire');

%-----point B : mapping -----

symbole = 2*data - 1 ;

subplot(312);

stem(symbole); grid on ;

xlim([0 20]);

title('mapping');

%------ sur echantillonnage-----

data_up=upsample(symbole,N);

%----point C : mise en forme rectangulaire ---

C = ones(1,N);

sig_C = filter (C,1,data_up);

subplot(313);

plot(t/T,sig_C); grid on ;

xlim([0 20]);

ylim([-2 2]);

title('mise en forme rect');

%----2)point D : la phase et FM-----

Delta = h/2/T;

y = 0 ;

for k = 2 : length(sig_C)

y(k) = y(k-1) + sig_C(k) ;

phi(k) = 2*pi*Delta*y(k)*dt;

end

figure;

subplot(211);

plot(t/T,phi); grid on ;

xlim([0 20]);

title('FM');

% -----signal modulé ---

sig_D = cos(2*pi*fp*t+phi) ;

subplot(212);

plot(t/T,sig_D); grid on ;

xlim([0 20]);

title('signal modulé');

%-----3)l'enveloppe complexe---

sig_Re=cos(phi);

figure;

subplot(211);

plot(t/T,sig_Re); grid on ;

xlim([0 20]);

title('envelope complexe :Re');

sig_Im=sin(phi);

subplot(212);

plot(t/T,sig_Im); grid on ;

xlim([0 20]);

title('envelope complexe :Im');

%-----4)constellation instantannée---

scatterplot(sig_Re+j*sig_Im,1,0,'b-');grid on; hold on;

title('constellation');

...