不同调制方式下误码率性能对比

包含2ASK、4ASK、8ASK等不同MASK调制方式下抗噪声性能对比和抗多径对比，单路多径、二路多径、三路多径，即锐利信道下性能； 包含2PSK、4PSK、8PSK等不同MPSK调制方式下抗噪声性能对比和抗多径对比，单路多径、二路多径、三路多径，即锐利信道下性能； 包含16QAM、64QAM、128QAM等不同MQAM调制方式下抗噪声性能对比和抗多径对比，单路多径、二路多径、三路多径，即锐利信道下性能； 包含OFDM-QPSK OFDM-QAM等不同OFDM调制方式下抗噪声性能对比和抗多径对比，单路多径、二路多径、三路多径，即锐利信道下性能；

close all;
clear all;
clc;
t0=0.049999;   %总仿真时间，单位为s
% 0.004999 50个码元 0.049999 500个码元  0.099999 1000个码元
% 码元宽度 0.0001s， 每个码元100个采样点
ts=1*10^-6;  %采样间隔1us
fc=50*10^3;    %载波频率50kHz
 snr=6; %信噪比，以dB来表示
fs=1/ts;  %采样频率 1M
df=0.3;  %频率分辨率
t=[0:ts:t0];
%计算仿真的数字序列长度
n=size（t）;
N=n（2）/100; % 码元个数

%载波信号生成，频率为fc
c=cos（2*pi*fc*t）;
%本地载波信号
c_local=cos（2*pi*fc.*t）;


m2=fix（2*rand（1N））; %生成消息信号（二进制数字序列）

%将二进制数字序列转换为时域信号，也就是对数字符号进行采样处理
for k=1:N
    for i=（100*k-99）:（100*k）
         m2_time（i）=m2（k）;
    end   
end

%清零处理
for i=（100*（N/2）+1）:（100*N）
    m2_time（i）=0;
end


ask2=m2_time.*c; %生成2ASK时域信号
%
figure（1）    
subplot（311）      %子图语句
plot（tm2_time）     %画出调制信号
axis（[0 0.001 -0.5 3.5]）
title（‘调制信号‘）;
grid on;

subplot（312）
plot（tc）          %画出载波信号
axis（[0 0.001 -3.5 3.5]）
title（‘载波信号‘）;
grid on;

subplot（313）
plot（task2）      %画出已调信号
axis（[0 0.001 -3.5 3.5]）
title（‘ASK已调信号‘）;
grid on;
%}


ask2_noise=awgn（ask2snr‘measured‘）;  %经过AWGN高斯白噪声

%绘制ASK信号波形（受噪声影响）
figure（2）
plot（task2_noise）
grid on;
title（‘ask信号波形（加噪）‘）

axis（[0 0.001 -2 2]）
%本地载波与接收已调信号相乘，实现相干解调
ask2_dem1=ask2_noise.*c_local;

%将相乘结果转换为频域
[ASK2_DEM1ask2_dem1df1]=fftseq（ask2_dem1tsdf）; %df1是频率采样间隔
ASK2_DEM1=ASK2_DEM1/fs; %乘以FFT系数
f=[0:df1:df1*（length（ask2_dem1）-1）]-fs/2; %f是FFT变换后的频域点
%
%将已调信号转换为频域
[ASK2_NOISEask2_noisedf1]=fftseq（ask2_noisetsdf）;
ASK2_NOISE=ASK2_NOISE/fs;
%}
%低通滤波器设计
f_cutoff=50000;          %截止频率
n_cutoff=floor（f_cutoff/df1）; %频率采样个数
H=zeros（size（f））;         %H是频域低通滤波器
H（1:n_cutoff）=2*ones（1n_cutoff）;
H（length（f）-n_cutoff+1:length（f））=2*ones（1n_cutoff）; %H的前n_cutoff个和后n_cutoff个值设为2，其他值为0

%接收信号通过低通滤波器
ASK2_DEM2=H.*ASK2_DEM1;
ask2_dem2=real（ifft（ASK2_DEM2））*fs; %重新转化为时域信号


figure（3）
%subplot（313）
plot（fabs（fftshift（ASK2_NOISE）））   %画出已调信号频率
axis（[-150000 150000 0 0.001]）
xlabel（‘频率‘）;
title（‘已调信号频谱‘）;
grid on;
%}

%抽样判决
k=1;
for i=50:100:（N*100-50）
    if ask2_dem2（i）<=0.5     %根据采样信号幅度来判决，大于0.5为1，小于0.5为0
        m2_sample（k）=0;
    else 
        m2_sample（k）=1;
    end
    k=k+1;
end

%将抽样结果转换为不归零码型
for k=1:（N/2）
    for i=（100*k-99）:（100*k）
         m2_sampletime（i）=m2_sample（k）;
    end   
end

%清零处理
for i=（100*（N/2）+1）:（100*N）
    m2_sampletime（i）=0;
end


%
figure（4）
%将调制信号与解调还原的信号作对比
subplot（311）
plot（tm2_time（1:length（t）））
axis（[0 0.001 -1 2]）
title（‘调制信号‘）
grid on;
subplot（312）
plot（task2_dem2（1:length（t）））
axis（[0 0.001 -1 2]）
title（‘通过LPF后的波形‘）
grid on;
subplot（313）
plot（tm2_sampletime（1:length（t）））
axis（[0 0.001 -1 2]）
title（‘抽样判决后波形‘）
grid on;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-04-18 20:49  MASK\
     目录           0  2017-04-18 20:49  MASK\2ASK\
     文件        3071  2014-04-20 22:45  MASK\2ASK\ask22.m
     文件         525  2014-04-16 22:58  MASK\2ASK\fftseq.m
     目录           0  2017-04-18 20:49  MASK\2ASK_4ASK_8ASK\
     文件         251  2014-04-23 00:11  MASK\2ASK_4ASK_8ASK\2askBER_0_15.mat
     文件         312  2014-04-23 00:11  MASK\2ASK_4ASK_8ASK\4askBER_0_15.mat
     文件         320  2014-04-23 00:12  MASK\2ASK_4ASK_8ASK\8askBER_0_15.mat
     文件        1366  2014-04-21 00:10  MASK\2ASK_4ASK_8ASK\ASK.asv
     文件        1484  2014-04-21 23:18  MASK\2ASK_4ASK_8ASK\ASK.m
     文件         146  2014-04-23 02:36  MASK\2ASK_4ASK_8ASK\BERplot.asv
     文件         146  2014-04-23 02:39  MASK\2ASK_4ASK_8ASK\BERplot.m
     文件        1687  2014-04-21 22:49  MASK\2ASK_4ASK_8ASK\Copy_of_ASK.asv
     文件        1998  2014-04-22 15:25  MASK\2ASK_4ASK_8ASK\Copy_of_ASK.m
     文件         351  2008-07-04 15:07  MASK\2ASK_4ASK_8ASK\F2T.m
     文件         451  2008-09-05 08:17  MASK\2ASK_4ASK_8ASK\jd.m
     文件         543  2014-04-23 01:06  MASK\2ASK_4ASK_8ASK\lpf.asv
     文件         547  2014-04-24 19:34  MASK\2ASK_4ASK_8ASK\lpf.m
     文件        2179  2014-04-22 23:16  MASK\2ASK_4ASK_8ASK\my_2ASK.m
     文件        2757  2014-04-23 01:16  MASK\2ASK_4ASK_8ASK\my_4ASK.asv
     文件        2765  2014-04-24 19:26  MASK\2ASK_4ASK_8ASK\my_4ASK.m
     文件        3847  2014-04-23 00:52  MASK\2ASK_4ASK_8ASK\my_8ASK.asv
     文件        4672  2014-04-24 19:34  MASK\2ASK_4ASK_8ASK\my_8ASK.m
     文件         493  2014-04-17 16:45  MASK\2ASK_4ASK_8ASK\T2F.asv
     文件         567  2014-04-22 00:22  MASK\2ASK_4ASK_8ASK\T2F.m
     文件        1994  2014-04-22 16:11  MASK\2ASK_4ASK_8ASK\test_my_ASK.m
     目录           0  2017-04-18 20:49  MFSK\
     目录           0  2017-04-18 20:49  MFSK\2FSK\
     文件         351  2014-04-24 19:51  MFSK\2FSK\F2T.m
     文件        1907  2014-04-24 19:51  MFSK\2FSK\FSK.m
     文件         543  2014-04-24 19:51  MFSK\2FSK\lpf.m
............此处省略55个文件信息