clear?all;
close?all;
clc;
EbN0dB=-10:10;%噪聲范圍
N_sample=8;%采樣數(shù)
N=1004;%產(chǎn)生的碼元數(shù)
A=1;%信號幅度
fc=2;%載波頻率
Ts=1;%信號持續(xù)周期
dt=Ts/N_sample;%抽樣間隔
gt=ones（1N_sample）;
L0=5;
m=1;p=2;
h=m/p;%調(diào)制指數(shù)
Qn=generateQn（mp）;
N0=10.^（-EbN0dB/10）;
sigma=sqrt（N0/2）;%噪聲
error=zeros（1length（N0））;
for?k=1:length（N0）
????num=0;
???while?num<100
????????d=randsrc（1N）;
????????dd=sigexpand（dN_sample）;???
????????d_NRZ=conv（ddgt）;
????????t=0:dt:N*Ts-dt;
????????s_xiangwei（1）=0;
????for?i=1:length（t）-1
??????s_xiangwei（i+1）=s_xiangwei（i）+d_NRZ（i）*dt;
????end
????s_cpfsk2=A*cos（2*pi*fc*t+pi*h*s_xiangwei（1:length（t）））;
????r=A*cos（2*pi*fc*t+pi*h*s_xiangwei（1:length（t）））+sigma（k）*（randn（1length（t）））;%信號加燥
????a=generate（L0）;%調(diào)用generate函數(shù)產(chǎn)生一個（8，3）的所有三符號組合
????a1=a（1:）;a2=a（2:）;a3=a（3:）;a4=a（4:）;a5=a（5:）;a6=a（6:）;a7=a（7:）;a8=a（8:）;
????a9=a（9:）;a10=a（10:）;a11=a（11:）;a12=a（12:）;a13=a（13:）;a14=a（14:）;a15=a（15:）;a16=a（16:）;
????a17=a（17:）;a18=a（18:）;a19=a（19:）;a20=a（20:）;a21=a（21:）;a22=a（22:）;a23=a（23:）;a24=a（24:）;
????a25=a（25:）;a26=a（26:）;a27=a（27:）;a28=a（28:）;a29=a（29:）;a30=a（30:）;a31=a（31:）;a32=a（32:）;
????dd1=sigexpand（a1N_sample）;dd2=sigexpand（a2N_sample）;dd3=sigexpand（a3N_sample）;dd4=sigexpand（a4N_sample）;
????dd5=sigexpand（a5N_sample）;dd6=sigexpand（a6N_sample）;dd7=sigexpand（a7N_sample）;dd8=sigexpand（a8N_sample）;
????dd9=sigexpand（a9N_sample）;dd10=sigexpand（a10N_sample）;dd11=sigexpand（a11N_sample）;dd12=sigexpand（a12N_sample）;
????dd13=sigexpand（a13N_sample）;dd14=sigexpand（a14N_sample）;dd15=sigexpand（a15N_sample）;dd16=sigexpand（a16N_sample）;
????dd17=sigexpand（a17N_sample）;dd18=sigexpand（a18N_sample）;dd19=sigexpand（a19N_sample）;dd20=sigexpand（a20N_sample）;
????dd21=sigexpand（a21N_sample）;dd22=sigexpand（a22N_sample）;dd23=sigexpand（a23N_sample）;dd24=sigexpand（a24N_sample）;
????dd25=sigexpand（a25N_sample）;dd26=sigexpand（a26N_sample）;dd27=sigexpand（a27N_sample）;dd28=sigexpand（a28N_sample）;
????dd29=sigexpand（a29N_sample）;dd30=sigexpand（a30N_sample）;dd31=sigexpand（a31N_sample）;dd32=sigexpand（a32N_sample）;
????d_NRZ1=conv（dd1gt）;d_NRZ2=conv（dd2gt）;d_NRZ3=conv（dd3gt）;d_NRZ4=conv（dd4gt）;
????d_NRZ5=conv（dd5gt）;d_NRZ6=conv（dd6gt）;d_NRZ7=conv（dd7gt）;d_NRZ8=conv（dd8gt）;
????d_NRZ9=conv（dd9gt）;d_NRZ10=conv（dd10gt）;d_NRZ11=conv（dd11gt）;d_NRZ12=conv（dd12gt）;
????d_NRZ13=conv（dd13gt）;d_NRZ14=conv（dd14gt）;d_NRZ15=conv（dd15gt）;d_NRZ16=conv（dd16gt）;
????d_NRZ17=conv（dd17gt）;d_NRZ18=conv（dd18gt）;d_NRZ19=conv（dd19gt）;d_NRZ20=conv（dd20gt）;
????d_NRZ21=conv（dd21gt）;d_NRZ22=conv（dd22gt）;d_NRZ23=conv（dd23gt）;d_NRZ24=conv（dd24gt）;
????d_NRZ25=conv（dd25gt）;d_NRZ26=conv（dd26gt）;d_NRZ27=conv（dd27gt）;d_NRZ28=conv（dd28gt）;
????d_NRZ29=conv（dd29gt）;d_NRZ30=conv（dd30gt）;d_NRZ31=conv（dd31gt）;d_NRZ32=conv（dd32gt）;
????for?L=1:N-4
???????for?st=1:4
?????????s_xiangwei1（1）=Qn（st）/（pi*h）;s_xiangwei2（1）=Qn（st）/（pi*h）;s_xiangwei3（1）=Qn（st）/（pi*h）;s_xiangwei4（1）=Qn（s

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????61979??2020-09-23?19:34??msd?╚2╖√║┼\3╖√║┼╧α╕╔╜Γ╡≈.jpg
?????文件???????55331??2020-09-23?19:34??msd?╚2╖√║┼\5.jpg
?????文件???????10350??2020-09-23?19:34??msd?╚2╖√║┼\FEIXIANGGAN.asv
?????文件???????10305??2020-09-23?19:34??msd?╚2╖√║┼\FEIXIANGGAN.m
?????文件?????????274??2020-09-23?19:34??msd?╚2╖√║┼\findmax.m
?????文件?????????271??2020-09-23?19:34??msd?╚2╖√║┼\fm.m
?????文件?????????419??2020-09-23?19:34??msd?╚2╖√║┼\generate.m
?????文件?????????236??2020-09-23?19:34??msd?╚2╖√║┼\generateQn.m
?????文件???????????6??2020-09-23?19:34??msd?╚2╖√║┼\read?me.txt
?????文件?????????136??2020-09-23?19:34??msd?╚2╖√║┼\sigexpand.m
?????文件????????2512??2020-09-23?19:34??msd?╚2╖√║┼\untitled.fig
?????文件???????63501??2020-09-23?19:34??msd?╚2╖√║┼\untitled.jpg