cdma的MATLAB仿真源程序

No = 2;

Eb = No * Eb_No; êlculate power p Tc = 1; Ts = N * Tc; p = Eb / Ts;

%generate BPSK symbols randomly with value +1 or -1 x = bingen(x_num);

%DS-SS modulate symbols with user code c = bingen(N); y = ds_mod(c(:),x);

%scale by appropriate power factor y = sqrt(p)*y;

-d Pulsed Noise Jammer to signal

y = [awgn(y(1:ro*length(y)),1/ro); y((ro*length(y)+1):length(y))];

%DS-SS demodulate symbols with user code x_de = ds_demod(c(:),y);

Tcision

x_de(find(x_de < 0)) = -1;

x_de(find(x_de >=0)) = 1;

Pe = length(find(x - x_de))/x_num; Plot_Pe = [Plot_Pe Pe]; end %end for EbNo

%--------------------------------------------- %return;

%---------------------------------------------

%display the calculated Pd and Pfa Plot_Pe

%plot Pe versus Eb/No %subplot(2,1,1)

semilogy(plot_EbNo,Plot_Pe,'m*-') xlabel('10log_{10}[(P/J)(W/R)] (dB)') ylabel('BER')

s=sprintf('BER versus 10log_{10}[(P/J)(W/R)] in pulsed noise jamming and AWGN'); title(s);

%********************************************************************* % This program computes the average BER of a DS-SS/BPSK %

communication system in barrage noise jamming and AWGN

%

% AUTHOR: Wenbin Luo

% DATE : 05/02/01 % % %

%********************************************************************

% function Plot_Pe = final12_extra() clear all; %close all;

format long;

%set up the threshold Vt Vt = 0;

Plot_Pe = []; N = 16; x_num = 10000;

%---------------------------------------------------------- %convert back from dB Eb_No = 5; % 2, 4, 6 dB Eb_No = 10.^(Eb_No/10); %assume No = 2;

No = 2;

Eb = No * Eb_No; êlculate power p Tc = 1; Ts = N * Tc; p = Eb / Ts;

%----------------------------------------------------------

plot_EbNj = 0:2:50; for EbNj = 0:2:50,

%convert back from dB Eb_Nj = EbNj; ? Eb_Nj = 10.^(Eb_Nj/10); Nj = Eb / Eb_Nj;

%generate BPSK symbols randomly with value +1 or -1 x = bingen(x_num);

%DS-SS modulate symbols with user code c = bingen(N); y = ds_mod(c(:),x);

%scale by appropriate power factor y = sqrt(p)*y;

-d barrage noise jamming and AWGN to signal y = awgn(y,(Nj/2)+1); %y = awgn(y,1);

%DS-SS demodulate symbols with user code x_de = ds_demod(c(:),y);

Tcision

x_de(find(x_de < 0)) = -1; x_de(find(x_de >=0)) = 1;

Pe = length(find(x - x_de))/x_num; Plot_Pe = [Plot_Pe Pe]; end %end for EbNo

%--------------------------------------------- %return;

%---------------------------------------------

%display the calculated Pd and Pfa Plot_Pe

%plot Pe versus Eb/No %subplot(2,1,1)