N_trials=1000;
N_number=100;
N_snr=10;
Es=1;
BER_m=[];
SER_m=[];
for trials=1:N_trials;
trials
s10=round(rand(1,N_number));
S=(s10*2-1)./sqrt(2);
S1=S(1:2:N_number);
S2=S(2:2:N_number);
Sc=S1+j.*S2;
niose=randn(1,N_number/2)
+j.*randn(1,N_number/2);
BER_v=[];
SER_v=[];
for snr_db=0:1:N_snr;
sgma=(1/2)*sqrt(10.^(-snr_db./10));
Y=Sc+sgma.*niose;
Y_r=sign(real(Y))./sqrt(2);
Y_i=sign(imag(Y))./sqrt(2);
Y_bit=[];
for k=1:length(Y_r);
Y_bit=[Y_bit,[Y_r(k),Y_i(k)]];
end;
Y_symbol=Y_r+j*Y_i;
X_b=S-Y_bit;
X_s=Sc-Y_symbol
ber_snr=0;
ser_snr=0;
for k=1:N_number
if X_b(k)~=0;
ber_snr=ber_snr+1;
end;
end;
for k=1:N_number/2;
if X_s(k)~=0;
ser_snr=ser_snr+1;
end;
end;
BER_v=[BER_v,ber_snr];
SER_v=[SER_v,ser_snr];
end;
BER=
sum(BER_m)./(N_trials*N_number);
SER=
sum(SER_m)./(N_trials*N_number/2);
BER_T=[];
SER_T=[];
for snr_db=0:1:N_snr;
snr=10.^(snr_db./10);
BER_THEORY=Qfunct(sqrt(2.*snr));
SER_THEORY=
1-(1-(1/2).*erfc(sqrt(snr))).^2;
BER_T=[BER_T,BER_THEORY];
SER_T=[SER_T,SER_THEORY];
end;
figure
i=0:1:N_snr;
semilogy(i,BER,'-r',i,SER,'-g',i,BER_T,'*b',i,SER_T,'*y');
legend('BER-simulation','SER-simulation','BER-theory','SER-theory');
xlabel('Eb/N0(db)');
ylabel('BER,SER');
QPSK.rar_qpsk误码率
版权申诉
133 浏览量
2022-09-23
09:17:53
上传
评论
收藏 688B RAR 举报
刘良运
- 粉丝: 68
- 资源: 1万+
最新资源
- 中医临床诊疗术语-证候同义词典
- NovAtelConvert-Setup
- MySql主备数据库配置
- BGP路由基本配置(可以互通)
- 《软件方法2024版》公开内容202405更新-epub版
- 适用于tensorflow-2.11.0 CUDA版本11.2的cuDNN8.1版本
- 5Y study学习平台2016计算机基础-综合测试(8)_哔哩哔哩_bilibili_2580252704.mp4
- (大赛作品)STM32实现的F072RB NUCLEO智能家居控制.zip
- STM32实现的数字示波器源码+数字信号处理教程、配套实例.zip
- 【cookie续续】【cookie续续】【cookie续续】
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈