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

while cntError<100 %最外層循環
%產生0、1隨機序列
source=randsrc(1,N,[0,1;.5,.5]);
%產生psk調制信號
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');
%計算誤碼率
%輸入接收信號，并判決

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=一個多進制符號的比特數