求解以下matlab程序！

这样改：

clc

clear

SNR_DB=0:1:12;

sum=1000000;

data= randsrc(sum,2,[0 1]);

[a1,b1]=find(data(:,1)==0&data(:,2)==0);

message(a1)=-1-1i;

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

message(a2)=-1+1i;

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

message(a3)=1-1i;

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

message(a4)=1+1i;

scatterplot(message)

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

A=1;

Tb=1;

Eb=A*A*Tb;

P_signal=Eb/Tb;

NO=Eb./(10.^(SNR_DB/10));

P_noise=P_signal*NO;

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*1i;

%resum=0;

%total=0;

m1=find(angle(receive)0);

remessage(1,m1)=1+1i;

redata(m1,1)=1;

redata(m1,2)=1;

m2= find( angle(receive)>pi/2&angle(receive)-pi&angle(receive)-pi/2&angle(receive)

百科狗 baikegou.com