OOK FSK PSK程序 PSK, FSK, OOK, 程序

xx7821096

lear all;
close all;
A=1;
fc=2;%2Hz
N_sample=8;
N=500;           %码元数
Ts=1;            %1Baud/s  码元宽度
dt=Ts/fc/N_sample;%波形采样间隔；
t=0:dt:N*Ts-dt;
Lt=length(t);
%产生2进制信源
d=sign(randn(1,N));
dd= sigexpand((d+1)/2,fc*N_sample);
gt=ones(1,fc*N_sample);
figure(1)
subplot(221);%NRZ波形
d_NRZ=conv(dd,gt);
plot(t,d_NRZ(1:length(t)));
axis([0 10 0 1.2]);
title('输入信号');
subplot(222)  %输入NRZ频谱
[f,d_NRZf]=T2F(t,d_NRZ(1:length(t)));
plot(f,10*log10(abs(d_NRZf).^2/Ts));
axis([-2 2 -50 50]);
title('输入信号功率谱密度（dB/Hz）');
%OOK信号
ht=A*cos(2*pi*fc*t);
s_2ask=d_NRZ(1t).*ht;
subplot(223);
plot(t,s_2ask);
axis([0 10 -1.2 1.2]);title('OOK信号');
%OOK信号的功率谱密度
[f,s_2askf]=T2F(t,s_2ask);
subplot(224);
plot(f,10*log10(abs(s_2askf).^2/Ts));
axis([-fc-4 fc+4 -50 50]);
title('OOK功率谱密度（dB/Hz）');
figure(2)
%2PSK信号
d_2psk=2*d_NRZ-1;
s_2psk=d_2psk(1t).*ht;
subplot(221)
plot(t,s_2psk);
axis([0 10 -1.2 1.2]);
title('2PSK信号');
%2PSK功率谱密度；
subplot(222);
[f,s_2pskf]=T2F(t,s_2psk);
plot(f,10*log10(abs(s_2pskf).^2/Ts));
axis([-fc-4 fc+4 -50 50]);
title('2PSK功率谱密度（dB/Hz）');
%2FSK信号
%s_2fsk=A*cos(2*pi*fc*t+int(2*d_NRZ-1));
sd_2fsk=2*d_NRZ-1;
s_2fsk=A*cos(2*pi*fc*t+...
    2*pi*sd_2fsk(1:length(t)).*t);
subplot(223);
plot(t,s_2fsk);
axis([0 10 -1.2 1.2]);
xlabel('t');
ylabel('2FSK');
title('2FSK信号');
%2FSK功率谱密度
subplot(224);
[f,s_2fskf]=T2F(t,s_2fsk);
plot(f,10*log10(abs(s_2fskf).^2/Ts));
axis([-fc-4 fc+4 -50 50]);
title('2FSK的功率谱密度（dB/Hz）');