clear
clc
t1=cputime;

Es=2;%功率約束值
N=1000000;%點數

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轉化為數值
????src=round（rand（1N））;%產生原始發送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?
????%-------產生BPSK調制信號---------------------
????M1=modem.pskmod（2）;
????s1=modulate（M1src）;
????pow1=norm（s1‘fro‘）^2/N;
????s1=s1*sqrt（Es/pow1）;
????%-------產生QPSK調制信號---------------------
????M2=modem.pskmod（4）;?
????s2=modulate（M2src1）;
????pow2=norm（s2‘fro‘）^2/N*2;
????s2=s2*sqrt（Es/pow2）;
????%-------產生噪聲----------------------------
????noise1=randn（1N）*（sqrt（Es/snr1））;
????noise2=（（randn（1N/2）+sqrt（-1）*randn（1N/2）））*（sqrt（（Es/snr1/2）））;
????%-------產生接收信號-------------------------
????y1=s1+noise1;
????y2=s2+noise2;
????%-------BPSK、QPSK解調----------------------
????M11=modem.pskdemod（2）;
????M22=modem.pskdemod（4）;
????sink