clear?all;
sn=0.1:0.01:100;??%定義信噪比序列
snlg=20*log10（sn）;%將信噪比轉化為dB表示
sdouble=sqrt（sn）;
ssingle=sqrt（sn/2）;
bdouble=Qxiu（sqrt（sn））;?%這里需要調用Qxiu函數?function?y=Qxiu（x）?y=erfc（x）/sqrt（2）;
bsingle=Qxiu（ssingle）;
semilogy（snlgbdouble）;
hold?on;
semilogy（snlgbsingle‘r:‘）;
axis（[-20?30?0.0000001?1]）;
i=1100;%通過觀察曲線，大致判斷出雙極性小于給定誤比特率的信噪比位置
while（i????if（bdouble（i）<10^（-6））
????????bsn=sn