clear
clc
t1=cputime;

Es=2;%功率約束值
N=1000000;%點(diǎn)數(shù)

src=zeros（1N）;
src1=zeros（1N/2）;
%sink1=zeros（1N）;
%sink2=zeros（1N）;
%sink22=zeros（1N）;%初始化

clear?ij;
SNR0=0;
SNR1=14;
for?j=SNR0:SNR1
????%error1=0;
????%error2=0;
????%error3=0;
????snr=j;
????snr1=10^（snr/10）;%將信噪比的值由dB轉(zhuǎn)化為數(shù)值
????src=round（rand（1N））;%產(chǎn)生原始發(fā)送0、1比特
????%--------將原始0、1比特映向射為0、1、2、3----
????k=1;
????for?i=1:2:N-1
????????src1（k）=2*src（i）+src（i+1）;
????????k=k+1;
????end?
????%-------產(chǎn)生BPSK調(diào)制信號(hào)---------------------
????M1=modem.pskmod（2）;
????s1=modulate（M1src）;
????pow1=norm（s1‘fro‘）^2/N;
????s1=s1*sqrt（Es/pow1）;
????%-------產(chǎn)生QPSK調(diào)制信號(hào)---------------------
????M2=modem.pskmod（4）;?
????s2=modulate（M2src1）;
????pow2=norm（s2‘fro‘）^2/N*2;
????s2=s2*sqrt（Es/pow2）;
????%-------產(chǎn)生噪聲----------------------------
????noise1=randn（1N）*（sqrt（Es/snr1））;
????noise2=（（randn（1N/2）+sqrt（-1）*randn（1N/2）））*（sqrt（（Es/snr1/2）））;
????%-------產(chǎn)生接收信號(hào)-------------------------
????y1=s1+noise1;
????y2=s2+noise2;
????%-------BPSK、QPSK解調(diào)----------------------
????M11=modem.pskdemod（2）;
????M22=modem.pskdemod（4）;
????sink