?? pskmoto.m
字號:
function [pb,ps]=pskmoto(snr_in_dB)
%pskmoto求出以db為單位的給定信噪比的比特誤碼率和符號誤碼率
N=1000;%符號數(shù)
E=1;
snr=10^(snr_in_dB/10); %計算信噪比的數(shù)值
sgma=sqrt(E/snr)/2;
%信號映射
n=[0 0];
s00=[1 0]; s01=[0 1]; s11=[-1 0]; s10=[0 -1];
for i=1:N
temp=rand
if(temp<0.25)
dsource1(i)=0;
dsource2(i)=0;
elseif(temp<0.5)
dsource1(i)=0;
dsource2(i)=0;
elseif(temp<0.75)
dsource1(i)=1;
dsource2(i)=0;
else
dsource1(i)=1;
dsource2(i)=1;
end;
end;
%判決、誤碼率的計算
numofsymbolerror=0; numofbiterror=0;
for i=1:N
%在判決器的接收端的信號,對于第i個符號為:
n=gngauss(sgma);
if((dsource1(i)==0)&(dsource2(i)==0))
r=s00+n;
elseif((dsource1(i)==0)&(dsource2(i)==1))
r=s01+n;
elseif((dsource1(i)==1)&(dsource2(i)==0))
r=s10+n;
else
r=s11+n;
end
%以下計算互相關量度
c00=dot(r,s00);
c01=dot(r,s01);
c10=dot(r,s10);
c11=dot(r,s11);
%第i個符號的判決如下進行
c_max=max([c00 c01 c10 c11]);
if(c00==c_max)
decis1=0; decis2=0;
elseif(c01==c_max)
decis1=0; decis2=1;
elseif(c10==c_max)
decis1=1; decis2=0;
elseif(c11==c_max)
decis1=1; decis2=1;
end;
%若判決不對,記錯器數(shù)加1
symbolerror=0;
if(decis1~=dsource1(i))
numofbiterror=numofbiterror+1; symbolerror=1;
end;
if(decis2~=dsource2(i))
numofbiterror=numofbiterror+1; symbolerror=1;
end;
if(symbolerror==1)
numofsymbolerror=numofsymbolerror+1;%誤符號數(shù)加1
end;
end;
ps=numofsymbolerror/N; %總共發(fā)出N個符號
pb=numofbiterror/(2*N); %總共發(fā)出2N個比特
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -