modulator.asv

Matlab 模拟 CDMA 直接扩频。m序列生成矩阵为p(x) = 1 + x^5 + x^7 + x^8 + x^9 + x^13 + x^15。能画出每阶段的时域、频域图。

%% Made by Zheng Xiao (052815 05b04)
%% Group Member: Ji Yue, Sun Xin, Zheng Xiao

%clear
% Generate Signal
%M = 16; %-the number of Sample per Period in input signal
N = 32; %-the number of Sample per Period in BPSK Modulated signal
%tk = 0:1/M:(1024/M);
%input = sin(2*pi*tk) + cos(pi*tk);
%ratio = max(input);
%input = input/ratio;
% A-law PCM encode

% Lengthen PCM code
%%pcmCode = [1];
original=[0 1 0];
pcmCode=expand_bit(original);
key=genM(64);
pcmCode=xor(pcmCode,[key,key,key]);


npcmCode = length(pcmCode); %-number of total symbol
pcmCode = [pcmCode;zeros((N-1),npcmCode)];
pcmCode = reshape(pcmCode,1,N*npcmCode);
pcmCode = conv(pcmCode,ones(1,N));
pcmCode = pcmCode(1:N*npcmCode);
pcmCode(find(pcmCode == 0)) = -1;
% BPSK Modulation
tk = 0:1/N:(1-1/N);
carrier = sin(4*pi*tk);
carrier = padarray(carrier,npcmCode-1,'circular','post');
carrier = reshapeSeq(carrier);
signal = carrier .* pcmCode; %<= THIS IS THE MODULATED SIGNAL
plot(signal);

%%%%% Here comes through the channel %%%%%
%=================[Media]================%

% BPSK Demodulation
signalIn = signal .*carrier; %<= INPUT SIGNAL OF DEMODULATOR
 %plot(signalIn)
judger = reshapeColumn(signalIn,N);
judger = sum(judger,2);
judger(find(judger>=0)) = 1;
judger(find(judger<0)) = 0;
judger = judger';

% A-law PCM decode
%output = pcm_decode(judger).*ratio;
%axis([0,102,-ratio,+ratio])
%input = input*ratio;
%err = output - input;
%plot(1:length(output),output,1:length(input),input, 1:length(err),err)