sys_comp.m

2psk的调制以及误码率等系统性能比较

% 主程序。典型系统、检测方法的AWGN误码率比较。
clear;

% 公共参数、数据、波形
N_sys=3;                % 系统的数量
SW=[1,1,1];             % 系统开关
num_dtct=50000;            % 检测总次数   
cnt_Err=zeros(1,N_sys); % 错误计数
cnt_dtct=0;             % 检测计数

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

% 各系统所使用的波形
% （1）系统I（2PSK，相干检测）：
g1=sin_wav(T,1,0,M);    % 正弦波。 E归一。参数：持续期，谐波次数，初相位，采样点数。

% 发送、传输、接收、检测。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相干检测）的误码率理论值