帮我写一个用匹配滤波器进行接收的QPSK系统MATLAB代码

以下是一个用匹配滤波器进行接收的QPSK系统的MATLAB代码：

%% 设置参数
fc = 10e6; % 载波频率
fs = 40e6; % 采样率
T = 1/fs; % 采样间隔
N = 1000; % 数据点数
SNR = 10; % 信噪比
M = 4; % 调制方式，这里为QPSK
Eb = 1; % 符号能量
Ts = 1/fs; % 符号间隔
fcut = 2*fc; % 截止频率

%% 生成QPSK信号
data = randi([0 M-1],N,1); % 生成随机数据
symbols = qammod(data,M,'gray'); % QPSK调制
txSig = reshape(symbols,N/2,2); % 两个符号为一组，组成复数信号
txSig = txSig(:,1) + 1j*txSig(:,2); % 转换为复数信号

%% 将信号通过带通滤波器
t = (0:N-1)*T; % 时间序列
carrier = cos(2*pi*fc*t); % 产生载波
txSig = real(txSig.*carrier'); % 将QPSK信号与载波相乘
[b,a] = butter(6,fcut/(fs/2)); % 产生带通滤波器
txSig = filter(b,a,txSig); % 将信号通过带通滤波器

%% 添加高斯噪声
rxSig = awgn(txSig,SNR); % 添加高斯噪声

%% 将信号通过匹配滤波器
h = fliplr(conj(carrier)); % 产生匹配滤波器
rxSig = filter(h,1,rxSig); % 将信号通过匹配滤波器

%% 对接收信号进行QPSK解调
rxSig = rxSig(2:end); % 去除匹配滤波器延迟
rxSymbols = rxSig(1:2:end) + 1j*rxSig(2:2:end); % 将信号分为两路
rxData = qamdemod(rxSymbols,M,'gray'); % 对解调后的信号进行解调

%% 计算误码率
BER = sum(rxData~=data)/N; % 计算误码率

%% 绘制信号图形
subplot(2,1,1);
plot(t,txSig);
title('发送信号');
xlabel('时间(s)');
ylabel('幅度(V)');

subplot(2,1,2);
plot(t,rxSig);
title('接收信号');
xlabel('时间(s)');
ylabel('幅度(V)');

%% 打印误码率
fprintf('误码率为：%f\n',BER);