SNR_DB=[0:1:12]; %Signal-to-noise ratio gradually improve

sum=1000000;

data= randsrc(sum,2,[0 1]); %等概产生1000000*2的矩阵：0和1

[a1,b1]=find(data(:,1)==0&data(:,2)==0);  %找到第一列和第二列为0

message(a1)=-1-j; % map [ 0 0] to 225°

[a2,b2]=find(data(:,1)==0&data(:,2)==1);

message(a2)=-1+j; % map [ 0 1] to 135°

[a3,b3]=find(data(:,1)==1&data(:,2)==0);

message(a3)=1-j; % map [ 1 0] to 275°

[a4,b4]=find(data(:,1)==1&data(:,2)==1);

message(a4)=1+j;% map [ 0 0] to 45°

scatterplot(message)

title('B点信号的星座图')

A=1;

Tb=1;

Eb=A*A*Tb;%带宽w=1hz

P_signal=Eb/Tb;

NO=Eb./(10.^(SNR_DB/10)); %SNR_DB=10.*log10(Eb./NO)

P_noise=NO; %noise power 单边功率谱密度(N0)主要用在复数信号中，双边功率谱密度(N0/2)主要用在实信号中。

sigma=sqrt(P_noise);           

for Eb_NO_id=1:length(sigma) 

    noise1=sigma(Eb_NO_id)*randn(1,sum);  

    noise2=sigma(Eb_NO_id)*randn(1,sum); 

    receive=message+noise1+noise2*j; %how to add gaussian noise to the original signal

倒数第三行为什么是乘号而不是加号