function?[ber]=BPSK（K）
fc=2;fs=16;fd=1;
snr=K-10*log10（fs/2）;???%SNR與Eb/N0之間的轉換，具體詳見實驗報告
alpha=1;
N=1000000;??????????????%取一百萬點，使計算更精確
ts=1/fs;Ts=1;
t=0:ts:N*Ts-ts;?
t=t‘;
B=randint（1N）;
sym=[-1;1];
A=sym（B+1）;
H=rcosine（fdfs‘sqrt‘alpha）;????????%根升余弦濾波器
A1=upsample（Afs/fd）;
A1=filter（H1A1）;
X=A1.*cos（2*pi*fc*t）;
X=awgn（Xsnr‘measured‘）;?????????????%信道中加入高斯白噪聲
R=X.*cos（2*pi*fc*t）;
R=filter（H1R）;
R2=downsample（Rfs/fd）;???????????????%下采樣，符號判決
R2=R2>0;
delay=（length（H）-1）/fs;???????????????%計算通過濾波器的延時
for?i=1:1:N-delay????????????????????????
????Sr（i）=R2（i+delay）;
end
%{
for?i=N-delay+1:N?????%由于延時，有一部分沒被接收，此處全置0，會有一定誤差
????Sr（i）=0;
end
%}
errbit=sum（Sr~=B（1:N-delay））;
ber=errbit/（N-delay）;??????????%誤比特率

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????847??2013-12-08?12:32??BPSK.m
?????文件?????????271??2013-12-01?14:09??mybpsk.m
?????文件?????????245??2013-12-01?13:34??myqpsk.m
?????文件????????1490??2013-12-01?14:15??QPSK.m