1. 首页
  2. 文章中心
  3. 通信系统

【数字基带传输】基础数字调制技术（Digital Modulation）的误码率分析和MATLAB实现简介一、基础数字调制技术框架二、调制技术的误码率和误比特率分析三、MATLAB实现

62 阅读 0 评论
我是靠谱客的博主 还单身可乐，最近开发中收集的这篇文章主要介绍【数字基带传输】基础数字调制技术（Digital Modulation）的误码率分析和MATLAB实现简介一、基础数字调制技术框架二、调制技术的误码率和误比特率分析三、MATLAB实现，觉得挺不错的，现在分享给大家，希望可以做个参考。

概述

简介

由于基础的数字调制技术比较基础，原理也简单，所以这篇文章不打算叙述各种基础地数字调制技术，只会给出一个大概的框架。但是会着重于误码率的分析。而数字调制技术在现代有很多新的发展，如最出名的OFDM等，那将会在本专栏别的文章中讲述。基础的误码率误比特率分析知识请看本专栏另一篇文章：

【数字基带传输】误码率的分析方法（BER Performance）_Tomatomatodo的博客-CSDN博客本篇文章将讲解分析各种调制技术的误码率（SER: SymbolError Rate）和误比特率（BER: Bit Error Rate）的基本方法。https://blog.csdn.net/Tomatomatodo/article/details/122031050?spm=1001.2014.3001.5501

目录

简介

一、基础数字调制技术框架

二、调制技术的误码率和误比特率分析

2.1 BPSK和BASK的误码率分析

2.2 QPSK和QAM的误码率分析

2.3 BFSK的误码率分析

三、MATLAB实现

一、基础数字调制技术框架

如李晓峰的 通信原理 一书所说：“许多重要的通信信道是带通型的，这些信道的传输频带位于某非零频率附近，例如，所有的无线电信道（微波与卫星通信）和很多的有线信道都是带通的。在带通信道上传输数字信号的方法统称为 数字频带传输技术，也称为数字调制技术。广义上讲，调制包括基带调制和带通调制，但在很多场合中，主要指带通调制”

而考虑正弦波：

s(t)=Acos(omega t+theta )

我们可以将信息封装进 相位（Phase），频率（Frequency），振幅（Amplitude）三种地方，从而对应诞生了 PSK. FSK. ASK。而进而诞生了：

1. BPSK. BFSK. BASK 等二元调制

2. QPSK, QAM 等调制

3. 差分调制： DPSK.DQPSK, 对应差分检测法

4. 多元传输调制： MPSK.MFSK,MASK,M-QAM

 

二、调制技术的误码率和误比特率分析

以下的算法知识如果不太清楚，可以看本专栏另一篇文章，讲述基本的误码率分析知识：

【数字基带传输】误码率的分析方法（BER Performance）_Tomatomatodo的博客-CSDN博客本篇文章将讲解分析各种调制技术的误码率（SER: SymbolError Rate）和误比特率（BER: Bit Error Rate）的基本方法。https://blog.csdn.net/Tomatomatodo/article/details/122031050?spm=1001.2014.3001.5501

以下误码率分析都考虑最佳接收，所以有：

P_e=Q(frac{y_{s1}-y_{s0}}{2sigma _n})=Q(sqrt{frac{E_d}{2N_0}})

且保证是控制变量，有以下：

1. 码元周期 Ts = 1s， 所以码率 Rs = 1 sym/s

2. 采用格雷编码（Gray Code）

一些命名规则：

1. 误码率为： P_e, 或 SER （Symbol Error Rate）

2. 误比特率为： P_b 或 BER (Bit Error Rate）

3. 码率为： R_s

4. 比特率为： R_b

2.1 BPSK和BASK的误码率分析

 

2.2 QPSK和QAM的误码率分析

 可以看到有趣的是，QPSK/QAM和BPSK/BASK的误比特率BER是一样的！

2.3 BFSK的误码率分析

注意了，MFSK使用的频率是彼此正交的，一种正交频率代表一个轴，于是我们推导：

 

更高阶的调制我就不写了，但是大家掌握这些分析方法，什么样的都可以算！

三、MATLAB实现

%%
clear;clc;
% Generation of Bits stream
N=12000; %Number of Bits
Bits_Sequence=randi([0 1],1,N); % original bits stream

%%
%BPSK Modulation
BPSK_Sequence=2.*Bits_Sequence-1; % uni-polar to bi-polar

% Adds Gaussian Noise to BPSK signal and Calculate the Bit Error Rate
SNR_dB = [0:2:20];
BPSK_BER_list=[];

% Received BPSK Signals
figure('name','Rx_BPSK','Position',[800,100,1100,900]);
subplot(3,4,1);
plot(real(BPSK_Sequence),imag(BPSK_Sequence),'o');
grid on;
hold on;
line([-3,3],[0,0],'Color','k');
line([0,0],[-3,3],'Color','k');
text(1.2,0.2,'1','FontSize',12,'Interpreter','latex');
text(-1.2,0.2,'0','FontSize',12,'Interpreter','latex');
title('Original Constellation Diagram');
xlim([-2,2]);
ylim([-2,2]);
hold off;
len = length(BPSK_Sequence);

for i=1:1:size(SNR_dB,2)
    N0 = 1/10^(SNR_dB(i)/10); % Noise variance
    Noise = sqrt(N0/2)*((randn(1,len)+1i*randn(1,len))); % AWGN Noise
    Rx_BPSK_sig = BPSK_Sequence + Noise; % Received signal
    subplot(3,4,i+1)
    plot(real(Rx_BPSK_sig(1:100)),imag(Rx_BPSK_sig(1:100)),'*');
    grid on
    hold on
    line([-5,5],[0,0],'Color','k');
    line([0,0],[-5,5],'Color','k');
    title(['When The SNR is',num2str(SNR_dB(i))]);
    xlim([-4,4]);
    ylim([-4,4]);
    hold off
    % Demodulation
    De_BPSK_Bits = (real(Rx_BPSK_sig)>0); %Demodulated Sequence
    Error_bits_BPSK = Bits_Sequence - De_BPSK_Bits; %Calculate Bit Errors
    BPSK_BER_list(end+1) = sum(abs(Error_bits_BPSK))/N; % BER
end

%%
%QPSK AWGN
% Split the bits into two streams
Bits1 = Bits_Sequence(1:2:end); % Odd order (shape 1*5000)
Bits2 = Bits_Sequence(2:2:end); % Even order

% QPSK pi/4 radians constellation
qpsk_sig =(((Bits1==0).*(Bits2==0)*(exp(1i*pi/4))+(Bits1==0).*(Bits2==1)...
*(exp(3*1i*pi/4))+(Bits1==1).*(Bits2==1)*(exp(5*1i*pi/4))...
+(Bits1==1).*(Bits2==0)*(exp(7*1i*pi/4))));

% Adds Gaussian Noise to QPSK signal
SNR_dB = [0:2:20];
QPSK_BER_list=[];
figure('name','Rx_QPSK','Position',[800,100,1100,900]);
subplot(3,4,1)
plot(real(qpsk_sig),imag(qpsk_sig),'o');
grid on
hold on
line([-3,3],[0,0],'Color','k');
line([0,0],[-3,3],'Color','k');
text(0.75,0.75,'00','FontSize',12,'Interpreter','latex');
text(-0.75,0.75,'01','FontSize',12,'Interpreter','latex');
text(-0.75,-0.75,'11','FontSize',12,'Interpreter','latex');
text(0.75,-0.75,'10','FontSize',12,'Interpreter','latex');
title('Original Constellation Diagram');
xlim([-2,2]);
ylim([-2,2]);
hold off
len = length(qpsk_sig);

for i=1:1:size(SNR_dB,2)
    N0 = 1/10^(SNR_dB(i)/10); % Noise variance
    Noise = sqrt(N0/2)*((randn(1,len)+1i*randn(1,len))); % AWGN Noise
    Rx_QPSK_sig = qpsk_sig + Noise; % Received signal
    subplot(3,4,i+1)
    plot(real(Rx_QPSK_sig(1:100)),imag(Rx_QPSK_sig(1:100)),'*'); % Select some of the points
    grid on
    hold on
    line([-5,5],[0,0],'Color','k');
    line([0,0],[-5,5],'Color','k');
    title(['When The SNR is',num2str(SNR_dB(i))]);
    xlim([-4,4]);
    ylim([-4,4]);
    hold off
    % Demodulation
    Bits4 = (real(Rx_QPSK_sig)<0);
    Bits3 = (imag(Rx_QPSK_sig)<0);
    Demod_qpsk_bits = zeros(1,2*length(Rx_QPSK_sig)); % Demodulated bits
    Demod_qpsk_bits(1:2:end) = Bits3;
    Demod_qpsk_bits(2:2:end) = Bits4;
    Error_bits_QPSK = Bits_Sequence - Demod_qpsk_bits; %Calculate Bit Errors
    QPSK_BER_list(end+1) = sum(abs(Error_bits_QPSK))/N; % BER
end

%%
%16-QAM AWGN
% Split the bits into four streams
%Source
%rand() generate 1*num_bits shape matrix with number randomly from (0-1)
QAM16_Bits = reshape(Bits_Sequence,4,length(Bits_Sequence)/4);
Bits1 = QAM16_Bits(1,:);
Bits2 = QAM16_Bits(2,:);
Bits3 = QAM16_Bits(3,:);
Bits4 = QAM16_Bits(4,:);
% normalizing factor
normalizer = sqrt(1/10);	% (2*4+10*8+18*4)/16=10			

% bit mapping        
qam16_sig = normalizer*(1j.*(2.*Bits2-1).*(-2.*Bits3+3)+(2.*Bits1-1).*(-2.*Bits4+3));

% Adds Gaussian Noise to QPSK signal
SNR_dB=[0:2:20];
QAM16_BER_list=[];
figure('name','Rx_QAM16','Position',[200,100,1600,900]);
subplot(3,4,1)
plot(real(qam16_sig),imag(qam16_sig),'ro'); 
grid on
hold on
line([-3,3],[0,0],'Color','k');
line([0,0],[-3,3],'Color','k');
text(0.3562,0.3562,'1111','FontSize',12,'Interpreter','latex');
text(-0.3562,0.3562,'0111','FontSize',12,'Interpreter','latex');
text(-0.3562,-0.3562,'0011','FontSize',12,'Interpreter','latex');
text(0.3562,-0.3562,'1011','FontSize',12,'Interpreter','latex');
text(0.9987,0.9987,'1100','FontSize',12,'Interpreter','latex');
text(-0.9987,0.9987,'0100','FontSize',12,'Interpreter','latex');
text(-0.9987,-0.9987,'0000','FontSize',12,'Interpreter','latex');
text(0.9987,-0.9987,'1000','FontSize',12,'Interpreter','latex');
text(0.3562,0.9987,'1101','FontSize',12,'Interpreter','latex');
text(-0.3562,0.9987,'0101','FontSize',12,'Interpreter','latex');
text(-0.3562,-0.9987,'0001','FontSize',12,'Interpreter','latex');
text(0.3562,-0.9987,'1001','FontSize',12,'Interpreter','latex');
text(0.9987,0.3562,'1110','FontSize',12,'Interpreter','latex');
text(-0.9987,0.3562,'0110','FontSize',12,'Interpreter','latex');
text(-0.9987,-0.3562,'0010','FontSize',12,'Interpreter','latex');
text(0.9987,-0.3562,'1010','FontSize',12,'Interpreter','latex');
title('16QAM Constellation');
xlim([-2,2]);
ylim([-2,2]);
hold off
len = length(qam16_sig);

for i=1:1:size(SNR_dB,2)
    N0 = 1/10^(SNR_dB(i)/10); % Noise variance
    Noise = sqrt(N0/2)*((randn(1,len)+1i*randn(1,len))); % AWGN Noise
    Rx_qam16_sig = qam16_sig + Noise; % Received signal
    subplot(3,4,i+1)
    plot(real(Rx_qam16_sig(1:100)),imag(Rx_qam16_sig(1:100)),'*'); % Select some of the points
    grid on
    hold on
    line([-5,5],[0,0],'Color','k');
    line([0,0],[-5,5],'Color','k');
    title(['When The SNR is',num2str(SNR_dB(i))]);
    xlim([-4,4]);
    ylim([-4,4]);
    hold off
    % Demodulation
    ab = 1/sqrt(10);
    Bits6 = imag(Rx_qam16_sig)>0;
    Bits7 = (imag(Rx_qam16_sig)<2*ab) & (imag(Rx_qam16_sig)>-2*ab);
    Bits5 = real(Rx_qam16_sig)>0;
    Bits8 = (real(Rx_qam16_sig)<2*ab) & (real(Rx_qam16_sig)>-2*ab);
    % Combine into single stream
    comb = [Bits5; Bits6; Bits7; Bits8];
    Demod_qam16_bits = reshape(comb,1,4*length(comb));
    Error_bits_qam16 = Bits_Sequence - Demod_qam16_bits;	  % Calculate Bit Errors
    QAM16_BER_list(end+1) = sum(abs(Error_bits_qam16))/N;	% BER
end

%%
% Raylaigh Channel - BPSK
SNR_dB=[0:2:20];
figure('name','Rx_BPSK','Position',[800,100,1100,900]);
subplot(3,4,1)
plot(real(BPSK_Sequence),imag(BPSK_Sequence),'o');
grid on
hold on
line([-3,3],[0,0],'Color','k');
line([0,0],[-3,3],'Color','k');
text(1.2,0.2,'1','FontSize',12,'Interpreter','latex');
text(-1.2,0.2,'0','FontSize',12,'Interpreter','latex');
title('Original Constellation Diagram');
xlim([-2,2]);
ylim([-2,2]);
hold off
len=length(BPSK_Sequence);
R_Chann=sqrt(0.5)*(randn(1,len)+1i*randn(1,len));
BPSK_R_BER_list=[];
for i=1:1:size(SNR_dB,2)
    N0 = 1/10^(SNR_dB(i)/10); % Noise variance
    Noise = sqrt(N0/2)*((randn(1,len)+1i*randn(1,len))); % AWGN Noise
    Rx_BPSK_sig = BPSK_Sequence.*R_Chann + Noise; % Received signal
    Es_BPSK_sig = Rx_BPSK_sig./R_Chann;
    subplot(3,4,i+1)
    plot(real(Es_BPSK_sig(1:100)),imag(Es_BPSK_sig(1:100)),'*');
    grid on
    hold on
    line([-5,5],[0,0],'Color','k');
    line([0,0],[-5,5],'Color','k');
    title(['When The SNR is',num2str(SNR_dB(i))]);
    xlim([-4,4]);
    ylim([-4,4]);
    hold off
    % Demodulation
    De_BPSK_Bits = (real(Es_BPSK_sig)>0); %Demodulated Sequence
    Error_bits_BPSK = Bits_Sequence - De_BPSK_Bits; %Calculate Bit Errors
    BPSK_R_BER_list(end+1) = sum(abs(Error_bits_BPSK))/N; % BER
end

%%
% Raylaigh to QPSK signal
SNR_dB = [0:2:20];
QPSK_R_BER_list=[];
figure('name','Rx_QPSK','Position',[800,100,1100,900]);
subplot(3,4,1)
plot(real(qpsk_sig),imag(qpsk_sig),'o');
grid on
hold on
line([-3,3],[0,0],'Color','k');
line([0,0],[-3,3],'Color','k');
title('Original Constellation Diagram');
xlim([-2,2]);
ylim([-2,2]);
hold off
len = length(qpsk_sig);
R_Chann=sqrt(0.5)*(randn(1,len)+1i*randn(1,len));
for i=1:1:size(SNR_dB,2)
    N0 = 1/10^(SNR_dB(i)/10); % Noise variance
    Noise = sqrt(N0/2)*((randn(1,len)+1i*randn(1,len))); % AWGN Noise
    Rx_QPSK_sig = qpsk_sig.*R_Chann + Noise; % Received signal
    Es_QPSK_sig = Rx_QPSK_sig./R_Chann;
    subplot(3,4,i+1)
    plot(real(Es_QPSK_sig(1:100)),imag(Es_QPSK_sig(1:100)),'*'); % Select some of the points
    grid on
    hold on
    line([-5,5],[0,0],'Color','k');
    line([0,0],[-5,5],'Color','k');
    title(['When The SNR is',num2str(SNR_dB(i))]);
    xlim([-4,4]);
    ylim([-4,4]);
    hold off
    
    % Demodulation
    Bits4 = (real(Es_QPSK_sig*sqrt(2))<0);
    Bits3 = (imag(Es_QPSK_sig*sqrt(2))<0);
    Demod_qpsk_bits = zeros(1,2*length(Es_QPSK_sig)); % Demodulated bits
    Demod_qpsk_bits(1:2:end) = Bits3;
    Demod_qpsk_bits(2:2:end) = Bits4;
    Error_bits_QPSK = Bits_Sequence - Demod_qpsk_bits; %Calculate Bit Errors
    QPSK_R_BER_list(end+1) = sum(abs(Error_bits_QPSK))/N; % BER
end

%%
% R Channel to QAM16
SNR_dB=[0:2:20];
QAM16_R_BER_list=[];
figure('name','Rx_QAM16','Position',[200,100,1600,900]);
subplot(3,4,1)
plot(real(qam16_sig),imag(qam16_sig),'ro'); 
grid on
hold on
line([-3,3],[0,0],'Color','k');
line([0,0],[-3,3],'Color','k');
text(0.3562,0.3562,'1111','FontSize',12,'Interpreter','latex');
text(-0.3562,0.3562,'0111','FontSize',12,'Interpreter','latex');
text(-0.3562,-0.3562,'0011','FontSize',12,'Interpreter','latex');
text(0.3562,-0.3562,'1011','FontSize',12,'Interpreter','latex');
text(0.9987,0.9987,'1100','FontSize',12,'Interpreter','latex');
text(-0.9987,0.9987,'0100','FontSize',12,'Interpreter','latex');
text(-0.9987,-0.9987,'0000','FontSize',12,'Interpreter','latex');
text(0.9987,-0.9987,'1000','FontSize',12,'Interpreter','latex');
text(0.3562,0.9987,'1101','FontSize',12,'Interpreter','latex');
text(-0.3562,0.9987,'0101','FontSize',12,'Interpreter','latex');
text(-0.3562,-0.9987,'0001','FontSize',12,'Interpreter','latex');
text(0.3562,-0.9987,'1001','FontSize',12,'Interpreter','latex');
text(0.9987,0.3562,'1110','FontSize',12,'Interpreter','latex');
text(-0.9987,0.3562,'0110','FontSize',12,'Interpreter','latex');
text(-0.9987,-0.3562,'0010','FontSize',12,'Interpreter','latex');
text(0.9987,-0.3562,'1010','FontSize',12,'Interpreter','latex');
title('16QAM Constellation');
xlim([-2,2]);
ylim([-2,2]);
hold off
len = length(qam16_sig);
R_Chann=sqrt(0.5)*(randn(1,len)+1i*randn(1,len));
for i=1:1:size(SNR_dB,2)
    N0 = 1/10^(SNR_dB(i)/10); % Noise variance
    Noise = sqrt(N0/2)*((randn(1,len)+1i*randn(1,len))); % AWGN Noise
    Rx_qam16_sig = qam16_sig.*R_Chann + Noise; % Received signal
    Es_qam16_sig = Rx_qam16_sig./R_Chann;
    subplot(3,4,i+1)
    plot(real(Es_qam16_sig(1:100)),imag(Es_qam16_sig(1:100)),'*'); % Select some of the points
    grid on
    hold on
    line([-5,5],[0,0],'Color','k');
    line([0,0],[-5,5],'Color','k');
    title(['When The SNR is',num2str(SNR_dB(i))]);
    xlim([-4,4]);
    ylim([-4,4]);
    hold off
    % Demodulation
    ab = 1/sqrt(10);
    Bits6 = imag(Es_qam16_sig)>0;
    Bits7 = (imag(Es_qam16_sig)<2*ab) & (imag(Es_qam16_sig)>-2*ab);
    Bits5 = real(Es_qam16_sig)>0;
    Bits8 = (real(Es_qam16_sig)<2*ab) & (real(Es_qam16_sig)>-2*ab);
    % Combine into single stream
    comb = [Bits5; Bits6; Bits7; Bits8];
    Demod_qam16_bits = reshape(comb,1,4*length(comb));
    Error_bits_qam16 = Bits_Sequence - Demod_qam16_bits;	  % Calculate Bit Errors
    QAM16_R_BER_list(end+1) = sum(abs(Error_bits_qam16))/N;	% BER
end

%%
%64-QAM AWGN
% Split the bits into four streams

%Source
QAM64_Bits = reshape(Bits_Sequence,6,length(Bits_Sequence)/6);
Bits1 = QAM64_Bits(1,:);
Bits2 = QAM64_Bits(2,:);
Bits3 = QAM64_Bits(3,:);
Bits4 = QAM64_Bits(4,:);
Bits5 = QAM64_Bits(5,:);
Bits6 = QAM64_Bits(6,:);

QAM64_Bits = zeros(1,length(Bits_Sequence)/6);
for i=1:1:length(Bits_Sequence)/6
    QAM64_Bits(i)=32*Bits1(i)+16*Bits2(i)+8*Bits3(i)+4*Bits4(i)+2*Bits5(i)+Bits6(i);
end

% normalizing factor
normalizer = sqrt(1/42);	% (2*4+10*8+18*4+26*8+34*8+50*4+50*8+58*8+74*8+98*4)/64=42			

% bit mapping        
qam64_sig = normalizer*qammod(QAM64_Bits,64);

% Adds Gaussian Noise to QPSK signal
SNR_dB=[0:2:20];
QAM64_BER_list=[];

% Constellation of 64 QAM
figure('name','Rx_QAM64','Position',[200,100,1600,900]);
subplot(3,4,1)
plot(real(qam64_sig),imag(qam64_sig),'ro'); grid on;hold on;
line([-10,10],[0,0],'Color','k');
line([0,0],[-10,10],'Color','k');
title('64QAM Constellation');
xlim([-2,2]);
ylim([-2,2]);
hold off;
len = length(qam64_sig);

for i=1:1:size(SNR_dB,2)
    N0 = 1/10^(SNR_dB(i)/10); % Noise variance
    Noise = sqrt(N0/2)*((randn(1,len)+1i*randn(1,len))); % AWGN Noise
    Rx_qam64_sig = qam64_sig + Noise; % Received signal
    subplot(3,4,i+1)
    plot(real(Rx_qam64_sig(1:100)),imag(Rx_qam64_sig(1:100)),'*'); % Select some of the points
    grid on;
    hold on;
    line([-5,5],[0,0],'Color','k');
    line([0,0],[-5,5],'Color','k');
    title(['When The SNR is',num2str(SNR_dB(i))]);
    xlim([-4,4]);
    ylim([-4,4]);
    hold off
    
    % Demodulation
    de_qam64_sym=qamdemod(Rx_qam64_sig/normalizer,64);
    de_Bits1=zeros(1,length(de_qam64_sym));
    de_Bits2=zeros(1,length(de_qam64_sym));
    de_Bits3=zeros(1,length(de_qam64_sym));
    de_Bits4=zeros(1,length(de_qam64_sym));
    de_Bits5=zeros(1,length(de_qam64_sym));
    de_Bits6=zeros(1,length(de_qam64_sym));
    for j=1:1:length(de_qam64_sym)
        buf=de_qam64_sym(j);
        de_Bits1(j)=floor(buf/32);
        de_Bits2(j)=floor((buf-32*de_Bits1(j))/16);
        de_Bits3(j)=floor((buf-32*de_Bits1(j)-16*de_Bits2(j))/8);
        de_Bits4(j)=floor((buf-32*de_Bits1(j)-16*de_Bits2(j)-8*de_Bits3(j))/4);
        de_Bits5(j)=floor((buf-32*de_Bits1(j)-16*de_Bits2(j)-8*de_Bits3(j)-4*de_Bits4(j))/2);
        de_Bits6(j)=floor((buf-32*de_Bits1(j)-16*de_Bits2(j)-8*de_Bits3(j)-4*de_Bits4(j)-2*de_Bits5(j)));
    end
    
    % Combine into single stream
    comb = [de_Bits1; de_Bits2; de_Bits3; de_Bits4;de_Bits5;de_Bits6];
    Demod_qam64_bits = reshape(comb,1,6*length(comb));
    Error_bits_qam64 = Bits_Sequence - Demod_qam64_bits;	  % Calculate Bit Errors
    QAM64_BER_list(end+1) = sum(abs(Error_bits_qam64))/N;	% BER
end

%%
% QAM64 Passing through Rayleigh Channel
% normalizing factor
normalizer = sqrt(1/42);	% (2*4+10*8+18*4+26*8+34*8+50*4+50*8+58*8+74*8+98*4)/64=42			

% Adds Gaussian Noise to QPSK signal
SNR_dB=[0:2:20];
QAM64_R_BER_list=[];

% Constellation of 64 QAM
figure('name','Rx_QAM64','Position',[200,100,1600,900]);
subplot(3,4,1)
plot(real(qam64_sig),imag(qam64_sig),'ro'); grid on;hold on;
line([-10,10],[0,0],'Color','k');
line([0,0],[-10,10],'Color','k');
title('64QAM Constellation');
xlim([-2,2]);
ylim([-2,2]);
hold off;
len = length(qam64_sig);
R_Chann=sqrt(0.5)*(randn(1,len)+1i*randn(1,len));

for i=1:1:size(SNR_dB,2)
    N0 = 1/10^(SNR_dB(i)/10); % Noise variance
    Noise = sqrt(N0/2)*((randn(1,len)+1i*randn(1,len))); % AWGN Noise
    Rx_qam64_sig = R_Chann.*qam64_sig + Noise; % Received signal
    E_qam64_sig=Rx_qam64_sig./R_Chann;
    subplot(3,4,i+1)
    plot(real(E_qam64_sig(1:100)),imag(E_qam64_sig(1:100)),'*'); % Select some of the points
    grid on;
    hold on;
    line([-5,5],[0,0],'Color','k');
    line([0,0],[-5,5],'Color','k');
    title(['When The SNR is',num2str(SNR_dB(i))]);
    xlim([-4,4]);
    ylim([-4,4]);
    hold off
    
    % Demodulation
    
    de_qam64_sym=qamdemod(E_qam64_sig/normalizer,64);
    de_Bits1=zeros(1,length(de_qam64_sym));
    de_Bits2=zeros(1,length(de_qam64_sym));
    de_Bits3=zeros(1,length(de_qam64_sym));
    de_Bits4=zeros(1,length(de_qam64_sym));
    de_Bits5=zeros(1,length(de_qam64_sym));
    de_Bits6=zeros(1,length(de_qam64_sym));
    for j=1:1:length(de_qam64_sym)
        buf=de_qam64_sym(j);
        de_Bits1(j)=floor(buf/32);
        de_Bits2(j)=floor((buf-32*de_Bits1(j))/16);
        de_Bits3(j)=floor((buf-32*de_Bits1(j)-16*de_Bits2(j))/8);
        de_Bits4(j)=floor((buf-32*de_Bits1(j)-16*de_Bits2(j)-8*de_Bits3(j))/4);
        de_Bits5(j)=floor((buf-32*de_Bits1(j)-16*de_Bits2(j)-8*de_Bits3(j)-4*de_Bits4(j))/2);
        de_Bits6(j)=floor((buf-32*de_Bits1(j)-16*de_Bits2(j)-8*de_Bits3(j)-4*de_Bits4(j)-2*de_Bits5(j)));
    end
    
    % Combine into single stream
    comb = [de_Bits1; de_Bits2; de_Bits3; de_Bits4;de_Bits5;de_Bits6];
    Demod_qam64_bits = reshape(comb,1,6*length(comb));
    Error_bits_qam64 = Bits_Sequence - Demod_qam64_bits;	  % Calculate Bit Errors
    QAM64_R_BER_list(end+1) = sum(abs(Error_bits_qam64))/N;	% BER
end

%%
figure('name','BER in different Channel','Position',[500,100,900,800]);
semilogy(SNR_dB,BPSK_BER_list,'r');text(4,0.002,'BPSK AWGN');
grid on;hold on;
semilogy(SNR_dB,QPSK_BER_list,'b');text(8,0.0005,'QPSK AWGN');
semilogy(SNR_dB,QAM16_BER_list,'g');text(14,0.0016,'QAM16 AWGN');
semilogy(SNR_dB,BPSK_R_BER_list);text(18,0.0025,'BPSK Rayleigh');
semilogy(SNR_dB,QPSK_R_BER_list);text(19,0.005,'QPSK Rayleigh');
semilogy(SNR_dB,QAM16_R_BER_list);text(19,0.021,'QAM16 Rayleigh');
semilogy(SNR_dB,QAM64_BER_list);text(19,0.0085,'QAM64 AWGN');
semilogy(SNR_dB,QAM64_R_BER_list);text(19,0.05,'QAM64 Rayleigh');
legend('BPSK-AWGN','QPSK-AWGN','QAM16-AWGN','BPSK-R','QPSK-R','QAM16-R','QAM64-AWGN','QAM64-R','Location','SouthWest');
title('BER in Rayleigh/AWGN Channel');xlabel('SNR in dB');ylabel('BER');

%%
Eb_N0_bpsk=10*log10(10.^(SNR_dB/10));
Eb_N0_qpsk=10*log10(10.^(SNR_dB/10)/2);
Eb_N0_16qam=10*log10(10.^(SNR_dB/10)/4);
Eb_N0_64qam=10*log10(10.^(SNR_dB/10)/6);
figure('name','BER in different Channel','Position',[500,100,900,800]);
semilogy(Eb_N0_bpsk,BPSK_BER_list,'r');text(4,0.0008,'BPSK AWGN');
grid on;hold on;
semilogy(Eb_N0_qpsk,QPSK_BER_list,'b');text(3,0.002,'QPSK AWGN');
semilogy(Eb_N0_16qam,QAM16_BER_list,'g');text(10,0.0016,'QAM16 AWGN');
semilogy(Eb_N0_bpsk,BPSK_R_BER_list);text(18,0.002,'BPSK Rayleigh');
semilogy(Eb_N0_qpsk,QPSK_R_BER_list);text(15,0.005,'QPSK Rayleigh');
semilogy(Eb_N0_16qam,QAM16_R_BER_list);text(12,0.02,'QAM16 Rayleigh');
semilogy(Eb_N0_64qam,QAM64_BER_list);text(10,0.008,'QAM64 AWGN');
semilogy(Eb_N0_64qam,QAM64_R_BER_list);text(11,0.0536,'QAM64 Rayleigh');
legend('BPSK-AWGN','QPSK-AWGN','QAM16-AWGN','BPSK-R','QPSK-R','QAM16-R','QAM64-AWGN','QAM64-R','Location','SouthWest');
title('BER in R/AWGN Channel');xlabel('Eb/N0 in dB');ylabel('BER');

 

最后

以上就是还单身可乐为你收集整理的【数字基带传输】基础数字调制技术（Digital Modulation）的误码率分析和MATLAB实现简介一、基础数字调制技术框架二、调制技术的误码率和误比特率分析三、MATLAB实现的全部内容，希望文章能够帮你解决【数字基带传输】基础数字调制技术（Digital Modulation）的误码率分析和MATLAB实现简介一、基础数字调制技术框架二、调制技术的误码率和误比特率分析三、MATLAB实现所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供，作为学习参考使用，或来自网络收集整理，版权属于原作者所有。
点赞(41)

相关文章

BPSK调制与解调
BPSK调制与解调
RS编码译码误码率性能matlab仿真1.算法描述2.仿真效果预览3.MATLAB部分代码预览4.完整MATLAB程序
RS编码译码误码率性能matlab仿真1.算法描述2.仿真效果预览3.MATLAB部分代码预览4.完整MATLAB程序
扩频bpsk误码率matlab代码,BPSK扩频误码率matlab仿真.docx
扩频bpsk误码率matlab代码,BPSK扩频误码率matlab仿真.docx
Matlab 编程实现(4，1)二元Repetition Code 基带仿真1 Repeat Code 的简单介绍2 编码与译码的原理3 仿真实现与结果4 总结5 参考资料1 Repeat Code 的简单介绍2 编码与译码的原理3 仿真实现与结果4 总结5 参考资料
Matlab 编程实现(4，1)二元Repetition Code 基带仿真1 Repeat Code 的简单介绍2 编码与译码的原理3 仿真实现与结果4 总结5 参考资料1 Repeat Code 的简单介绍2 编码与译码的原理3 仿真实现与结果4 总结5 参考资料
【数字基带传输】基础数字调制技术（Digital Modulation）的误码率分析和MATLAB实现简介一、基础数字调制技术框架二、调制技术的误码率和误比特率分析三、MATLAB实现
【数字基带传输】基础数字调制技术（Digital Modulation）的误码率分析和MATLAB实现简介一、基础数字调制技术框架二、调制技术的误码率和误比特率分析三、MATLAB实现
matlab中 直扩系统误码率_直接序列扩频系统MATLAB仿真(BPSK调制)
matlab中 直扩系统误码率_直接序列扩频系统MATLAB仿真(BPSK调制)
MIMO大规模天线阵列原理与matlab仿真（含GUI）一 设计实现流程 三 结果分析MATLAB代码展示
MIMO大规模天线阵列原理与matlab仿真（含GUI）一 设计实现流程 三 结果分析MATLAB代码展示
matlab误码率的代码,matlab实现BPSK误码率仿真
matlab误码率的代码,matlab实现BPSK误码率仿真

评论列表共有 0 条评论

立即
投稿
返回
顶部