clear,clc;
L=10000;
t=0:0.01*pi:2*pi;
SNR=-4:0.5:10;
snr=SNR/10;
for ii=1:length(snr)
Sigma(ii)=1/sqrt(2*(10^snr(ii)));%标准差
psk_error1=0; fsk_error1=0; ask_error1=0; %误码数初始值为0
psk_error2=0; fsk_error2=0; ask_error2=0;
S=randi(L);
for k=1:L
n=Sigma(ii)*randn(1,length(t)); %一定信噪比下的高斯白噪声
for i=1:length(t)
psk1(i)=sin(t(i)); psk2(i)=-sin(t(i)); %2PSK的S0(t)和S1(t)
fsk1(i)=sin(3*t(i)); fsk2(i)=sin(5*t(i)); %2FSK的S0(t)和S1(t)
ask1(i)=sin(t(i)); ask2(i)=0; %2ASK的S0(t)和S1(t)
end
%------------------2PSK
overlay11=sum(psk1.^2); overlay12=sum(psk2.^2);
psk1=psk1/sqrt(overlay11); psk2=psk2/sqrt(overlay12);%进行归一化
xp1=psk1+n; %接收信号(含有噪声的)
xp2=psk2+n