4qam、16qam、64qam、256qam理论仿真曲线

本博文给出了4qam、16qam、64qam、256qam理论仿真曲线，画出了EbN0 vs BER的曲线图，可以作为大家学习的一个参考。

仿真结果:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Theoretical ber curves of different orders of QAM modulation
% Author: 苏大axing
% Version: 1.2
% Date: 2022/4/3
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
close all;
clear all;
clc%% Parameters
EbN0_dB                         = -5:20;
bits_per_symbol_for_4qam        = 2;
bits_per_symbol_for_16qam       = 4;
bits_per_symbol_for_64qam       = 6;
bits_per_symbol_for_256qam      = 8;
Nsym                            = 2^18;
SNR_for_4qam                    = EbN0_dB+10*log10(bits_per_symbol_for_4qam);
SNR_for_16qam                   = EbN0_dB+10*log10(bits_per_symbol_for_16qam);
SNR_for_64qam                   = EbN0_dB+10*log10(bits_per_symbol_for_64qam);
SNR_for_256qam                  = EbN0_dB+10*log10(bits_per_symbol_for_256qam);disp('           ');disp('           ');
disp('***********************************');
disp('generating binary bits');
disp('***********************************');
% Generating seeds
randn('state',12013)
rand('state',3693)
binary_bits_forqpsk             = randsrc(1,Nsym*bits_per_symbol_for_4qam,[0 1]);
binary_bits_for16qam            = randsrc(1,Nsym*bits_per_symbol_for_16qam,[0 1]);
binary_bits_for64qam            = randsrc(1,Nsym*bits_per_symbol_for_64qam,[0 1]);
binary_bits_for256qam           = randsrc(1,Nsym*bits_per_symbol_for_256qam,[0 1]);% generating MQAM signals
disp('           ');disp('           ');
disp('***********************************');
disp('generating MQAM signal');
disp('***********************************');%%%%%%%%%% 4QAM
paraller_for_4qam               = reshape(binary_bits_forqpsk,bits_per_symbol_for_4qam,Nsym);
k_4qam                          = 2^bits_per_symbol_for_4qam;
[Output]                        = BitsToInteger(paraller_for_4qam);
qam4_signal                     = qammod(Output,k_4qam,'gray','UnitAveragePower', true);% figure(1)
% subplot(2,2,1);
% plot(qam4_signal,'.');
% title('4QAM signal');%%%%%%%%%% 16QAM
paraller_for_16qam              = reshape(binary_bits_for16qam,bits_per_symbol_for_16qam,Nsym);
k_16qam                         = 2^bits_per_symbol_for_16qam;
[Output]                        = BitsToInteger(paraller_for_16qam);
qam16_signal                    = qammod(Output,k_16qam,'gray','UnitAveragePower', true);% figure(1)
% subplot(2,2,2);
% plot(qam16_signal,'.');
% title('16QAM signal');%%%%%%%%%% 64QAM
paraller_for_64qam              = reshape(binary_bits_for64qam,bits_per_symbol_for_64qam,Nsym);
k_64qam                         = 2^bits_per_symbol_for_64qam;
[Output]                        = BitsToInteger(paraller_for_64qam);
qam64_signal                    = qammod(Output,k_64qam,'gray','UnitAveragePower', true);% figure(1)
% subplot(2,2,3);
% plot(qam64_signal,'.');
% title('64QAM signal');%%%%%%%%%% 256QAM
paraller_for_256qam             = reshape(binary_bits_for256qam,bits_per_symbol_for_256qam,Nsym);
k_256qam                        = 2^bits_per_symbol_for_256qam;
[Output]                        = BitsToInteger(paraller_for_256qam);
qam256_signal                   = qammod(Output,k_256qam,'gray','UnitAveragePower', true);% figure(1)
% subplot(2,2,4);
% plot(qam256_signal,'.');
% title('256QAM signal');for i=1:length(EbN0_dB)
received_qam4                   = awgn(qam4_signal,SNR_for_4qam(i),'measured');
received_qam16                  = awgn(qam16_signal,SNR_for_16qam(i),'measured');
received_qam64                  = awgn(qam64_signal,SNR_for_64qam(i),'measured');
received_qam256                 = awgn(qam256_signal,SNR_for_256qam(i),'measured');%% demodulation
%%%%%%%%%% 4QAM
demod_signal_for_4qam           = qamdemod(received_qam4,k_4qam,'gray','UnitAveragePower', true);
[Output]                        = IntegerToBits(demod_signal_for_4qam,bits_per_symbol_for_4qam);
received_Bits_for_4qam          = reshape(Output,1,Nsym*bits_per_symbol_for_4qam);
ErrorNumber_for_4qam            = size(find(received_Bits_for_4qam-binary_bits_forqpsk),2);
BER_for_4qam(i)                 = ErrorNumber_for_4qam/(Nsym*bits_per_symbol_for_4qam)%%%%%%%%%% 16QAM
demod_signal_for_16qam          = qamdemod(received_qam16,k_16qam,'gray','UnitAveragePower', true);
[Output]                        = IntegerToBits(demod_signal_for_16qam,bits_per_symbol_for_16qam);
received_Bits_for_16qam         = reshape(Output,1,Nsym*bits_per_symbol_for_16qam);
ErrorNumber_for_16qam           = size(find(received_Bits_for_16qam-binary_bits_for16qam),2);
BER_for_16qam(i)                = ErrorNumber_for_16qam/(Nsym*bits_per_symbol_for_16qam)%%%%%%%%%% 64QAM
demod_signal_for_64qam          = qamdemod(received_qam64,k_64qam,'gray','UnitAveragePower', true);
[Output]                        = IntegerToBits(demod_signal_for_64qam,bits_per_symbol_for_64qam);
received_Bits_for_64qam         = reshape(Output,1,Nsym*bits_per_symbol_for_64qam);
ErrorNumber_for_64qam           = size(find(received_Bits_for_64qam-binary_bits_for64qam),2);
BER_for_64qam(i)                = ErrorNumber_for_64qam/(Nsym*bits_per_symbol_for_64qam)%%%%%%%%%% 256QAM
demod_signal_for_256qam         = qamdemod(received_qam256,k_256qam,'gray','UnitAveragePower', true);
[Output]                        = IntegerToBits(demod_signal_for_256qam,bits_per_symbol_for_256qam);
received_Bits_for_256qam        = reshape(Output,1,Nsym*bits_per_symbol_for_256qam);
ErrorNumber_for_256qam          = size(find(received_Bits_for_256qam-binary_bits_for256qam),2);
BER_for_256qam(i)               = ErrorNumber_for_256qam/(Nsym*bits_per_symbol_for_256qam)
endfigure(2)
semilogy(EbN0_dB,BER_for_4qam,'g-<','MarkerFaceColor','g','LineWidth',1);hold on;
semilogy(EbN0_dB,BER_for_16qam,'m->','MarkerFaceColor','m','LineWidth',1);hold on;
semilogy(EbN0_dB,BER_for_64qam,'b-^','MarkerFaceColor','b','LineWidth',1);hold on;
semilogy(EbN0_dB,BER_for_256qam,'r-<','MarkerFaceColor','r','LineWidth',1);hold on;
axis([-5 20 1e-5 1]);
plot([-5 20],[3.8e-3,3.8e-3],'r','LineWidth',1);hold on;
grid on;
xlabel('E_{b}N_{0}(dB)');ylabel('BER');
legend('4qam','16qam','64qam ','256qam',...'FEC(3.8x10^-^3)','Location','SouthWest');function [Output] = BitsToInteger(input)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Convert Bits to Integer
% Author: Zhongxing Tian
% Version: 1.1
% Date: 2022/2/18
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%[M,N] = size(input);Value = zeros(1,N);for i = 2:MValue = Value+input(i,:).*2^(i-1); % 阶数从上到下（第一行最小
endValue = Value+input(1,:);Output = Value;
endfunction [Output] = IntegerToBits(input,bits_per_symbol)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Convert Integer to Bits
% Author: Zhongxing Tian
% Version: 1.1
% Date: 2022/2/18
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%[~,N] = size(input);Value = zeros(bits_per_symbol,N);for i = 1:Nk = 2^bits_per_symbol;for j = 1:bits_per_symbolif input(1,i)<k/2Value(bits_per_symbol-j+1,i) = 0;% 阶数从上到下（第一行最小）k = k-2^(bits_per_symbol-j);                elseValue(bits_per_symbol-j+1,i) = 1;k = k+2^(bits_per_symbol-j);endend
endOutput = Value;
end

4qam、16qam、64qam、256qam理论仿真曲线相关推荐

自适应均衡matlab仿真,对比RLS,LMS以及NLMS的均衡前后星座图效果,调制采用4QAM,16QAM,64QAM
目录 1.算法描述 2.仿真效果预览 3.MATLAB核心程序 4.完整MATLAB 1.算法描述在无线通信系统中,由于多径效应及码间干扰的存在,信号误码率会升高.均衡技术是一种对抗码间干扰的重要技 ...
基于MATLAB的自适应调制解调通信系统的误码率仿真,对比BPSK,QPSK,16QAM,64QAM
目录 1.算法仿真效果 2.MATLAB核心程序 3.算法涉及理论知识概要 4.完整MATLAB 1.算法仿真效果 matlab2022a仿真结果如下: 2.MATLAB核心程序 .......... ...
BPSK,QPSK,2FSK,16QAM,64QAM信号在高斯信道与瑞利信道下的误码率性能仿真
BPSK,QPSK,2FSK,16QAM,64QAM信号在高斯信道与瑞利信道下的误码率性能仿真_南大小王-CSDN博客 16QAM调制解调仿真(matlab,详细介绍仿真方案的设计.结果及结论.完整代 ...
基于ADS的c语言程序设计实验,实验一：基于ADS软件传输线理论仿真设计与分析.docx...
龙 lerrnl'' 龙 lerrnl'' $ Num=l 1Z=50 Ohm h|」M -- TL1 Z=50 0.0hm E=50 F=5GHz ggj I SrPARAM ETERS . | Z ...
matlab两个曲线的名称怎么显示出来的,求MATLAB的高手，怎么把两张仿真曲线合并显示在一张图上？...
求MATLAB的高手,怎么把两张仿真曲线合并显示在一张图上? 关注:194 答案:5 手机版解决时间 2021-02-10 14:48 提问者终究是陌生了 2021-02-09 16:38 程序 ...
凸轮设计c语言程序,凸轮理论轮廓曲线设计c语言程序.doc
/************习题4-1凸轮轮廓曲线绘制程序**********/ #include #include #define PI 3.1415926 void main() { double ...
R语言使用econocharts包创建微观经济或宏观经济图：需求曲线、供给曲线、供需曲线、无差异曲线、Laffer曲线、前景理论价值函数曲线、税收图、生产可能性边界线、劳动力供给曲线、计算曲线交叉点等
R语言使用econocharts包创建微观经济或宏观经济图:需求曲线.供给曲线.供需曲线.无差异曲线.Laffer曲线.前景理论价值函数曲线.税收图.生产可能性边界线.劳动力供给曲线.计算曲线交叉点等 ...
c语言凸轮轮廓曲线设计程序,凸轮理论轮廓曲线设计c语言程序.doc
凸轮理论轮廓曲线设计c语言程序.doc /习题 4-1 凸轮轮廓曲线绘制程序/includeincludedefine PI 3.1415926void maindouble e10.0,ro40.0 ...
4QAM、16QAM 调制与解调仿真电路，观察并分析QAM星座图和误码率曲线【matlab代码】
源码: https://download.csdn.net/download/qq_44394952/86236776 要求完成两种调制方式的调制解调:4QAM.16QAM (1)画出系统框图,搞清 ...

4qam、16qam、64qam、256qam理论仿真曲线

4qam、16qam、64qam、256qam理论仿真曲线相关推荐

最新文章

热门文章