可以写一下QPSK调制解调MATLAB仿真代码

当然可以，以下是QPSK调制解调MATLAB仿真代码：

调制端代码：

clear all;
close all;
clc;

%定义调制参数
M = 4; %符号数目
k = log2(M); %比特数目
n = 10000; %符号数
num_sample = 16; %每个符号的采样数
SNR_dB = 15; %信噪比

%生成随机比特流
data = randi([0,1],n,k);

%将比特流转换为符号序列
data_symbols = bi2de(data,'left-msb');

%将符号序列进行QPSK调制
constellation = qammod(data_symbols,M,'gray');

%将符号序列并行转换为串行
tx = reshape(constellation,[],1);

%将串行信号进行升采样
tx_upsampled = upsample(tx,num_sample);

%定义发送端滤波器
rolloff = 0.25;
span = 10;
sps = num_sample;
filter_tx = rcosdesign(rolloff,span,sps);

%将升采样后的串行信号通过发送端滤波器
tx_filtered = filter(filter_tx,1,tx_upsampled);

%定义通道传输过程
SNR = 10^(SNR_dB/10);
n0 = 1/SNR;
noise = sqrt(n0/2)*randn(length(tx_filtered),1);
rx_signal = tx_filtered + noise;

%接收端滤波器
filter_rx = filter_tx;

%将接收到的信号进行接收端滤波
rx_filtered = filter(filter_rx,1,rx_signal);

%将接收端滤波后的信号进行降采样
rx_downsampled = downsample(rx_filtered,num_sample);

%将接收端降采样后的信号进行QPSK解调
rx_symbols = qamdemod(rx_downsampled,M,'gray');

%将解调后的符号进行并行转换为比特流
rx_data = de2bi(rx_symbols,'left-msb');

%计算误码率
[number,ber] = biterr(data,rx_data);

%输出结果
fprintf('错误比特数：%d\n',number);
fprintf('误码率：%d\n',ber);

解调端代码：

clear all;
close all;
clc;

%定义调制参数
M = 4; %符号数目
k = log2(M); %比特数目
n = 10000; %符号数
num_sample = 16; %每个符号的采样数
SNR_dB = 15; %信噪比

%生成随机比特流
data = randi([0,1],n,k);

%将比特流转换为符号序列
data_symbols = bi2de(data,'left-msb');

%将符号序列进行QPSK调制
constellation = qammod(data_symbols,M,'gray');

%将符号序列并行转换为串行
tx = reshape(constellation,[],1);

%将串行信号进行升采样
tx_upsampled = upsample(tx,num_sample);

%定义发送端滤波器
rolloff = 0.25;
span = 10;
sps = num_sample;
filter_tx = rcosdesign(rolloff,span,sps);

%将升采样后的串行信号通过发送端滤波器
tx_filtered = filter(filter_tx,1,tx_upsampled);

%定义通道传输过程
SNR = 10^(SNR_dB/10);
n0 = 1/SNR;
noise = sqrt(n0/2)*randn(length(tx_filtered),1);
rx_signal = tx_filtered + noise;

%接收端滤波器
filter_rx = filter_tx;

%将接收到的信号进行接收端滤波
rx_filtered = filter(filter_rx,1,rx_signal);

%将接收端滤波后的信号进行降采样
rx_downsampled = downsample(rx_filtered,num_sample);

%将接收端降采样后的信号进行QPSK解调
rx_symbols = qamdemod(rx_downsampled,M,'gray');

%将解调后的符号进行并行转换为比特流
rx_data = de2bi(rx_symbols,'left-msb');

%计算误码率
[number,ber] = biterr(data,rx_data);

%输出结果
fprintf('错误比特数：%d\n',number);
fprintf('误码率：%d\n',ber);

注意事项：

1. 本代码只适用于理解QPSK调制解调的基本原理，具体应用中需要根据实际情况进行修改。
2. 本代码中使用的是QPSK调制解调的基本流程，没有做任何通道编码和解码操作。在实际应用中，需要根据通信系统的需要进行编解码操作。
3. 本代码中使用的是理想信道模型，没有考虑多径效应、噪声等现实因素，实际应用中需要根据实际情况进行修改。