4qam调制与解调代码MATLAB，画出星座图，点图，计算误码率

以下是4QAM调制解调的MATLAB代码：

调制：

% 4QAM调制
M = 4; % 星座大小
N = 10000; % 数据点数
dataIn = randi([0 M-1], N, 1); % 随机数据生成
dataSym = qammod(dataIn, M); % QAM调制
scatterplot(dataSym); % 星座图

解调：

% 4QAM解调
dataOut = qamdemod(dataSym, M); % QAM解调
[num, ratio] = biterr(dataIn, dataOut); % 计算误码率
fprintf("误码率为 %d\n", ratio); % 输出误码率

其中，`qammod` 和 `qamdemod` 分别是4QAM调制和解调函数，`scatterplot` 可以画出星座图，`biterr` 可以计算误码率。

运行代码后，将会得到星座图和误码率输出。

点图可以通过以下代码实现：

% 4QAM点图
figure;
plot(dataSym, '-o');
title('4QAM点图');
xlabel('实部');
ylabel('虚部');

相关问题

4qam调制与解调代码MATLAB，画出星座图，眼图，计算误码率

以下是4QAM调制与解调的MATLAB代码：

调制：

% 4QAM调制
clear all; close all; clc;
N = 1000; % 信号长度
M = 4; % 码元数
% 生成随机整数序列
data = randi([0 M-1],1,N);
% 显示原始数据
figure;
stem(data);
title('原始数据');
% QAM调制
x = qammod(data,M); % 调制
% 显示星座图
figure;
plot(x,'o');
title('星座图');

解调：

% 4QAM解调
% QAM解调
data_hat = qamdemod(x,M); % 解调
% 显示解调数据
figure;
stem(data_hat);
title('解调数据');

画出星座图：

% 显示星座图
figure;
plot(x,'o');
title('星座图');

计算误码率：

% 计算误码率
error_rate = sum(data~=data_hat)/N;
fprintf('误码率为：%f\n',error_rate);

画出眼图：

% 画出眼图
eyediagram(x,4);
title('眼图');

完整代码如下：

% 4QAM调制与解调代码MATLAB，画出星座图，眼图，计算误码率
clear all; close all; clc;
N = 1000; % 信号长度
M = 4; % 码元数
% 生成随机整数序列
data = randi([0 M-1],1,N);
% 显示原始数据
figure;
stem(data);
title('原始数据');
% QAM调制
x = qammod(data,M); % 调制
% 显示星座图
figure;
plot(x,'o');
title('星座图');
% QAM解调
data_hat = qamdemod(x,M); % 解调
% 显示解调数据
figure;
stem(data_hat);
title('解调数据');
% 计算误码率
error_rate = sum(data~=data_hat)/N;
fprintf('误码率为：%f\n',error_rate);
% 画出眼图
eyediagram(x,4);
title('眼图');

运行以上代码可以得到4QAM调制与解调的结果，包括星座图、眼图和误码率等。

16qam 调制解调 并画出16qam星座图，误码率曲线，眼图 matlab代码仿真

16QAM调制解调的MATLAB代码如下：

调制：

% 设置参数
M = 16; % 星座点数
fc = 5000; % 载波频率
fs = 100000; % 采样频率
T = 1/fs; % 采样时间间隔
numBits = 1000; % 待调制的比特数
nSamp = 4; % 每个符号的采样数
EbNo = 10; % 信噪比（dB）
EsNo = 10*log10(M) + EbNo; % 符号噪声比（dB）
sigma = 1/sqrt(2*EsNo); % 噪声标准差

% 生成随机比特流
bitsIn = randi([0 1],1,numBits);

% 将比特流分组成符号
M = 16;
k = log2(M);
symbolsIn = bi2de(reshape(bitsIn,length(bitsIn)/k,k),'left-msb');

% 生成16QAM星座图
constellation = qammod(0:M-1,M);

% 16QAM调制
txSig = qammod(symbolsIn,M,'gray');

% 将符号映射到星座图上
txSigMapped = constellation(txSig+1);

% 为星座图添加噪声
noise = sigma*randn(size(txSigMapped));
rxSig = txSigMapped + noise;

% 绘制星座图
scatterplot(constellation); hold on;
plot(real(txSigMapped),imag(txSigMapped),'ko');
plot(real(rxSig),imag(rxSig),'r.');
title('16QAM Constellation');
xlabel('In-Phase');
ylabel('Quadrature');
grid on;

解调：

% 16QAM解调
rxSigMapped = rxSig ./ constellation;

% 将解调后的符号映射到星座图上
[~,rxSymbols] = min(abs(rxSigMapped.' - constellation),[],2);
rxBits = de2bi(rxSymbols-1,k,'left-msb').';
rxBits = rxBits(:).';

% 统计误码率
numErrs = sum(rxBits ~= bitsIn);
ber = numErrs / length(bitsIn);

% 绘制误码率曲线
EbN0 = 0:2:20;
theoryBer = berawgn(EbN0,'qam',M);
semilogy(EbN0,theoryBer,'b-','LineWidth',2); hold on;
semilogy(EbNo,ber,'ro','LineWidth',2);
title('16QAM BER');
xlabel('Eb/No (dB)');
ylabel('Bit Error Rate');
legend('Theory','Simulation');
grid on;

% 绘制眼图
eyediagram(rxSigMapped,nSamp*2);
title('16QAM Eye Diagram');
xlabel('Time');
ylabel('Amplitude');
grid on;

注意，在以上代码中，`qammod`和`qamdemod`函数可以分别用于16QAM调制和解调。眼图可以使用`eyediagram`函数来绘制。