- 经验
- 8
- 分贝
- 0
- 家园分
- 35
- 在线时间:
- 16 小时
- 最后登录:
- 2016-7-25
- 帖子:
- 5
- 精华:
- 0
- 注册时间:
- 2014-5-26
- UID:
- 1020149
注册:2014-5-261
|
发表于 2014-11-19 21:45:51
|显示全部楼层
%%等效基带调制
clc;
clear all;
close all;
J=sqrt(-1);%%虚数符号
nsymbol = 1e6;%%每种信噪比下符号数的发送符号数2^12 10^6
data=randi([0,1],1,nsymbol);%%产生1行,nsymbol列均匀分布的随机数0,1
bpsk=2*data-1;%%调制,0转化为-1;1转化为1
spow = norm(bpsk_mod).^2/nsymbol;%%求每个符号的平均值,其中norm是求向量2范数函数
SNR_dB = 1:10;%%%信噪比dB形式
SNR = 10.^(SNR_dB/10);%%信噪比转化为线性值
for loop= 1:length(SNR)
sigma = sqrt(spow/(2*SNR(loop)));%%%根据符号功率求噪声功率
s_receive = bpsk_mod+sigma*(randn(1,length(bpsk_mod))+J*randn(1,length(bpsk_mod)));%%添加复高斯白噪声
%%%解调
[err,ser(loop)] = symerr(data,data_receive);%误码率
end
ser_theory = qfunc(sqrt(2*SNR));%理论误码率,注意Q函数和误差函数的对应关系
semilogy(SNR_dB,ser,'-k*',SNR_dB,ser_theory,'-bo');
title('BPSK信号在AWGN信道下的性能');
xlabel('信噪比/dB');ylabel('误码率');
legend('误码率','理论误码率');
grid on;
|
|