sim_2psk.asv

关于BPSK的调制和误码的分析。matlab源程序。

% 主。2PSK相干检测系统仿真。AWGN。信噪比。误码率测试。理论误码率。
clear;

% 公共参数、数据、波形
num_dtct=5000;          % 检测总次数   
cnt_Err=0;              % 错误计数
cnt_dtct=0;             % 检测计数

T=1;                    % 符号时宽
M=2*10^3;                 % 符号时宽内采样点数
h=T/M;                  % 采样间隔，步长

% 各系统所使用的波形
% （1）系统I（2PSK，相干检测）：
g1=sin_wav(T,1,0,M);    % 正弦波。 E归一。参数：持续期，谐波次数，初相位，采样点数。
subplot(2,1,1);
plot(g1);
Fg1=fft(g1)/(M/2);      % 谱分析 
AFg1=abs(Fg1);
subplot(2,1,2);
stem(AFg1(1:10));
a1=AFg1(2)

% 发送、传输、接收、检测。num_dtct次。
SNR1=4.32                  % dB。 Eb/N0，S/N。
N0=1/10^(SNR1/10);      % 每谐波成分所具有的功率，W。
Pn=N0*(M/2);
Pn_dB=10*log10(Pn);

for cnt_dtct=1:num_dtct
Di=randsrc(1,1);        % 发送数据, 1行1位伪随机码，码元取值为+1,-1
n=wgn(1,M,Pn_dB);        % 白噪声波形 n(t)，1行，M点，dBW

% （1）系统I（2PSK，相干检测）：
x1=g1;                   % 发送。发送用波形
s1=Di*x1;                % 发送信号
r=s1+n;                  % 接收信号
y1=g1;                   % 检测。相关用波形
R1=projct(r,y1,T,M);      % 接收信号r在y1上的投影
if R1>0
    Dr1=1;
else
    Dr1=-1;
end    
if Dr1~=Di
    cnt_Err(1)=cnt_Err(1)+1;
end

%cnt_dtct
end                   % 一次检测循环的end
cnt_Err;
Pe_test=cnt_Err/num_dtct

sn_r=10^(SNR1/10);             
Pe1=(1/2)*erfc(sqrt(sn_r))    % 系统I（2PSK相干检测）的误码率理论值