% Matlab version: 7.0.4.365 (R14) Service Pack 2
% Date: 2009/01/15
% Rs = 1000; % rate is 1kbps
% fc = 4000; % carrier is 4k Hz
% '1' employs carrier
% '0' no carrier
% t0 = 1; % time duration is 1 second
N = 1000; % Sample number, equal to rate
fc = 4000; % Carrier frequency
fs = 100; % Smapling rate per bit
ts = 1/fs; % Sampling interval
t0 = 1; % Time duration
tc = t0/N; % Each bit lasts 1 millisecond
t = [0:tc:t0-tc]; % Time axis of message signal
ct = [0:ts/N:tc-ts/N]; % Time axis of one bit duration
ca = sin(2*pi*fc*ct); % Carrier signal
test = ca.^2;
limitation = sum(test);
snr = -5; % SNR in dB
bercnt=0; % bit error counter
% Message signal
xn = round(rand(1,1000)); % 1000