信号检测试验三.docVIP

实验三、信号检测的基本准则 实验目的： 掌握最大似然准则的基本原理； 掌握贝叶斯准则的基本原理； 掌握最小均方误差准则的基本原理； 掌握带限系统的等效基带仿真方案； 学习Simulink工具的使用； 掌握蒙特卡罗的仿真方法和思想。 实验内容： 对如下图所示的二元数字通信系统进行仿真 假设系统中的噪声是高斯白噪声，若对每个符号进行1次采样，仿真采用最大似然准则时（此时每信号帧包含100个信息比特）系统的误码性能，画出比特信噪比（亦即）与误码率的关系曲线（要求的信噪比区间为0~15dB）。答：利用matlab软件建立simulink仿真模型如图1： 图1.二元数字通信系统simulink仿真模型 结果讨论：最大似然中，先验等概§=0.5，判决门限为XT=0.5,此时误比特信噪比与误码率的关系曲线如图2所示： 图2.仿真结果与理论曲线对比matlab代码如下： clear; s0=0;s1=5; %定义信号幅度 p0=0.5;p1=1-p0; A1=0:0.5:15;%信噪比范围（DB） A2=10.^(A1/10);%信噪比 Sigma2=(s1^2)./A2;%定义方差Sigma2 N=1e5;%序列长度 for k=1:length(Sigma2) X=(rand(1,N)p0);%信源信号输入 n=sqrt(Sigma2(k))*randn(1,N);%定义噪声 Xi=s1*X+n;%接收机判决输入信号 B=(s0+s1)/2+Sigma2(k)/(s1-s0)*log(p0/p1);%定义最佳判决门限B Y=(XiB);%判决输出 err(k)=(sum(X-Y~=0))/N;%误码率统计 end semilogy(A1,err,*r); hold on %仿真输出 for k=1:length(Sigma2) B=(s0+s1)/2+Sigma2(k)/(s1-s0)*log(p0/p1);%定义最佳判决门限B pe0=0.5-0.5*erf((B-s0)/(sqrt(2*Sigma2(k))));%发0出错率 pe1=0.5+0.5*erf((B-s1)/(sqrt(2*Sigma2(k))));%发1出错率 pe(k)=p0*pe0+p1*pe1;%平均出错率 end semilogy(A1,pe); xlabel(信噪比A^2/Sigma^2（dB）); ylabel(错误率pe); legend(仿真结果,理论曲线); 假设系统中的噪声是高斯白噪声，若对每个符号进行1次采样，仿真采用贝叶斯准则时（此时每信号帧包含100个信息比特）系统的误码性能，画出比特信噪比（亦即）与误码率的关系曲线（要求的信噪比区间为0~15dB）。 答：根据上一问题的模型，再加入代价因子,并修改门限值得结果如图3：  图3.贝叶斯准则下的仿真结果matlab代码如下：clear; s0=0;s1=5; p0=0.5;p1=1-p0; C00=0;C11=0;C01=1;C10=1; A1=0:0.5:15;%信噪比范围（DB） A2=10.^(A1/10);%信噪比 Sigma2=(s1^2)./A2;%定义方差Sigma2 N=1e5;%序列长度 for k=1:length(Sigma2) X=(rand(1,N)p0);%信源信号输入 n=sqrt(Sigma2(k))*randn(1,N);%定义噪声 Xi=s1*X+n;%接收机判决输入信号 B=(p0*(C10-C00))./(p1*(C01-C11));%定义最佳判决门限B Y=(XiB);%判决输出 err(k)=(sum(X-Y~=0))/N;%误码率统计 end semilogy(A1,err,+r); hold on %仿真输出 for k=1:length(Sigma2) B=(p0*(C10-C00))./(p1*(C01-C11));%定义最佳判决门限B pe0=0.5-0.5*erf((B-s0)/(sqrt(2*Sigma2(k))));%发0出错的概率 PE0=C00*(1-pe0)+C10*pe0;%发0出错的风险 pe1=0.5+0.5*erf((B-s1)/(sqrt(2*Sigma2(k))));%发1出错的概率 PE1=C01*pe1+(1-pe1)*C11;%发1出错的风险 pe(k)=p0*PE0+p1*PE1;%平均出错率 end semilogy(A1,pe); xlabel(信噪比A^2/Sigma^2（dB）); ylabel(错误率pe); leg