%OOK,2PSK,2FSK,文件名binarymod.m
clear all;
close all;

A=1;
fc=2;                    %2Hz
N_sample=8;
N=500;                   %碼元數
Ts=1;                    %1 Baud/s
T=1;

dt=Ts/fc/N_sample;       %波形采樣間隔
t=0:dt:N*Ts-dt;
Lt=length(t);

%產生二進制源
d=sign(randn(1,N));
dd=sigexpand((d+1)/2,fc*N_sample);
gt=ones(1,fc*N_sample);  %NRZ波形

figure(1)
subplot(221);%輸入NRZ信號波形（單極性）
d_NRZ=conv(dd,gt);
plot(t,d_NRZ(1:length(t)));
axis([0 10 0 1.2]);   ylabel('輸入信號');

subplot(222);%輸入NRZ頻譜
[f,d_NRZf]=T2F(t,d_NRZ(1:length(t)));
plot(f,10*log10(abs(d_NRZf).^2/T));
axis([-2 2 -50 50]);ylabel('輸入信號功率譜密度（dB/Hz）');
ht=A*cos(2*pi*fc*t);
s_2ask=d_NRZ(1:Lt).*ht;
subplot(223)
plot(t,s_2ask);
axis([0 10 -1.2 1.2]);ylabel('OOK');

[f,s_2askf]=T2F(t,s_2ask);
subplot(224)
plot(f,10*log10(abs(s_2askf).^2/T));
axis([-fc-4 fc+4 -50 50]);ylabel('OOK功率譜密度（dB/Hz）');

figure(2)
%2PSK信號
d_2psk=2*d_NRZ-1;
s_2psk=d_2psk(1:Lt).*ht;
subplot(221)
plot(t,s_2psk);
axis([0 10 -1.2 1.2]);ylabel('2PSK');

subplot(222)
[f,s_2pskf]=T2F(t,s_2psk);
plot(f,10*log10(abs(s_2pskf).^2/T));
axis([-fc-4 fc+4 -50 50]);ylabel('2PSK功率譜密度（dB/Hz）');

%2FSK
%s_2fsk=Acos(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')
subplot(224)
[f,s_2fskf]=T2F(t,s_2fsk);
plot(f,10*log10(abs(s_2fskf).^2/T));
axis([-fc-4 fc+4 -50 50]);xlabel('f');ylabel('2FSK功率譜密度(dB/Hz)');