% 主。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相干檢測）的誤碼率理論值