%?計算2FSK的信噪比
function?[p]=BFSK（snr_in_dB）
N=10000;
Eb=1;				??
snr=10^（snr_in_dB/10）;	?	??	
sgma=sqrt（Eb/（2*snr））;	??	??	
phi=0;??????????????????????????%?相移
numoferr=0;?????????????????????%?誤碼率

%?產生信號
for?i=1:N	
??temp=rand;			??	
??if?（temp<0.5）
????dsource（i）=0;
??else
????dsource（i）=1;
??end
end


for?i=1:N
??if?（dsource（i）==0）
????r0c=sqrt（Eb）*cos（2*pi+phi）+Gngauss（sgma）;
????r0s=sqrt（Eb）*sin（2*pi+phi）+Gngauss（sgma）;
????r1c=Gngauss（sgma）;
????r1s=Gngauss（sgma）;
??else
????r0c=Gngauss（sgma）;
????r0s=Gngauss（sgma）;
????r1c=sqrt（Eb）*cos（4*pi+phi）+Gngauss（sgma）;
????r1s=sqrt（Eb）*sin（4*pi+phi）+Gngauss（sgma）;
??end
??
??r0=r0c^2+r0s^2;
??r1=r1c^2+r1s^2;
??
??%?解調
??if?（r0>r1）
????decis=0;
??else
????decis=1;
??end
??
??%?計算誤碼率
??if?（decis~=dsource（i））
????numoferr=numoferr+1;
??end
end

p=numoferr/N;	

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2020-07-20?19:01??1\
?????文件?????????890??2018-12-29?20:47??1\BFSK.m
?????文件????????1177??2018-12-30?17:14??1\BFSK3way.m
?????文件????????1127??2018-12-30?15:04??1\BFSKrayleigh.m
?????文件?????????627??2018-12-29?20:50??1\BPSK.m
?????文件????????1118??2018-12-30?17:15??1\BPSK3way.m
?????文件????????1059??2018-12-30?15:24??1\BPSKrayleigh.m
?????文件?????????345??2018-12-29?20:37??1\Gngauss.m
?????文件????????1164??2018-12-30?09:02??1\QAM_16.m
?????文件????????1384??2018-12-30?17:15??1\QAM_16_3way.m
?????文件????????1338??2018-12-30?16:26??1\QAM_16_rayleigh.m
?????文件????????2185??2018-12-30?09:08??1\QAM_64.m
?????文件????????1529??2018-12-30?17:15??1\QAM_64_3way.m
?????文件????????1483??2018-12-30?16:26??1\QAM_64_rayleigh.m
?????文件????????1501??2018-12-30?09:32??1\QPSK.m
?????文件????????1236??2018-12-30?17:14??1\QPSK3way.m
?????文件????????1193??2018-12-30?15:50??1\QPSKrayleigh.m
?????文件????????1274??2020-07-20?19:01??1\Untitled.m
?????文件?????????901??2018-12-30?17:13??1\simulation_3way.m
?????文件?????????831??2018-12-30?09:11??1\simulation_awgn.m
?????文件????????1011??2018-12-30?16:21??1\simulation_rayleigh.m