wumalv.m

以2PSK为例仿真说明数字频带传输的整个过程。

clear all
close all
A=1;    %载波振幅
fc=2;   %载波频率（Hz）
SNRindB1=-5:1:12;     % 信噪比取值向量，dB为单位
SNRindB2=-5:0.2:12;   % 信噪比取值向量，dB为单位
N_sample=3;%每个码元的的采样点数
N=10000;  %  码元数
Ts=1;   % 码元宽度
d=sign(randn(1,N));
df=0.01;
simu_err_prb=zeros(1,length(SNRindB1));
%理论误码率
for i=1:length(SNRindB2)
    % 计算信噪比值
    SNR=10^(SNRindB2(i)/10);
 % 计算普通接收机的理论误码率。
theo_err_prb(i)=0.5*erfc(sqrt(SNR));%一般蒙特卡洛仿真是针对最佳接收机的，对于来自同一信道的接收信号，信道高斯噪声的功率谱密度是不变的，最佳接收机性能之所以比一般接收机好，是由于其输入噪声比是一般接收机的2倍（一般接收机带通滤波器带宽按2/Ts计算）；故信道高斯白噪声的单边带功率谱密度即可按最佳接收机设计，也可按普通接收机设计。
%互补误差函数
theo_err_prb1(i)=0.5*erfc(sqrt(2*SNR));%最佳接收机理论误码率曲线          
end
%计算普通接收机实际误码率
for i=1:length(SNRindB1)
[numoferr,panjue,desingal,t]=bpskberr(A,fc,SNRindB1(i),N_sample,N,Ts,d,df);
simu_err_prb(i)=numoferr;
end
% 误码率曲线图：估计值和理论值曲线对比图
semilogy(SNRindB2,theo_err_prb,SNRindB1,simu_err_prb,'o',SNRindB2,theo_err_prb1,'*');
axis([-5 12 0.00000001 1]);
xlabel('SNR in dB');
ylabel('Prb of Err');
legend('普通接收机理论误码率','普通接收机实际误码率','最佳接收机理论误码率');