cdma_demod_test1.asv

CDMA下行解调的MatLab程序

clear all;
close all;
clc;
%AD采样率为20倍符号速率；
fs=1228800*20;
%读取数据，16量化，
datain=read_binary ('PILOTSIG.dat',50000);
% datain=read_binary ('SKYFSIG.dat',50000);
datain=datain';
%中频采样后剩余载波频率为 fc=3.7280e6Hz
fc=3.7280e6;%-100*4;
%fc=3.72755e6;%锁相环跟踪后的真实的载波频率;
length_wave=length(datain);
t=[0:length_wave-1]/fs;
phi0=0;%初始相位
lo_cos=cos(2*pi*fc*t+phi0);
lo_sin=sin(2*pi*fc*t+phi0);
hi_data=datain.*lo_cos';	% down converting
hq_data=datain.*lo_sin';	% down converting

%低通滤波器，
b=fir1(30,0.7e6/fs);
ii_data=filter(b,1,hi_data);
qq_data=filter(b,1,hq_data);

%采样频率将为4倍符号速率
fs=fs/5;
i_data=ii_data(100:5:end)';
q_data=qq_data(100:5:end)';
clear datain lo_cos lo_sin t ii_data qq_data hi_data hq_data;

offset_n=3; % 1~4之间
n_ratio=4;

N1=length(i_data)/n_ratio-10;
for kk=1:N1
    baseband_data(kk)=i_data(offset_n+n_ratio*kk)+1i*q_data(offset_n+n_ratio*kk);
end;
figure,
plot(baseband_data,'*');
I_data=real(baseband_data(500:end));
Q_data=imag(baseband_data(500:end));

%读取短PN序列，该序列已仿真产生并验证
load PNDATA_align.mat M_PN_I M_PN_Q;

%任意选64个符号数据做滑动相关
[PeakValue1, PeakPos1 ,corrvalue1] = Pilot_PN_Search(I_data(1:64),Q_data(1:64),M_PN_I,M_PN_Q);
figure,
plot(corrvalue1);