nocoop.m

来自「这个代码是matlab仿真的在没有协作用户的情况下」· M 代码 · 共 38 行

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %NOCOOP
 %
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
echo off;clear all;close all;clc;
tic
N=100;
L=65;    %一帧长度
BerSnrTable=zeros(20,5);
for snr=0:25
    BerSnrTable(snr+1,1) = snr;
    sig=1/sqrt(10^(snr/10));
    for i=1:N
        BitsTx = floor(rand(1,L)*2);
        BitsTxcnv=cnv(BitsTx);
        Mod8Tx=mod_8psk(BitsTxcnv);
        M=length(Mod8Tx);
        %以下为假设信道模型和噪声模型
        H1d=randn(1,M)+j*randn(1,M);
        Z1d=randn(1,M)+j*randn(1,M);
        Y1d=H1d.*Mod8Tx+sig*Z1d;
            R1d=conj(H1d).*Y1d; 
            BitR1d=demod_8psk(R1d);
            BitR1dviterbi=viterbi(BitR1d);
            BitR1dviterbi=BitR1dviterbi(1:length(BitR1dviterbi)-1);
            [Num,Ber] = symerr(BitR1dviterbi,BitsTx);
            BerSnrTable(snr+1,2)=BerSnrTable(snr+1,2)+Num;    
    end
    BerSnrTable(snr+1,3)=BerSnrTable(snr+1,2)/(L*N);  %此处将M改成N 
end   
semilogy(BerSnrTable(:,1),BerSnrTable(:,3),'r*-');
figure
semilogy(BerSnrTable(:,1),BerSnrTable(:,4),'g*-');
time_of_sim = toc
echo on;