【数字信号调制】16QAM数字信号调制解调(含误码率 星座图 时域 频谱)【含Matlab源码 4762期】

在这里插入图片描述

⛄一、获取代码方式

获取代码方式1
完整代码已上传我的资源:【数字信号调制】基于matlab 16QAM数字信号调制解调(含误码率 星座图 时域 频谱)【含Matlab源码 4762期】
点击上面蓝色字体,直接付费下载,即可。

获取代码方式2:
付费专栏Matlab信号处理(初级版)

备注:
点击上面蓝色字体付费专栏Matlab信号处理(初级版),扫描上面二维码,付费29.9元订阅海神之光博客付费专栏Matlab信号处理(初级版),凭支付凭证,私信博主,可免费获得1份本博客上传CSDN资源代码(有效期为订阅日起,三天内有效);
点击CSDN资源下载链接:1份本博客上传CSDN资源代码

⛄二、部分源代码

clear;
clc;
echo off;
close all;
%==========================================================================
% 参数设置
f=64e3; %调制后的信号符号率
fb=4f; %基带信号比特率
fs=32
fb; %抽样频率.一个点由32个点表示
fc=4*fb; %载波频率,为便于观察已调信号,我们把载波频率设的较低

Kbase=2; % Kbase=1,不经基带成形滤波,直接调制;
% Kbase=2,基带经成形滤波器滤波后,再进行调制
N=10000; %设定码元数量

%==========================================================================
%生成二进制信号序列,并作图
info=random_binary(N); %产生二进制信号序列
x=(1:N)/fb;
stem(x,info);
axis([0 100/fb 0 1.5])

%==========================================================================
% 基带信号调制,并做出已调信号频谱和星座图
[y,I,Q]=qam(info,Kbase,fs,fb,fc); %对基带信号进行16QAM调制
y1=y; y2=y; %备份信号,供后续仿真用
T=length(info)/fb;
m=fs/fb;
nn=length(info);
dt=1/fs;
t=0:dt:T-dt;

%画出已调信号波形
figure(2);
%便于观察,这里显示的已调信号及其频谱均为无噪声干扰的理想情况
%由于测试信号码元数量为10000个,在这里我们只显示其总数的1/10
plot(t(1:1000),y(1:1000),t(1:1000),I(1:1000),t(1:1000),Q(1:1000),[0 1000*dt],[0 0],‘b:’);
title(‘已调信号(In:red,Qn:green)’);

⛄三、运行结果

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 沈再阳.精通MATLAB信号处理[M].清华大学出版社,2015.
[2]高宝建,彭进业,王琳,潘建寿.信号与系统——使用MATLAB分析与实现[M].清华大学出版社,2020.
[3]王文光,魏少明,任欣.信号处理与系统分析的MATLAB实现[M].电子工业出版社,2018.
[4]李树锋.基于完全互补序列的MIMO雷达与5G MIMO通信[M].清华大学出版社.2021
[5]何友,关键.雷达目标检测与恒虚警处理(第二版)[M].清华大学出版社.2011

3 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除

  • 7
    点赞
  • 21
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
以下是16QAM调制解调误码率曲线、星座图、眼图的Matlab代码以及注释: %% 16QAM调制解调 % 假设我们要传输的数字信号为:[1 0 1 1 0 1 0 0] data = [1 0 1 1 0 1 0 0]; % 将数字信号转换为符号(即将2进制转换为16QAM星座图中的点) M = 16; % 星座图的大小 symbols = bi2de(reshape(data, [], log2(M)).','left-msb'); % 将二进制转换为十进制作为符号 % 显示16QAM星座图 scatterplot(qammod(0:M-1, M)); % 对符号进行16QAM调制 Fs = 1000; % 采样率 fc = 100; % 载波频率 t = 0:1/Fs:(length(symbols)-1)/Fs; % 时间轴 modulated = qammod(symbols, M); % 16QAM调制 transmitted = real(modulated .* exp(1j*2*pi*fc*t)); % 按照载波频率进行移频 % 在时域和频域中显示调制后的信号 figure; subplot(2,1,1); plot(t, transmitted); title('Modulated signal in time domain'); xlabel('Time (s)'); ylabel('Amplitude'); subplot(2,1,2); f = -Fs/2:Fs/length(transmitted):Fs/2-Fs/length(transmitted); plot(f, fftshift(abs(fft(transmitted)))); title('Modulated signal in frequency domain'); xlabel('Frequency (Hz)'); ylabel('Magnitude'); % 对调制后的信号进行16QAM解调 received = transmitted .* exp(-1j*2*pi*fc*t); % 进行移频还原 demodulated = qamdemod(received, M); % 16QAM解调 % 将解调后的符号转换为二进制 binary = reshape(de2bi(demodulated, log2(M)).', [], 1).'; disp(['Original data:', num2str(data)]); disp(['Demodulated data:', num2str(binary)]); %% 16QAM误码率曲线 % 生成随机的数字信号 data = randi([0 1], 1, 1000); % 将数字信号转换为符号 M = 16; % 星座图的大小 symbols = bi2de(reshape(data, [], log2(M)).','left-msb'); % 将二进制转换为十进制作为符号 % 对符号进行16QAM调制 Fs = 1000; % 采样率 fc = 100; % 载波频率 t = 0:1/Fs:(length(symbols)-1)/Fs; % 时间轴 modulated = qammod(symbols, M); % 16QAM调制 transmitted = real(modulated .* exp(1j*2*pi*fc*t)); % 按照载波频率进行移频 % 添加高斯白噪声 EbNo = 0:2:20; % 信噪比范围 ber = zeros(size(EbNo)); % 误码率 for i = 1:length(EbNo) snr = EbNo(i) + 10*log10(log2(M)); % 转换为信噪比 noisy = awgn(transmitted, snr, 'measured'); % 添加高斯白噪声 received = noisy .* exp(-1j*2*pi*fc*t); % 进行移频还原 demodulated = qamdemod(received, M); % 16QAM解调 binary = reshape(de2bi(demodulated, log2(M)).', [], 1).'; % 将解调后的符号转换为二进制 [~, ber(i)] = biterr(binary, data); % 计算误码率 end % 绘制误码率曲线 figure; semilogy(EbNo, ber); title('16QAM Bit Error Rate'); xlabel('Eb/No (dB)'); ylabel('Bit Error Rate'); %% 16QAM星座图 % 显示16QAM星座图 scatterplot(qammod(0:15, 16)); %% 16QAM眼图 % 生成随机的数字信号 data = randi([0 1], 1, 1000); % 将数字信号转换为符号 M = 16; % 星座图的大小 symbols = bi2de(reshape(data, [], log2(M)).','left-msb'); % 将二进制转换为十进制作为符号 % 对符号进行16QAM调制 Fs = 1000; % 采样率 fc = 100; % 载波频率 t = 0:1/Fs:(length(symbols)-1)/Fs; % 时间轴 modulated = qammod(symbols, M); % 16QAM调制 transmitted = real(modulated .* exp(1j*2*pi*fc*t)); % 按照载波频率进行移频 % 绘制眼图 eyediagram(transmitted, Fs/16); % 每个符号的采样点数为16

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Matlab领域

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值