2PSK.rar_2PSK matlab_2psk_matlab 2psk_信道检测

% main_2psk，主程序。2PSK调制解调（相干检测）系统仿真。AWGN信道。误码率测试。理论误码率。 clear; % 公共参数： num_dtct=5000; % 检测总次数（发送、接收的二进制符号数） cnt_Err=0; % 错误计数 cnt_dtct=0; % 检测计数 T=1; % 符号时宽 M=2*10^3; % 符号时宽内采样点数 h=T/M; % 采样间隔，步长 % 2PSK调制所使用的波形： g1=sub_sin_wav(T,1,0,M); % 正弦波。 能量归一。参数：时宽，谐波次数，初相位，采样点数。 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。 N0=1/10^(SNR1/10); % 每谐波成份所具有的功率，W。噪声功率谱密度。 Pn=N0*(M/2); % 噪声功率。 Pn_dB=10*log10(Pn); % 噪声功率的dBW数。 for cnt_dtct=1:num_dtct % 循环，调制、解调、检测num_dtct次。 Di=randsrc(1,1); % 发送数据, 1行1位伪随机码，码元取值为+1,-1 n=wgn(1,M,Pn_dB); % 白噪声波形 n(t)，1行，M点。 x1=g1; % 发送。发送用波形x1。 s1=Di*x1; % 调制。产生已调2PSK信号。 r=s1+n; % 接收信号 y1=g1; % 接收端相关使用的波形y1。（恢复的载波） R1=sub_projct(r,y1,T,M); % 相关，r,y1相乘、[0，T]内积分。（接收信号r在y1上的投影） if R1>0 % 判决，得到恢复的数据 Dr1。 Dr1=1; else Dr1=-1; end if Dr1~=Di % 错误计数。 cnt_Err=cnt_Err+1; end end % 一个符号检测循环的end % num_dtct个符号检测完毕后，求实际误码率： Pe_test=cnt_Err/num_dtct % 2PSK相干检测的误码率理论值： sn_r=10^(SNR1/10); Pe1=(1/2)*erfc(sqrt(sn_r)) % 结束。