clc
a=round(rand(1,1000))%生成消息码
f=2;%假设载波频率为2
for SNR=0:15%给出信噪比范围,求误比特率
cp=[];%存储BPSK映射后的消息序列
sigma=1/sqrt(2*SNR);%方差
for n=1:num
die=-1;
die=1;
cp=[cp,die];
b=[b,cp+Noise];%*ones(1,10)%传输过程中的加上了噪声后的消息码
%输出消息码
out=[];%检测器输出的消息码
if b(m)>0
else
end
end
%比较检测器输出的消息码和原消息码,并统计错误码元个数
count=0;%码元错误个数
if a(k)~=out(k)
end
P=[P,count/n]
semilogy(x,P)