msd CPM多符号检测算法调制解调

CPM多符号检测算法调制解调，译码性能比viterbi译码优秀2~3dB CPM multi-symbol detection algorithm modulation and demodulation, decoding performance even better than viterbi decoder 2 ~ 3dB)

clear all;
close all;
clc;
EbN0dB=-10:10;%噪声范围
N_sample=8;%采样数
N=1004;%产生的码元数
A=1;%信号幅度
fc=2;%载波频率
Ts=1;%信号持续周期
dt=Ts/N_sample;%抽样间隔
gt=ones（1N_sample）;
L0=5;
m=1;p=2;
h=m/p;%调制指数
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）;%调用generate函数产生一个（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 ╚²╖√║┼\3╖√║┼╧α╕╔╜Γ╡≈.jpg
     文件       55331  2020-09-23 19:34  msd ╚²╖√║┼\5.jpg
     文件       10350  2020-09-23 19:34  msd ╚²╖√║┼\FEIXIANGGAN.asv
     文件       10305  2020-09-23 19:34  msd ╚²╖√║┼\FEIXIANGGAN.m
     文件         274  2020-09-23 19:34  msd ╚²╖√║┼\findmax.m
     文件         271  2020-09-23 19:34  msd ╚²╖√║┼\fm.m
     文件         419  2020-09-23 19:34  msd ╚²╖√║┼\generate.m
     文件         236  2020-09-23 19:34  msd ╚²╖√║┼\generateQn.m
     文件           6  2020-09-23 19:34  msd ╚²╖√║┼\read me.txt
     文件         136  2020-09-23 19:34  msd ╚²╖√║┼\sigexpand.m
     文件        2512  2020-09-23 19:34  msd ╚²╖√║┼\untitled.fig
     文件       63501  2020-09-23 19:34  msd ╚²╖√║┼\untitled.jpg