%分別測(cè)試SNR從0到10時(shí)的系統(tǒng)性能
%初始化
SNR=-2; %信噪比
N=100; %隨機(jī)序列長(zhǎng)度
Fdc=10; %歸一化的采樣頻率
n=0; %總考察信號(hào)數(shù)
cntError=0; %誤碼個(gè)數(shù)
%生成積分算子，解調(diào)使用
m=1:1:10
temp2=sin(2*pi*m/Fdc)
temp2=temp2'; 

while cntError<100 %最外層循環(huán)
%產(chǎn)生0、1隨機(jī)序列
source=randsrc(1,N,[0,1;.5,.5]);
%產(chǎn)生psk調(diào)制信號(hào)
for i=1:N
for j=1:Fdc
modsrc(j+(i-1)*Fdc)=sin(pi*source(i)+2*pi*j/Fdc);
end
end
%加入高斯白噪聲
modsrcnoise=awgn(modsrc,SNR,'measured','db');
%計(jì)算誤碼率
%輸入接收信號(hào)，并判決

for i=1:N
temp1=modsrcnoise((i-1)*Fdc+m);
if temp1*temp2>0
receive(i)=0
else
receive(i)=1
end
end
for i=1:N
if source(i)==receive(i)
else
cntError=cntError+1;
end
end

n=n+N;
ErrorRate=cntError/n;
end
snr=10log（E*m/No）dB m=一個(gè)多進(jìn)制符號(hào)的比特?cái)?shù)
信噪比的意思是信號(hào)的功率/噪聲的功率譜密度
Eb/N0中的Eb代表每傳輸1bit信息所需要的能量,N0代表高斯白噪聲信道的單邊功率譜密度